Research and application of a novel grey multivariable model in port scale prediction under the impact of Free Trade Zone

Yuyu Sun (School of Economics, Ocean University of China, Qingdao, China)
Yuchen Zhang (School of Economics, Ocean University of China, Qingdao, China)
Zhiguo Zhao (School of Economics, Ocean University of China, Qingdao, China)

Marine Economics and Management

ISSN: 2516-158X

Article publication date: 14 May 2024

36

Abstract

Purpose

Considering the impact of the Free Trade Zone (FTZ) policy on forecasting the port cargo throughput, this paper constructs a fractional grey multivariate forecasting model to improve the prediction accuracy of port cargo throughput and realize the coordinated development of FTZ policymaking and port construction.

Design/methodology/approach

Considering the effects of data randomization, this paper proposes a novel self-adaptive grey multivariate prediction model, namely FDCGM(1,N). First, fractional-order accumulative generation operation (AGO) is introduced, which integrates the policy impact effect. Second, the heuristic grey wolf optimization (GWO) algorithm is used to determine the optimal nonlinear parameters. Finally, the novel model is then applied to port scale simulation and forecasting in Tianjin and Fujian where FTZs are situated and compared with three other grey models and two machine learning models.

Findings

In the Tianjin and Fujian cases, the new model outperforms the other comparison models, with the least mean absolute percentage error (MAPE) values of 6.07% and 4.16% in the simulation phase, and 6.70% and 1.63% in the forecasting phase, respectively. The results of the comparative analysis find that after the constitution of the FTZs, Tianjin’s port cargo throughput has shown a slow growth trend, and Fujian’s port cargo throughput has exhibited rapid growth. Further, the port cargo throughput of Tianjin and Fujian will maintain a growing trend in the next four years.

Practical implications

The new multivariable grey model can effectively reduce the impact of data randomness on forecasting. Meanwhile, FTZ policy has regional heterogeneity in port development, and the government can take different measures to improve the development of ports.

Originality/value

Under the background of FTZ policy, the new multivariable model can be used to achieve accurate prediction, which is conducive to determining the direction of port development and planning the port layout.

Keywords

Citation

Sun, Y., Zhang, Y. and Zhao, Z. (2024), "Research and application of a novel grey multivariable model in port scale prediction under the impact of Free Trade Zone", Marine Economics and Management, Vol. ahead-of-print No. ahead-of-print. https://doi.org/10.1108/MAEM-03-2024-0005

Publisher

:

Emerald Publishing Limited

Copyright © 2024, Yuyu Sun, Yuchen Zhang and Zhiguo Zhao

License

Published in Marine Economics and Management. Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) licence. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this licence may be seen at http://creativecommons.org/licences/by/4.0/legalcode


1. Introduction

China has been deepening its basic state policy of open-door to the outside world and exploring establishing a Free Trade Zone (FTZ) for the past few years. It aims to facilitate investment and foreign trade through the construction of multifunctional special economic zones, probe new ways and accumulate new experiences for the comprehensive deepening of reform (Yao and Whalley, 2016). Since 2013, China has laid out 21 FTZs across the country through multiple batches, initially forming a basic pattern of “1 + 3+7 + 1+6 + 3.” Each FTZ is situated in a significant economic region, serving specific economic zones and having distinct development positioning. With the gradual improvement of the construction layout of the FTZ, a reform and opening-up innovation pattern covering the east, west, north, south and centre has been formed. Bold exploration and remarkable results have been achieved in areas such as investment, foreign trade and financial services for the real economy (Chen et al., 2022; Wang et al., 2022a).

Ports are important transportation hubs and open gateways and play an essential role in accelerating international trade and creating economic benefits (Gonzalez Aregall et al., 2018). Its development is of great significance to the construction of FTZs while establishing FTZs will be conducive to port improvement. Port development studies indicate port cargo throughput as an essential indicator of port production capacity and operation. It is also the basis for port planning and construction, which can reflect the modernization level of ports and the degree of foreign trade (Wiegmans et al., 2015). However, due to the persistent global economic downturn, the size of Chinese ports has been saturated, and their development prospects face many challenges and risks (Chen et al., 2020; Wan et al., 2021). Then, in the context of FTZ policy, improving the prediction exactness of port cargo throughput and realizing precise forecasting are conducive to determining port development direction, planning the port layout and better planning of the port logistics industry. Thus, policy formulation and port construction in the FTZs can be realized through synergistic development.

In summary, considering that the port cargo throughput data is difficult to obtain and affected by many factors, this paper constructs a fractional grey multivariate forecasting model considering the impact of the FTZ policy to forecast the port cargo throughput. The three main contributions of this study are as follows:

  1. Although the traditional grey multivariate model has good forecasting performance for time series data, the problem of forecasting time series with randomness has not been perfectly solved. In this paper, a new grey multivariate forecasting model (FDCGM(1,N)) is constructed based on the GM(1,N) model by introducing the fractional-order AGO and a generalized time response function with a forecasting function is derived. Meanwhile, optimal nonlinear parameters of the new model are obtained by the GWO algorithm, which significantly enhances the description of the new model.

  2. FDCGM(1,N) model is used to simulate and forecast port cargo throughput of Tianjin and Fujian Province implementing FTZ policy. Importantly, compared with three other grey models and two machine learning models, the new model has the smallest mean absolute percentage errors (MAPEs) in the in-sample prediction interval, which are 6.70% and 1.63%, respectively.

  3. This paper applies the FDCGM(1,N) model to forecast and analyse the port cargo throughput of Tianjin and Fujian Province from 2022 to 2025. It is found that port development under the FTZ policy is regionally heterogeneous. It is estimated that in 2025, the port cargo throughput of Tianjin will reach 762,348,800 tons and Fujian province’s port cargo throughput will reach 1,017,638,300 tons.

The rest of the paper is organized as follows. Section 2 details previous studies on FTZS and ports as well as grey models. Section 3 introduces the model evaluation method, the parameter solution algorithm and the research framework. Section 4 applies the FDCGM(1,N) model to the port cargo throughput forecasts of Tianjin and Fujian province implementing FTZ policy and discusses the port cargo growth rates and the forecast results. Section 5, then, presents the conclusions and the future outlook.

2. Literature review

As a functional node for the flow of goods trade, ports are the ligament of the entire supply chain of FTZs and the centre of the flow of goods, capital and information (Wan et al., 2014). Consequently, many scholars have explored the impact of FTZ policy on port development. For example, Chen et al. (2018) constructed a system to assess the development performance of typical FTZs. They conducted an empirical study on six selected distinctive FTZs. Liu et al. (2021) evaluated port efficiency in the context of FTZ policy, identified its influencing factors and proved that the scale of port operation significantly positively impacts port efficiency. Li et al. (2021) studied the impact of FTZ policy on the evolution of port-listed companies in FTZs. They found that port-listed companies’ development performance continuously grew under the FTZ policy. Considering the port throughput and the scale of urban import and export trade, Fan et al. (2022) selected four ports established earlier in FTZs and used the allometric growth model and linear scale factor to empirically test the FTZ policy’s impact on the ports.

In addition, a lot of scholars have also conducted in-depth research on the prediction of port cargo throughput. At present, the prediction methods of port cargo throughput mainly include traditional econometric model, artificial intelligence model and grey model. In terms of traditional econometric models, Sanguri et al. (2022) proposed an intertemporal forecasting model based on exponential smoothing to forecast container throughput at the Port of Los Angeles. In terms of the artificial intelligence model, Cuong et al. (2022) utilized a neural network predictive controller and adaptive fractional-order supervision sliding mode control to handle throughput under external disturbances. In the meantime, some scholars have also used the decomposition integration method (Du et al., 2019; Jin et al., 2023), SARIMA and machine learning hybrid method (Huang et al., 2022; Mo et al., 2018) for forecasting container throughput prediction. In fact, port cargo throughput data is usually provided by port authorities, customs, shipping companies, etc., and its availability is limited especially when some ports may be unwilling or unable to disclose it. As a result, the accuracy of the model predictions mentioned above, which require a large amount of data support, may be compromised. For this reason, some scholars have used grey forecasting models to predict port cargo throughput. For example, Hsu et al. (2020) proposed a new model combining grey share prediction with Markov chain and grey residual correction, and proved that the model has good prediction performance.

Deng (1982) first proposed the grey theory. Based on this theory, later scholars optimized the grey univariate prediction model from various perspectives, such as model structure (Yang and Wang, 2022; Zeng et al., 2020c), initial value optimization (Zeng et al., 2020b), background value (Wu et al., 2020a) and cumulative generation operator (Wang et al., 2022c; Zeng et al., 2023a, b), and achieved better prediction results. Currently, grey prediction models have been widely used to deal with the prediction problems in the fields of energy price (Duan et al., 2022; Duan and Liu, 2021), air quality (Du et al., 2021; Shi and Wu, 2020; Zeng et al., 2021), food production (Zeng et al., 2020a, b, c), new energy vehicle sales (Liu et al., 2022), energy consumption (Liu and Wu, 2021; Moonchai and Chutsagulprom, 2020; Wu et al., 2020b), output of high-technology industries (Ding et al., 2022) and epidemic disease transmission (Saxena, 2021). However, considering the influence of external factors, the multivariable grey model is gradually being used more and more to compensate for the shortcomings of univariate models. As the most basic multivariate grey model, the GM(1,N) model is widely used (Duan and Luo, 2022; Zeng et al., 2023a, b; Zhang et al., 2022). Unfortunately, the model also has some problems in the modelling process (Ren et al., 2023; Zeng et al., 2019). Therefore, many scholars optimize the GM(1,N) model from different perspectives to enhance the applicability of the model.

In terms of dummy variables, Ding et al. (2018) added dummy variables into the GM (1,N) model, gave a specific modelling approach, and verified the availability of the neoteric model by examples. Wan et al. (2022) introduced dummy variables’ interactions into Ding’s model to establish a new grey multivariate model, which improved the model’s prediction accuracy. Meanwhile, considering the impact of fractional-order accumulative generation, they introduced fractional-order into the new model to enhance the model’s simulation and forecasting performance for time series with randomness. For instance, Wang et al. (2022b) introduced a new kind of fractional order, namely Grünwald-Letnikov fractional-order calculus, to enhance model adaptability. Yan et al. (2022) introduce fractional-order cumulative generation into the multivariable time-delayed grey model to reduce the impact of randomness of online public opinion data on the prediction results and optimize the modelling parameters using particle swarm optimization.

3. Methods

3.1 The existing GM(1,N) model

In this paper, the GM(1,N) model is selected as the benchmark model. In order to fully demonstrate the new model proposed, the GM(1,N) model is briefly introduced at first.

Definition 1.

Zeng et al. (2016) Let the original sequences X1(0) and Xi(0)(i=2,3,,N) denote the system behaviour sequences and driver sequences, respectively, and their 1-order accumulative generation sequences X1(1) and Xi(1), respectively, and Z1(1) is called the background value and is generated by X1(1) where Z1(1)(k)=12×[X1(1)(k)+X1(1)(k1)],k=2,3,,n.

Then,

(1)x1(0)(k)+az1(0)(k)=i=2Nbixi(1)(k)
is the GM(1,N) model. Where a is the system development factor, bixi(1)(k) is the driving term, and bi is the driving coefficient.
Definition 2.

Let aˆ=[a,b2,,bN]T be the parameter column of the model. Then,

(2)dx(1)(t)dt+ax(1)(t)=i=2Nbixi(1)(k)
is the whitening equation of the GM(1,N) model.

3.2 The proposed FDCGM(1,N) model

In this section, a new grey multivariate model is proposed, referred to as the FDCGM(1,N) model. Before the new model is introduced, the fractional-order cumulative generating series and the cumulative-decreasing generating series are first introduced.

Definition 3.

Let sequence Xi(0) be as described in Definition 1. Then, the sequence Xi(r)=(xi(r)(1),xi(r)(2),,xi(r)(n)) is known as the r-order accumulation of Xi(0), and the sequence Xi(r)=(xi(r)(1),xi(r)(2),,xi(r)(n)) is known as the r-order inverse accumulation of Xi(0). Specifically,

(3)xi(r)(k)=s=1kΓ(r+ks)Γ(ks+1)Γ(r)xi(0)(s),k=1,2,,n.
(4)xi(r)(k)=s=0k1(1)sΓ(r+1)Γ(s+1)Γ(rs+1)xi(0)(ks),k=1,2,,n.
Definition 4.

Let the sequences Xi(0) and Xi(r) be as described in Definition 1 and Definition 3, respectively, and the sequence of dummy variables Dj(0)=(dj(0)(1),dj(0)(2),,dj(0)(n)),dj(0)(k)=0 or 1, which the r-order accumulative generation sequence is Dj(r)(j=M+1,M+2,,N). Z1(r) is the background value sequence of X(r), which is Z1(r)=(z1(r)(1),z1(r)(2),,z1(r)(n)), z1(r)(k)=12×[x1(r)(k)+x1(r)(k1)].

Then,

(5)x1(r)(k)x1(r)(k1)+az1(r)(k)=i=2Mbixi(r)(k)+j=M+1Nbjdj(r)(k)
is basic form of FDCGM(1,N) model.

i=2Mbixi(r)(k) is the quantified variable driver; j=M+1Nbjdj(r)(k) is the dummy variable driver, which takes into account the influence of qualitative factors on the dependent variable; bj is the driving coefficient of the accumulative generation term of the dummy variable.

Theorem 1.

Let Xi(0), Xi(r), Dj(0) and Dj(r) be as described in Definition 1, Definition 3 and Definition 4. The model parameters are listed as bˆ=[a,b2,,bM,,bN]T are known. Where,

(6)Y=(x1(r)(2)x1(r)(1)x1(r)(3)x1(r)(2)x1(r)(n)x1(r)(n1))B=(z1(r)(2)x2(r)(2)xM(r)(2)dM+1(r)(2)dN(r)(2)z1(r)(3)x2(r)(3)xM(r)(3)dM+1(r)(3)dN(r)(3)z1(r)(n)x2(r)(n)xM(r)(n)dM+1(r)(n)dN(r)(n))

Then, the least squares of the parameter columns satisfy:

  • (1) When n=N+1, bˆ=B1Y,|B|0;

  • (2) When n>N+1, bˆ=(BTB)1BTY,|BTB|0;

  • (3) When n<N+1, bˆ=BT(BTB)1Y,|BTB|0.

Proof. Bringing k=2,3,,n into the FDCGM(1,N) model, we get

(7)x1(r1)(2)=ax1(r1)(2)+i=2Mbixi(r)(2)+j=M+1Nbjdj(r)(2)x1(r1)(3)=ax1(r1)(3)+i=2Mbixi(r)(3)+j=M+1Nbjdj(r)(3)x1(r1)(n)=ax1(r1)(n)+i=2Mbixi(r)(n)+j=M+1Nbjdj(r)(n)

That is, from the least squares method, we have Y=Bbˆ.

  1. When n=N+1 and |B|0, the inverse matrix of B exists and the system of equations has a unique solution, we can obtain bˆ=B1Y.

  2. When n>N+1 and B is column-full rank, there is a column-full rank solution of B to B=DC. In turn, the generalized matrix B+ of B can be obtained as

B+=CT(CCT)1(DTD)1DT,βˆ=CT(CCT)1(DTD)1DTY
Since B is a full-rank matrix, C can be taken as a unit matrix, B=DIN, B=D, we get
bˆ=CT(CCT)1(DTD)1DTY=(DTD)1DTY=(BTB)1BTY
  1. When n<N+1 and B is a row full-rank matrix, D can be taken as a unit matrix, B=In1C, B=C, then we gain

bˆ=CT(CCT)1(DTD)1DTY=CT(CCT)1Y=BT(BBT)1Y
Theorem 2.

Let Xi(0), Xi(r), Dj(0) and Dj(r) be as described in Definition 1, Definition 3 and Definition 4, we have:

  • (1) The solution of the whitening differential equation of the FDCGM(1,N) model is

(8) x1(r)(t)=eat[x1(r)(0)t(i=2Mbixi(r)(0)+j=M+1Nbjdj(r)(0))+i=2N(i=2Mbixi(r)(t)+j=M+1Nbjdj(r)(t))eatdt]
  • (2) When the magnitude of change in the driver sequence Xi(1)(i=2,3,,N) is small, the model drivers i=2Mbixi(1)(t) and j=M+1Nbjdj(1)(t) can be considered as grey constants. The approximate time corresponding function sequence of the grey differential equation of the FDCGM(1,N) model is:

(9) xˆ1(r)(k)=ea(k1)[x1(r)(1)(i=2Mbixi(r)(k)+j=M+1Nbjdj(r)(k))1a]+(i=2Mbixi(r)(k)+j=M+1Nbjdj(r)(k))1a
  • (3) The time response equation of xˆ1(r)(k) can be obtained by the r-order accumulative generation of xˆ1(0)(k):

(10) xˆ1(0)(k)={x1(0)(1)k=1(xˆ1(r))(r)(k)=j=0k1(1)jΓ(r+1)Γ(j+1)Γ(rj+1)xˆ1(r)(kj)k=2,3,,n

Proof.

  1. From the whitening equation, we can obtain the general solution equation:

(11)x1(r)(t)=eat[i=2N(i=2Mbixi(r)(t)+j=M+1Nbjdj(r)(t))eatdt+e]
Where e is a constant to be determined. Bringing xi(r)(0) into the above equation yields, we can get
(12)e=x1(r)(0)i=2N(i=2Mbixi(r)(t)+j=M+1Nbjdj(r)(t))eatdt

Then, Eq. (8) is proved.

  • (2) Let the model drivers i=2Mbixi(1)(t) and j=M+1Nbjdj(1)(t) be considered as grey constants, and then the approximate time corresponding function sequence of the grey differential equation of the FDCGM(1,N) model is

xˆ1(r)(k)=ea(k1)[x1(r)(1)(i=2Mbixi(r)(k)+j=M+1Nbjdj(r)(k))1a]+(i=2Mbixi(r)(k)+j=M+1Nbjdj(r)(k))1a
  • (3) Based on (2),

When k=1, xˆ1(0)(1)=x1(0)(1).

When k=2,3,,n, by the r-order inverse accumulative generation operator (I-AGO)

xi(r)(k)=s=0k1(1)sΓ(r+1)Γ(s+1)Γ(rs+1)xi(0)(ks),k=1,2,,n.

We can get

(13)xˆ1(0)(k)=(xˆ1(r))(r)(k)=s=0k1(1)sΓ(r+1)Γ(s+1)Γ(rs+1)xˆ1(r)(ks)

In summary, Eq. (10) is proved.

Theorem 3.

When bj and r in the FDCGM(1,N) model take disparate values, it can be transformed into existing grey models, including GM(1,N) model, FGM(1,N) model and DVCGM(1,N) model.

Proof.

  1. When bj=0 and r=0, the new model can be transferred to the GM(1,N) model (Zeng et al., 2016).

(14)dx1dt+ax1(t)=i=2Mbixi(t)
  1. When bj=0 and r0, the new model can be turned to the FGM(1,N) model (Wang and Li, 2020).

(15)dx1(r)dt+ax1(r)(t)=i=2Mbixi(r)(t)
  1. When bj0 and r=0, the new model can be the DVCGM(1,N) model (Ding et al., 2018).

(16)dx1dt+ax1(t)=i=2Mbixi(t)+j=M+1Nbjdj(t)

3.3 Optimization of the fractional order based on GWO

In the FDCGM(1,N) model, the fractional order r needs to be optimized. It can be solved by constructing a nonlinear optimization problem and choosing MAPEs as the objective function. The objective function and each constraint condition can be expressed as:

(17)minrMAPES=1m1k=1m|x1(0)(k)xˆ1(0)(k)x1(0)(k)|×100%
s.t.{bˆ=[a,b2,,bM,,bN]TB=(z1(r)(2)x2(r)(2)xM(r)(2)dM+1(r)(2)dN(r)(2)z1(r)(3)x2(r)(3)xM(r)(3)dM+1(r)(3)dN(r)(3)z1(r)(n)x2(r)(n)xM(r)(n)dM+1(r)(n)dN(r)(n))Y=(x1(r)(2)x1(r)(1)x1(r)(3)x1(r)(2)x1(r)(n)x1(r)(n1))xˆ1(r)(k)=ea(k1)[x1(r)(1)(i=2Mbixi(r)(k)+j=M+1Nbjdj(r)(k))1a]+(i=2Mbixi(r)(k)+j=M+1Nbjdj(r)(k))1axˆ1(0)(k)=(xˆ1(r))r(k)

The above objective function and constraints suggest that parameter optimization of FDCGM(1,N) model is a nonlinear optimization problem with multiple nonlinear constraints. In this paper, we will apply the grey wolf optimization (GWO) algorithm to find new model’s optimal fractional order r. According to previous studies (Mirjalili et al., 2014), the steps of GWO include social stratification, follow the prey, and find prey and attack it.

3.4 Performance metrics

This section will define the following two metrics to estimate the exactness of the FDCGM(1,N) model.

The first one is absolute percentage error (APE), which is set to be represented concretely:

(18)APE(k)=|x1(0)(k)xˆ1(0)(k)x1(0)(k)|×100%

The second one is the MAPE. It can reflect the difference between fitting and original values that APE cannot measure. Suppose there is a total of raw data used as modelling data to judge the model’s simulation accuracy. The remaining data are applied as prediction data to evaluate the model’s prediction veracity. The following two equations can express the average percentage error of simulation and prediction:

(19)MAPES=1m1k=1m|x1(0)(k)xˆ1(0)(k)x1(0)(k)|×100%
(20)MAPEP=1nmk=m+1n|x1(0)(k)xˆ1(0)(k)x1(0)(k)|×100%

In the above equation, x(0)(k) and xˆ(0)(k) represent the actual and predicted values of the original data, respectively.

3.5 Modelling process

According to the modelling idea of the FDCGM(1,N) model, this section elaborates on the modelling steps (see Figure 1).

  • Step 1: Selection of model variables. Based on practical application cases, policy dummy variables are identified, and three main influencing factors of port cargo throughput are decided by grey correlation analysis.

  • Step 2: Based on the selected variables, the original data sequence Xi(0) is created and the corresponding r-order accumulating generation sequence is Xi(r).

  • Step 3: The FDCGM(1,N) model is developed according to Theorem 1, and the model parameters are estimated. The time response equation is calculated according to Theorem 2 and used for subsequent simulations and predictions.

  • Step 4: The nonlinear parameter optimization problem is established, and GWO obtains the optimal fractional order with the minimization of the mean percentage error of the simulation as the objective function.

  • Step 5: Sequence X1(r) is generated by r-order I-AGO to obtain the original data sequence’s simulated and predicted sequence X1(0). If the new model’s MAPEs and MAPEp meet the prediction accuracy requirements, they can be used for out-of-sample prediction for practical applications. Otherwise, the model needs to be improved.

  • Step 6: The new model through evaluation is applied to predict the port cargo throughput of Tianjin and Fujian, which implements the policy of FTZ from 2022 to 2025, and at the same time prediction results and the growth rate are compared and analysed, and corresponding policy recommendations are put forward.

4. Empirical analysis

This article selects the Tianjin and Fujian provinces that have established FTZs and have relatively unabridged index data as objects of study. Given that more and more factors influence port cargo throughput, grey correlation analysis is applied to ascertain port cargo throughput’s main influencing factors in two coastal provinces (cities). Specifically, the main influencing factors of Tianjin’s port cargo throughput are the secondary industry-added value, road freight volume and total import and export volume, among which the secondary industry-added value has the most significant influence, with a correlation value of 0.6798. The correlation between port throughput and total import and export volume, GDP and water freight volume is strong in Fujian province, and total import and export volume is most closely related to port cargo throughput. The data on port cargo throughput (measured in millions of tons) and the influencing factors are from statistical yearbooks published by the National Bureau of Statistics and local statistical bureaus from 2005 to 2021.

4.1 Case 1: simulation and prediction of Tianjin’s port cargo throughput

4.1.1 Model comparison and analysis of port cargo throughput in Tianjin

Tianjin FTZ was established in 2015. Benefiting from the FTZ policy dividend, Tianjin’s ports actively carry out offshore trade, which helps port development. Through grey correlation analysis, the secondary industry-added value, road freight volume, and total import and export volume are identified as the three main influencing factors of port cargo throughput in Tianjin. The data can be seen in Table 1.

The optimal order of the new model is −0.6195 after optimization by the GWO. According to the modelling steps, the data in Table 1 are modelled, and the calculation formula of the FDCGM(1,N) model can be obtained as follows:

xˆ1(0.6195)(k)=e0.6687*(k1)*x1(0.6195)(1)+10.6687*(1e0.6687*(k1))(7.3537x2(0.6195)(k)0.0315x3(0.6195)(k)3.5034x4(0.6195)(k)2018.7736d5(0.6195)(k))

From Table 2 and Figure 2, the port cargo throughput of Tianjin in 2005–2021 fluctuates relatively smoothly and shows a stable upward trend overall. Compared with the other three grey and two machine learning models, the new model best fits the original data. Especially for the simulations and forecasts of recent years, the simulated and forecasted values are approximate to raw values. Meanwhile, other comparison models have large deviations from the actual values. Figure 3 shows that the new model has the least APE fluctuation and is more stable, while the Backpropagation Neural Network (BPNN) and Long Short-Term Memory (LSTM) models will have a larger APE and poor stability at a certain time.

In terms of model analysis, DVCGM(1,N) model’s MAPEs and MAPEp are smaller than those of the GM(1,N) and DGM(1,N) model. It implies that FTZ policy does have an impact on Tianjin’s port cargo throughput. Meanwhile, the FDCGM(1,N) model’s MAPEs and MAPEp are smaller than those of the DVCGM(1,N) model. It means that the impact of port cargo throughput in previous years needs to be considered when forecasting port cargo throughput in Tianjin. Figure 4 shows that the FDCGM(1,N) model performs well in simulation and forecasting. Its simulation and in-sample prediction errors are the smallest, at 6.07% and 6.70%, respectively. These indicate that the new model effectively identifies serial trends of port cargo throughput in Tianjin. It may be due to the inclusion of policy dummy variables, and fractional-order AGO in the new model, significantly improving the time series’ estimated performance and stability.

4.1.2 Forecast Tianjin’s port cargo throughput in 2022–2025

The FDCGM(1,N) model performs well in the Tianjin simulation and forecasting port cargo throughput. Thus, using all the data from 2005 to 2021, the new model is used to predict Tianjin’s port cargo throughput in 2022–2025. As seen from Table 3, the port cargo throughput of Tianjin shows a stable growth trend from 2022 to 2025. In recent years, Tianjin has taken advantage of FTZ’s innovation and used its strategic port resources and hard-core benefits to refine the port. Tianjin’s port carrying capacity has been significantly improved; port cargo throughput and infrastructure capacity have been substantially improved; and the construction of intelligent, green and safe ports has been practical. Tianjin will further enhance the port collection and distribution network, strengthen the role of shipping hubs and accelerate the formation of Tianjin port as the centre of the northern international shipping hub. At the same time, the port must implement the “smart supervision” mode of operation, build an innovative, green world-class port, upgrade the port economy, better serve the Jing-Jin-Ji region and strengthen Tianjin FTZ’s functional positioning. In 2025, the port cargo throughput of Tianjin will reach 762,348,800 tons.

4.2 Case 2: simulation and prediction of Fujian province’s port cargo throughput

4.2.1 Model comparison and analysis of port cargo throughput in Fujian province

Fujian FTZ is based on promoting cross-strait economic and trade cooperation. The correlation between port cargo throughput and its influencing factors in Fujian province is obtained through grey correlation analysis, and then three main influencing factors are obtained, namely total import and export volume, GDP and water freight volume. Table 4 manifests the required data.

FDCGM(1,N) model’s optimal order is −0.9396, obtained by GWO. According to the modelling steps, the data in Table 4 are modelled, and the calculation formula of the FDCGM(1,N) model can be obtained as follows:

xˆ1(0.9396)(k)=e1.1345*(k1)*x1(0.9396)(1)+11.1345*(1e1.1345*(k1))(7.4598x2(0.9396)(k)+0.1488x3(0.9396)(k)+0.7167x4(0.9396)(k)778.0284d5(0.9396)(k))

Table 5 shows the raw data and six models’ simulation and forecast results. Fujian province’s port cargo throughput of achieved more than three times growth from 2005 to 2021. Figure 5 shows that the new model is not only closer to actual values in terms of simulation and forecasting results but also very consistent with the actual situation in terms of change trends, which all show an upward trend overall and can reflect the actual state of port throughput in Fujian province. Especially for the in-sample prediction stage, the predicted values of other models have a significant deviation from the true value.

Figure 6 exhibits that the APE of the FDCGM(1,N) model is within 5% except for a few years, and the degree of fluctuation is the most stable among all the comparison models. It indicates that new model has more durable and accurate performance and more robust adaptability than the comparison models. Meanwhile, Figure 7 shows that the MAPEs of the FDCGM(1,N) model is 4.16%, a little bit higher than that of the DGM(1,N) model. The LSTM model has the largest MAPEs of 14.62%. The FDCGM(1,N) model performs well with the smallest MAPEp among the six models at 1.63%. The new model considering the fractional-order AGO can overcome the limitations of other models in the simulation and forecasting process and has better forecasting accuracy and stable performance. It can perform out-of-sample forecasting of cargo throughput in Fujian province.

4.2.2 Forecast Fujian’s port cargo throughput in 2022–2025

A new model is used to forecast Fujian province’s port cargo throughput in 2022–2025. From Table 6, in the next four years, Fujian province’s port cargo throughput will continue to grow and reach 1,017,638,300 tons in 2025. The construction of Fujian FTZ provides an endogenous impetus for developing ports. With a series of preferential policies brought by the FTZ, Fujian provincial ports have been significantly strengthened, the disadvantages of “one bay and two ports” management have been effectively solved, and the level of port intensification, scale and modernization have been significantly improved. The free trade and convenient investment brought by Fujian FTZ have accelerated port development and enhanced port enterprises’ competitiveness. In the future, Fujian province will continue to promote port integration actively, clarify the scope of port areas, functional positioning and development goals of each port, strengthen cross-strait industrial docking and promote advantageous industrial cooperation.

4.3 Comprehensive contrastive analysis

Figures 8 and 9 show the growth rates of port cargo throughput from 2005 to 2025, and the port cargo throughput prediction results from 2022 to 2025 for Tianjin and Fujian provinces. From them, the effects of the FTZ policy on port cargo throughput differ from region to region.

Tianjin’s port cargo throughput has shown a growth trend after the constitution of the Tianjin FTZ, but the growth rate has slowed. It may be because Tianjin FTZ strives to develop its transit transport business, which diverts part of the import and export cargo flow. Meanwhile, the cargo throughput mainly comes from the Beijing-Tianjin-Hebei region. After establishing the FTZ, Tianjin port faces more and more competition from the Bohai Seaport group. The growth rate will be reduced as a result. Overall, the FTZ policy has brought stable development to Tianjin’s ports. In 2022–2025, the cargo throughput of Tianjin ports will continue to grow, and its growth rate will remain stable. It indicates that the Tianjin FTZ policy will continue to fuel the development of Tianjin’s ports.

Meanwhile, the Fujian FTZ policy has contributed significantly and continuously to the port cargo throughput. With the preferential policies of the Fujian FTZ, easier market access has attracted many enterprises to move in, driving import and export business. It has directly caused the rapid growth of port cargo throughput in Fujian province, thus maintaining a stable growth rate. The port cargo throughput continues its growth trend in 2022–2025, reaching 1,017,638,300 tons in 2025. Its growth rate increases and then decreases, indicating that Fujian needs to strengthen port infrastructure construction, improve port productivity and realize the boosting effect of the FTZ on port development in the long term.

5. Conclusions

The establishment of the FTZs has advanced port production and development. As an essential indicator of port production capacity and operation, accurate prediction of port cargo throughput is vital for policymakers in formulating port development plans. This paper forecasts and analyses the port cargo throughput under the FTZ policy by constructing a new model and obtains the following conclusions:

  1. By introducing the fractional-order AGO and considering the policy effect, this paper proposes a new grey multivariate forecasting model, which enhances the adaptability and forecasting performance of the new model and compensates for the shortcomings of the traditional multivariate model. Further, the GWO algorithm is employed to track the optimal introduced fractional-order AGO, which weakens the randomness of the original data series and improves the ability of the FDCGM(1,N) model to mine the information of the original data series.

  2. By simulating and forecasting port cargo throughput from 2005 to 2021 for Tianjin where FTZs are located, the results show that the new model curve can fit the raw data curve well. Compared with five comparison models, GM(1,N), DGM(1,N), DVCGM(1,N), BPNN and LSTM, the new model has good performance with MAPE values of 6.07% in the simulation phase and 6.70% in the forecasting phase, respectively. It indicates that the FDCGM(1,N) model is a practical tool for forecasting port cargo throughput. It can be used for out-of-sample forecasting and port cargo throughput of Tianjin will reach 762,348,800 tons in 2025. Meanwhile, comparative analysis reveals that Tianjin’s port cargo throughput has shown a slow growth trend under the influence of FTA policy.

  3. The port cargo throughput of Fujian province where FTZs are located is simulated from 2005 to 2021, and the results show that the new model curve can better reflect the changing trend of the original data. Compared with the five comparison models GM(1,N), DGM(1,N), DVCGM(1,N), BPNN and LSTM, the MAPE value of the new model is 4.16% in the simulation stage and 1.63% in the prediction stage. It shows that the FDCGM(1,N) model can be used to predict the port cargo throughput in Fujian province out of sample. Through the forecast, it is found that the cargo throughput of ports in Fujian province will reach 1,017,638,300 tons in 2025. In addition, the results of the comparative analysis show that, thanks to the FTZ policy, Fujian’s port cargo throughput has exhibited rapid growth.

The grey multivariate model proposed in this paper, which considers the effects of both policy and fractional order, has better simulation and forecasting performance than the traditional grey forecasting and machine learning models. In the future, other intelligent optimization algorithms can be considered to improve the new model’s performance. The autoregressive time-lag term can be also introduced to improves the new model’s performance. In addition, the FDCGM(1,N) model can also be used to predict the port cargo throughput of other provinces (cities) implementing FTZ policy, providing a reliable basis for decision-makers to make port development plans.

Figures

Modelling step diagram of the FDCGM(1,N) model

Figure 1

Modelling step diagram of the FDCGM(1,N) model

Fitting curves of five models to Tianjin’s port cargo throughput

Figure 2

Fitting curves of five models to Tianjin’s port cargo throughput

APE(%) of Tianjin’s port cargo throughput amongst five models

Figure 3

APE(%) of Tianjin’s port cargo throughput amongst five models

MAPE(%) of six models to Tianjin’s port cargo throughput

Figure 4

MAPE(%) of six models to Tianjin’s port cargo throughput

Fitting curves of five models to Fujian province’s port cargo throughput

Figure 5

Fitting curves of five models to Fujian province’s port cargo throughput

APE(%) of Tianjin’s port cargo Fujian province’s amongst six models

Figure 6

APE(%) of Tianjin’s port cargo Fujian province’s amongst six models

MAPE(%) of six models to Fujian province’s port cargo throughput

Figure 7

MAPE(%) of six models to Fujian province’s port cargo throughput

The growth rate of port cargo throughput of Tianjin and Fujian from 2006 to 2025

Figure 8

The growth rate of port cargo throughput of Tianjin and Fujian from 2006 to 2025

Port cargo throughput of Tianjin and Fujian from 2022 to 2025

Figure 9

Port cargo throughput of Tianjin and Fujian from 2022 to 2025

Observations of Tianjin’s port cargo throughput and the related factors from 2005 to 2021

YearPort cargo throughput (10000 tons)Secondary industry-added value (100 million yuan)Road freight volume (10000 tons)Total import and export volume (USD 100 million)D
200524069.001630.5319850.00533.870
200625760.001834.5420290.00645.730
200730946.002123.6323500.00715.500
200835593.002659.7118160.00805.390
200938111.002808.7419800.00639.440
201041325.003259.7420855.00822.010
201145338.003756.2623505.001033.910
201247697.004134.0327735.001156.230
201350063.004407.1028206.001285.280
201454002.004615.5031130.001339.120
201554051.004489.5930551.001143.471
201655056.004367.9732841.001026.511
201750056.004564.0634720.001129.451
201850774.004835.3034711.001225.111
201949220.004947.1831250.001066.451
202050290.004911.7732261.001059.311
202152954.005854.2734527.001325.651

Source(s): Table created by authors

Fitting and forecasting results of six models to Tianjin’s port cargo throughput

YearReal dataFDCGM(1,N)GM(1,N)DGM(1,N)
FittingAPE (%)FittingAPE (%)FittingAPE (%)
200524069.0024069.000.0024069.000.0024069.000.00
200625760.0030675.8719.0820491.4420.4530160.8217.08
200730946.0033332.437.7142782.8838.2534948.5912.93
200835593.0037530.325.4432569.688.4940914.0214.95
200938111.0038849.721.94−20409.11153.5543078.5213.03
201041325.0041735.300.99−72633.17275.7646508.4912.54
201145338.0045085.770.56−131536.77390.1251600.0313.81
201247697.0047705.360.02−210393.26541.1057274.7020.08
201350063.0049447.651.23−325222.87749.6366167.1132.17
201454002.0050853.865.83−482056.20992.6677471.5243.46
201554051.0047116.7612.83−962511.821880.7589805.5966.15
201655056.0045950.0016.54−1483889.662795.24104767.5990.29
201750056.0047188.055.73−2043942.044183.31126557.57152.83
201850774.0049415.612.68−3039098.036085.54157853.43210.89
201949220.0051414.894.46−5658615.4111596.58198566.41303.43
MAPE(%) 2122.25 71.69
202050290.0050813.971.04−8051467.5616110.08256814.57410.67
202152954.0059494.8812.35−12530458.3823762.91337972.18538.24
MAPE(%) 19936.49 474.45
YearReal dataDVCGM(1,N)BPNNLSTM
FittingAPE(%)FittingAPE(%)FittingAPE(%)
200524069.0024069.000.0027713.6415.1440860.3169.76
200625760.0022200.9313.8229163.7213.2138510.6149.50
200730946.0045838.9948.1333256.097.4639740.8728.42
200835593.0045361.1527.4435690.490.2740716.2414.39
200938111.0055502.6845.6338374.810.6941472.668.82
201041325.0057147.6938.2940777.121.3342813.843.60
201145338.0061507.7535.6643406.074.2644634.281.55
201247697.0074595.3456.3948446.281.5746526.282.45
201350063.0073482.4346.7851094.682.0648528.153.07
201454002.0084352.8556.2056204.624.0850402.476.67
201554051.0047996.1311.2044611.8617.4650791.546.03
201655056.0061230.6611.2242131.8123.4750789.087.75
201750056.0065221.3130.3035914.6228.2550490.430.87
201850774.0058988.1816.1835697.1329.6950080.161.37
201949220.0041848.4114.9852226.566.1149286.400.13
MAPE(%) 32.30 10.34 13.63
202050290.0046496.597.5449623.491.3348365.363.83
202152954.0039175.3326.0243590.6017.6847917.659.51
MAPE(%) 16.78 9.50 6.67

Source(s): Table created by authors

Tianjin's port cargo throughput predicted by the FDCGM(1,N) model from 2022 to 2025

2022202320242025
Port cargo throughput (10000 tons)64516.2868172.3672074.3576234.88

Source(s): Table created by authors

Observations of Fujian province’s port cargo throughput and the related factors in 2005–2021

YearPort cargo throughput (10000 tons)Total import and export volume (USD 100 million)GDP (100 million yuan)Water freight volume (10000 tons)D
200519809.25544.119210.006415.500
200623865.61626.5910841.007468.600
200723880.90744.5112130.009325.600
200827422.06848.2115193.0010931.800
200930831.81796.4914272.0012418.100
201033069.011087.8016803.0015002.500
201137695.951435.2218871.0017917.700
201241359.231559.3821100.0020190.700
201345911.191693.2223162.0022503.800
201449541.241774.0825782.0024942.100
201550652.091688.4628419.0026819.501
201651140.091568.1931664.0029609.401
201751995.491710.3533453.0033842.401
201856130.881875.7636854.0038687.801
201959483.991930.8642263.0042326.601
202062132.472033.1745018.0043608.601
202169190.282852.5050224.0048810.401

Source(s): Table created by authors

Fitting and forecasting results of six models to Fujian province’s port cargo throughput

YearReal dataFDCGM(1,N)GM(1,N)DGM(1,N)
FittingAPE(%)FittingAPE(%)FittingAPE(%)
200519809.2519809.250.0019809.250.0019809.250.00
200623865.6126309.4510.2419321.5019.0423143.583.03
200723880.9029560.2323.7829117.7821.9325146.885.30
200827422.0631762.8015.8334213.0124.7629705.888.33
200930831.8132104.634.1327472.8610.8928152.038.69
201033069.0135783.118.2133601.991.6131917.443.48
201137695.9539996.336.1039281.884.2137505.140.51
201241359.2342382.292.4742829.683.5641939.521.40
201345911.1944856.432.3046087.660.3845775.670.30
201449541.2447151.694.8249619.810.1649354.310.38
201550652.0947328.136.5651379.511.4451382.081.44
201651140.0948629.834.9152447.812.5652286.242.24
201751995.4952386.950.7551794.850.3952563.631.09
201856130.8856902.791.3854437.783.0254430.193.03
201959483.9960200.701.2062314.384.7659466.070.03
MAPE(%) 7.05 2.80
202062132.4761977.450.2568586.8210.3965107.764.79
202169190.2871271.523.0186089.3324.4279119.1114.35
MAPE(%) 17.41 9.57
YearReal dataDVCGM(1,N)BPNNLSTM
FittingAPE(%)FittingAPE(%)FittingAPE(%)
200519809.2519809.250.0024863.6025.5240870.00106.32
200623865.6119180.8619.6325521.596.9430766.2528.91
200723880.9029221.1622.3627081.1113.4028160.5617.92
200827422.0634453.9625.6428256.753.0430734.8612.08
200930831.8127817.369.7829134.125.5134638.1212.35
201033069.0133806.582.2333494.811.2937534.5413.50
201137695.9539311.854.2939021.023.5240799.688.23
201241359.2342915.813.7641955.431.4443458.645.08
201345911.1946221.490.6845083.911.8046087.970.39
201449541.2449887.450.7047922.433.2748717.441.66
201550652.0950835.200.3649129.093.0150991.790.67
201651140.0952307.972.2850113.412.0153065.703.77
201751995.4951808.850.3651992.740.0154630.185.07
201856130.8854663.142.6154233.133.3856350.720.39
201959483.9962858.315.6757019.614.1457737.662.94
MAPE(%) 7.17 5.22 14.62
202062132.4769172.2311.3357628.067.2558461.785.91
202169190.2886211.7824.6058242.1215.8260340.6312.79
MAPE(%) 17.97 11.54 9.35

Source(s): Table created by authors

Fujian province’s port cargo throughput predicted by the FDCGM(1,N) model from 2022 to 2025

2022202320242025
Port cargo throughput (10000 tons)78720.2885580.9293231.36101763.83

Source(s): Table created by authors

References

Chen, J., Wan, Z., Zhang, F., Park, N., Zheng, A. and Zhao, J. (2018), “Evaluation and comparison of the development performances of typical free trade port zones in China”, Transportation Research Part A: Policy and Practice, Vol. 118, pp. 506-526, doi: 10.1016/j.tra.2018.09.009.

Chen, J., Fei, Y., Wan, Z., Yang, Z., Li, H., Choi, K.-S. and Xie, X. (2020), “Allometric relationship and development potential comparison of ports in a regional cluster: a case study of ports in the Pearl River Delta in China”, Transport Policy, Vol. 85, pp. 80-90, doi: 10.1016/j.tranpol.2019.10.018.

Chen, W., Hu, Y., Liu, B., Wang, H. and Zheng, M. (2022), “Does the establishment of Pilot Free Trade Test Zones promote the transformation and upgradation of trade patterns?”, Economic Analysis and Policy, Vol. 76, pp. 114-128, doi: 10.1016/j.eap.2022.07.012.

Cuong, T.N., Kim, H.-S., You, S.-S. and Nguyen, D.A. (2022), “Seaport throughput forecasting and post COVID-19 recovery policy by using effective decision-making strategy: a case study of Vietnam ports”, Computers and Industrial Engineering, Vol. 168, 108102, doi: 10.1016/j.cie.2022.108102.

Deng, J.L. (1982), “Control problems of grey systems”, Systems and Control Letters, Vol. 1 No. 5, pp. 288-294, doi: 10.1016/S0167-6911(82)80025-X.

Ding, S., Dang, Y.G. and Xu, N. (2018), “Construction and application of GM(1,N) based on control of dummy variables”, Control and Decision, Vol. 33 No. 2, pp. 309-315, doi: 10.13195/j.kzyjc.2016.1613.

Ding, S., Tao, Z. and Hu, J. (2022), “Forecasting the economic indices of the high-tech industries in China using the grey multivariable convolution model”, Applied Soft Computing, Vol. 126, 109301, doi: 10.1016/j.asoc.2022.109301.

Du, P., Wang, J., Yang, W. and Niu, T. (2019), “Container throughput forecasting using a novel hybrid learning method with error correction strategy”, Knowledge-Based Systems, Vol. 182, 104853, doi: 10.1016/j.knosys.2019.07.024.

Du, P., Wang, J., Niu, T. and Yang, W. (2021), “PM2.5 prediction and related health effects and economic cost assessments in 2020 and 2021: case studies in Jing-Jin-Ji, China”, Knowledge-Based Systems, Vol. 233, 107487, doi: 10.1016/j.knosys.2021.107487.

Duan, H. and Liu, Y. (2021), “Research on a grey prediction model based on energy prices and its applications”, Computers and Industrial Engineering, Vol. 162, 107729, doi: 10.1016/j.cie.2021.107729.

Duan, H. and Luo, X. (2022), “A novel multivariable grey prediction model and its application in forecasting coal consumption”, ISA Transactions, Vol. 120, pp. 110-127, doi: 10.1016/j.isatra.2021.03.024.

Duan, H., Liu, Y. and Wang, G. (2022), “A novel dynamic time-delay grey model of energy prices and its application in crude oil price forecasting”, Energy, Vol. 251, 123968, doi: 10.1016/j.energy.2022.123968.

Fan, G., Xie, X., Chen, J., Wan, Z., Yu, M. and Shi, J. (2022), “Has China's Free Trade Zone policy expedited port production and development?”, Marine Policy, Vol. 137, 104951, doi: 10.1016/j.marpol.2021.104951.

Gonzalez Aregall, M., Bergqvist, R. and Monios, J. (2018), “A global review of the hinterland dimension of green port strategies”, Transportation Research Part D: Transport and Environment, Vol. 59, pp. 23-34, doi: 10.1016/j.trd.2017.12.013.

Hsu, C.-I., Huang, Y.-C. and Wong, K.I. (2020), “A Grey hybrid model with industry share for the forecasting of cargo volumes and dynamic industrial changes”, Transportation Letters, Vol. 12 No. 1, pp. 25-36, doi: 10.1080/19427867.2018.1502505.

Huang, D., Grifoll, M., Sanchez-Espigares, J.A., Zheng, P. and Feng, H. (2022), “Hybrid approaches for container traffic forecasting in the context of anomalous events: the case of the Yangtze River Delta region in the COVID-19 pandemic”, Transport Policy, Vol. 128, pp. 1-12, doi: 10.1016/j.tranpol.2022.08.019.

Jin, J., Ma, M., Jin, H., Cui, T. and Bai, R. (2023), “Container terminal daily gate in and gate out forecasting using machine learning methods”, Transport Policy, Vol. 132, pp. 163-174, doi: 10.1016/j.tranpol.2022.11.010.

Li, S., Liu, J. and Kong, Y. (2021), “Pilot free trade zones and Chinese port-listed companies performance: an empirical research based on quasi-natural experiment”, Transport Policy, Vol. 111, pp. 125-137, doi: 10.1016/j.tranpol.2021.07.022.

Liu, L. and Wu, L. (2021), “Forecasting the renewable energy consumption of the European countries by an adjacent non-homogeneous grey model”, Applied Mathematical Modelling, Vol. 89, pp. 1932-1948, doi: 10.1016/j.apm.2020.08.080.

Liu, J., Wang, X. and Guo, J. (2021), “Port efficiency and its influencing factors in the context of pilot free trade zones”, Transport Policy, Vol. 105, pp. 67-79, doi: 10.1016/j.tranpol.2021.02.011.

Liu, L., Liu, S., Wu, L., Zhu, J. and Shang, G. (2022), “Forecasting the development trend of new energy vehicles in China by an optimized fractional discrete grey power model”, Journal of Cleaner Production, Vol. 372, 133708, doi: 10.1016/j.jclepro.2022.133708.

Mirjalili, S., Mirjalili, S.M. and Lewis, A. (2014), “Grey wolf optimizer”, Advances in Engineering Software, Vol. 69, pp. 46-61, doi: 10.1016/j.advengsoft.2013.12.007.

Mo, L., Xie, L., Jiang, X., Teng, G., Xu, L. and Xiao, J. (2018), “GMDH-based hybrid model for container throughput forecasting: selective combination forecasting in nonlinear subseries”, Applied Soft Computing, Vol. 62, pp. 478-490, doi: 10.1016/j.asoc.2017.10.033.

Moonchai, S. and Chutsagulprom, N. (2020), “Short-term forecasting of renewable energy consumption: augmentation of a modified grey model with a Kalman filter”, Applied Soft Computing, Vol. 87, 105994, doi: 10.1016/j.asoc.2019.105994.

Ren, Y., Xia, L. and Wang, Y. (2023), “Forecasting China's hydropower generation using a novel seasonal optimized multivariate grey model”, Technological Forecasting and Social Change, Vol. 194, 122677, doi: 10.1016/j.techfore.2023.122677.

Sanguri, K., Shankar, S., Punia, S. and Patra, S. (2022), “Hierarchical container throughput forecasting: the value of coherent forecasts in the management of ports operations”, Computers and Industrial Engineering, Vol. 173, 108651, doi: 10.1016/j.cie.2022.108651.

Saxena, A. (2021), “Grey forecasting models based on internal optimization for Novel Corona virus (COVID-19)”, Applied Soft Computing, Vol. 111, 107735, doi: 10.1016/j.asoc.2021.107735.

Shi, K. and Wu, L. (2020), “Forecasting air quality considering the socio-economic development in Xingtai”, Sustainable Cities and Society, Vol. 61, 102337, doi: 10.1016/j.scs.2020.102337.

Wan, Z., Zhang, Y., Wang, X. and Chen, J. (2014), “Policy and politics behind shanghai's free trade zone program”, Journal of Transport Geography, Vol. 34, pp. 1-6, doi: 10.1016/j.jtrangeo.2013.10.015.

Wan, C., Zhao, Y., Zhang, D. and Yip, T.L. (2021), “Identifying important ports in maritime container shipping networks along the Maritime Silk Road”, Ocean and Coastal Management, Vol. 211, 105738, doi: 10.1016/j.ocecoaman.2021.105738.

Wan, G., Li, X., Yin, K. and Zhao, Y. (2022), “Forecasting carbon emissions from energy consumption in Guangdong Province, China with a novel grey multivariate model”, Environmental Science and Pollution Research, Vol. 29 No. 39, pp. 59534-59546, doi: 10.1007/s11356-022-19805-8.

Wang, J. and Li, N. (2020), “Influencing factors and future trends of natural gas demand in the eastern, central and western areas of China based on the grey model”, Natural Gas Industry B, Vol. 7 No. 5, pp. 473-483, doi: 10.1016/j.ngib.2020.09.005.

Wang, L., Liu, Z. and Shi, H. (2022a), “The impact of the pilot free trade zone on regional financial development”, China and World Economy, Vol. 30 No. 5, pp. 154-184, doi: 10.1111/cwe.12441.

Wang, Y., Wang, L., Ye, L., Ma, X., Wu, W., Yang, Z., He, X., Zhang, L., Zhang, Y., Zhou, Y. and Luo, Y. (2022b), “A novel self-adaptive fractional multivariable grey model and its application in forecasting energy production and conversion of China”, Engineering Applications of Artificial Intelligence, Vol. 115, 105319, doi: 10.1016/j.engappai.2022.105319.

Wang, Y., Yang, Z., Wang, L., Ma, X., Wu, W., Ye, L., Zhou, Y. and Luo, Y. (2022c), “Forecasting China’s energy production and consumption based on a novel structural adaptive Caputo fractional grey prediction model”, Energy, Vol. 259, 124935, doi: 10.1016/j.energy.2022.124935.

Wiegmans, B., Witte, P. and Spit, T. (2015), “Characteristics of European inland ports: a statistical analysis of inland waterway port development in Dutch municipalities”, Transportation Research Part A: Policy and Practice, Vol. 78, pp. 566-577, doi: 10.1016/j.tra.2015.07.004.

Wu, L.Z., Li, S.H., Huang, R.Q. and Xu, Q. (2020a), “A new grey prediction model and its application to predicting landslide displacement”, Applied Soft Computing, Vol. 95, 106543, doi: 10.1016/j.asoc.2020.106543.

Wu, W., Ma, X., Wang, Y., Cai, W. and Zeng, B. (2020b), “Predicting China’s energy consumption using a novel grey Riccati model”, Applied Soft Computing, Vol. 95, 106555, doi: 10.1016/j.asoc.2020.106555.

Yan, S., Su, Q., Gong, Z. and Zeng, X. (2022), “Fractional order time-delay multivariable discrete grey model for short-term online public opinion prediction”, Expert Systems with Applications, Vol. 197, 116691, doi: 10.1016/j.eswa.2022.116691.

Yang, Y. and Wang, X. (2022), “A novel modified conformable fractional grey time-delay model for power generation prediction”, Chaos, Solitons and Fractals, Vol. 158, 112004, doi: 10.1016/j.chaos.2022.112004.

Yao, D. and Whalley, J. (2016), “The China (shanghai) pilot free trade zone: background, developments and preliminary assessment of initial impacts”, The World Economy, Vol. 39 No. 1, pp. 2-15, doi: 10.1111/twec.12364.

Zeng, B., Luo, C., Liu, S., Bai, Y. and Li, C. (2016), “Development of an optimization method for the GM(1,N) model”, Engineering Applications of Artificial Intelligence, Vol. 55, pp. 353-362, doi: 10.1016/j.engappai.2016.08.007.

Zeng, B., Duan, H. and Zhou, Y. (2019), “A new multivariable grey prediction model with structure compatibility”, Applied Mathematical Modelling, Vol. 75, pp. 385-397, doi: 10.1016/j.apm.2019.05.044.

Zeng, B., Li, H. and Ma, X. (2020a), “A novel multi-variable grey forecasting model and its application in forecasting the grain production in China”, Computers and Industrial Engineering, Vol. 150, 106915, doi: 10.1016/j.cie.2020.106915.

Zeng, B., Ma, X. and Zhou, M. (2020b), “A new-structure grey Verhulst model for China's tight gas production forecasting”, Applied Soft Computing, Vol. 96, 106600, doi: 10.1016/j.asoc.2020.106600.

Zeng, B., Tong, M. and Ma, X. (2020c), “A new-structure grey Verhulst model: development and performance comparison”, Applied Mathematical Modelling, Vol. 81, pp. 522-537, doi: 10.1016/j.apm.2020.01.014.

Zeng, B., Zhou, W. and Zhou, M. (2021), “Forecasting the concentration of sulfur dioxide in Beijing using a novel grey interval model with oscillation sequence”, Journal of Cleaner Production, Vol. 311, 127500, doi: 10.1016/j.jclepro.2021.127500.

Zeng, B., Li, H., Mao, C. and Wu, Y. (2023a), “Modeling, prediction and analysis of new energy vehicle sales in China using a variable-structure grey model”, Expert Systems with Applications, Vol. 213, 118879, doi: 10.1016/j.eswa.2022.118879.

Zeng, B., Yin, F., Yang, Y., Wu, Y. and Mao, C. (2023b), “Application of the novel-structured multivariable grey model with various orders to forecast the bending strength of concrete”, Chaos, Solitons and Fractals, Vol. 168, 113200, doi: 10.1016/j.chaos.2023.113200.

Zhang, M., Guo, H., Sun, M., Liu, S. and Forrest, J. (2022), “A novel flexible grey multivariable model and its application in forecasting energy consumption in China”, Energy, Vol. 239, 122441, doi: 10.1016/j.energy.2021.122441.

Corresponding author

Zhiguo Zhao can be contacted at: zzhiguo2021@stu.ouc.edu.cn

Related articles