A Simplified Network Model for Travel Time Reliability Analysis in a Road Network

This paper proposes a simplified network model which analyzes travel time reliability in a road network. A risk-averse driver is assumed in the simplified model. The risk-averse driver chooses a path by taking into account both a path travel time variance and a mean path travel time. The uncertainty addressed in this model is that of traffic flows (i.e., stochastic demand flows). In the simplified network model, the path travel time variance is not calculated by considering all travel time covariance between two links in the network. The path travel time variance is calculated by considering all travel time covariance between two adjacent links in the network. Numerical experiments are carried out to illustrate the applicability and validity of the proposed model. The experiments introduce the path choice behavior of a risk-neutral driver and several types of risk-averse drivers. It is shown that the mean link flows calculated by introducing the risk-neutral driver differ as a whole from those calculated by introducing several types of risk-averse drivers. It is also shown that the mean link flows calculated by the simplified network model are almost the same as the flows calculated by using the exact path travel time variance.


Introduction
Conventional frameworks for analyzing and modeling transportation systems have been confined to average representations of the network state (e.g., average link flow or average travel flow).For instance, in the traditional traffic assignment model, one can obtain a deterministic prediction of a future flow on a certain link in the network based on average origin-destination (O-D) flows, link capacities, and a form of proportional path choice model (either deterministic user equilibrium (DUE) or stochastic user equilibrium (SUE)).This represents a deterministic view of the environment and the modeler's postulation that the variability or uncertainty in the system is not influential in system design and evaluation.
There has been a growing concern over the uncertainty of travel time in transport systems and its effect on the reliability of transport services [1].Research on network reliability has begun to address this problem [2][3][4][5].From the traveler's perspective, the issue of travel time reliability has been a major concern.Travelers may experience excessive variability of travel time from day to day on the same trips [6][7][8].
In transport modeling, some advances have been made toward incorporating traffic flow uncertainties into the network modeling framework (i.e., developing a stochastic network model).Watling [9] proposed a second-order network equilibrium model that explicitly considers random path choice behavior.His model, in contrast to the conventional SUE model, uses path choice probability, as predicted by a SUE model, to define stochastic path flows that follow a multinomial distribution.The path flows derived using the traditional SUE model are, in fact, the expected flows of this multinomial distribution.In this model, the drivers choose their paths so as to minimize their perceived long-run expected travel costs.With a nonlinear travel cost function, this long-run expected travel cost will differ from the equilibrium cost computed by the conventional SUE model [10].
Clark and Watling [6] extended this stochastic network model to the case with a Poisson distribution of O-D flows.Similarly, Nakayama and Takayama [11] proposed a stochastic network model with random path choice behavior but using a binomial distribution of O-D flows.These models fully represent stochastic/uncertain path choice behavior travel time variance as a measure of travel time variability in this study.Then, link and path travel time under stochastic demand flows are formulated in Section 3. In Section 4, two network equilibrium models under stochastic demand flows are formulated considering a risk-averse driver's path choice behavior in a road network.The first model introduces path travel time variance which is calculated considering all travel time covariances between two links in the network.The generalized travel time in this model is not additive since the generalized path travel time is not equal to the sum of the generalized link travel times related to that path.The second model, which we propose in this study, is a simplified version of the first model.In this model, the path travel time variance is not calculated by considering all travel time covariance between two links in the network.The path travel time variance is calculated by considering all travel time covariance between two adjacent links in the network.The generalized path travel time in this model is additive.It is shown that a unique solution is provided by the simplified network model.Numerical experiments are carried out to illustrate the applicability and validity of the proposed model.Finally, concluding remarks are provided in Section 6.

Measures of Travel Time Variability
We will briefly review how to obtain a measure of travel time variability based on the expected utility maximization principle.Vickrey [21] considered a separable or additive utility function which is a sum of utilities obtained from time spent at an origin and time spent at a destination of a trip.Using such formulation of utility, it is possible to consider a driver who chooses a departure time optimally in order to maximize expected utility when facing uncertain travel time.Noland and Small [22], Bates et al. [23], Fosgerau and Karlström [24], Fosgerau and Engelson [25], and Engelson [26] have shown how measures of travel time variability can be derived from the drivers' scheduling preferences.
A popular formulation of scheduling preferences is the -- preference in which the marginal utility of time (MUT) at the origin is constant and that at the destination is a step function [27,28].By assuming an exponential travel time distribution or a uniform travel time distribution, Noland and Small [22] derived the scheduling utility which is linear in (, ), where  and  are the mean and standard deviation (SD) of the stochastic travel time.Fosgerau and Karlström [24] generalized this result to any travel time distributions.Fosgerau and Engelson [25] considered the value of travel time reliability under scheduling preferences that were defined in terms of linear MUTs being at the origin and at the destination.They found that the scheduling utility was linear in (,  2 ,  2 ) and that this result was independent of the shape of a travel time distribution.
Engelson [26] derived the scheduling utility for the two cases when the MUTs at both the origin and the destination are either quadratic or exponential in form, and demonstrated special cases when the scheduling utility is additive.The necessary condition when the scheduling utility is additive is that the MUT at the origin is a positive constant.
Engelson and Fosgerau [29] derived a measure of travel time variability for travelers equipped with scheduling preferences defined in terms of MUT and who chose departure time optimality.In the case of ℎ() = ℎ 0 and () = ℎ 0 + (/) ⋅ ( ⋅(− 0 ) − 1) which are, respectively, MUT at the origin and that at the destination evaluated at clock time  (the parameters in the MUTs satisfy the following conditions: ℎ 0 > 0,  ≥ 0,  > 0, and ℎ( 0 ) = ( 0 )), the scheduling utility is where  is an independently distributed random travel time with a mean value of  and an SD of  and  0 is the time such that the individual prefers being at the origin before this time and at the destination after this time.From the first-order condition, the scheduling utility is derived as where ln [ ⋅(−) ] is the cumulant-generating function (CGF) of the centralized travel time distribution.A limiting case of  → 0 yields lim If the travel time distribution has compact support, then the scheduling utility is finite for any  and can be presented as the convergent Taylor series where   is the cumulant of order  of the travel time distribution.If the travel time follows a normal distribution, however, the normal distribution does not have compact support, and for any  the scheduling utility is A special case of this problem with constant ℎ() ≡  and the two-valued function () =  −  for  <  0 and () =  +  for  ≥  0 leads to the -- preferences model.It was reported that (2) and ( 3) have advantages over the -- preferences model.First, they do not depend on the shape of the travel time distribution.The second is additivity with respect to parts of trip with independent travel time, which is an important property to analyze travel time reliability in a road network.Travel time variance as a measure of travel time variability may be criticized for not taking into account the skewness of the travel time distribution [30].The CGF depends on the skewness ( 3 / 3 = [( − ) 3 ]/ 3 ) of the travel time distribution for nonzero .However, the travel time covariance between two links is not taken into account in a CGF in which independent link travel time is assumed.The effect of link travel time covariance terms on the path travel time variance becomes larger than that of the skewness as the number of links in a path increases.For example, if a path is comprised of  links, the number of link travel time covariance terms taken into account in calculating the path travel time variance is   2 .We recognize that the travel time covariance between two links is a more important factor in analyzing travel time reliability in a road network than the skewness of the travel time distribution.Hjorth et al. [31] analyzed the stated preference data by applying the scheduling preferences model that assumes MUTs at the origin and at the destination.They have shown that the value of travel time variability can be proportional to the variance of travel time.This result can partially support the use of travel time variance as a measure of travel time variability.
The additivity of the scheduling utility is a convenient property for a network equilibrium model from the practical viewpoint (calculation cost efficiency, ease of handling the network equilibrium model, etc.).If we employ an SD related measure of the travel time variability (SD, percentile value, etc.), unrealistic drivers' path choice behavior as shown next may be generated.
The following example is cited from Cominetti and Torrico [32].The generalized travel time of a random travel time  is given by where we assume ω = 1 without loss of generality.We consider then the traffic situation shown in Figure 1 in which a road network consists of three nodes and three links, and the stochastic link travel time is shown.In the network, link travel time is denoted by   ( = 1, . . ., 3), where  is link number, which follows the normal distribution (,  2 ) with a mean of  and a variance of  2 .From the link travel times shown in the figure, we obtain ( 1 ) = 12, ( 2 ) = 10 + √ 5, ( 1 +  3 ) = 21 + √ 3, and ( 2 +  3 ) = 20 + √ 7. From these four generalized travel times, the following results can be obtained.The minimum generalized travel time between nodes 1 and 2 is 12, and the minimum path consists of link 1.
The minimum generalized travel time between nodes 1 and 3 is min(( 1 +  3 ) ( 2 +  3 )) = 20 + √ 7, and the minimum path consists of links 2 and 3.This example shows that a driver in a network with an origin node 1 and a destination node 2 will choose the path comprising link 1 if he/she prefers smaller generalized travel time.However, the same driver in the network whose origin and destination nodes are, respectively, 1 and 3 will choose the path comprising two links 2 and 3.However, the path choice criterion is clear, and the path choice behavior of the driver is somehow unrealistic.Even though the mean deviation is used instead of SD, such an unrealistic case can occur since the generalized path travel time is not equal to the sum of the generalized travel times of the links that comprise that path in the case of mean deviation.
We examined the convenience and importance of the additivity of the generalized path travel time when addressing it in a network problem.In the following, we formulate some network equilibrium models in which travel time variance is employed as a measure of the travel time variability in reference to (5).According to (5), the generalized travel time is given by where

Link and Path Travel Times under Stochastic Flows
The variance of   is given by The conservation of the path flow variance in relation to the O-D flow variance holds (Appendix A).The stochastic flow of link ,   , is given by The mean and covariance of the stochastic link flow are then cov where ∑ ∈ ∑ ∈    ⋅   ⋅   is the sum of all stochastic path flows that pass through both links  and .

Stochastic Link Travel Time and Stochastic Path Travel
Time.In this study, link travel time is represented by the following BPR function [34]: where  0  is free flow travel time of link  and (≥ 0) and (≥ 1) are calibration parameters.By substituting V  in (15) with   , we obtain where κ =  0  ⋅ /(  )  .Next, we will show how to calculate both a mean value and variance of the stochastic link travel time shown by (16).By performing an th-order ( ≥ 1) Taylor expansion to (16) at   = V  , we obtain where   is the coefficient of the th term of the Taylor expansion given by The mean link travel time is then calculated as The travel time covariance between two links is As shown in Clark and Watling [6], ( 19) and ( 20) can be calculated by applying a method proposed by Isserlis [35] given the moments of   by assuming that the link flow follows normal distribution [17,18,33].This assumption was supported by Rakha et al. [36], which demonstrated that the normality assumption may be sufficient from a practical standpoint given its computational simplicity.Equation ( 19) can be calculated as follows: The results of  = 4 are provided in Appendix B. We now assume that the coefficient of each O-D flow takes a specific value.By applying this assumption to (14) where V  in ( 23 By substituting ( 22) into ( 19) and ( 20), we obtain cov where in (25) is the coefficient for the th term.The results of  = 4 are provided in Appendix D. Note that it is shown from ( 24) and ( 25) that the mean, variance, or covariance of link travel time is expressed by using only mean link flow(s) with some given parameters, and it will be shown that both [  (  )] and cov[  (  ),   (  )] are increasing functions with respect to V  and V  , respectively.It will be shown that these two mathematical properties are convenient for developing a network equilibrium model.The most dominant reason for these two properties is that the coefficients of variation of all O-D flows are assumed to be the same.Thanks to this assumption, any moments for the stochastic link flow can be calculated by using its mean value.Also, as far as the Taylor series expansion provides good approximation, the method presented in this study can be applied to any functional forms.Since the Taylor series expansion can approximate well the function of () = 1/(1 − ) for 0 <  ≤ 1, the proposed method can be applied to the Davidson type link cost function.From (25), cov[  (  ),   (  )] > 0 if and only if V  > 0. From (23), if V  > 0, then V  > 0 and V  > 0; however, the inverse relationship does not always hold for any two links in the network (i.e., even though V  > 0 and V  > 0, V  can be zero).These two mathematical properties show that the travel time covariance of two links, cov[  (  ),   (  )], is greater than zero if and only if V  is greater than zero and that V  and V  can influence the travel time covariance of two links, cov[  (  ),   (  )], if and only if V  is greater than zero.Therefore, a calculation of V  is important to calculate cov[  (  ),   (  )] in the network.As discussed in the next section, the calculation of V  for two adjacent links in the network is easily implemented since there is no need to enumerate a path set in the network.In contrast, the calculation of V  for two unconnected links in the network becomes more difficult than that for two adjacent links since that may need to enumerate a path set in the network.
For notational simplicity, in the rest of the paper, [  (  )], var[  (  )], and cov[  (  ),   (  )] are denoted by t (V  ),  2   (V  ), and   (V  ; V  , V  ), respectively.The travel time of path  which serves O-D pair  (Ξ  ) is given by The mean path travel time and path travel time variance are, respectively, given by The path travel time covariance is For calculation methods of the mean travel time and travel time variance when each link flow in the network follows a lognormal distribution, the reader is referred to Tani and Uchida [37] in which each link capacity in the network is also assumed to follow a lognormal distribution.

Network Equilibrium Model under Stochastic Flows
where  ≥ 0 is a relative weight assigned to var[Ξ  ].The risk-averse driver assumed in this study chooses the path with lower path travel time variance if the mean path travel times of all the alternative paths are the same.Such risk-averse driver's path choice problem based on the DUE principle can be formulated as follows: subject to (10), (13), and (23), where   is the minimum generalized travel time of O-D pair .The superscript * is used to denote the variables that are obtained at equilibrium.It is known that this problem is equivalent to the following nonlinear complementary problem (NCP): subject to (10), (13), and (23).The equilibrium path flows can be obtained by solving the following variational inequality (VI) problem [38]. where There are efficient solution algorithms for solving the DUE traffic assignment problem.However, most of such algorithms cannot be applied to solve the VI problem shown above in which the path travel time variance is nonadditive due to the link travel time covariance between two links in the network.Therefore, path-based solution algorithms [39] which, in general, require enumeration of a path set need to be applied in order to solve the VI problem.However, the enumeration of all possible paths is almost impossible for the case of a large network.Therefore, in the next section, we will propose a simplified network model in which only the covariance terms between two adjacent links in the network are taken into account in calculating the path travel variance by considering practicality.Thus, an efficient linkbased algorithm for DUE traffic assignment problem can be applied to the simplified network model.

Simplification.
Paths enumeration may be required for solving the VI problem presented in the previous section.Enumeration of all noncyclic paths in a large road network is impossible.From a practical standpoint, it may be reasonable to enumerate several paths for each O-D pair (e.g., two or three paths for each O-D pair).However, different solutions can be obtained depending on paths enumerated, and that may be a troublesome issue in estimating the benefit of a project.
If we assume that the link travel time follows an independent distribution, then path travel time variance is the sum of the link travel time variance related to that path.In this case, the path choice problem can be formulated as the following convex programming problem: subject to (10), (13), and (23), where . This problem has the same mathematical structure as the standard DUE traffic assignment model and thus can be solved easily by applying standard link-based algorithms (MSA (Method of Successive Averages), the Frank-Wolfe algorithm, etc.) [40].However, ignoring all travel time covariance in (28) may bring about unrealistic solutions.
Rakha et al. [36] presented extensive evidence of a significant correlation between travel times of two adjacent links in the network.Although they analyzed travel time variability over vehicles, this evidence may support travel time variability over days.By utilizing this evidence, we now take into account all travel time covariance between two adjacent links in the network when calculating the path travel time variance.The corresponding path choice problem can be then formulated as the following link-based VI problem: simplified network model (SNM). where () is the set of links which are adjacent to link  in front of it.SNM includes no stochastic variable, although SNM analyzes travel time reliability in the network.
To solve SNM, one can apply a network representation which may be used when addressing intersection delays, in which dummy links connecting between all adjacent links in the network are added to the original network.Consider an original network that consists of a set of links and a set of nodes.It is assumed that each node in the network also has an identical number.Consider then a directed link in the original network.We can find each link of which origin node number is the same as the destination node of the link.By using this relationship, we can construct the augmented network (e.g., right-hand side of Figure 2), that corresponds to the original network (e.g., left-hand side of Figure 2), by inserting a directed dummy link between these two links.In the augmented network, the generalized travel time of link  ∀ ∈  is   (V  ) and that of dummy link  ∀ ∈ , ∀ ∈ () is   (V  , V  , V  ).Once the augmented network is constructed, V  is easily calculated as the flow of link .By applying this network representation, SNM can be solved by a diagonalization method in which V  and V  in   (V  , V  , V  ) are regarded as constant terms [41][42][43][44][45]. Also, MSA can be applied for solving SNM.Due to the mathematical structure of   (V  , V  , V  ), however, SNM has the same mathematical property as an asymmetric DUE traffic assignment problem which can have multiple solutions.Therefore, SNM may have multiple solutions.Next, we will examine the uniqueness of the solution of SNM.
If both generalized link travel time functions,   (V  ) ∀ ∈  and   (V  , V  , V  ) ∀ ∈ , ∀ ∈ (), in the augmented network are monotone functions, SNM has a unique solution.Or equivalently, if the Jacobian matrix ∇g(k, k), where g = ((  ) ∈ , (  ) ∈,∈() ), is positive-definite, SNM has a unique solution.The Jacobian matrix ∇g(k, k) is given by the following lower triangle matrix: If every eigenvalue of a Jacobian matrix is positive, the Jacobian matrix is positive-definite.The eigenvalues of a lower triangle matrix are equal to the values of diagonal elements of the matrix.In SNM, the following two conditions hold for each link in the augmented network (see Appendix E for the proofs): Therefore, ∇g(k, k) is positive-definite and thus SNM has a unique solution.If the link travel time covariance for any two links in the network is taken into consideration in calculating path travel time variance, the uniqueness of the solution is not obtained.
Since SNM calculates the stochastic variables at the equilibrium, the stability as well as the uniqueness of the solution is guaranteed in SNM.Even in a problem in a static context, the stability and uniqueness of the solution can be analyzed by applying theory in a dynamic context.This is true for our static problem.

Settings.
In this section, we carry out numerical experiments for illustrating the application and validity of SNM.We adopt the network of Nguyen and Dupuis [46] with 4 O-D pairs, 25 paths, and 19 directed links (Figure 3).The link sequences of the paths are shown in Table 1.
We employed the mean link travel time and travel time covariance shown in (24) and ( 25) that are calculated by assuming  = 4, respectively.Parameters for the BPR function,  and , are 2.62 and 5, respectively.The coefficient of variation of total O-D flow, V, is 0.1.The other parameters used in (24) and ( 25) are shown in Table 2.
In this study, we prepared four experimental cases by changing the effect of the path travel time variance on the driver's path choice behavior in the network.In the first case, we assume a risk-neutral driver who chooses a path based only on mean path travel time in the network.Therefore, the relative weight assigned to path travel time variance, , in (30) is 0 in this case.The other three cases assume three types of risk-averse drivers.The relative weight assigned to var[Ξ  ] is assumed as 0.3 in these three cases.In the second case, we assume that there is no travel time correlation between two different links in the network.In the third case, SNM is employed (i.e., all travel time covariance between two adjacent links in the network is introduced to calculate the path travel time variance).In the fourth case, all travel time covariance between two different links in the network is introduced to calculate the path travel time variance.In the latter three cases, the path travel time variance is calculated as (39) by using both the relative weight assigned to path travel time variance, , and the set of link(s), , shown in Table 3.  (1) 3-5-8-14-16 (25) 3-6-12-14-16 For solving the first three cases, we employed MSA as a link-based solution algorithm.The final case can be solved by minimizing the gap function for NCP [17,47,48].In fact, the paths set shown in Table 1 is prepared only for the final case.Only case 4 cannot be solved by MSA which does not require paths enumeration and thus is applicable to a large network.On the other hand, it is difficult to solve the problem of a large network using the method based on the gap function that requires paths enumeration.

Results
. Table 4 shows the mean path flows for each case.For the first three cases, since we applied MSA to solve corresponding path choice problems, the path flows for each case were not uniquely determined.Even so, presenting path flows for the first three cases may be useful to understand roughly how path flows change according to different expressions of path travel time variance.Also, by presenting the path flows for cases 1-3 to which MSA was applied, it is easy to understand that the generalized travel times for the paths between each O-D pair that are used by the drivers are the same although both the mean travel time and travel time variance of a path between the O-D pair can be different from the others.In contrast, the path flows for case 4 were uniquely determined, since we applied a path-based algorithm when solving its path choice problem.In all cases, if a path flow is zero, such path flow is denoted by bold figures in Table 4 so that we can know that the corresponding path is not chosen by the drivers in the network.Table 5 shows generalized path travel times for all cases which are calculated by using both the mean path travel times shown in Table 6 and the path travel time variance shown in Table 7.It is observed that the paths chosen by the drivers have the minimum generalized travel time and that the paths which are not chosen by the drivers have generalized travel times which are equal to or greater than the minimum generalized travel time.The generalized path travel times for unused paths are denoted by bold figures in Table 5.It is shown in Tables 4, 6, and 8 that although the mean path flows of each case are different from the other cases, mean link flows are similar among all cases.Since mean travel time of a path is calculated by using mean link flows, therefore similar mean path travel times are obtained among the four cases.
For O-D pair 1, the flows of paths 4, 5, and 8 in cases 2-4 are smaller than those in case 1 whereas the flows of path 1 in cases 2-4 are larger than that in case 1.These differences can be explained as follows.The travel time variances of path 1 in cases 2-4, which are denoted by bold figures in Table 7, are much smaller than those of paths 2-8, although the mean travel time of path 1 is longer than those of paths 2-8, which are denoted by bold figures in Table 6.However, path 1 has longer mean travel time than paths 2-8, and that path is more reliable in terms of travel time variance than paths 2-8 in cases 2-4.In total, the path choice probabilities of path 1 in cases 2-4 are larger than that in case 1.In contrast, since the drivers in case 1 choose their paths based only on mean travel times, the flow of path 1 is larger than those in cases 2-4.This path choice switch can be observed in the other O-D pairs (e.g., from path 10 to path 9 for O-D pair 2, from paths 17 and 19 to paths 16 and 18 for O-D pair 3, and from paths 20 and 24 to path 25 for O-D pair 4.
We can find from Table 6 a tendency of the mean path travel time in case 1 to be greater than those in cases 2-4.Surprisingly, this tendency holds for all paths in the network except for path 1. Obviously, this tendency was derived from the introduction of travel time variance to the drivers' path choice behavior.We will then check how total mean travel time in the network is shortened by the introduction of travel time variance.An index for the total mean travel time for cases  ∈ {1, . . ., 4}, TTT  , can be given by By using both the mean link flows shown in Table 8 and the mean link travel time shown in Table 9, TTT  ∀ ∈ {1, . . ., 4} are calculated as 2.847×10 5 , 2.789×10 5 , 2.797×10 5 ,  an efficient use of the road network in terms of the mean travel time.This effect is similar to the one obtained by applying the system optimal principle (i.e., Wardrop's second principle).From the mean link flows shown in Table 8, it is observed that the mean flows of links 5, 9, 11, 12, 16, and 18 in cases 2-4 are larger than those in case 1, whereas those of links 4, 7, 8, 13, 14, 15, and 17 in cases 2-4 are smaller than those in case 1.These changes happened by the path flow changes in cases 1-4 which we have discussed.
As discussed in the previous section, if we do not consider the effect of travel time reliability on the driver's path choice behavior in the network, unrealistic results may be obtained.However, the introduction of all travel time covariance between two different links in the network to the calculation of the path travel time variance is difficult to implement.Therefore, we proposed SNM, in which all travel time covariance between two adjacent links in the network is considered in calculating the path travel time variance.Our concern is now how far the results obtained in cases 1-3 are from the results obtained in case 4. The results obtained in case 4 can be regarded as the exact solution, since it is obtained by adopting the path travel time variance denoted by (28).Table 10 shows the coefficients of correlation of mean link flows among four cases.The bold figures in the table show the coefficients of correlation between case 4 and the other three cases.One can see that SNM (case 3) reproduces almost the same mean link flows as case 4, since the coefficient of correlation between case 4 and case 3 is 1 whereas the other coefficients of correlation are less than 1.Note that the coefficient of correlation between case 4 and case 1 is the smallest.Almost the same results were obtained by calculating the coefficients of correlation of mean link travel times.These results may bring about two practical implications that all travel time covariance between two different links in the network may not be required in calculating mean link flows.Instead, all travel time covariance between two adjacent links in the network is required in calculating mean link flows, and that calculation of all travel time covariance between two adjacent links in the network is easy even in a large network.

Conclusions
In this study, we proposed a simplified network model, that is, SNM, for travel time reliability analysis.The uncertainty addressed in this model is that of O-D flows.In this model, the generalized path travel time is a linear combination of mean path travel time and path travel time variance.In calculating the path travel time variance, we consider all travel time covariance between two adjacent links in the network in SNM.A risk-averse driver in the network is assumed.By applying a network representation used for addressing intersection delays, SNM can be solved by applying a standard link-based algorithm.The other property of SNM which needs to be emphasized here is that its formulation requires only mean network flows.This property may be important for practitioners, since once the coefficient of variation of total O-D flow is determined, one can apply SNM to real road network analysis for which network data sets for a conventional network model (e.g., a deterministic/stochastic user equilibrium traffic assignment model) are already prepared.Numerical experiments are carried out for illustrating the applications and validity of SNM.The experiments assumed four types of drivers in the network.The first type of driver is a risk-neutral driver who chooses a path based only on mean path travel time.The other three types of drivers are risk-averse drivers who choose their paths based on both mean path travel time and path travel time variance.The second type of driver's path travel time variance is calculated by assuming statistically independent link travel time.The third type of driver's path travel time variance is calculated by considering all travel time covariance between two adjacent links in the network.The fourth type of driver's path travel time variances is calculated by considering all travel time covariance between two different links in the network.The path travel time variance of the fourth type of driver is the exact one in the network.It is shown that mean network flows obtained by assuming the risk-neutral driver differ as a whole from those obtained by assuming the risk-averse drivers.These differences in mean network flows are generated by the effects of travel time reliability on path choice behavior by the driver.It is also shown that mean link flows obtained by assuming the third type of driver, that is, mean link flows calculated by SNM, are almost the same as the mean link flows calculated by assuming the fourth type of driver.In a practical sense, it may be difficult to calculate network flows in a large road network by assuming the fourth type of driver.In contrast, SNM can be easily applied to a large road network in calculating network flows.
In this study, we recognize that the introduction of path travel time variance considering all travel time covariance between two links in a road network to the generalized path travel time is more important in expressing the driver's route choice behavior than the introduction of the skewness of link travel time to the generalized path travel time.Therefore, SNM introduces the path travel time variance considering all travel time covariance between two adjacent links in the network to the generalized path travel time.It is interesting to see how the path choice probabilities which are calculated by assuming both the statistically independent link travel time and the skewness of link travel time differ from the probabilities calculated by SNM.In light of this, there is a need for a network equilibrium model that introduces a driver's path choice preference based on (2).A deterministic path choice model based on Wardrop's first principle is employed in this study in order to express the driver's path choice behavior in the network.The introduction of stochastic models, for example, logit-based models or probitbased models, to SNM is needed, in order to express the driver's perception error on the generalized travel time.These two challenges are our future tasks.
) is the mean flow that passes through both links  and .If  =  in (22), then we obtain var[  ] = (V⋅V  ) 2 ∀ ∈ A. In fact, this assumption can be justified if we regard total O-D flow in the network,  = ∑ ∈   , as a random variable with a mean of [] = (= ∑ ∈   ) and a variance of var[] = (V⋅) 2 , where   =   ⋅.  =   / is the proportion of the O-D flow,   , to total O-D flow, .In this case, all O-D flows are statistically dependent on each other.The conservation of the O-D flow variance in relation to the total O-D flow variance holds (Appendix C).

Figure 2 :
Figure 2: Network representation for addressing intersection delays.

Figure 3 :
Figure 3: The Nguyen and Dupuis network.

Table 1 :
Paths and link sequences.
Conservation of the O-D Flow Variance in Relation to the Total O-D Flow Variance under the Stochastic Total O-D Flow Total O-D flow variance is calculated as the sum of O-D flow variance as follows: var [∑ ∈   ] = ∑ ∈ var [  ] [33] assumed to be a random variable with a mean of [  ] =   and a variance of var[  ] = (V  ⋅   )2, where V  is the coefficient of variation of the random flow   .Following Lam et al.[33], the stochastic flow on path  ∈   ,   , is then given by  =   ⋅   ∀ ∈ , ∀ ∈   . is a random variable with a mean of   =   ⋅   ≥ 0 and a covariance of cov[  ,   ] =   ⋅   ⋅ var[  ], where   ( ∈   ) is path choice probability which can be determined by a path choice model (DUE, SUE, etc.).The following flow conservation law holds for each O-D pair:

Table 2 :
Link travel time parameters.

Table 3 :
Assumptions of each case.

Table 5 :
Generalized path travel time.

Table 6 :
Mean path travel time.

Table 7 :
Path travel time variance.

Table 9 :
Mean link travel time.

Table 10 :
Coefficients of correlation of mean link flows.