Traffic Flow Modeling of Freeway Variable Speed Limit Control Based on the Big Data of Driving Behavior

Variable speed limit (VSL) control is a flexible restriction on the rate at which motorists can drive on a given stretch of road. Effective VSL control can increase safety and provide clear guidance for motorists. Previous traffic flow models of VSL control were mostly based on the influence of VSL on average speed (macro) or driver’s expected speed (micro). Few models considered the influence of VSL on driver’s actual driving behavior. In this paper, we first briefly introduce the big traffic data involved in this study and explain the mapping relationship between the data and driving behavior. (en, we analyze the driver’s actual driving behavior under the VSL control. (en, an improved single-lane cellular automaton model is established based on the driving behavior characteristics under VSL control. After that, we calibrate the parameters of the single-lane cellular automaton model with the left lane as the calibration object. Finally, this paper uses the proposed single-lane cellular automaton model to simulate the traffic flow characteristics under VSL control. (e numerical simulation results show that the simulation of the variable speed limit in different density intervals presents different results, but these results are consistent with the actual situation of variable speed limit control, which verifies the validity of the proposed model.


Introduction
e variable speed limit (VSL) control system adjusts the road speed limit in real time based on the traffic flow state detected by advanced traffic flow and environment detection technology. As one of the advanced highway management control methods, VSL has been applied in Western countries such as Germany, the United States, Netherlands, the United Kingdom, Finland, Israel, New Zealand, and Australia [1]. e implementation effect indicates that VSL has a great impact on improving road traffic safety and alleviating traffic congestion [2][3][4][5].
e essential basic research content of VSL is traffic flow modeling. Macroscopic models are an important component of VSL modeling, which can be divided into two categories in terms of the characteristics of different VSL traffic flows [6][7][8][9]. One category is well studied, considering the influence of VSL on the equilibrium speed relationship [10][11][12][13][14]. e equilibrium speed function is replaced by a new one with the VSL control parameter according to the modified fundamental diagram under VSL control. Another category establishes the macroscopic model using the correspondence between macroscopic and microscopic parameters when simulating the driver behavior under VSL control.
e macroscopic model has the advantage of high operational efficiency in numerical simulation, but most of the existing macroscopic models lack the behavior description of individual vehicles, which makes it impossible to simulate the safety benefits of VSL. Hence, the optimization goal of VSL is limited to efficiency. e microscopic model is preferred because both the efficiency benefits and security benefits of VSL strategies can be evaluated. Lee et al. [15] used PARAMICS to simulate VSL control; Torday and Bierlaire [16] modeled the VSL microscopic model based on the MITSIM simulator; and Park and Yadlapati [17] used VISSIM to model the VSL with different compliance rates. Other microscopic traffic simulations of Cruise Control (CC) and Adaptive Control (ACC) are used [18][19][20]; such models were used for freeway traffic for the evaluation of the effectiveness of VSL [21,22]. However, in most of the microscopic models, the phenomenon that VSL affects the vehicle's desired speed is pondered, but the changes in drivers' behavior are ignored. In fact, VSL performance depended greatly on compliance rates [23]; the VSL not only changes the desired speed but also changes the following behavior of drivers. is is reflected in a previous paper where the proportion of the smaller headway was reduced [24]. Besides, models on considering the desired speed are not reasonable. In most of the studies, vehicles' speed was set equal to the VSL value, which does not correspond with the real estimation. When simulating large or small value of VSL control, the desired speed would be slightly above or beneath the limit value, which is worth exploring but had not been investigated in existing models. erefore, it is necessary to explore the impact regularity of VSL control under different traffic conditions and compliance rates.
Traffic data always plays an important role in the process of traffic flow modeling; especially in the research of microscopic traffic flow modeling, it is necessary to describe driving behavior characteristics through massive microscopic traffic data. Recently, many scholars have tried to use driving behavior data as a research entry point to solve a series of related transportation problems. Li et al. [25] analyzed drivers' deceleration behavior based on naturalistic driving data and, finally, a brake control system was designed. In order to evaluate whether a driving behavior is fuel-efficient, Hao et al. [26] mined the traffic operating data and fuel consumption data and finally built a corresponding analysis model. Deng et al. [27] established a curve safe speed model on the basis of a thorough analysis of driving behavior data. e use of empirical data to investigate driving behaviors with VSL control was well-studied [10,12,17]. e consensus emerging from these studies is that VSL systems can improve safety by decreasing the mean speed value, the speed difference, and the percentage of small-time headway and can increase the mean time headway. Studies using empirical traffic data as input can objectively analyze the impacts of VSL on driving behaviors and derive credible conclusions. However, with the constraint of predefined VSL control strategies, the majority of these studies only have access to data from a limited range of traffic conditions and limited VSL values. erefore, most of the researchers merely performed qualitative comparisons between driving behavior parameters whether they used VSL control or not while missing the quantitative examination of the relationship between traffic variables and VSL values. In a realworld implementation, VSL values vary in accordance with the prevailing traffic conditions and weather conditions. However, few scholars evaluated the potential outcomes of each VSL value. In the meantime, only a limited number of dissertations controlled traffic states when examining VSL impacts; therefore, it is difficult to determine whether the obtained results are caused by VSL control or by the difference in traffic states within the collected empirical data.
is paper summarized the literature on traffic flow models and driving behavior under VSL control and proposed an improved cellular automaton model that represented individual driving behavior and reproduced corresponding traffic flow characteristics under VSL control. e rest of this paper is organized as follows. Driving behavior under VSL control is analyzed with data from a motorway in Section 2. Section 3 establishes the single cellular automaton model taking account of the compliance parameter. In Section 4, empirical data of a highway segment is used to identify the distribution of driver's individual radical degree and calibrate the parameters of the proposed model. In Section 5, numerical simulation is conducted to validate the superiority of the proposed model in reproducing traffic flow characteristics under VSL control.

Data Preparation.
e data used in this study was collected from a two-way European motorway segment with three lanes in each direction. e traffic flow state changes obviously, including recurrent and nonrecurrent congestions. e traffic flow information, such as the vehicular speed, the headway, and the vehicle length, was collected by loop detectors installed every 500-600 m in each lane. In order to alleviate congestion and improve driving safety, the expressway was equipped with VSL control facilities that adopt a control strategy based on mitigating traffic shocks. A wide variety of speed limits including 50 km/h, 60 km/h, 80 km/h, 100 km/h, and 120 km/h were implemented in the system. is research collected traffic flow data and corresponding VSL control speed data from 7:00 AM to 6:00 PM for two weeks at a nonweaving area site. Only data in good weather condition (no rain or fog) and visibility conditions (at daytime) was selected to avoid external interference. Missing and wrong data was removed. Finally, the dataset including 4,266 minutes of data samples with 355,599 vehicles was constituted. e traffic state is a crucial factor when investigating the VSL impacts on driving behavior. Under various traffic states, driving behavior will be differently affected by the VSL system. Hence, it is difficult to determine whether the obtained results are caused by the VSL control or by the difference in traffic states through the raw empirical data. To eliminate the impact of traffic states, the classification of the samples according to traffic states is indeed critical. In this paper, the collected samples were firstly classified into different traffic density intervals to measure traffic congestion. en the traffic data in the same density interval with different VSL control was compared to that without VSL control to investigate the effect of VSL control on driving behavior.

Desired Speed Distribution.
From the previous study, it can be inferred that the VSL has the ability to reduce the desired speed, since the average speed in each density interval under free-flow state was descended with different VSL values. Gartner's study mentioned that when the space headway is greater than 125 m (corresponding to the headway of 4∼6 s at 20∼30 m/s vehicle speed), the vehicle will not be in the following state [28]. According to the conclusion, vehicle speed in the above-mentioned state reflects the driver's desired speed, which is not affected by the preceding vehicle. A statistical analysis of the vehicle speed with a headway greater than 5 s in the free-flow state was performed. Figure 1 shows the vehicle speed distribution of the left lane in the uncontrolled state and VSL control state including 120 km/h, 100 km/h, and 80 km/h. Figure 1 demonstrates that the effect of the VSL on the driver's desired speed is holistic. After setting VSL, the driver basically responded and reduced the speed. e exponential regression line in the figure shows that the driving speed gradually increases as the headway increases, which is in line with reality. At the same time, it can be seen from the figure that under different VSL values, a few drivers drive at a speed obviously higher than the VSL value. It is believed that drivers who are significantly speeding are relatively aggressive drivers who may not follow VSL. Besides, at the high-speed limit, these drivers will be speeding by more than 30%, while the percentage will rise to 50% at the low VSL value. We believe that some drivers might challenge the rationality of VSL control when a low-speed limit is implemented, whereby the reduced VSL obedience rate leads to the phenomenon of higher mean speeds under lower VSL values.

Average Desired Speed.
e average desired speed reflects the compliance rate of drivers to some extent. In order to quantify the change in the average expected speed, the arithmetic average speed is calculated in the case of different speed limit values on different lanes. Table 1 shows the average desired speed with the headway time greater than 5 s under different VSL values in each lane ("-" indicates that the number of samples is less than 50; no statistics were performed). It can be observed that the driver's desired speed decreased significantly with VSL in all lanes compared with the uncontrolled state. In addition, from the speed limit value of 120 km/h to 60 km/h, the average speed decreases as the variable speed limit decreases. However, the reduction is not linear, since the difference between the driver's desired speed and the VSL value is bigger under the control of the lower VSL value like 60 km/h.

Small Headway Distribution.
e distribution of the headway, especially the proportion of the smaller headway, is critical to the stability of the traffic flow. is paper mainly focuses on the effect of the variable speed limit on the percentage of vehicles whose headway is less than 1 second. In order to quantify the effects of different VSL values on the headway, we calculated the percentage of small headways within different density intervals under different speed limits on different lanes. Table 2 summarizes the percentage of vehicles with smaller than 1 s headway under all combinations of density intervals and speed limits on the left, the middle, and the right lane ("-" indicates that less than 300 vehicles were observed within the associated density interval).
It can be observed from the table that the ratio of less than 1 s headway in the uncontrolled state increases with the increase of the density, while the ratio under the VSL control does not change significantly with the increase of the density. In general, the percentage of less than 1 s headway in the VSL control is obviously lower than that in the uncontrolled state. e percentage at the speed limit of 100 km/h and 80 km/h is the least, not even less than 50% of the uncontrolled state. e results show that some drivers become more cautious under VSL control and tend to maintain a large headway to ensure safety.

Single-Lane Cellular Automaton Model
is section presents a single-lane cellular automaton model for variable speed limit (VSL) control. On the basis of improved cellular automaton model incorporating impaired driver's radical feature (RF) raised by Qu et al. [29], this paper introduces the VSL compliance parameter to describe the driver's compliance characteristics, setting the desired speed according to the compliance rate of each vehicle, and then modify the moderation probability of the RF model combined with compliance rate; that is, if a driver with VSLs compliance and within the range of the VSLs control keeps headway less than 1 second, the driver will possess higher decelerating probability.
In the model, the driver will possess a higher decelerating probability while being within the range of the VSL control with headway less than 1 second. e model was established under two assumptions: (1) e VSL control changes the desired speed of most drivers, but the desired speed and the value of VSL are not equal and their relationship is not linear.
(2) e ratio of the smaller headway (less than 1 s) is significantly reduced under the VSL control [24], which reveals that some drivers become more cautious in driving and tend to maintain a larger headway to ensure safety.
e random probability of deceleration expression is Journal of Advanced Transportation if c n � 1 and a n (t) � 1 and t h < 1 p d : in all other cases, where c n is the compliance parameter of n th vehicle, c n � 1 means compliance, c n � 0 indicates noncompliance, the compliance of vehicle is affected by the driver's personal characteristics, and the compliance property of each vehicle is initially determined and will not change along with the simulation or varies after each update circle. e compliance parameter of each vehicle is determined by the following formula: c n � 1: if r and(0, 1) < C 0: else. C is the VSL compliance rate, defined as the percentage of the drivers affected by VSL control, abbreviated as "compliance rate" afterward. "r and (0,1)" indicates a random number between 0∼1. a n (t) is a parameter to evaluate whether the vehicle is within the VSL control range, a n (t) � 1 indicates that the n th vehicle is within the VSL control range at time t, while a n (t) � 0 indicates that it is out of the range. Other parameters are the same as the RF model defined.
e update rules of the proposed model are as follows: (1) Determining the desired speed as if c n � 1 and Xvsl i ≤ x n (t) < Xvsl i+1 , (2) Determining the probability of randomization as , t h , t s , c n , a n (t) , b n (t + 1) � 0, a n (t + 1) � 0.
(6) Location update is Determining whether it is in the range of VSL control as if Xvsl i ≤ x n (t + 1) < Xvsl i+1 and Vvsl i ≠ 0 , then: a n (t + 1) � 1, where Xvsl i is the starting position of the i th VSL control area. v denotes the average value of vehicle's desired speed (AVDS) under the VSL control, determined by the speed limit value in the i th VSL control area, which is called the AVDS of VSL, described as v � f(Vvsli). Vvsl i indicates the  On the basis of a number of samples, the average speed of the vehicle with headway larger than 5 seconds can be selected as the value of AVDS under a certain speed limit control. If it is not possible to obtain a sufficient sample size, the average speed in the lower density state may be selected as the value of AVDS. In this study, the AMS is equivalent to the AVDS under the condition of no VSL control, and the calculation method is the same as in calculating AVDS. For the speed limit values of 120 km/h, 100 km/h, 80 km/h, and 60 km/h, we selected the average speed of the vehicle with headway greater than 5 seconds as the AVDS of the VSL control. For the 50 km/h speed limit value, the average speed in the lower density state was selected as its AVDS. e specific values are shown in Table 3. e result of analyzing exiting data reveals that the formation of regression function is consistent with quadratic polynomial; in order to supplement the missing values to calculate the corresponding AVDS, the quadratic polynomial regression can be performed through the existing data, and the values are listed in Table 3. Herein, only the speed limit values, which can be integers divided by 10 and greater than 50 km/h, are listed. e regression function is modeled as (10), where x is cell speed value and y is the corresponding AVDS. e coefficient of determination is 0.97:

Value of Compliance
Rate. e value of the compliance rate is affected by many factors, among which the VSLs control traffic regulations and the driver's recognition are the most significant. In general, the compliance rate under the mandatory regulations is superior to recommended regulations and is more prominent along with the penalty rate rising or driver's recognition improvement.
ere are some drawbacks of analyzing compliance through a survey. Firstly, it depends on whether the driver answers the questionnaire honestly. What is more, driving behavior may be different between humans. erefore, in this paper, the statistics method was conducted. rough observing the distribution of vehicle speeds with its headway larger than 5 seconds, we found that some vehicles' speeds are significantly higher than others in any speed limit control. ese vehicles have noncompliance characteristics because they were not affected by VSL control cases.
It can be concluded that, among all vehicles with headway larger than 5 seconds, except for some vehicles with speed beyond the critical value, the percentage of the remaining vehicles is the compliance rate of VSL control. e critical speed value should be kept in a reasonable scope to decrease the errors in statistical analysis; in this study, the average speed over 30% of the limit value was selected. In the case of a 100 km/h limit control situation, according to statistics results, 5% of vehicles were beyond 130 km/h. As a consequence, the compliance rate is 95%, drafted as C � 95%.

Calibration of Random Probability of Deceleration p c .
e random probability of deceleration p c is to simulate the driver's cautious driving characteristics under the VSL control; its value will directly affect the proportion of the small headway.
erefore, we will calibrate the parameter p c by determining the ratio of the headway less than 1 second in the actual VSL control.
To calibrate the p c values corresponding to different speed limit values, we performed simulations under speed limit values of 120 km/h, 100 km/h, 80 km/h, 60 km/h, and 50 km/h, respectively. For a certain speed limit value, the value of the p c is gradually changed, and other parameters were fixed. e road length was set to 4000 cell units, the time step was set to 1 second, and the density was set to 15 veh/km, respectively. Periodic boundary conditions were adopted. e initial position and speed of the vehicle were randomly generated. Each simulation runs 10,600 time steps, where the first 10,000 time steps were not counted in order to eliminate the effects of transients, while in the next 600 time steps, the ratio of the headway less than 1 second was calculated. For the same p c value, the simulation was taken 10 times to obtain an average value. Figure 2 shows the ratio of the headway less than 1 second at different p c values under the 100 km/h speed limit control.
It can be seen from the figure that the ratio of the headway less than 1 second is relatively sensitive to the parameter p c . Qu et al. [24] mentioned that, in practice, the ratio of the headway less than 1 second under the control of 100 km/h speed limit in the free-flow state is about 15%∼20%; therefore, the value of p c ranges from approximately 0.5 to 0.8 under the limit speed value of 100 km/h. In the same way, we calibrate the p c values of other speed limit values. e results are shown in Table 4.
As can be seen from the results, the p c values corresponding to the excessively high-and low-speed limit values are small, while under the moderate speed limit, values (80 km/h, 100 km/h) are large. is reveals that the driver's cautious driving behavior under the control of the moderate speed limit is the most obvious. For the sake of simplicity, this study divides the value of the parameter p c into two parts: when the speed limit value is between 80 km/h and 100 km/h, p c � 0.6; under other speed limit values, p c � 0.4.

Numerical Simulation of the Model
In this section, the simulation of a fundamental diagram, time-space diagram, and speed differences characteristics was conducted to evaluate the performance of the singlelane cellular automaton model. In the simulation, the length of the road was set to 4000 cell units, and the time step was set to 1 second with periodic boundary condition. e initial position and speed of the vehicle were randomly generated. Each simulation ran 10,600 time steps, of which the first 10,000 time steps did no statistics to eliminate the impact of the transient. e value for the new model's parameters is the same as that in the calibration session.

Fundamental Diagram.
is study carried out a pseudonumerical simulation of VSL control with values of 120 km/ h, 100 km/h, and 80 km/h with no VSL control, which has been run 10 times for each density value; then the average traffic flow was calculated. e fundamental diagram under different speed limit values is shown in Figure 3. It can be seen from the figure that, in the low-density region, the average speed (curve slope) under the VSL control is minor compared to that of the uncontrolled situation, with the tendency of the speed limit value and the average speed being the same. In the high-density region, it can be found that the fundamental diagram at different speed limit values does not significantly differ from the uncontrolled situation. In the moderate density region, the critical density (the density region where the traffic flow changes phase) under the VSL control is larger than that in an uncontrolled situation, and the corresponding key density value is growing with the descent of the speed limit value. Accompanied by the decrease of speed limit value (80 km/h), the characteristics of the traffic flow phase transition have become vague. In addition, under the VSL control, the road capacity declined with its scope being weightier as the speed limit value decreased. e features of the fundamental diagram in the above simulations are basically consistent with those in actual observation. Figure 4 shows the time-space diagram of the 100 km/h speed limit and no VSL control cases at 0.1, 0.15, and 0.3 density, respectively. It can be seen from the figure that, in the small density (k � 0.1) and large density (k � 0.3) regions, the VSL control does not significantly change the basic form of the time-space diagram. In the small density region, the traffic flow is in the free-flow state as a form of the vehicle fleet, while the large density region exhibits a wide motion-blocking feature with a largescale blockage group. In the moderate density (k � 0.15) region, the VSL control significantly changes the shape of the time-space diagram, because it does not exhibit a phase separation phenomenon consisting of a distinct blocking phase and a moving phase compared to the uncontrolled situation, but possesses free-flow characteristics. It can be confirmed in Figure 4 that when k � 0.15, the free-flow state changes to the blocking phase without VSL control and maintains stability while being under control. is corresponds with a previous study reporting that VSL can increase the critical density [12].

Time-Space Diagram.
It can be found that, for a certain speed limit value, the spatiotemporal characteristics changed every time, but some     Table 5.
e result reveals that the speed limit value decreased as the density increased, and the probability of generating a free-flow state expanded. is indicates that the VSL control could boost the critical density of traffic flow, which is consistent with the fundamental diagram and the actual observation of traffic flow characteristics.

Average Speed Difference of Two Adjacent Vehicles.
e average speed difference (ASD) of two adjacent vehicles was selected to measure the speed difference characteristics of traffic flow. It is defined as the mean value of speed difference between the adjacent vehicles passing a given point at a specified time interval [31], as shown as follows: where v n represents the speed of the n th vehicle when passing the fixed point and N represents the number of vehicles passing a given point. In this study, simulations were carried out 10 times for each density condition following the singlelane cellular automaton model. rough simulations, the ASD of two adjacent vehicles in different speed limit values and different density conditions is shown in Figure 5.
It can be seen from the figure that the ASD of two adjacent vehicles in different VSL control cases is basically similar to the trend of the density varying. Except for excessive small density conditions, the ASD of two adjacent vehicles decreased with the density growth. In the small density region, the value under the VSL control is less than that in an uncontrolled state; this tendency is similar to actual situations [24]. In the large density region, different speed limit values are not significantly different.

Conclusions
is paper studies the real driving behavior data under VSL control and finds the following conclusions: (1) rough the processing and analysis of traffic big data, the traffic data in the same density interval with different VSL control was compared with that without VSL control to investigate the effect of VSL control on driving behavior (2) Under the control of the VSL, there always exists a certain proportion of drivers who fail to comply with the VSL control (3) e VSL control changes the driver's desired speed, but the average desire speed and the value of VSL are not equal and show a nonlinear relationship (4) e proportion of smaller headway time (less than 1 s) under VSL control is significantly reduced, which shows that some drivers become more cautious under the VSL control Based on the above-mentioned driving behavior, this study introduces the compliance rate parameter during the modeling process. In addition, it improves the slowing probability of the RF model [26], so that the vehicle slows down with a higher probability when the headway time is less than 1 s, in order to simulate the prudent driving behavior of the vehicle under the VSL control. In this study, the parameters of the model are calibrated using real vehicle data under VSL control. e numerical simulation results show that the model can not only fit the macroscopic traffic phenomena such as the fundamental diagrams and the time-space diagrams under the VSL control, but also fit the influence of the VSL on the speed dispersion.
It should be noted that this study assumed that the driver's radical feature is subject to the normal distribution when defining the function of radical feature parameter. Whether this assumption is reasonable or not needs to be discussed further in future research. Besides, this study is just early stage work to incorporate impaired driving behavior. More work should be done to develop a more accurate calibration process as well as multiple lane models. In our future research, we will focus on collecting real data to help discuss the driver's radical degree and its distribution under the car-following state.

Data Availability
e basic data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.