Speed-Density Model of Interrupted Traffic Flow Based on Coil Data

As a fundamental traffic diagram, the speed-density relationship can provide a solid foundation for traffic flow analysis and efficient traffic management. Because of the change in modern travel modes, the dramatic increase in the number of vehicles and traffic density, and the impact of traffic signals and other factors, vehicles change velocity frequently, which means that a speed-density model based on uninterrupted traffic flow is not suitable for interrupted traffic flow. Based on the coil data of urban roads in Wuhan, China, a new method which can accurately describe the speed-density relation of interrupted traffic flow is proposed for speed fluctuation characteristics. The model of upper and lower bounds of critical values obtained by fitting the data of the coils on urban roads can accurately and intuitively describe the state of urban road traffic, and the physical meaning of each parameter plays an important role in the prediction and analysis of such traffic.


Introduction
Flow, speed, and density are known as the basic elements of traffic flow theory.Flow can measure the number of vehicles and the demand for traffic infrastructure.Speed is an important control index in road planning, and it is also an evaluation index of vehicle operation efficiency.Density reflects the intensity of the vehicles on the road and determines traffic management and control measures.The relationships between flow, speed, and density called fundamental diagrams play a very important role in traffic flow theory and traffic engineering.For example, the speed-flow relationship can be used in highway capacity analysis in order to determine the highway service quality, and the speed-density relationship can reflect dynamic change in traffic flow, which can be used to study the disturbance propagation between vehicles.Therefore, sound mathematical models provide a solid foundation for traffic flow analysis and efficient traffic management.The relationship between speed and density which can reflect the quality of service received from the road is attracting considerable research attention.
The earliest speed-density model was a linear model proposed by Greenshields et al. [1] in 1935.The linear model overlaps and classifies the observed data groups, which is proved to be unreasonable, and observation time is a holiday, with a narrow range of representations, so there are some deviations between the derived speed-density relation and the actual situation.Later, the relationship between speed and density was studied in greater depth, and the Greenberg logarithmic model, Edie model, Underwood exponent model, Pipes-Munjal model, modified Greenshields model, Newell model, and so forth, emerged in turn [2,3].Heydecker and Addison [4] studied the relationship between speed and density under various speed limits and found that zero speed induces traffic jams, not the other way around.Ma et al. [5] derived a general logistic model of traffic flow characteristics, which includes several traffic flow parameters with clear physical meanings and analyzed the effects of the parameters on speed-density logistic curves.The experimental results showed that this model can well describe the traffic flow characteristics in different states.Shao et al. [6] proposed a speed-density model 2 Mobile Information Systems under congested traffic conditions combined with the minimum safety spacing constraint, and the experimental results showed that the absolute error of this model was smaller than that of other models fitting the traffic data of two freeways.Wang et al. [7] proposed a family of speed-density models with different numbers of parameters with important physical significance and got good performance in the final experiment.
All of the above studies are based on continuous traffic flow data.These data, also called uninterrupted traffic flow, are traffic flow with no effect of external fixation factors, such as freeway, urban expressway, and so forth.Discontinuous traffic flow, referred to as interrupted traffic flow, is periodically influenced by external fixation factors.The most common interrupted traffic flow is originated by signal lamps of urban intersections.Because of the variety of vehicle types, the periodic effect of signal lamps, shunts in the canal section, and other factors, the characteristics of interrupted traffic flow are very complex compared with uninterrupted traffic flow.In addition, the city is still in a rapid increase in population and, with the development of economy, people are more inclined to self-driving travel, thus more and more vehicles and more and more congestion in the city, which leads to the increase of travel time, the growth of fuel consumption [8], the aggravation of environmental pollution, and other awful issues [9,10].Compared with the highway, the urban road has a strong influence on the individual, society, and the environment.Therefore, further study of the characteristics of interrupted traffic flow to provide support for management decisions is particularly important.
Research on interrupted traffic flow has attracted a lot of attention [11][12][13][14][15].Many scholars see traffic flow located at a certain distance from the intersection as continuous traffic flow, believing that it can be described by continuous traffic flow models.Some of the literature [16,17] suggests, however, that because of the short distance between intersections in the city and the influence of signal lamps, there are differences between traffic flow located at a certain distance from the intersection and the traffic flow of freeways.Because traffic data are difficult to obtain and for other objective reasons, only a few scholars focus on the speed-density model of discontinuous traffic flow.Wang et al. [18] introduced a fourparameter logit model for complete data fitting and established a speed-density logit model for left-turning, straight, and right-turning traffic flow.However, the experimental data were obtained by VISSIM simulation, and the simulation parameters were not accurate enough to depict the complex city road environment, so the experimental results have certain limitations.Wang et al. [19] thought that the stochastic model would contain more traffic information and put forward the stochastic speed-density model.This stochastic model can generate a probabilistic traffic flow model and can achieve real-time traffic prediction.
In order to provide favorable data analysis and presentation for city traffic, thus to provide decision support for intelligent transportation, characterizing the speed-density relationship of interrupted traffic flow more accurately is full of importance.By analyzing a large amount of data, we propose a description method for a speed-density relationship model which is suitable for discontinuous traffic flow, using the upper and lower curves to describe the upper and lower bounds of velocity values.Because of the discrepant characteristics of the traffic flow in the outer and inner lanes, the coil data of the outer and inner lanes are analyzed and verified.

Speed-Density Model
Three basic parameters (flow , speed , and density ) are the core content of the traffic flow model.The three have the following relationship: that is, flow is the product of density and speed.The relationship between two parameters of the three is of great significance in traffic flow, and the relationship between speed and density has received a lot of research attention.Greenshields et al. was an early researcher, who proposed the speed-density linear relationship [1]: where   is the speed of free flow, that is, the speed of vehicles unimpeded when the traffic density tends to zero, and   is the density of block flow, that is, the density when the traffic flow is blocked and cannot move.As shown in Figure 1, when  = 0, the speed can reach the theoretical maximum value, namely, the free flow velocity   .The area surrounded by the abscissa, the ordinate of any point on the line, and the coordinate origin is the traffic flow.Equation (2) can change to Respectively, introduce (2) and ( 3) into (1), and we get Equations ( 4) illustrate that - and - are quadratic function relations, as shown in Figure 1.
The linear model is too simple, and there are many deficiencies.In order to improve the model, scholars have proposed models based on the linear model but with a higher degree of accuracy.simulation software to set up and change six parameters of road traffic, including section length , stretch section length , cart rate , signal period , the ratio of the time span of left-turn green signal to signal period   , and the ratio of the time span of the straight green signal to signal period   , and established 22 groups of parameters.The simulation results showed that the relationship between speed and density presents an inverse S curve.Therefore, a four-parameter logit model is proposed here to describe the speed-density inverse S curve, and its expression is as follows: where  min is the mean value of the minimum speed,  max is the mean value of the maximum speed,  is the flow value of a section,   is the flow value at the inflection point of the curve, and  is a parameter determining curve shape.Then, the data obtained from the 22 groups of simulation parameters were fitted.The four parameters ( min ,  max ,   , and ) were calculated for each simulation environment.  and  were, respectively, fitted in left-turn, straight, and rightturn cases, and the fitting results are as follows:

The Description Method of the Speed-Density Model for Interrupted Traffic Flow
3.1.The Characteristics of the Data of Interrupted Traffic.Coil data for one day, three days, seven days, and fourteen days were selected to compare and analyze the discontinuous flow data and the existing six speed-density models, as shown in Figure 2. We found the following: (1) The six models' performance was poor when the coil data of interrupted traffic flow were fitted, illustrating that although suitable for uninterrupted traffic flow (2) The interval value of critical densities   of one-day, three-day, seven-day, and fourteen-day data sets was [62.56 pcu/km, 71.23 pcu/km], and most of the data were located in  <   range, meaning unimpeded flow data accounted for the absolute proportion, so the traffic flow of the location coil was in a state of flow most of the time.(3) When  <   , with the increase of density, the velocity decreased sharply; when  >   , as the density increased, the velocity decreased slowly, and the speed variation amplitude was very small.(4) When the density was small, the speed had a large range of values, of which the largest was [23 km/h, 72 km/h].We filtered out the small-density data to obtain a scatter diagram of flow and occupancy which were directly collected by a loop detector, as shown in Figure 3.In Figure 3 it is obvious that the loop detector acquires large-range flow values for the same occupancy value, and the largest range can reach 100 pcu/h.Hence, after calculating speed and density by density formula and velocity formula, speed accordingly has a large range of values for the same density in speed-density diagram.(5) In addition, the density values were found to be near a number of points, and the difference between adjacent points was approximately equal to a certain value.
From the above analysis, we found that, because of the big differences between uninterrupted and interrupted traffic flow, existing models suitable for uninterrupted traffic flow are unsuited for describing the speed-density relation of interrupted traffic flow.What is more, the flow collected by a loop detector has a large range of values.Therefore, for the speed-density relation of interrupted traffic flow, we must find a new descriptive method.

Description Method of Speed-Density Relationship for
Interrupted Traffic Flow.Because of the difference between the uninterrupted and interrupted traffic flow and the volatility of speed, the speed-density relationship cannot be adequately described by a single model, so we use two curves,  upper and  lower , to describe the supremum and infimum of velocity values: where  upper and  lower are, respectively, the upper and lower bounds of velocity and  upper and  lower are fitting functions.
Divide the density interval [ min ,  max ] into  connected intervals  1 ,  2 , . . .,   .Partition data  as  1 ,  2 , . . .,   by density intervals, and correspondingly get speed sets  1 ,  2 , . . .,   , causing that, for any  ∈ (1, 2, . . ., ), we have where (  ) is used test for   with the Shapiro-Wilk normal test method.Sort  independent observations in   by nondescending order, recorded as  1 ,  2 , . . .,   , and construct the -test statistic where   is the coefficient when sample size is .When the population distribution is normal distribution, the value of  should be close to one. quantile   of statistic  can be obtained by the look-up table method.When  ≤   , the original hypothesis should be rejected at the significant level, indicating that   does not obey normal distribution; when  >   , the original hypothesis cannot be rejected, and   satisfies normal distribution.Under the conditions of ( 8), for every  ∈ (1, 2, . . ., ), extract the upper quantile Fit  upper and  lower using the nonlinear least square method.The tabulated function   = (  ),  = 1, 2, . . ., , is Mobile Information Systems available by (10).Then we need to obtain the fitting function, () =  0 +  1 ×  1 () + ⋅ ⋅ ⋅ +   ×   (), making the sum of squared deviations Take the minimum, of which  1 (),  2 (), . . .,   () are  nonmergeable monomials of variable , and  0 ,  1 , . . .,   are the coefficients of monomials. is a nonnegative polynomial of  0 ,  1 , . . .,   , so there must be a minimum value.Respectively, calculate partial derivatives of  for  0 ,  1 , . . .,   , and make them equal to zero.

Experiment and Analysis
The experiment data is collected by the coil detectors underground closed to Optical Valley Walking Street in Wuhan, China.Coil detectors collect data every 15 minutes, recording time, flow, occupancy, and so forth, as shown in Table 2.
Use the method in [20,21] to calculate speed and density, and the ratio of the amount of data between two model curves to the total amount of experiment data is used to describe the performance of model.The loop detector in the outer lane measures the traffic flow of straight and right-turning lanes, and the loop detector in the inner lane measures the traffic flow of the left-turning lane.The traffic flow characteristics of two loop detectors must have certain differences.Therefore, analyze the coil data of both the outer lane and the inner lane to find the diversity of their speed-density relationship.

Coil Data Analysis of the Outer Lane.
The experimental steps are as follows.
Step 1. Analyze coil data of the outer lane and find that density values are clustered at a number of points  1 ,  2 , . . .,   , where the mean value of the difference between the adjacent points is about 2.5 pcu/km.Divide density  into a number of intervals with length 2.5 pcu/km by  1 ,  2 , . . .,   .
Step 3. Execute a distribution test for   where the result shows that one data set is too small to meet the requirements of the test.Merge the adjacent density segments in Step 2 to enlarge the amount of the small data set.Redo the distribution test for the new data set, more than 80% of which meets the normal distribution, with totally 95% of the total data satisfying the normal distribution, which makes it reasonable to consider all the small data set satisfying the normal distribution. Step Figure 4 shows the validation result of the speed-density logarithmic model of the outer lane when upper value = 0.95 and lower value = 0.05.Equations ( 18) correspondingly are the green and blue curves in Figure 4, which is the speed-density model of interrupted traffic flow created by the new description method.Significant test results indicate that  values of two regression coefficients of two curves are minima ( < 2 − 16), which means that coefficients are significant and two log models constructed with density as the independent variable are applied to estimate velocity as the dependent variable.
The coil data of the outer lane for two weeks, four weeks, six weeks, and eight weeks are, respectively, selected  3 gives the ratio of the data between two logarithmic curves to the total amount of data in each case.Make a longitudinal observation; it is obvious that, with upper value increasing and lower value decreasing, the proportion increases accordingly, where amplitudes are obvious, respectively, 7.2%, 6.2%, and 6.9%.On the other hand, the main transverse trend is that the proportion increases along with the increase of experiment data loosely, where, however, sixweek data has the best performance.The above suggests that the two logarithmic models are able to describe the speeddensity relation of the outer lane.Figure 5 shows the four groups' validation results when upper value = 0.95 and lower value = 0.05.

Coil Data Analysis of the Inner Lane.
We select coil data of the inner lane and follow Steps 1 to 5 as for the outer lane.When fitting sets  upper and  lower at Step 6, we find that the speed-density models proposed by scholars all have poor performance with goodness of fit of less than 0.5, which suggests that a single model cannot accurately describe the quantile set of the coil data.Thus we consider using a segmentation model.
In the density-flow curve there is a critical density   , which is the density of maximum traffic flow, as shown in Figure 1.When the density  <   , the traffic is in a state of flow; when  >   , the traffic flow gradually becomes crowded.Therefore, consider using   as the critical value of the subsection. A Figure 6 shows the fitting result of a segmentation model of the outer lane when upper value = 0.95 and lower value = 0.05, and ( 19) are the models corresponding to the green curve and blue curve in Figure 6, which is the speed-density model of interrupted traffic flow via the new description method. value of each parameter is very small, suggesting the coefficient is very significant.
The coil data of the inner lane for two weeks, four weeks, six weeks, and eight weeks are, respectively, selected and four groups of parameters are established for the model validation, the same as that for the outer lane.Table 4 gives the ratio of the data between two logarithmic curves to the total amount of data in each case.Comparing the result with that of the outer lane, we find that the validation results of the model of the inner lane are better with greater ratio.Take a longitudinal observation; similarly, it is obvious that with upper value increasing and lower value decreasing, the proportion increases accordingly, where amplitudes are smaller than that of outer lane, respectively 3.7%, 1.9%, and 4.3%.The main transverse trend is the same as outer lane except the case of upper value = 0.95 and lower value = 0.05.The result indicates that the two segmentation models are suitable for describing the speed-density relation of the inner lane.Figure 7 shows the four groups' validation results when upper value = 0.95 and lower value = 0.05.4.3.1.Difference between the Models of the Outer Lane and the Inner Lane.The loop detector of the outer lane measures right-turning and straight lanes, and the coil is located in a road adjacent to a commercial pedestrian street with a heavy flow of people and traffic.A logarithmic model is applied to describe traffic flow with large density, and therefore it is accepted that the coil data of the outer lane satisfy the logarithmic model.

Experimental Result Analysis
The loop detector of the inner lane measures the leftturning lane which also has heavy traffic flow.The speeddensity relation of the inner lane does not satisfy the single log model but is suitable for the segmentation model.The logarithmic model and the inner lane meets the segment model.What is more, the road segment where loop detectors located is unblocked at most of the time; the inner lane is more unimpeded than the outer lane.

Conclusion
In this paper, the characteristics of urban interrupted flow data were analyzed, and it was found that they differ from the data of uninterrupted flow.Since the existing classical models cannot describe them very well, a description method of speed-density relation for interrupted traffic flow was proposed where the upper and lower curves were used as the upper and lower bounds of the predicted speed.In this method, the speed was divided into small data sets which satisfied the normal distribution, and two quantiles of normal distribution were obtained as the predicted values.Then two quantile sets were fitted to get two curves as the speed-density relation model of the interrupted traffic flow.Finally, the coil data of the outer and inner lanes were applied for model validation.The results showed that the new method can give a good description of the speed-density relationship of interrupted traffic flow and get different model results for the outer lane and inner lane, whereby the speed-density relation of the outer lane satisfies the logarithmic model and the inner lane satisfies the segment model instead of the single model, where when the density is less than critical density, it conforms to the exponential model and otherwise the logarithmic model.The fitting results of the internal and external lanes were analyzed in combination with the actual local road environment and traffic flow theory.So this model can provide favorable data analysis and presentation for city traffic, thus to provide decision support for intelligent transportation.

Figure 2 :
Figure 2: Comparison of six speed-density models.

Figure 4 :
Figure 4: Speed-density logarithmic model of the outside lane.

Figure 6 :
Figure 6: Speed-density multisession model of the inside lane.

Table 1 :
Table 1 lists results for the speed-density model, including the Greenberg model, Underwood model, Northwestern model, Newell's model, Pipes-Munjal model, Drew model, Modified Greenshields model, Del Castillo and Benitez model, Van Aerde model, MacNicholas model.These models with the parameters of important physical meaning provide good results.Speed-density models.
Wang et al. [19] established a speed-density logit probability model with four parameters.Wang et al. used VISSIM

Table 2 :
The data example.
4. Get two quantiles    as upper and lower critical values of velocity for density   .(fit  upper and  lower ).Because the loop detector is located near commercial street which has heavy traffic, we use the logarithmic model to formulize the data.

Table 4 :
Validation results of the model of the inner lane.
density-flow curve is obtained by local polynomial regression fitting, and the density value at the curve vertex is just   .Take   as the critical value and piecewise analyze