Next Article in Journal
Endotoxin as a Marker for Water Quality
Next Article in Special Issue
Modeling and Simulation of Crowd Pre-Evacuation Decision-Making in Complex Traffic Environments
Previous Article in Journal
Acculturative Stress, Perceived Social Support, and Mental Health: The Mediating Effect of Negative Emotions Associated with Discrimination
Previous Article in Special Issue
CEEMDAN-IPSO-LSTM: A Novel Model for Short-Term Passenger Flow Prediction in Urban Rail Transit Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Environment Estimation Method of Whole Sample Traffic Flows and Emissions Based on Multifactor MFD

1
Beijing Advanced Innovation Center for Future Urban Design, Beijing University of Civil Engineering and Architecture, Beijing 100044, China
2
Key Laboratory of Transport Industry of Big Data Application Technologies for Comprehensive Transport, Beijing Jiaotong University, Beijing 100044, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(24), 16524; https://doi.org/10.3390/ijerph192416524
Submission received: 1 November 2022 / Revised: 27 November 2022 / Accepted: 7 December 2022 / Published: 9 December 2022

Abstract

:
Vehicle emissions seriously affect the air environment and public health. The dynamic estimation method of vehicle emissions changing over time on the road network has always been the bottleneck of air quality simulation. The dynamic traffic volume is one of the important parameters to estimate vehicle emission, which is difficult to obtain effectively. A novel estimation method of whole sample traffic volumes and emissions on the entire road network based on multifactor Macroscopic Fundamental Diagram (MFD) is proposed in this paper. First, the intelligent clustering and recognition methods of traffic flow patterns are constructed based on neural network and deep-learning algorithms. Then, multifactor MFD models are developed considering different road types, traffic flow patterns and weekday peak hours. Finally, the high spatiotemporal resolution estimation method of whole sample traffic volumes and emissions are constructed based on MFD models. The results show that traffic flow patterns are clustered efficiently by the Self-Organizing Maps (SOM) algorithm combined with the direct time-varying speed index, which describe 91.7% traffic flow states of urban roads. The Deep Belief Network (DBN) algorithm precisely recognizes 92.1% of the traffic patterns based on the speeds of peak hours. Multifactor MFD models estimate the whole sample traffic volumes with a high accuracy of 91.6%. The case study shows that the vehicle emissions are evaluated dynamically based on the novel estimation method proposed in this paper, which is conducive to the coordinated treatment of air pollution.

1. Introduction

Vehicle exhaust pollution is an important factor threatening public health [1]. Dynamic estimation of mobile pollution sources changing over time on the road network is the fundamental data for air pollution control and public health promotion [2]. Pollution sources are required to be controlled in coordination in air quality management [3]. Fixed pollution sources such as steel plants [4], power plants [5] and meteorological parameters [6] have been detected with high space-time resolution. However, the emissions of mobile vehicle sources still remain at the level of annual total amount measurement due to the limitation of emission detection technology [7]. The time and space granularity of vehicle emissions evaluation makes it difficult to measure it consistently in comparison with other pollution sources.
Traffic volume covering every link of the entire road network is one of the important parameters for dynamic emission estimation, which is difficult to collect based on the existing traffic detection with low coverage and is not estimated effectively [8]. The short-term traffic flow prediction method based on machine learning, which depends on huge historical traffic flow data, has been well studied [9]. Historical flow data of high-grade roads are sufficient, but the flow data of low-grade roads are extremely lacking [10]. Traffic flow characteristics extracted from historical laws based on machine learning reflect with difficulty the real-time traffic state on the road network, especially when accidental traffic events occur.
Although the whole sample traffic volumes of the entire road network are difficult to collect, the detection technology of speed data from the large-scale road network is well developed [11,12]. The traffic fundamental diagram model describes the relationship between speed and volume [13]. It is an effective method to estimate volume by the traffic fundamental diagram model based on speed data. The prediction accuracy for traffic volume based on the traffic fundamental diagram is improved with the increase of time granularity, which is relatively low when the time granularity is at the minute level and is relatively high when the time granularity is at the hour level [14]. This method is not widely used for short-term traffic flow prediction. The hourly granularity volume has met the demand for emission data. The traffic fundamental diagram is used to estimate the dynamic real-time flow based on the timely speed data in this paper.
The influencing factors of MFD are analyzed quantitatively, including road types, weather conditions, traffic patterns and travel periods. A novel estimation method of whole sample traffic flows and emissions based on multifactor MFD is proposed. The research results of this paper can be implemented for the preparation of a dynamic vehicle emission inventory on the road network, which is helpful for transportation and environmental departments to formulate policies of energy conservation and emission reduction. It plays a major role in promoting public health.

2. Literature Review

The quantitative estimation of dynamic emissions is fundamental for pollution control. Mobile sources have become an important part of air pollution, accounting for 55.6% of the total nitrogen oxide emissions in China in 2020 [15]. The Ministry of Ecology and Environment of China plans to evaluate emission sources with a uniform time granularity. Industrial, agricultural and living sources are fixed sources, which make it easy to obtain continuous dynamic data due to fewer pollution sources and simple monitoring technology [4,5,6]. Mobile sources have complex detection technologies and high testing costs, which remain at the level of static annual total quantity evaluation [16]. The evaluation granularity of vehicle emissions cannot be consistent with other pollution sources in terms of time and space, which limits the real-time detection and overall evaluation of air pollution. Therefore, dynamic estimation of vehicle emissions has always been a bottleneck problem in dynamic air quality simulation.
One of the important reasons why it is difficult to estimate dynamic vehicle emissions is the lack of an effective method to estimate whole sample dynamic traffic volumes of urban road networks. Vehicle emissions can be calculated based on single vehicle emission factors and vehicle kilometers traveled (VKT) [14]. The single vehicle emission factor is the emissions per kilometer of a single vehicle, which has a good research foundation and can be effectively obtained [17]. The whole sample dynamic traffic volume is the traffic volume that dynamically changes over time on every road link of a road network, which is an important parameter to calculate the dynamic VKT [16]. The main sources of traffic volume in existing studies include field traffic surveys, traffic simulations and traffic flow models. Field traffic surveys include fixed detector detection and manual surveys, which are costly and difficult. The road network coverage of field traffic survey data is limited, and it is difficult to obtain large-scale real-time dynamic data [18]. The traffic simulation method [19] has difficulty ensuring the consistency of the simulation environment and the actual traffic conditions. The traffic flow model method uses speed data to estimate volume based on traffic fundamental diagrams [20]. Dynamic speed data of large-scale road networks can be easily obtained based on massive data, such as floating car data (FCD) [11] and travel trajectories collected from mobile applications [21], which have high temporal and spatial continuity. Dynamic traffic control, such as signal timing optimization, requires small time granularity (five minutes) and high estimation accuracy of traffic volume [22]. The estimation accuracy of small time granularity traffic volume calculated using a fundamental diagram model cannot meet the requirement of the traffic control field; therefore, a fundamental diagram is not widely used in the traffic control field [16]. The hourly time granularity volume meets the emission estimation requirement and can be estimated with high accuracy based on a traffic flow fundamental diagram [14]. This paper proposes a method to calculate whole sample dynamic volumes of urban road networks based on traffic flow fundamental diagrams.
The traffic fundamental diagram was originally proposed by Greenshields in the form of a linear model in the 1930s [23]. Then, the logarithmic, exponential and vehicle-following models were separately introduced by Greenberg [24], Underwood [25] and Pipes [26]. In 1995, the Van Aerde model was proposed [27] with four parameters, a single structure, continuous feature and flexible calibration. It was found that the Van Aerde model fits the traffic flow of highways and expressways well [28]. The Underwood model has good performance in describing major arterial road and minor arterial road traffic flows [20]. However, the data used to develop these models are mostly collected from a single observation point. With the development of urban road networks and increasing traffic flow complexity, traffic fundamental diagrams have gradually expanded from single-section to multi-section networks. The concept of a macroscopic fundamental diagram (MFD) was proposed [29], and the existence of MFDs was verified based on field survey data and simulation data [30,31]. Some scholars [32,33] proposed new fundamental traffic diagrams based on a MFD. Existing research on MFDs is mainly based on simple structure road networks, excellent traffic conditions or simulation methods. MFD models based on real traffic conditions and complex road networks need to be further studied.
Traffic flow change patterns are divided to improve the effectiveness and accuracy of traffic volume prediction. It was pointed out that the travel behaviors of residents are cyclical and repetitive [18], which led to multiple but limited patterns of traffic flow patterns. Representative traffic flow patterns were constructed based on traffic volume and speed data in existing studies as follows. Jia et al. [34] proposed a spatiotemporal neural network model to predict traffic flow for each road segment. The traffic flow was divided into recent, daily and weekly parts. A multimode traffic flow prediction method with clustering-based attention convolution LSTM was proposed by Huang et al. [35] to model spatial-temporal data of traffic flow by integrating weather, wind speed, holidays and other factors to improve the prediction accuracy. Ma et al. [36] proposed a dynamic time warping method to select the appropriate historical data for daily traffic flow forecasting based on daily traffic flow pattern influence factors, such as season, day of the week, weather and holiday. The prediction effectiveness and accuracy of traffic volume are improved by dividing traffic flow patterns. Field survey traffic flows on highways or expressways are used as clustering indexes in most studies, which do not take into account the influence of the geographic location, road type and land use function on the division of traffic flow patterns. It is difficult to discover the differences among various traffic flow change patterns. A traffic flow clustering method considering various factors is proposed based on various road types of traffic flow data in this paper.
Above all, traffic flows can be clustered into limited typical patterns due to the regularity of travel and the grading characteristics of roads. Predicting flows for different patterns separately is an effective method to improve the prediction accuracy. Although whole sample dynamic traffic volumes of urban road networks are difficult to obtain, the large-scale speed data of a road network can be obtained based on a floating car system. The traffic flow model is used in this paper to predict traffic volume based on speed data for different patterns.

3. Methodology

3.1. Research Steps

The general methodology of this study is divided into the following steps, as shown in Figure 1.
First, the data preparation of traffic speed and volume for developing traffic flow patterns and traffic fundamental diagram models is introduced. The speed and volume data are aggregated into a 1-h granularity. Then, the time and location information of different data sources are matched.
Second, a clustering method of traffic flow patterns based on multidimensional features is constructed. The optimal traffic flow pattern clustering method is designed using different clustering methods. A library of typical traffic flow patterns under various traffic conditions is constructed.
Third, a traffic flow pattern recognition method based on speed indexes in different periods is constructed. Indexes suitable for different traffic application scenarios are designed to recognize traffic flow patterns rapidly.
Finally, a method is proposed to estimate whole sample dynamic traffic volumes of urban road networks based on clustering and recognition of traffic flow patterns. Multifactor MFD models are constructed based on road types, traffic flow patterns and morning and evening peak hours on weekdays. The dynamic traffic volume is estimated based on real-time speed data. The dynamic emissions of a road network before and after the traffic restriction policy are calculated as a case study.

3.2. Data Source and Preparation

The data used in this study include three parts.
(1) The average spatial speed data are from FCD in 5-min intervals of various road types in Beijing. The number of days for highway, expressway, major arterial road and minor arterial road data were 164, 175, 265 and 213, respectively, in 2018. Data on weekdays, weekends and holidays under different weather conditions are included. The attributes of average spatial speed data include link ID, road name, date, time, road type and average spatial speed.
(2) The traffic volume data from remote traffic microwave sensors (RTMSs) in 2-min intervals of various road types in Beijing are used. The number of days for highway, expressway, major arterial road and minor arterial road data were 132, 144, 205 and 153, respectively, in 2018. The attributes of average spatial speed data include detector number, road name, date, time and traffic volume.
(3) The traffic volume from traffic survey data at 5-min intervals under various road types in Beijing is used. The RTMS data on highways and expressways are relatively abundant, whereas the data on major arterial roads and minor arterial roads are relatively lacking. The traffic volume survey data are used to supplement the RTMS volume data. The attributes of traffic survey data include survey date, road name, vehicle type and traffic volume.
After data quality control, the speed data are integrated according to Equation (1).
V j = i n v i n
where V j is the speed of hour j in 1-h intervals, vi is the original (60/n) min speed in each 1-h interval and n is the number of the original (60/n) min data in each 1-h interval. The traffic volume data are integrated according to Equation (2).
Q j = i n q i
where is the traffic volume of hour j in 1-h intervals and q i is the original (60/n) min traffic volume in each 1-h interval.
Then, the speed and traffic data are normalized using Equation (3) to avoid the effect of the quantity of the data.
x i * = x i x m i n x m a x x m i n
where x i is the speed (km/h) or traffic volume (pcu/h/lane), x i * is the normalized speed or traffic volume and x m i n and x m a x are the minimum and maximum values of all records in a day, respectively.

4. Speed Clustering Method

The pattern clustering method of speeds includes the following processes: (1) clustering algorithm design; (2) clustering index selection; and (3) clustering effect evaluation.

4.1. Traffic Speed Pattern Definition

Affected by differences in the road function and travel demand, speeds on different roads and dates show different change patterns. The travel purpose on weekdays is mostly commuting, and the travel time is relatively fixed. The speed presents a bimodal distribution of morning and evening peaks with fixed time. The commuting demand decreases, whereas the entertainment demand increases on weekends and holidays. The traffic speed changing patterns of weekends and holidays differ from those of weekdays. The representative form of the traffic flow change curves over time is defined as the traffic flow pattern in this paper.

4.2. Determination of the Optimal Number of Clusters

A method to automatically determine the optimal number of clusters is proposed. The silhouette coefficient (SC) is adopted to determine the number of categories. The silhouette coefficients of different clustering number of speed curves are automatically calculated based on the K-means clustering algorithm. The optimal clustering number is determined when the profile coefficient reaches the maximum.
The SC is used to evaluate the tightness of samples within the categories and the separation of samples between categories, as shown in Equation (4).
S i l ( a ) = a a b a max { a a , b a } = min K K ( a ) { mean A c K | A a A c | } mean A b K ( a ) | A a A b | max { min K K ( a ) { mean A c K | A a A c | } , mean A b K ( a ) | A a A b | }
where Sil(a) is the silhouette of Aa, aa is the minimum value of the average distance between sample Aa and the samples in other categories, ba is the average distance between sample Aa and other samples in the same category, Ab is a sample in the same category of Aa, Ac is a sample of different categories from Aa, K(a) is the category of Aa and K’ is a different category from K(a).

4.3. Clustering Algorithms

The hierarchical clustering method, K-means clustering method and SOM neural network clustering method are used to classify the speed time-varying pattern, which proves to be effective in clustering traffic flows.

4.3.1. Hierarchical Clustering Method

The solution of hierarchical clustering method is carried out hierarchically, which is simple and intuitive. The pseudo code of hierarchical clustering algorithm is described as follows [37].
Step 1: Enter N cluster samples and the number of expected clusters K. Each sample is regarded as a cluster, and the distance di j between every two clusters is calculated to obtain the initial clustering matrix, as shown in Equation (5).
d i j = k = 1 p ( x i k x j k ) 2
where, di j is the Euclidean distance between sample point i and sample point j; p is the dimensions of clustering index; x i k and x j k are the coordinates of point i and point j respectively.
Step 2: Merge the two clusters with the smallest di j into a new cluster.
Step 3: Recalculate the distance dij between the new cluster and other clusters;
Step 4: Repeat Step 2 and Step 3 until the number of clusters becomes K.

4.3.2. K-Means Clustering Method

K-means clustering method converges fast and has a better clustering effect when the sample gap between classes is obvious. The pseudo code of K-means clustering algorithm is described as follows [38].
Step 1: Enter N cluster samples and the number of expected clusters K. K samples are randomly selected from all samples as the initial cluster center.
Step 2: Calculate the Euclidean distance from each sample to every cluster center and classify it with the nearest center.
Step 3: Recalculate the average coordinate of each cluster as the new cluster center.
Step 4: Calculate the square error E(t) of all sample points and compare it with the previous error E(t − 1).
Step 5: If E(t) − E(t − 1) < 0, turn to step 2, otherwise the algorithm ends.

4.3.3. SOM Clustering Method

The SOM neural network was proposed by Professor Kohonen [39]. The SOM network includes an input layer and output layer, which are connected by weight vectors. The input layer is a high-dimensional vector, and the output layer is a series of ordered nodes organized on a two-dimensional grid. The clustering results of SOM are highly visual and interpretable. SOM clustering requires the following steps [40]:
Step 1: Initialize the network and assign the initial value to each node weight of the output layer.
Step 2: Randomly select the input vector from the input samples to find the weight vector with the minimum distance from the input vector.
Step 3: Define the winning unit and adjust the weight in the adjacent area of the winning unit to make it close to the input vector.
Step 4: Provide new samples and conduct training.
Step 5: Shrink the neighborhood radius, reduce the learning rate, repeat until it is zero and output the clustering results.

4.4. Clustering Index Design

Clustering index refers to the feature vector that is related to the clustering analysis target, reflects the characteristics of the classification object and has obvious differences in different objects. Three indexes are designed in this study: (1) direct index (DI) based on time-varying speed differences, (2) comprehensive index (CI) based on speed statistical characteristics and curve geometric characteristics, and (3) simplified index (SI) based on speed statistical characteristics and curve geometric characteristics.

4.4.1. Direct Index Based on Time-Varying Speed Differences

The speeds of 0:00–24:00 at 1-h granularity are used as a clustering index to describe the difference in the traffic flow at different times directly. The speed of different time periods has different abilities to describe the characteristics of traffic flow. The speeds of peak hours can reflect the difference in traffic flow better than speeds of early morning. The coefficient of variation of each hour speed is used as the weight in this study. A larger weight is given to speeds with more different characteristics to improve the discrimination of traffic flows on different days.
The coefficient of variation calculation method is shown in Equation (6):
C V j = σ j ( x ) E j ( x ) = x i j / m [ x i j E j ( x ) ] 2 / m
where C V j is the coefficient of variation of the initial index of j, σ j ( x ) is the population standard deviation of the initial index of all samples j to be clustered, Ej(x) is the mean value of j initial index values of all samples to be clustered, xij is the initial index value of j in sample i and mj is the number of samples to be clustered.
The calculation method of the clustering index is shown in Equation (7):
a i j = x i j × C V j
where a i j is the final clustering index value of j in sample i. The final input cluster sample index matrix is shown in Equation (8):
[ a 11 a 12 a 1 n a 21 a 22 a 2 n a m 1 a m 2 a m n ]
where n is the number of indexes. In the DI, n is 24.

4.4.2. Comprehensive Index Based on Speed Statistical Characteristics and Curve Geometric Characteristics

The CI includes statistical characteristics of speeds and geometric characteristics of the speed time-varying curve. Statistical characteristics of speeds include the maximum speed, the minimum speed, the ratio of the maximum speed and the minimum speed, the average speed of different time periods, such as 0:00 to 24:00, 0:00 to 6:00, 20:00 to 23:00, 10:00 to 12:00, 2:00 to 4:00, morning peak hours, evening peak hours, the ratio of the morning peak hours speeds and the evening peak hours speeds and the variance of speeds of 0:00–24:00. The time-varying speed curve is divided into concave and convex shapes: (1) concave shape from 6:00 to 12:00; (2) concave shape from 12:00 to 22:00; and (3) convex shape from 10:00 to 14:00, as shown in Figure 2. The geometric characteristics of the speed time-varying curve include the kurtosis and skewness of each concave and convex shape. The comprehensive indexes are shown in Table 1.

4.4.3. Simplified Index Based on Speed Statistical Characteristics and Curve Geometric Characteristics

There are problems in the CI, such as the excessive number of indexes, correlation between indexes and different capabilities to describe traffic flow characteristics. The SI is constructed using the principal component analysis method to extract the most representative key indexes from the CI, which can effectively classify traffic flow, simplify the feature matrix and improve the calculation efficiency. The SI are shown in Table 2.

4.5. Clustering Effect Evaluation Indexes

The clustering effect evaluation indexes used in this paper are the silhouette coefficient (SC), consistency within a category, differences between categories and the coverage ratio (CR).

4.5.1. Consistency within a Category

The consistency within a category is the coefficient of variation (CV) of samples within a category, which is used to measure the consistency of speeds within a category, as shown in Equation (9).
C V in = 1 n k = 1 n C V in k = 1 n k = 1 n ( 1 l i = 1 l C V i , in k )
where C V in is the average CV within a category, C V in k is the CV within category k, n is the number of clusters, C V i , in k is the CV within a category at time i in category k and l is the number of time intervals in a day.

4.5.2. Difference between Categories

The difference between categories is the coefficient of variation of samples between categories, which is used to measure the differences between categories, as shown in Equation (10).
C V out = 1 l i = 1 l C V i , out
where C V out is the CV of the speed of all categories and C V i , out iis the CV s between categories at time i.

4.5.3. Coverage Ratio

The coverage ratio is the proportion of reasonably classified samples to the total samples in a category, which reflects the rationality of the clustering results and the representativeness of the traffic flow pattern. A reasonably categorized sample is defined as the average Pearson correlation coefficient between the sample and other samples within a category greater than 0.95. The coverage ratio is calculated by Equation (11).
C R = D 1 D × 100 %
where D1 is the number of correctly classified days and D is the total number of days.

4.5.4. Clustering Results Evaluation

The silhouette coefficients, consistencies within a category, differences between categories and coverage ratios of traffic flow patterns under different indexes are shown in Figure 3.
As shown in Figure 3, the categories of the DI have the highest C V out and C R and have higher SC. The differences between samples of different categories are high, and the proportion of samples reasonably classified reaches 91.7%. The categories of the CI have higher CR, but lower SC and C V out . The difference between samples of different categories is lower than DI and SI. The categories of the SI have higher SC and C V out . The difference between samples of different categories is great. However, the categories of the SI have the lowest CR. The proportion of samples reasonably classified is lower than DI and CI. Above all, DI has the best clustering effect, which has the strongest ability to classify the samples reasonably and distinguish the differences between different categories.
The categories of DI based on different algorithms are shown in Figure 4.
The hierarchical clustering algorithm based on DI clusters speed curves into ordinary weekdays and weekends, rainy days, non-congestion weekdays, days in lieu, The Dragon Boat Festival, festivals, prominent evening peak and prominent evening peak patterns. There is no distinction between weekdays and weekends, and the samples of some patterns are too few to be practical based on hierarchical clustering algorithm and the DI method.
The K-means clustering algorithm based on DI clusters speed curves into Mondays to Thursdays, Fridays, weekends, rainy days, festivals, weekends and festivals, prominent evening peak and prominent evening peak patterns. The categories of weekends and festivals are not clear. The common features of the samples in some categories are not obvious.
The SOM clustering algorithm based on DI clusters speed curves into Mondays to Thursdays, Fridays, Saturdays, Sundays, rainy days, festivals, prominent evening peak and prominent morning peak patterns. The traffic flow patterns are reasonable and practical, which have the following characteristics: (1) the differences between categories are the highest; (2) the traffic patterns are more detailed, and factors such as weather conditions are taken into account to better reflect the characteristics of different traffic flow patterns.
Above all, the SOM neural network clustering algorithm based on the DI is the most effective clustering method, which is used to cluster the traffic flow. The speed patterns and volume patterns on an expressway are shown in Figure 5 and Figure 6.

5. Traffic Flow Pattern Recognition Method

5.1. Traffic Flow Pattern Recognition Algorithms

Pattern recognition algorithms based on intelligent algorithms have good performance in the traffic flow field. A back propagation (BP) neural network is widely used because of the simplicity and flexibility of the structure. An optimization algorithm is used to improve the slow convergence and easily fall into the extreme value of a BP neural network, such as a genetic algorithm (GA) and a simulated annealing genetic algorithm (SAGA). A DBN has the ability to extract high-dimensional abstract features of traffic flows through a hidden multilayer structure, which recognizes the pattern of large-scale data efficiently. Therefore, the BP, GA-BP, SAGA-BP and DBN algorithms are used to recognize the traffic flow patterns [41,42,43,44]. The flow chart of each algorithm is shown in Figure 7. The parameters of the algorithms are shown in Table 3.

5.2. Traffic Flow Pattern Recognition Indexes

A variety of traffic flow pattern recognition indexes are designed according to traffic flow characteristics to adapt to different collection conditions of traffic data and traffic application scenarios.
(1)
Index I: Speeds from 0:00 to 24:00
Speeds from 0:00 to 24:00 contain the complete traffic flow characteristics over time, which are suitable for pattern recognition of historical traffic flow data, such as the evaluation of effect on traffic policy.
(2)
Index II: Speeds from 0:00 to 6:00
Speeds from 0:00 to 6:00 describe the characteristics of the early morning traffic flow, which are suitable for real-time traffic scenarios.
(3)
Index III: Speeds from 0:00 to 12:00
Speeds from 0:00 to 12:00 describe the traffic flow characteristics in the early morning and morning peak hours, which have a degree of real-time performance.
(4)
Index IV: Speeds from 6:00 to 10:00
Speeds from 6:00 to 10:00 describe the traffic flow characteristics of morning peak hours, which have a degree of real-time performance. Index IV has fewer dimensions, which is useful to improve the efficiency of pattern recognition.
(5)
Index V: Speeds from 6:00 to 10:00 and 17:00 to 21:00
Speeds from 6:00 to 10:00 and 17:00 to 21:00 describe the traffic flow characteristics of morning and evening peak hours, which have a better ability to reflect the difference in traffic flow characteristics of different categories.
The pattern recognition results of different algorithms based on each index are analyzed on different road types, as shown by taking an expressway as an example in Figure 8.

5.3. Pattern Recognition Evaluation Indexes

The evaluation indexes of the pattern recognition effect include the recognition rate based on a confusion matrix (CM) and the pattern recognition simulation time. The recognition rates and simulation times of different algorithms based on each index are shown in Table 4.
The results of pattern recognition of the different algorithms based on each index are analyzed on different road types, and the conclusions are as follows.
(1) The accuracy and efficiency of pattern recognition of the DBN are the highest under various conditions. The GA-BP and SAGA-BP have better performance than the BP algorithm.
(2) The index of speeds from 0:00 to 24:00 has the best recognition effect, and the average correct recognition rate of all road types is 94.87%. The index of speeds from 6:00 to 10:00 and 17:00 to 20:00 has a good recognition effect, and the average recognition of all road types is 93.26%.
(3) The average recognition rates of the index of speeds from 0:00–12:00 and the index of speeds from 6:00–10:00 exceed 80%, and the highest average recognition rate is 89.39%.
(4) The index of speeds from 0:00 to 6:00 has a low recognition rate, and the correct recognition rate of all road types is 59.36%. The early-morning data in each category have relatively small differences, which are not suitable as indexes of traffic flow patterns.

6. Traffic Flow Estimation Method

The traffic flow estimation method is proposed based on the traffic flow fundamental diagram. The influencing factors on the fundamental diagram are analyzed, and then the fundamental diagrams are constructed for different conditions.

6.1. Analysis of the Factors Influencing the Fundamental Diagram

The traffic flow characteristics are used to develop the traffic flow fundamental diagram based on traffic flow theory, such as the traffic composition, the morning and evening peak and the weather conditions.

6.1.1. Traffic Composition

Traffic composition is the proportion of different vehicle types in the traffic flow. Large vehicles occupy more road space than small passenger vehicles. The ability of large vehicles to accelerate and decelerate is lower than that of small passenger cars. Large vehicles have a larger headway [45]. The proportion of large vehicles is an important factor affecting traffic capacity.
The large-vehicle proportions in different traffic flow patterns are calculated based on the traffic flow data of the North Second Ring Road, North Third Ring Road and North Fourth Ring Road, as shown in Figure 9. The large-vehicle proportion from 0:00 to 6:00 in each pattern is higher, and it is above 20% in each hour except on Sundays. The large-vehicle proportion is below 15% each hour after 7:00 with low volatility. There are differences in the large-vehicle proportion among different patterns, and the large-vehicle proportion on weekends is lower than that on weekdays. Due to the difference in the large-vehicle proportion, there are differences among the parameters of the traffic pattern.
To analyze the influence of the large-vehicle proportion on traffic flow parameters, fundamental diagrams of different large-vehicle proportions at intervals of 10% are developed in this study, as shown in Figure 9. It was found that the traffic capacity decreases as the large-vehicle proportion increases. It is necessary to establish the fundamental diagrams for different patterns and different time periods.

6.1.2. Morning and Evening Peak Hours

The fundamental diagrams in the morning and evening peak hours are developed, taking expressways as an example. The RTMS data and FCD data of 15 sections of the Beijing Expressway from February to October 2018 are used to develop fundamental diagrams. To ensure that there are enough data, the data from 6:00 to 9:00 are selected as the morning peak hours, and the data from 17:00 to 20:00 are selected as the evening peak hours to build fundamental diagrams, as shown in Figure 9. The traffic capacity of morning peak hours is 1600 pcu/h, which is 3.75% larger than that in the evening peak hours. The travel purpose of travelers in the morning peak hours is commuting, which has high driving behavior consistency and traffic flow stability.
To analyze whether the differences in traffic capacity during morning and evening peak hours are consistent on different weekdays, a week of fundamental diagrams of morning and evening peak hours are developed on expressways, as shown in Figure 10. The traffic capacity in the morning peak hours is on average 3.24% higher than that in the evening peak hours, which is consistent and repetitive on different weekdays.

6.1.3. Weather Conditions

It is pointed out in the Highway Capacity Manual [41] that bad weather conditions reduce road speed and capacity compared with normal weather conditions. The data of rainy days on 17 July 2018, and normal days of the same month are used to analyze the influence of weather conditions on the speed and fundamental diagrams, as shown in Figure 11.
In summary, the road type, traffic composition, travel purpose and weather conditions are all factors influencing traffic flow fundamental diagrams, which are different under different traffic flow patterns. Therefore, it is necessary to construct fundamental diagrams for different traffic flow patterns and weekday peak hours.

6.2. Fundamental Diagram Models

The traffic flow fundamental diagrams of highways and expressways are constructed based on the Van Aerde model. The traffic flow fundamental diagrams of major arterial roads and minor arterial roads are constructed based on the Underwood model. The expressway is taken as an example to introduce the MFD model construction method. There are 142 days of expressway flow and speed data selected in this study, among which 85 days of data are used to construct MFD models, and 57 days of data are used to analyze the model uncertainty. The fundamental diagrams for different road types, traffic flow patterns and peak hours of workdays are developed for three proposals, as shown in Figure 12. Traffic flow characteristic parameters of different proposals on expressways are in Table 5.

6.2.1. Proposal I: Fundamental Diagrams for Different Road Types

Proposal I is to develop fundamental diagrams for different road types. The fundamental diagrams on expressways are shown in Figure 12a. The characteristic parameters on expressways for Protocol I are shown in Table 5.

6.2.2. Proposal II: Fundamental Diagrams for Different Traffic Flow Patterns

Proposal II is to develop fundamental diagrams for different traffic flow patterns and road types based on the traffic flow patterns constructed in this study, as shown in Figure 12b. The characteristic parameters for Protocol II are shown in Table 5.
The fundamental diagram of each pattern is different. There are lower large-vehicle proportions and higher traffic capacities on weekends and holidays than on weekdays. The traffic capacity on rainy days is lower than that on days of normal weather. The traffic capacity on Saturday, Sunday and holidays patterns is 2.61%, 3.58%, 5.21% higher than that in the Monday to Thursday weekday pattern, respectively, and 3.28%, 4.26%, 5.90% higher than that in Friday weekday pattern. The capacity of the rainy-day pattern is reduced by 6.84% and 6.23% compared with the Monday to Thursday weekday pattern and the Friday weekday pattern, respectively.

6.2.3. Proposal III: Fundamental Diagrams for Morning/Evening Peak Hours of Weekday Patterns

Proposal III is to develop the fundamental diagrams for peak hours of weekday patterns and road types, as shown in Figure 12c. The characteristic parameters for Protocol III are shown in Table 5.
The capacity of morning peak hours is higher than that of the evening peak hours in each pattern. The capacity of morning peak hours of the Monday to Thursday weekday pattern and Friday weekday pattern is 3.66% and 4.11% higher than that of evening peak hours, respectively.
The fundamental diagrams of Proposal I, Proposal II and Proposal III on the Beijing highways, expressways, major arterial roads and minor arterial roads are developed in this study.

6.2.4. Error Analysis

The RTMS and FCD data that have not been used for the model construction of 52, 57, 93 and 74 days on Beijing highways, expressways, major arterial roads and minor arterial roads are selected to estimate the accuracy of the MFD models under the three proposals. The FCD speed data are used to estimate the volume data based on the fundamental diagrams. The errors from 6:00 to 24:00 between the fundamental diagrams estimated volumes and the field measured volumes are calculated, as shown in Table 6.
The prediction accuracy of the fundamental diagrams for different weekday peak hours and road types from 6:00 to 24:00 are improved. The average relative errors of Proposal I, Proposal II and Proposal III on the highway are 15.76%, 11.82% and 9.76%, respectively, and the average relative errors of Proposal III and Proposal II are 6.06% and 3.94% lower than that of Proposal I, respectively. The average relative errors of Proposal I, Proposal II and Proposal III on the expressway are 12.65%, 8.88% and 7.64%, respectively, and the average relative errors of Proposal III and Proposal II are 5.01% and 3.77% lower than that of Proposal I, respectively. The average relative errors of Proposal I, Proposal II and Proposal III on the major arterial road are 14.45%, 11.09% and 8.58%, respectively, and the average relative errors of Proposal III and Proposal II are 5.87% and 3.36% lower than that of Proposal I, respectively. The average relative errors of Proposal I, Proposal II and Proposal III on the minor arterial road are 14.35%, 8.74% and 7.81%, respectively, and the average relative errors of Proposal III and Proposal II are 6.54% and 5.61% lower than that of Proposal I, respectively.
Above all, Proposal III has the highest accuracy. Representative roads of different road types on 3 April 2018 are selected to further evaluate the accuracy of the fundamental diagrams of Protocol III. The estimated volume of Protocol III is compared with the measured RTMS volume, as shown in Figure 13.
The average relative error of the Beijing-Kaisha highway is 7.19%, and the average absolute error of that road is 48 pcu/h. The average relative error of the North Fourth Ring Expressway is 7.47%, and the average absolute error of that road is 49 pcu/h. The average relative error of the Fangzhuang major arterial road is 6.83%, and the average absolute error of that road is 32 pcu/h. The average relative error of the Guang’an minor arterial road is 7.25%, and the average absolute error of that road is 16 pcu/h. The accuracy of volumes calculated by the fundamental diagrams of Protocol III constructed in this paper is relatively high.

6.3. Traffic Flow Estimation Method of the Whole Sample on the Urban Road Network

Fundamental diagrams for road types, traffic flow patterns, and weekday peak hours are developed in this study. Whole sample dynamic traffic volumes of the urban road network are calculated based on the real-time FCD data. The calculation steps for whole sample dynamic traffic volumes of urban road network are as follows:
(1)
The speed data input and quality control.
(2)
The map data input: including the road type, the number of lanes, the road length, etc.
(3)
The speed data and map data matching: the speed is matched with the map data.
(4)
Pattern recognition: the input speed is used for traffic flow pattern recognition based on the typical traffic flow pattern library and traffic flow pattern recognition method established in this study.
(5)
Fundamental diagram selection and flow calculation: the fundamental diagram selection based on traffic flow pattern and time information, flow data calculation and quality control.
The traffic flow of each road link can be calculated based on the traffic flow estimation method, which can be further used to calculate the dynamic vehicle emissions of the urban road network.

7. Case Study

The traffic restriction implementation effect of diesel trucks of the national stage III emission standard is evaluated as a case study. The changes in dynamic vehicle emissions before and after the implementation of the traffic restriction policy are calculated.
The traffic restriction policy of diesel trucks of the national stage III emission standard was published on 1 December 2018. A diesel truck of the national stage III emission standard is forbidden to drive within the Fifth Ring Road (excluding) from 6:00 to 23:00 every day. Trucks are forbidden to drive over 8 tons (including) on the Fifth Ring Road. Data from weekdays on 5 November, 6 November and 7 November 2018, before the implementation of the policy, and 3 December, 4 December and 5 December 2018, after the implementation of the policy, are used to calculate the emissions within the Fifth Ring Road.
The calculation method of the road link emissions is shown in Equation (12):
E i = E F j R , v × Q i × L i × P j R
where Ei is the emission of road link i (ton); E F j R , v is the emission factor (g/km) of vehicle type j when its speed is v on road type R, Q i is the traffic volume of road link i (pcu/h); is the length of road link i (km); and L i is the proportion of vehicle type j on road type R. The vehicle emissions of the road network were the total emissions of each road link.
The emissions on the road network of each pollutant before and after the implementation of the traffic restriction policy in different periods is calculated, such as from 0:00 to 24:00, 6:00 to 23:00, 7:00 to 9:00, 17:00 to 19:00 and 20 to 22:00, as shown in Table 7.
The emissions of each pollutant on the whole road network within the Fifth Ring Road all decrease after the implementation of the traffic restriction policy from 6:00 to 23:00. The CO2 emissions decrease from 36,817.12 tons to 34,144.84 tons, with a decrease rate of 7.26%. The NOx emission decreases from 135.96 tons to 118.05 tons, with a decreasing ratio of 13.18%. The CO emission decreases from 522.25 tons to 495.58 tons, with a decreasing ratio of 5.11%. The PM emissions decrease from 1.95 tons to 1.57 tons, with a decreasing ratio of 19.78%. The HC emission decreases from 96.56 tons to 91.56 tons, with a decreasing ratio of 5.18%.
The emission intensity on the road network of each pollutant before and after the implementation of the traffic restriction policy in different periods is calculated, such as from 0:00 to 24:00, 6:00 to 23:00, 7:00 to 9:00, 17:00 to 19:00 and 20 to 22:00, as shown in Table 8.
The emission intensity on the road network of each pollutant before and after the implementation of the traffic restriction policy in different periods is calculated, such as from 0:00 to 24:00, 7:00 to 9:00, 17:00 to 19:00 and 20 to 22:00 before and after the implementation of the traffic restriction policy is shown on the GIS map in Figure 14.
The emission intensities of CO2, NOx, CO, HC and PM in the road network all decrease after the traffic restriction. The emission intensities of CO2, NOx, CO, HC and PM decreased by 5.50%, 11.48%, 3.28%, 3.33% and 19.05%, respectively, from 0:00 to 24:00. The emission intensities of CO2, NOx, CO, HC and PM decreased by 6.03%, 13.08%, 3.48%, 3.55% and 21.44%, respectively, from 6:00 to 23:00. The emission intensities of CO2, NOx, CO, HC and PM decreased by 27.91%, 59.52%, 14.03%, 15.68% and 73.76%, respectively, from 2:00 to 4:00. The emission intensities of CO2, NOx, CO, HC and PM decreased by 3.12%, 7.84%, 1.20%, 1.43% and 14.89%, respectively, from 7:00 to 9:00. The emission intensities of CO2, NOx, CO, HC and PM decreased by 5.53%, 10.66%, 3.26%, 3.31% and 18.16%, respectively, from 10:00 to 12:00. The emission intensities of CO2, NOx, CO, HC and PM decreased by 6.52%, 11.37%, 4.57%, 4.43% and 18.27%, respectively, from 17:00 to 19:00. The emission intensities of CO2, NOx, CO, HC and PM decreased by 6.53%, 14.46%, 3.89%, 3.93% and 22.50%, respectively, from 20:00 to 22:00.

8. Conclusions

The main research results and conclusions of this paper are as follows:
(1) The SOM clustering method based on the direct indexes of time-varying speed differences effectively clusters urban traffic flow. Traffic flow patterns describe more than 90% of the traffic flows on urban road networks, such as Monday to Thursday, Friday, Saturday, Sunday, rainy days, holidays, prominent evening peak and prominent morning peak patterns.
(2) The recognition accuracy and efficiency of the DBN are the highest under various conditions. The GA-BP and SAGA-BP have better performance than the BP algorithm. The average correct recognition rates of different road types of speeds from 0:00 to 24:00, 0:00 to 12:00, 6:00 to 10:00, 6:00 to 10:00 and 17:00 to 21:00 are 94.87%, 82.45%, 80.96%, 93.26%, respectively. The early-morning data in each category are not suitable for traffic flow pattern recognition.
(3) In addition to the large-vehicle proportion and weather conditions identified by existing studies, it was found that the difference in driving behavior caused by different travel purposes significantly affected the fundamental diagrams, such as morning and evening peak hours, weekdays and holidays. The traffic capacity of the same road link in the morning peak hours is 3.47% higher than that in the evening peak hours. The traffic capacity of the same road link on holidays is 4.73% higher than that on weekdays. MFD models considering the road type, morning and evening peak hours, weekdays and holidays improve the flow measurement accuracy by 6.51% compared with that considering only road type.
(4) The case study shows that based on the traffic flow pattern clustering and the recognition method proposed in this paper, the traffic flow pattern is identified quickly, and the whole sample traffic flow on the urban road network is calculated quickly by speed data. The dynamic emissions on the road network are further acquired.

Author Contributions

This work is the result of collaboration among J.Z., P.J., G.S., Z.L. and T.P. All authors have equally contributed, reviewed, and improved the manuscript. All authors have revised the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Fundamental Research Funds for the Universities of Beijing (X21063) and the Natural Science Foundation of China (NSFC) (No. 52172301, 51678045), Major Program of the National Social Science Foundation of China (21ZDA029).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Some or all data and models that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ahmad, M.; Manjantrarat, T.; Rattanawongsa, W.; Muensri, P.; Saenmuangchin, R.; Klamchuen, A.; Aueviriyavit, S.; Sukrak, K.; Kangwansupamonkon, W.; Panyametheekul, S. Chemical Composition, Sources, and Health Risk Assessment of PM2.5 and PM10 in Urban Sites of Bangkok, Thailand. Int. J. Environ. Res. Public Health 2022, 19, 14281. [Google Scholar] [CrossRef]
  2. Gupta, P.; Jangid, A.; Kumar, R. COVID-19-associated 2020 lockdown: A study on atmospheric black carbon fall impact on human health. Environ. Geochem. Health 2022, 1–14. [Google Scholar] [CrossRef]
  3. Lin, H.; Zhu, J.; Jiang, P.; Cai, Z.; Yang, X.; Yang, X.; Zhou, Z.; Wei, J. Assessing drivers of coordinated control of ozone and fine particulate pollution: Evidence from Yangtze River Delta in China. Environ. Impact Assess. Rev. 2022, 96, 106840. [Google Scholar] [CrossRef]
  4. Kakiuchi, T.; Yajima, T.; Shigaki, N.; Kawajiri, Y. Modeling and optimal design of multicomponent vacuum pressure swing adsorber for simultaneous separation of carbon dioxide and hydrogen from industrial waste gas. Adsorption 2022, 1–19. [Google Scholar] [CrossRef]
  5. Smaisim, G.F.; Abed, A.M.; Alavi, H. Analysis of pollutant emission reduction in a coal power plant using renewable energy. Int. J. Low-Carbon Technol. 2022. [Google Scholar] [CrossRef]
  6. Tao, H.; Xing, J.; Pan, G.; Pleim, J.; Ran, L.; Wang, S.; Chang, X.; Li, G.; Chen, F.; Li, J. Impact of anthropogenic heat emissions on meteorological parameters and air quality in Beijing using a high-resolution model simulation. Front. Environ. Sci. Eng. 2022, 16, 44. [Google Scholar] [CrossRef]
  7. Davison, J.; Rose, R.A.; Farren, N.J.; Wagner, R.L.; Murrells, T.P.; Carslaw, D.C. Verification of a National Emission Inventory and Influence of On-road Vehicle Manufacturer-Level Emissions. Environ. Sci. Technol. 2021, 55, 4452–4461. [Google Scholar] [CrossRef]
  8. Zou, C.; Wu, L.; Wang, Y.; Sun, S.; Wei, N.; Sun, B.; Ni, J.; He, J.; Zhang, Q.; Peng, J.; et al. Evaluating traffic emission control policies based on large-scale and real-time data: A case study in central China. Sci. Total Environ. 2022, 160435. [Google Scholar] [CrossRef]
  9. Zhu, Y.; Wu, Q.; Xiao, N. Research on highway traffic flow prediction model and decision-making method. Sci. Rep. 2022, 12, 19919. [Google Scholar] [CrossRef]
  10. Niu, J.; He, J.; Li, Y.; Zhang, S. Highway Temporal-Spatial Traffic Flow Performance Estimation by Using Gantry Toll Collection Samples: A Deep Learning Method. Math. Probl. Eng. 2022, 2022, 8711567. [Google Scholar] [CrossRef]
  11. Ambros, J.; Usami, D.S.; Valentová, V. Developing speed-related safety performance indicators from floating car data. IET Intell. Transp. Syst. 2022. [Google Scholar] [CrossRef]
  12. Wang, H.-W.; Peng, Z.-R.; Wang, D.; Meng, Y.; Wu, T.; Sun, W.; Lu, Q.-C. Evaluation and prediction of transportation resilience under extreme weather events: A diffusion graph convolutional approach. Transp. Res. Part C Emerg. Technol. 2020, 115, 102619. [Google Scholar] [CrossRef]
  13. Li, Y.; Ramezani, M. Quasi revenue-neutral congestion pricing in cities: Crediting drivers to avoid city centers. Transp. Res. Part C Emerg. Technol. 2022, 145, 103932. [Google Scholar] [CrossRef]
  14. Jiang, Y.; Song, G.; Zhang, Z.; Zhai, Z.; Yu, L. Estimation of Hourly Traffic Flows from Floating Car Data for Vehicle Emission Estimation. J. Adv. Transp. 2021, 2021, 6628335. [Google Scholar] [CrossRef]
  15. Ministry of Ecological Environment, PRC. China Mobile Source Environmental Management Annual Report in 2021; Ministry of Ecological Environment: Beijing, China, 2021. [Google Scholar]
  16. Zang, J. Method of Traffic Flow Estimation of the Whole Sample of Urban Road Network for Emission Estimation; Beijing Jiaotong University: Beijing, China, 2020. [Google Scholar]
  17. Meng, D.; Song, G.; Wu, Y.; Zhai, Z.; Yu, L.; Zhang, J. Modification of Newell’s car-following model incorporating multidimensional stochastic parameters for emission estimation. Transp. Res. Part D Transp. Environ. 2021, 91, 102692. [Google Scholar] [CrossRef]
  18. Ge, M.; Song, G.; Zang, J.; Wu, Y.; Yu, L. Powered Link-based traffic volume forecasting for dynamic emission estimation based on pattern clustering and recognition. In Proceedings of the 100th Transportation Research Board meeting, Washington, DC, USA, 5–29 January 2021. [Google Scholar]
  19. Hu, Z.; Huang, J.; Yang, D.; Zhong, Z. Constraint-tree-driven modeling and distributed robust control for multi-vehicle cooperation at unsignalized intersections. Transp. Res. Part C Emerg. Technol. 2021, 131, 103353. [Google Scholar] [CrossRef]
  20. Fei, W. Prediction Model of Congestion Propagation and Dissipation for Dynamic Route Guidance under Incidents; Beijing Jiaotong University: Beijing, China, 2017. [Google Scholar]
  21. Xia, D.; Jiang, S.; Yang, N.; Hu, Y.; Li, Y.; Li, H.; Wang, L. Discovering spatiotemporal characteristics of passenger travel with mobile trajectory big data. Phys. A Stat. Mech. Its Appl. 2021, 578, 126056. [Google Scholar] [CrossRef]
  22. Ma, D.; Xiao, J.; Song, X.; Ma, X.; Shen, J. A Back-Pressure-Based Model with Fixed Phase Sequences for Traffic Signal Optimization Under Oversaturated Networks. IEEE Trans. Intell. Transp. Syst. 2020, 22, 5577–5588. [Google Scholar] [CrossRef]
  23. Greenshields, B. A study of traffic capacity. Highw. Res. Board Proc. 1935, 14, 448–477. [Google Scholar]
  24. Greenberg, H. An analysis of traffc flow. Oper. Res. 1959, 7, 79–85. [Google Scholar] [CrossRef]
  25. Underwood, R. Speed, volume, and density relationships: Quality and theory of traffic flow. Yale Bur. Highw. Traffic 1961, 141–188. [Google Scholar]
  26. Pipes, L. Car following models and the fundamental diagram of road traffic. Transp. Res. 1967, 1, 21–29. [Google Scholar] [CrossRef]
  27. Van Aerde, M. Single regime speed-flow-density relationship for congested and uncontested highways. In Proceedings of the 74th TRB Annual Meeting, Washington, DC, USA, 22–28 January 1995. [Google Scholar]
  28. Wang, J.; Yu, L.; Gao, Y.; Zhang, J.; Song, G. Analysis of Speed Differences in Data from Remote Traffic Microwave Sensors and Floating Car Data Systems on Expressways. Transp. Res. Rec. J. Transp. Res. Board 2017, 2643, 112–120. [Google Scholar] [CrossRef]
  29. Daganzo, C.F. Urban gridlock: Macroscopic modeling and mitigation approaches. Transp. Res. Part B Methodol. 2007, 41, 49–62. [Google Scholar] [CrossRef]
  30. Guo, Y.; Yang, L.; Hao, S.; Gu, X. Perimeter traffic control for single urban congested region with macroscopic fundamental diagram and boundary conditions. Phys. A Stat. Mech. Its Appl. 2021, 562, 125401. [Google Scholar] [CrossRef]
  31. Chen, Z.; Wu, W.-X.; Huang, H.-J.; Shang, H.-Y. Modeling traffic dynamics in periphery-downtown urban networks combining Vickrey’s theory with Macroscopic Fundamental Diagram: User equilibrium, system optimum, and cordon pricing. Transp. Res. Part B Methodol. 2022, 155, 278–303. [Google Scholar] [CrossRef]
  32. Cheng, Q.; Liu, Z.; Lin, Y.; Zhou, X. An s-shaped three-parameter (S3) traffic stream model with consistent car following relationship. Transp. Res. Part B Methodol. 2021, 153, 246–271. [Google Scholar] [CrossRef]
  33. Cheng, Q.; Liu, Z.; Guo, J.; Wu, X.; Pendyala, R.; Belezamo, B.; Zhou, X. Estimating key traffic state parameters through parsimonious spatial queue models. Transp. Res. Part C Emerg. Technol. 2022, 137, 103596. [Google Scholar] [CrossRef]
  34. Jia, T.; Yan, P. Predicting Citywide Road Traffic Flow Using Deep Spatiotemporal Neural Networks. IEEE Trans. Intell. Transp. Syst. 2020, 22, 3101–3111. [Google Scholar] [CrossRef]
  35. Huang, X.; Ye, Y.; Wang, C.; Yang, X.; Xiong, L. A multi-mode traffic flow prediction method with clustering based attention convolution LSTM. Appl. Intell. 2021, 52, 14773–14786. [Google Scholar] [CrossRef]
  36. Ma, D.; Ben Song, X.; Zhu, J.; Ma, W. Input data selection for daily traffic flow forecasting through contextual mining and intra-day pattern recognition. Expert Syst. Appl. 2021, 176, 114902. [Google Scholar] [CrossRef]
  37. Lin, X.; Sun, R.; Wang, Y. Improved Key Performance Indicator-Partial Least Squares Method for Nonlinear Process Fault Detection Based on Just-In-Time Learning. J. Frankl. Inst. 2022. [Google Scholar] [CrossRef]
  38. Banerjee, P.; Chattopadhyay, T.; Chattopadhyay, K. Comparison among different Clustering and Classification Techniques: Astronomical data-dependent study. New Astron. 2022, 100, 101973. [Google Scholar] [CrossRef]
  39. Kohonen, T. The self-organizing map. Proc. IEEE 1990, 78, 1464–1480. [Google Scholar] [CrossRef]
  40. Yang, S.; Yuan, H.; Dong, L. Offshore wind resource assessment by characterizing weather regimes based on self-organizing map. Environ. Res. Lett. 2022, 17, 124009. [Google Scholar] [CrossRef]
  41. Li, J.; Cheng, J.H.; Shi, J.Y.; Huang, F. Brief introduction of back propagation (BP) neural network algorithm and its improvement. In Advances in Computer Science and Information Engineering; Springer: Berlin/Heidelberg, Germany, 2012; pp. 553–558. [Google Scholar] [CrossRef]
  42. Yi, H.; Zhang, K.; Ma, K.; Zhou, L.; Tang, F. Prediction of Natural Rubber Customs Declaration Price Based on Wavelet Decomposition and GA-BP Neural Network Group. Mathematics 2022, 10, 4264. [Google Scholar] [CrossRef]
  43. Wang, W.; Zhu, Q.; Wang, Z.; Zhao, X.; Yang, Y. Research on Indoor Positioning Algorithm Based on SAGA-BP Neural Network. IEEE Sens. J. 2021, 22, 3736–3744. [Google Scholar] [CrossRef]
  44. Cai, Y.; Fan, P.; Lang, S.; Li, M.; Muhammad, Y.; Liu, A. Downscaling of SMAP soil moisture data by using a Deep Belief Network. Remote Sens. 2022, 14, 5681. [Google Scholar] [CrossRef]
  45. National Research Council (U.S.). HCM2010: Highway Capacity Manual, 5th ed.; Transportation Research Board: Washington, DC, USA, 2010. [Google Scholar]
Figure 1. Methodology flow chart.
Figure 1. Methodology flow chart.
Ijerph 19 16524 g001
Figure 2. Geometric characteristics of traffic flow speed curves.
Figure 2. Geometric characteristics of traffic flow speed curves.
Ijerph 19 16524 g002
Figure 3. Evaluation of clustering methods under different indexes with different colors.
Figure 3. Evaluation of clustering methods under different indexes with different colors.
Ijerph 19 16524 g003
Figure 4. The categories of DI based on different algorithms (curves with different colors mean different days).
Figure 4. The categories of DI based on different algorithms (curves with different colors mean different days).
Ijerph 19 16524 g004
Figure 5. Clustering of expressway speed (curves with different colors mean different days).
Figure 5. Clustering of expressway speed (curves with different colors mean different days).
Ijerph 19 16524 g005
Figure 6. Clustering of expressway volume (curves with different colors mean different days).
Figure 6. Clustering of expressway volume (curves with different colors mean different days).
Ijerph 19 16524 g006
Figure 7. The flow chart of each algorithm.
Figure 7. The flow chart of each algorithm.
Ijerph 19 16524 g007
Figure 8. Pattern recognition results of different algorithms based on each index on an expressway.
Figure 8. Pattern recognition results of different algorithms based on each index on an expressway.
Ijerph 19 16524 g008
Figure 9. Average speed-flow model under different large-vehicle proportions.
Figure 9. Average speed-flow model under different large-vehicle proportions.
Ijerph 19 16524 g009
Figure 10. Average speed-flow models of peak hours on different weekdays on expressways.
Figure 10. Average speed-flow models of peak hours on different weekdays on expressways.
Ijerph 19 16524 g010
Figure 11. Fundamental diagrams of rainy days and normal days.
Figure 11. Fundamental diagrams of rainy days and normal days.
Ijerph 19 16524 g011
Figure 12. Fundamental diagram models of different proposals on expressways: (a) Proposal I, (b) Proposal II, (c) Proposal III.
Figure 12. Fundamental diagram models of different proposals on expressways: (a) Proposal I, (b) Proposal II, (c) Proposal III.
Ijerph 19 16524 g012
Figure 13. Measured volume and calculated volume curves on different road types.
Figure 13. Measured volume and calculated volume curves on different road types.
Ijerph 19 16524 g013
Figure 14. Emission intensity of the road network in different periods.
Figure 14. Emission intensity of the road network in different periods.
Ijerph 19 16524 g014
Table 1. Meaning of CI.
Table 1. Meaning of CI.
IndexCIIndexCI
X1maximum speedX11average speed from 0:00–24:00
X2minimum speedX12speed variance from 0:00–24:00
X3ratio of the maximum and the minimum speedX13speed standard deviation from 0:00–24:00
X4average speed from 0:00–6:00X14kurtosis from 6:00–12:00
X5average speed from 20:00–23:00X15skewness from 6:00–12:00
X6average speed from 10:00–12:00X16kurtosis from 12:00–22:00
X7average speed from 2:00–4:00X17skewness from 12:00–22:00
X8average speed of morning peakX18kurtosis from 10:00–14:00
X9average speed of evening peakX19skewness from 10:00–14:00
X10ratio of the speed of morning peak and evening peak
Table 2. Meaning of SI.
Table 2. Meaning of SI.
CompositionFirstSecondThirdFourthFifth
FeatureAverage FeatureHigh Speed FeatureConcave FeatureConvex FeatureTraffic Flow and Geometric Feature
IndexX11X10X16X18X1
X9X4X14X10X4
X7X1X10X8X18
X2X15X6X19X3
X6X14X13X15X12
Table 3. Parameters of algorithms.
Table 3. Parameters of algorithms.
AlgorithmParameter Design
Parameter of BPNameHidden layer nodeMaximum epochsLearning rateAccuracy
Value m = n + l + α 1000.110.001
Parameter of GANamePopulation sizeMutation probabilityCrossover probabilityEvolutionary algebra
Value500.050.4200
Parameter of SANameInitial temperatureAnnealing speedIterations
Value10000.85100
Table 4. Recognition rates and simulation times of different algorithms based on each index.
Table 4. Recognition rates and simulation times of different algorithms based on each index.
Road TypeAlgorithmRecognition Rate of Indexes (%)Simulation Time (s) of Indexes
IIIIIIIVVIIIIIIIVV
HighwayBP78.7962.1266.6781.8277.275.923.565.003.953.80
GA-BP89.3950.0066.6789.3989.390.441.111.341.271.16
SAGA-BP83.3362.1272.7387.8889.390.391.090.771.221.06
DBN93.9466.6777.2765.1580.300.070.070.070.080.07
ExpresswayBP82.8641.4372.8672.8688.572.984.594.344.144.75
GA-BP90.0047.1475.7172.8684.290.411.031.271.081.23
SAGA-BP90.0052.8672.8661.4385.710.341.311.311.190.38
DBN94.2944.2980.0064.2994.290.070.070.070.070.10
Major Arterial RoadBP87.6249.5282.8678.1092.383.273.863.583.582.91
GA-BP89.5250.4882.8683.8188.570.441.030.310.310.22
SAGA-BP90.4856.1988.5782.8692.380.421.060.330.330.27
DBN96.1955.2488.5777.1494.290.070.110.100.100.09
Minor Arterial RoadBP77.7861.7374.0771.6188.895.974.394.144.272.94
GA-BP87.6554.3277.7874.0792.590.421.191.531.251.27
SAGA-BP95.0659.2666.6772.8482.720.421.051.641.170.52
DBN87.6559.2683.9577.7895.060.080.080.080.080.07
Optimal Algorithm Recognition Rate94.8759.3682.4580.9693.260.080.080.080.080.08
Table 5. Traffic flow characteristic parameters of different proposals on expressways.
Table 5. Traffic flow characteristic parameters of different proposals on expressways.
ParameterCapacity (pcu/h)Free Flow Speed (km/h)Critical Speed (km/h)Blocking Density (pcu/km)
Parameters of Proposal I
Value158072.531.2148.5
Parameters of Proposal II
Pattern 1153571.232.6149.4
Pattern 2152571.131.8150.6
Pattern 3157571.433.2170.2
Pattern 4159071.332.8178.2
Pattern 5143070.232.1148
Pattern 6161571.833.8179
Pattern 7155079.833.2148
Pattern 8160078.533.2156.5
Parameters of Proposal III
Mon. to Thur. peakAM164073.532.9150.2
PM158072.931.9148.6
Fri. peakAM157070.833.6152.1
PM151569.830.2136.2
Table 6. Error between the measured flow and calculated flow of the three proposals.
Table 6. Error between the measured flow and calculated flow of the three proposals.
Road TypeProposal IProposal IIProposal IIIProposal III Lower than Proposal IProposal II Lower than Proposal I
Highway15.76%11.82%9.70%6.06%3.94%
Expressway12.65%8.88%7.64%5.01%3.77%
Major arterial road14.45%11.09%8.58%5.87%3.36%
Minor arterial road14.35%8.74%7.81%6.54%5.61%
Table 7. The emissions of each pollutant before and after the implementation of the traffic restriction policy in different periods.
Table 7. The emissions of each pollutant before and after the implementation of the traffic restriction policy in different periods.
Pollutant (ton)Time0:00–24:006:00–23:002:00–4:007:00–9:0010:00–12:0017:00–19:0020:00–22:00
CO2Before38,780.9836,817.12527.754726.344061.094340.772683.81
After35,734.0934,144.84428.194309.213774.884030.482315.86
Change−7.86%−7.26%−18.87%−8.83%−7.05%−7.15%−13.71%
COBefore546.71522.256.6167.3156.6361.2735.06
After517.55463.655.9762.7754.0158.1132.29
Change−5.33%−5.11%−9.77%−6.74%−4.63%−5.15%−7.91%
NOXBefore144.50135.962.9417.6415.3116.078.97
After122.22118.051.7615.2413.3513.965.28
Change−15.42%−13.18%−40.17%−13.64%−12.79%−13.12%−41.15%
HCBefore100.9796.561.1712.3710.4711.516.49
After95.5091.561.0511.529.9810.935.97
Change−5.42%−5.18%−10.34%−6.92%−4.64%−5.10%−7.97%
PMBefore2.091.950.050.250.220.240.14
After1.611.570.020.200.180.190.06
Change−22.84%−19.78%−53.44%−19.80%−19.71%−18.97%−56.13%
Table 8. The emissions intensity before and after the implementation of the traffic restriction policy in different periods.
Table 8. The emissions intensity before and after the implementation of the traffic restriction policy in different periods.
Pollutant (kg/km/h)Time0:00–24:006:00–23:002:00–4:007:00–9:0010:00–12:0017:00–19:0020:00–22:00
CO2Before163,130.88134,208.3714,588.51 172,638.51 165,367.55 198,057.24 145,382.73
After154,186.954126,115.5810,516.21 167,256.72 156,215.96 185,152.65 136,193.86
Change−5.50%−6.03%−27.91%−3.12%−5.53%−6.52%−6.32%
COBefore624.292518.8382.20 665.68 647.79 743.56 545.60
After553.772450.9433.27 613.47 578.75 659.00 466.70
Change−11.48%−13.08%−59.52%−7.84%−10.66%−11.37%−14.46%
NOXBefore2378.9361945.46169.02 2528.96 2396.91 2927.59 2095.76
After2300.6421877.75145.31 2498.63 2318.82 2793.70 2014.31
Change−3.28%−3.48%−14.03%−1.20%−3.26%−4.57%−3.89%
HCBefore434.278355.0129.92 462.06 436.31 528.26 389.75
After419.802342.4025.23 455.43 421.85 504.88 374.45
Change−3.33%−3.55%−15.68%−1.43%−3.31%−4.43%−3.93%
PMBefore9.247.751.52 9.65 9.63 10.75 8.42
After7.56.090.40 8.21 7.88 8.79 6.53
Change−19.05%−21.44%−73.76%−14.89%−18.16%−18.27%−22.50%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zang, J.; Jiao, P.; Song, G.; Li, Z.; Peng, T. A Novel Environment Estimation Method of Whole Sample Traffic Flows and Emissions Based on Multifactor MFD. Int. J. Environ. Res. Public Health 2022, 19, 16524. https://doi.org/10.3390/ijerph192416524

AMA Style

Zang J, Jiao P, Song G, Li Z, Peng T. A Novel Environment Estimation Method of Whole Sample Traffic Flows and Emissions Based on Multifactor MFD. International Journal of Environmental Research and Public Health. 2022; 19(24):16524. https://doi.org/10.3390/ijerph192416524

Chicago/Turabian Style

Zang, Jinrui, Pengpeng Jiao, Guohua Song, Zhihong Li, and Tingyi Peng. 2022. "A Novel Environment Estimation Method of Whole Sample Traffic Flows and Emissions Based on Multifactor MFD" International Journal of Environmental Research and Public Health 19, no. 24: 16524. https://doi.org/10.3390/ijerph192416524

APA Style

Zang, J., Jiao, P., Song, G., Li, Z., & Peng, T. (2022). A Novel Environment Estimation Method of Whole Sample Traffic Flows and Emissions Based on Multifactor MFD. International Journal of Environmental Research and Public Health, 19(24), 16524. https://doi.org/10.3390/ijerph192416524

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop