Next Article in Journal
Actual and Perceived Motor Competence in Relation to Body Mass Index in Primary School-Aged Children: A Systematic Review
Previous Article in Journal
Built Environment Characteristics, Daily Travel, and Biometric Readings: Creation of an Experimental Tool Based on a Smartwatch Platform
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Target-Oriented User Equilibrium Considering Travel Time, Late Arrival Penalty, and Travel Cost on the Stochastic Tolled Traffic Network

Department of Management Science and Engineering, School of Business, Qingdao University, Qingdao 266071, China
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(17), 9992; https://0-doi-org.brum.beds.ac.uk/10.3390/su13179992
Submission received: 24 July 2021 / Revised: 25 August 2021 / Accepted: 1 September 2021 / Published: 6 September 2021

Abstract

:
In this paper, we employ a target-oriented approach to analyze the multi-attribute route choice decision of travelers in the stochastic tolled traffic network, considering the influence of three attributes, which are (stochastic) travel time, (stochastic) late arrival penalty, and (deterministic) travel cost. We introduce a target-oriented multi-attribute travel utility model for this analysis, where each attribute is assigned a target by travelers, and travelers’ objective is to maximize their travel utility that is determined by the achieved targets. Moreover, the interaction between targets is interpreted as complementarity relationship between them, which can further affect their travel utility. In addition, based on this travel utility model, a target-oriented multi-attribute user equilibrium model is proposed, which is formulated as a variational inequality problem and solved with the method of successive average. Target for travel time is determined via travelers’ on-time arrival probability, while targets for late arrival penalty and travel cost are given exogenously. Lastly, we apply the proposed model on the Braess and Nguyen–Dupuis traffic networks, and conduct sensitivity analysis of the parameters, including these three targets and the target interaction between them. The study in this paper can provide a new perspective for travelers’ multi-attribute route choice decision, which can further show some implications for the policy design.

1. Introduction

Traffic assignment models play a fundamental role in transportation planning and real-time applications, e.g., traffic prediction and optimal routing, and this is important for the desirable development of cities, e.g., a low-carbon city [1]. The work of [2] proposed two widely used principles for assigning vehicular flows to the traffic network, and his first principle is closely related to our study, i.e., the user equilibrium principle, which can be stated as: No one can decrease his (or her) route travel time by unilaterally changing his (or her) route choice decision. Although Wardrop’s user equilibrium principle is widely used, it has some unrealistic assumptions, e.g., the traffic network is deterministic, only travel time is considered in travelers’ route choice decision, travelers always grasp all information about the network conditions, and always incline to adopt the route with the shortest travel time, i.e., they are naturally assumed to be “rational people” [3]. One of the aims of this paper is to partially overcome the unrealistic assumptions of this principle.
Uncertainty is inevitable in the real-world traffic network, which could come from the demand side (e.g., peak-hour demand surge), supply side (e.g., traffic accident), or both [4,5]. Uncertainty is categorized into two kinds by [6], which are risk and ambiguity. The fundamental difference between these two kinds is whether the full distribution information of outcomes is known or not. Particularly, it is known in the risk modeling, while unknown in the ambiguity modeling. Scholars have conducted numerous studies on travelers’ route choice behavior under risk (e.g., [4,7,8]) and ambiguity (e.g., [5,9,10,11]), and the resultant models are fundamental to the traffic assignment, e.g., in the formulation of variational inequality (e.g., [4,7,8,12]). For example, the authors in [4] proposed the mean-excess travel time (METT) model based on the travel time budget (TTB) model in [12]. The authors in [7] assumed travelers’ confidence level belongs to an interval due to day-to-day travel time variations, and proposed the subjective-utility TTB model. The work of [8] introduced a new travelers’ non-expected disutility model, which can incorporate METT and TTB as the special cases. The authors in [9] formulated the shortest path problem with distributional uncertainty as a distributionally robust mixed integer linear programming. The authors in [10] discussed the robust shortest path problem with distributional uncertainty. The authors in [5] assessed the travel time reliability of the traffic network with partial information of the distribution. The authors in [11] formulated Wasserstein’s distributionally robust shortest path as a mixed 0–1 convex problem. Among all the studies, risk could come from the distributions of free-flow travel time, e.g., gamma distribution and normal distribution, distribution of the link capacity degradation, e.g., uniform distribution, and the distribution of travel demand, e.g., normal distribution [4], while the ambiguity set of probability distributions could be described by the moment information, e.g., [9,10], or by the statistical distance [11]. The study in this paper focuses on the risk modeling, i.e., route choice decision making by travelers under risk.
Substantial studies have demonstrated that there are multiple attributes affecting the route choice decision made by travelers (e.g., [13,14,15,16]). For example, the authors in [13] conducted the statistical analysis to examine the route attributes which travelers consider to be important during their route choice, and the results show that three most important attributes are shorter travel time (40% of respondents choose this), travel time reliability (32% of respondents choose this), and shorter distance (31% of respondents choose this). Besides the empirical studies, scholars also make substantial progress on travelers’ multi-attribute route choice behavior form the theoretical side, and two main methodologies are widely used in the current studies. One of these two is to aggregate different criteria into the generalized travel cost (e.g., [4,8,12,17,18]), and the other one is to study different criteria separately, resulting in a multi-objective route choice model (e.g., [19,20,21]). For example, in the TTB model proposed by the authors in [12], travel time and travel time reliability are considered, while in the METT model proposed by the authors in [4], travel time, travel time reliability, and travel time unreliability are discussed. The work of [19] treated travel time and travel time reliability separately, resulting in the bi-objective route choice model, and the authors in [21] discussed TTB and distance separately. Research in this paper also discusses travelers’ multi-attribute route choice behavior, and combines all the attributes aggregately.
Another noteworthy topic closely related to this paper is to use behavioral economics theory in the study of traffic and transportation problems, e.g., prospect theory (PT) proposed in [22], and cumulative prospect theory (CPT) proposed in [23]. (C)PT can be used to explain a variety of phenomena that cannot be explained with the expected utility theory (e.g., [24]). The authors in [25] firstly proposed to incorporate CPT in the equilibrium analysis on the stochastic network, where the reference is exogenous. Following this paper, the authors in [3] proposed an endogenous reference when applying the CPT in the equilibrium analysis, and discussed a behaviorally consistent congestion pricing problem based on the equilibrium principle. In the recent decade, (C)PT are also widely used in other aspects of traffic and transportation study, e.g., choice in the unreliable transport networks [26], day-to-day learning [27], mode choice [28], and electric vehicle drivers’ charging behavior of battery [29], to name a few. (C)PT and the methodology used in this paper both belong to the descriptive paradigm, i.e., they have similarities, while they also have differences following the discussion in [30], which are summarized as follows: (1) The target under this methodology serves as a specific reference under (C)PT; (2) (C)PT has two typical characteristics, namely, probability distortion function and reference-dependent value function, while stochastic correlation among these attributes and the interaction between targets, namely the complementarity relationship between targets, are the two typical characteristics of this methodology.
The closest paper related to our study here is [31], where the author proposed a new general methodology to study travelers’ multi-attribute route choice behavior on the traffic network (in Section 2.1, we review this methodology to make this paper self-contained). Based on this methodology, the author considered stochastic travel cost and travel time simultaneously in a tolled traffic network, and further discussed a target-oriented bi-attribute travel route utility model and bi-attribute user equilibrium based on aforementioned utility model. Nevertheless, although the authors in [31] proposed a general methodology for the multi-attribute route choice modeling, they only discuss a bi-attribute application, which motivates the multi-attribute study in this paper. To be specific, we apply the methodology proposed in [31] in the multi-attribute analysis, and extend the application scope of this methodology.
Additionally, the study in this paper focuses on travelers’ a priori route choice behavior, which means that they will not change their route choice decision after this is determined. In contrast, there is another kind of research called adaptive route choice, where travelers can change their route choice behavior based on the information of a traffic network. For example, the authors in [32] firstly adopted (C)PT utility function to capture the impact of risk attitude and online information on route decision adjustments. The authors in [33] formulated information location models for minimizing travel cost of travelers, which can be classified into three categories according to congestion effect and vehicle number. The authors in [34] proposed the adaptive routing policy problem and its necessary conditions for optimality in the stochastic network, considering three types of partial online information, namely, delayed global, global pre-trip, and up-to-date radio. The author in [35] developed a label-setting-based algorithm to identify the adaptive least-expected time hyperpaths. The work of [36] proposed user equilibrium with recourse (UER) models in the stochastic network, where the link states is defined by probability mass function. Furthermore, the work of [37] formulated the continuous network design problem under UER as bi-level programming. In the studies of experimental economics, the authors in [38] compared travelers’ route choice decision under no online information (expected user equilibrium) and perfect online information (UER), and the results demonstrate that the provision of online information brings heavy travel costs. The authors in [39] concluded the provision of information can reduce travelers’ risk aversion behavior, and further reduce their valuation on information and reliability.
With all the above discussions, we adopt the target-oriented methodology to analyze the multi-attribute route choice decision of travelers in the scenario of tolled traffic network, considering three attributes, which are (stochastic) travel time, (stochastic) late arrival penalty (LAP), and (deterministic) travel cost, following the work of [18] (in the work, the authors used travel distance, which can also be reinterpreted as travel cost following our definitions, as our paper assumes travel cost is given in advance). That is, we extend the application scope of the general methodology in [31]. LAP can be used to model travel time unreliability, which brings the benefit for the practical use of our model. Moreover, travel cost denotes the congestion tolls on the traffic network. We introduce a target-oriented multi-attribute travel utility model for this analysis, where each attribute is assigned a target by travelers and travelers’ objective is to maximize their travel utility that is determined by the achieved targets. Moreover, the interaction between targets is interpreted as a complementarity relationship between them, which can further affect their travel utility. The interaction between targets means that travelers are willing to acquire more targets for additional utility. Particularly, we propose a rule to determine the utility values based on the essentiality of the interaction between targets.
Next, the user equilibrium model based on aforementioned travel utility model is proposed, which is formulated as a variational inequality problem and solved with method of successive average. Target for travel time is determined via travelers’ on-time arrival probability, while targets for late arrival penalty and travel cost are given exogenously. Lastly, we apply the proposed model to the Braess and Nguyen–Dupuis traffic networks, focusing on the impact of these three targets and the target interaction among them, i.e., the sensitivity analysis.
The rest of the sections and contents of the paper are shown As follows. Review of the general methodology and target-oriented multi-attribute travel utility model considering the travel time, LAP, and travel cost are presented in Section 2. Target-oriented multi-attribute user equilibrium based on this new route choice model and its solution algorithm are listed in Section 3. In Section 4, we apply the proposed model on the Braess and Nguyen–Dupuis traffic networks, and conduct sensitivity analysis of these three targets and the target interaction between them. Finally, the major conclusions and findings are shown in Section 5.

2. Target-Oriented Multi-Attribute Travel Utility Model

The main content of this section is to study the target-oriented route choice decision of travelers based on the methodology described in [31], where travelers need to consider travel time, LAP, and travel cost simultaneously. Next, we review this general methodology considering N attributes first, although we only study the impact of three attributes.

2.1. Review of the Target-Oriented Methodology

Definition 1. 
Considering N attributes for a specific route, denoted by X = { X 1 , X 2 , , X N } , which are assumed to be stochastic and correlated, a traveler can be seen as target-oriented if his (or her) utility for an outcome, denoted by x = { x 1 , x 2 , , x N } , is determined by the targets acquired from this outcome, where each attribute is assigned one target by travelers, denoted by d = { d 1 , d 2 , , d N } .
Here, we assume all the attributes are stochastic and correlated. This assumption demonstrates the generality of this methodology, and the methodology for deterministic attributes or stochastic attributes without correlation can be obtained based on our methodology. That is, this methodology can handle the deterministic attributes as shown in the following application, as the deterministic attribute is a special case of the stochastic one. We introduce an indicator function for each attribute X i , ( i = 1 , 2 , , N ) , which is formulated as
I i = { 1 ,   if   x i d i 0 ,   otherwise  
When I i = 1 , we say this target is achieved, which means that travelers are satisfied with this route with outcome x i if only this attribute is considered, while it is not achieved when I i = 0 . Considering there are N attributes, travelers’ satisfaction depends on the targets that are achieved as shown in Definition (1). We introduce a set of indices { i | I i = 1 } , denoted by A , i.e., we use A to denote the targets which are achieved. As each I i can be 0 or 1, there are 2 N different A . For example, when all the I i are 0, A is an empty set; when x 1 d 1 , x 3 d 3 , and x i > d i ( i = 2 , 4 , 5 , , N ) , A is { 1 , 3 } .
Accordingly, the occurrence probability for each A is written as P ( x i d i , i A ; x j > d j , j A ) considering that all the outcomes are stochastic. The utility when different targets are achieved, termed as target-oriented utility, is denoted as ξ ( I i = 1 , i A ; I j = 0 , j A ) , which is rewritten as ξ ( A ) or ξ i A for brevity hereinafter, depending on the context. Moreover, we assume that the value of ξ ( A ) is larger when more targets are achieved, and ξ ( A ) = 1 when all the targets are achieved and ξ ( A ) = 0 when no target is achieved. Therefore, only 2 N 1 situations need to be considered, and with the consideration that they are mutually exclusive, the target-oriented route utility ω is formulated as the sum of 2 N 1 values, where each value is the product of the target-oriented utility and the corresponding occurrence probability, i.e.,
ω = A N , A ξ A P A
If the joint probability P is known, all the occurrence probabilities used in Equation (2) can be evaluated. However, it is rarely known in practice, as it takes a huge cost to grasp perfect information. We adhere to the inclusion-exclusion rule proposed in [31] and reformulate Equation (2) as
ω = A N , A ω A P ( X i d i , i A )
where ω A = B A 1 | A | | B | ξ B .
This reformulation brings two benefits. Firstly, the marginal distributions, which might be estimated with the data-driven methods, can be used to evaluate the joint probabilities. The second benefit is that ω A can be employed to grasp the interaction between targets, which can also be interpreted as the complementarity relationship between them. Both of these benefits will be shown in detail in the next subsections.

2.2. Network Representation and Attributes

We use G ( V , A ) to denote a connected and directed traffic network, where V denotes the node set and A denotes the link set. Let R and S denote the set of origins and destinations, respectively, and r s denote O D pair from an origin r R to a destination s S . Let q r s denote the travel demand between this O D pair. We use p r s to denote all the routes between the O D pair, and the route flow vector ( , f p r s , ) is denoted by f , where f p r s denotes the flow on route p P r s . The link travel time is stochastic as we discuss the stochastic traffic network, which is denoted by T a , a A , and the link flow is denoted by v a , a A . We use Δ = [ δ p a r s ] to denote the link-route incidence matrix, where δ p a r s = 1 if link a is on the route p , and 0, otherwise. Finally, we use c a to denote the toll on link a , a A .
As aforementioned, we study traveler’s target-oriented route choice decision considering (stochastic) travel time, (stochastic) LAP, and (deterministic) travel cost. In this section, we discuss the details on these three attributes, e.g., the probability distribution functions and the corresponding targets we used.

2.2.1. Stochastic Route Travel Time and Its Corresponding Target

Several studies have discussed how to derive the probability distribution for route travel time (e.g., [40,41]), and we adopt the method in [40]. The link travel time throughout this paper is the BPR performance function, which is written as
t a ( x a ) = t a 0 [ 1 + β ( v a C a ) n ]
where for each link, a , t a 0 is the free-flow travel time, C a is the capacity, v a is the flow as shown before, and t a ( x a ) is the travel time with flow v a . β and n are deterministic parameters. The value of β is 0.15 and the value of n is 4. In [40], stochasticity comes from the link capacity degradation, and route travel time is shown to follow the normal distribution N ( E ( T p r s ) , σ ( T p r s ) ) , where E ( T p r s ) denotes the mean of the stochastic travel time T p r s for route p between O D pair r s , and σ ( T p r s ) denotes its standard deviation. The detailed derivation is relegated into Appendix A.
Next, we denote the route travel time as γ t and show how to determine the target for it, which is motivated by the method proposed in [3]. Travelers are assumed to own the same on-time arrival probability, denoted by θ ( θ 0.5 ) . θ = 0.5 means that travelers are risk-neutral, and larger values of θ means that travelers are more risk-averse. Target for travel time is represented as the travel time budget specified by the travelers to ensure their desired on-time arrival probability. Therefore, when travelers choose a route p from multiple routes of O D pair r s , target on this route, denoted by γ p r s , is formulated as
γ p r s = min { γ Pr ( T p r s γ ) θ }
As O D pair r s is connected by multiple routes, target for travel time of a route belonging to this O D  pair is represented as the minimum value of the budget of all travel times on all routes, i.e., γ t r s = min p P r s γ p r s . When we employ this method to obtain the target for travel time of a route, it seems that travelers specify the exogenous on-time arrival probability first, and then the target is endogenously determined with this exogenously given on-time arrival probability.

2.2.2. Stochastic Route LAP and Its Corresponding Target

Given the target for travel time of a route, denoted by γ t , following the idea in [18], route LAP is defined as L p r s = max { 0 , T p r s γ t } , p P r s , r R , and s S . In the work of [18], route LAP is defined, given a longest possible route travel time. While we assume this longest possible route travel time is the target for route travel time. One can also define the longest possible route travel time as the sum of target for route travel time and some buffer time, and our arguments can be amplified accordingly.
In this paper, target for LAP is assumed to be exogenously given, denoted by γ l , which can be interpreted as the allowable delay by the company. Next, we construct the probability distribution function for route LAP. We see that the range of route LAP can be written as [ 0 , L max ] , where L max is the maximum value of route LAP. Given the target for travel time of a route (denoted as γ t ), and the above-derived probability distribution function of route travel time, we construct the following probability distribution function for route LAP.
P ( L p r s γ l ) = { 0   if   γ l < 0 F T ( γ t + γ l )   if   γ l 0
where F T denotes the cumulative distribution function of route travel time, and γ l denotes the target for route LAP. Equation (6) can be verified as follows: when γ l < 0 , P ( L p r s γ l ) = 0 can be obtained, as 0 is minimum value of the range of route LAP; when γ l 0 , we see that P ( L p r s γ l ) = P ( T p r s γ t γ l ) = P ( T p r s γ l + γ t ) = F T ( γ l + γ t ) . Especially when γ l = 0 , P ( L p r s = γ l ) = F T ( γ t ) .

2.2.3. Deterministic Route Travel Cost and Its Corresponding Target

The route travel cost C p r s , p P r s , r R , and s S  can be obtained according to the link travel cost, and is formulated as C p r s = a A δ p a r s c a . In this paper, target for travel cost is also assumed to be exogenously given, denoted by γ c , which can be interpreted as the acceptable expense by the travelers. Therefore, the probability that the target for travel cost can be reached on some particular routes is 1, while it cannot be achieved on other routes for the O D pair r s .
Considering the attributes and the corresponding targets, Equation (3) can be rewritten as
ω p r s = ξ 1 P ( T p r s γ t ) + ξ 2 P ( L p r s γ l ) + ξ 3 P ( C p r s γ c ) + ( ξ 13 ξ 1 ξ 3 ) P ( T p r s γ t , C p r s γ c ) + ( ξ 23 ξ 2 ξ 3 ) P ( L p r s γ l , C p r s γ c ) + ( ξ 12 ξ 1 ξ 2 ) P ( T p r s γ t , L p r s γ l ) + ( 1 ξ 13 ξ 23 ξ 12 + ξ 1 + ξ 2 + ξ 3 ) P ( T p r s γ t , L p r s γ l , C p r s γ c )
where ω p r s denotes the target-oriented multi-attribute travel utility (ToMaTU) on route p from an origin r R to a destination s S . That is, Equation (7) is the specific form of Equation (3) considering the attributes and targets used in this study. The target-oriented utility ξ 1 represents the situation where only the first target is achieved, the target-oriented utility ξ 13 represents the situation where the first and the third targets are achieved as discussed before, and others have the similar meanings. Moreover, hereinafter P ( T p r s γ t ) and P ( L p r s γ l ) are called target achievement probability for travel time (TAPt), and LAP (TAPl), respectively.

2.3. Joint Probability Evaluation Derived from the Marginal Distributions

In this section, we show how to evaluate the joint probabilities in Equation (7) derived from the marginal distributions, namely the marginal distribution of LAP, travel cost, and travel time mentioned in last section.
First, we discuss the effect of the marginal distribution of route travel cost. We know that route travel cost is deterministic, and thus, determined by the scenario of whether the target of travel cost is finally achieved, i.e., P ( C p r s γ c ) is 0 or 1, we have the following two situations. When target for travel cost is achieved, ToMaTU considering LAP, travel cost, and travel time is rewritten as
ω p r s = ξ 1 P ( T p r s γ t ) + ξ 2 P ( L p r s γ l ) + ξ 3 + ( ξ 13 ξ 1 ξ 3 ) P ( T p r s γ t ) + ( ξ 23 ξ 2 ξ 3 ) P ( L p r s γ l ) +   ( ξ 12 ξ 1 ξ 2 ) P ( T p r s γ t , L p r s γ l ) + ( 1 ξ 13 ξ 23 ξ 12 + ξ 1 + ξ 2 + ξ 3 ) P ( T p r s γ t , L p r s γ l )   = ( ξ 13 ξ 3 ) P ( T p r s γ t ) + ( ξ 23 ξ 3 ) P ( L p r s γ l ) + ξ 3 + ( 1 ξ 13 ξ 23 + ξ 3 ) P ( T p r s γ t , L p r s γ l )
When the target for travel cost fails to be achieved (denoted with tilde symbol), ToMaTU considering LAP, travel cost, and travel time is rewritten as
ω ˜ p r s = ξ 1 P ( T ˜ p r s γ t ) + ξ 2 P ( L ˜ p r s γ l ) + ( ξ 12 ξ 1 ξ 2 ) P ( T ˜ p r s γ t , L ˜ p r s γ l )
From Equation (8) and Equation (9), we see the effect of marginal distribution of route travel cost on the ToMaTU, and we also see that joint probability P ( T p r s γ t , L p r s γ l ) needs to be evaluated further. For this probability, we follow the method in [42], and obtain that the random vector composed of route travel time and LAP is comonotonic. Thus, the joint probability P ( T p r s γ t , L p r s γ l ) equals to P ( T p r s γ t ) .
Therefore, we can further simplify the ToMaTU, i.e., Equations (8) and (9) in our application. In the situation where the target for travel cost is successfully achieved, ToMaTU considering LAP, travel cost, and travel time is further rewritten as
ω p r s = ( ξ 13 ξ 3 ) P ( T p r s γ t ) + ( ξ 23 ξ 3 ) P ( L p r s γ l ) + ξ 3 + ( 1 ξ 13 ξ 23 + ξ 3 ) P ( T p r s γ t , L p r s γ l )   = ( ξ 13 ξ 3 ) P ( T p r s γ t ) + ( ξ 23 ξ 3 ) P ( L p r s γ l ) + ξ 3 + ( 1 ξ 13 ξ 23 + ξ 3 ) P ( T p r s γ t )   = ( 1 ξ 23 ) P ( T p r s γ t ) + ( ξ 23 ξ 3 ) P ( L p r s γ l ) + ξ 3
In the situation where the target for travel cost fails to be achieved (denoted with tilde symbol), ToMaTU considering LAP, travel cost, and travel time is further rewritten as
ω ˜ p r s = ξ 1 P ( T ˜ p r s γ t ) + ξ 2 P ( L ˜ p r s γ l ) + ( ξ 12 ξ 1 ξ 2 ) P ( T ˜ p r s γ t , L ˜ p r s γ l )   = ξ 1 P ( T ˜ p r s γ t ) + ξ 2 P ( L ˜ p r s γ l ) + ( ξ 12 ξ 1 ξ 2 ) P ( T ˜ p r s γ t )   = ( ξ 12 ξ 2 ) P ( T ˜ p r s γ t ) + ξ 2 P ( L ˜ p r s γ l )
From Equations (10) and (11), we see the hybrid effect of marginal distribution of route LAP and marginal distribution of route travel time on the ToMaRU. After this, we see that all the probabilities in ToMaTU can be obtained.

2.4. Target Interaction

In this paper, we follow the idea in [31] to model the interaction between targets, which can also be interpreted as the complementarity relationship between them. One of the essential meanings of target interaction is that achieving more targets could result in receiving more utility to the travelers, which is further reflected in the supermodularity of aforementioned utility function.
If, for all given subset A an B of N , the utility function always satisfies the inequation ξ ( A ) + ξ ( B ) ξ ( A B ) + ξ ( A B ) , then target interaction is called complementarity in the multi-attribute modeling. Here, we eliminate all the identical equations, e.g., B = { 1 } and A = { 1 , 2 , 3 } . The complementarity relationship is motivated by the idea in economics, which means that combining two or more goods can bring more values, even though they have less or no value by themselves [43]. The complementarity relationship is further divided into two additional cases, namely, imperfect complementarity relationship and perfect complementarity relationship, respectively. In the former case, the inequation ξ ( A ) + ξ ( B ) ξ ( A B ) + ξ ( A B ) is satisfied for all subsets A and B of N , while in the latter case, the inequation ξ ( A ) + ξ ( B ) < ξ ( A B ) + ξ ( A B ) is satisfied for all subsets A and B of N . We see that the imperfect complementarity corresponds to the modularity of the utility function. Moreover, comparing the imperfect complementarity relationship with the perfect complementarity relationship, it is clear that travelers in the perfect complementarity relationship have greater willingness to achieve more targets, as achieving more targets strictly represents that travelers can receive more utility; we also know that mixed situations, namely perfect complementarity among some targets and, meanwhile, imperfect complementarity among other targets, are not discussed, e.g., ξ 1 + ξ 2 = ξ 12 and ξ 1 + ξ 3 < ξ 13 . That is, if the target interaction is a perfect complementarity relationship, achieving more targets will bring more utilities to travelers strictly.
In the following, we propose a rule to determine the values of different target-oriented utilities. We assume that ξ 1 / ξ 2 = α 1 , and ξ 1 / ξ 3 = α 2 . Here α i ( i = 1   or   2 ) measures the relative importance between different target achievements, e.g., larger value of α 1 means that the target achievement of travel time on route deserves more attention, in contrast to the achievement of target set for LAP on route. Hereinafter, we term α 1 and α 2 as utility ratios. The relative importance between target achievement of LAP on route and that of travel cost on route can be obtained as α 2 / α 1 , and thus we only need two utility ratios. Starting from the imperfect complementarity relationship, we solve the subsequent system of equations to obtain different target-oriented utilities.
{ ξ 1 + ξ 2 + ξ 3 = 1 ξ 1 + ξ 2 = ξ 12 ,   ξ 1 + ξ 3 = ξ 13 , ξ 2 + ξ 3 = ξ 23 ξ 1 / ξ 2 = α 1 ,   ξ 1 / ξ 3 = α 2
While in the perfect complementarity relationship, we introduce two additional parameters, β b and β s , and use prime symbol to denote the values in this relationship. We formulate ξ 12 (termed as bi-attribute target-oriented utility) as ξ 12 / β b , and formulate ξ 1 (termed as single-attribute target-oriented utility) as ξ 1 / ( β b β s ) in the perfect complementarity relationship. Other target-oriented utilities in this relationship can be defined similarly. According to the definition of perfect complementarity relationship, we have β b > 1 , and β s > 1 , meanwhile, β b > 2 ( 1 / β s ) , where 2 is equal to ξ 12 + ξ 13 + ξ 23 , and 1 is equal to ξ 1 + ξ 2 + ξ 3 , based on the system of equations (12). The values of β b and β s , and the relationship between them, are determined according to the definition of perfect complementarity relationship. In particular, if we choose β b = 1 and β s = 1 , the imperfect complementarity relationship is reproduced. Hereinafter, we term β b and β s as the complementarity ratios. Furthermore, β b is called bi-attribute complementarity ratio, as it can be used with bi-attribute target-oriented utility, while β s is called single-attribute complementarity ratio as it can only be used with single-attribute target-oriented utility. The bi-attribute complementarity ratio β b shows the interaction among three targets, while the single-attribute complementarity ratio β s shows the interaction between two targets. Moreover, we see that larger values of β b and β s implicates travelers have more motivation to reach these targets, as achieving more targets strictly brings extra utility. With this parameter setting rule, all the feasible results of travelers’ target-oriented utilities could be achieved if we change the values of the four parameters properly.
From Equations (10) and (11), we see that it seems that some target-oriented utilities, e.g., ξ 1 , have no impact on travelers’ route choice decision. However, according to the proposed parameter setting rule, these target-oriented utilities also have their impact, not in a direct manner, but in an indirect manner, as a change their values will cause the change of other target-oriented utilities shown in Equations (10) and (11). This demonstrates the validity of our parameter setting rule, which means that this rule might be consistent with travelers’ behavior.

3. Target-Oriented Multi-Attribute User Equilibrium

The main content in this section is to examine the long-run impact of the route choice model proposed in last section, and propose the target-oriented multi-attribute user equilibrium (ToMaUE) considering LAP, travel cost, and travel time, which can be seen as an expansion of Wardrop user equilibrium.

3.1. Equilibrium Condition

With the above consideration, a user equilibrium will be reached in the long run where no traveler can unilaterally change route to increase their ToMaTU, termed as ToMaUE. With the notations and definitions described in Section 2.2, ToMaUE can be stated as follows, where π r s denote the maximum ToMaTU for O D pair r s .
Definition 2. 
(ToMaUE) It is assumed that all travelers choose the route with the purpose of maximizing their ToMaTU on the traffic network, i.e., considering any O D pair, ToMaTUs of the routes with flows are equal and maximum, which are more than or equal to those of the routes without flows.
Under ToMaUE, the following conditions need to be satisfied:
π r s ω p r s ( f ) = { = 0   if   ( f p r s ) > 0 0   if   ( f p r s ) = 0 s S ,   r R ,   p P r s  
where the equilibrium route flows are denoted by the tag *.
Furthermore, we reformulate the aforementioned ToMaUE conditions as the nonlinear complementarity problem (NCP), as shown in the following:
q r s = p P r s f p r s , r R , s S
v a = r R s S p P r s δ p a r s f p r s , a A
f p r s 0 , p P r s , r R , s S
f p r s π r s ω p r s f * = 0 , p P r s , r R , s S
π r s ω p r s ( f ) 0 ,   p P r s ,   r R ,   s S
In the above NCP system, Equation (14) describes that for each O D pair r s , the total travel demand is equal to the sum of flows on the routes that connect this O D pair, i.e., the travel demand conservation. Equation (15) describes that the link flow is the sum of flows on the routes that use this link. Equation (16) describes that all the route flows are non-negative. Equations (14)–(16) show the feasible flow set of given traffic network, denoted as Ω.
Therefore, a variational inequality (VI) problem V I ( f , Ω ) is obtained by reformulating the ToMaUE conditions, i.e., the equilibrium route flow vector f Ω needs to satisfy the inequation ω ^ ( f ) T ( f f ) 0 , f Ω , where ω ^ ( f ) = ω ( f ) . In addition, the properties of the VI formulation have been well studied in the literature, and one can refer to [44] for more information. Here, we discuss three properties of our VI problem, namely the equivalence between our VI problem and the ToMaUE conditions, the solution existence and solution uniqueness.
Th equivalence between our VI problem and the ToMaUE conditions (14)–(18) can be validated by examining the optimality condition of the linear programming problem min f Ω ω ^ ( f * ) T f . As aforementioned, we learn that the probability distribution function about the travel time on route is parameterized by f, and we assume that this function is continuous with respect to f. Moreover, the probability distribution function of route LAP is also continuous with respect to f, based on Equation (6). Therefore, we obtain that the ToMaTU is continuous with respect to f. Considering that Ω is compact and convex, there exists at least one solution to the above VI problem [45]. The solution uniqueness requires that the mapping in the VI problem is strictly monotone. Due to the complex formulation of the ToMaTU, ω ^ ( f ) might not be strictly monotone on Ω. In particular, most user equilibrium models may not meet this requirement in the presence of uncertainty (e.g., [4,46]). Thus, to determine a reasonable equilibrium solution from the set (Ω) of the feasible flow is a promising research direction, and it is also an open issue that has received sustained attention to the best of our knowledge.

3.2. Solution Algorithm

As shown in the obtained VI problem, the ToMaTU condition is non-additive, i.e., the VI problem follows the route-based form. In order to address it, we intend to utilize the approach of successive average (e.g., [47]), which is shown as follows. During the solution process, route enumeration is needed due to this non-additive issue, which causes high computational cost, especially for the large-scale complex network, whereas the objective of our study here is to provide insights into target-oriented route choice decision of travelers, considering these three attributes, and to design a computationally efficient algorithm is beyond our scope, which is research reserved for the future.
  • Stage 0: Initialization. Set the value of the convergence tolerance as ϵ 0 , which means the maximum number of iterations, as well as k = 1 .
  • Stage 1: All-Or-Nothing assignment. Conduct all-or-nothing assignment in terms of the present vacant link flow v a 1 , as well as obtain the flow f p r s , 1 of route connecting each O D pair.
  • Stage 2: Link flow, route travel time distribution and route LAP distribution updating. According to the present flow f p r s , k of route, orderly update the flow v a of its link, the route travel time distribution, and the route LAP distribution.
  • Stage 3: Calculation for ToMaTU. Acquire the value of ToMaTU ω p k , r s based on Equation (10) or (11), and obtain the maximum route utility π k , r s .
  • Stage 4: Check the Convergence. Set ϵ = ϵ + q r s π k , r s q f q r s ω p k , r s q r s π k , r s . Stop the algorithm if ϵ ϵ 0 or the maximum number of iterations is reached, if not, go to Stage 5.
  • Stage 5: Route flow updating. Detect the direction of search ( d p k , r s ) and the step size ( s k ) judging from the present values, and the flow of route is updated as the rule f p r s , k + 1 = f p k , r s + s k d p k , r s , and go the Stage 2. Set the direction of search as d p k , r s = f ˜ p k , r s f p k , r s , where f ˜ p k , r s is called as auxiliary flow. If ω p k , r s = π k , r s , f ˜ p k , r s = q r s / m r s with m r s being the number of routes that have the maximum value of ToMaRUs in step k for O D pair r s , and otherwise, f ˜ p k , r s is 0. The step size is set as s k = 1 / ( k + 1 ) .

4. Numerical Analysis

The main content in this section is running proposed model on Braess traffic network, as well as the more general Nguyen and Dupuis’s network, examining the models’ performance and conducting sensitivity analysis of the parameters including these three targets and the target inter-action between them. The convergence tolerance is 10 6 , and the total number of iterations is limited to 10 6 .

4.1. Test on the Braess Traffic Network

Firstly, we apply our proposed model to a stochastic tolled Braess network, and its performance is shown in Figure 1. There are four nodes and five links on this network, and one O D pair ( 1 , 4 ) is connected by three routes. The total units of demand between this O D pair is 1500. Other characteristics of this network are summarized in Table 1, where ϕ reflects the link capacity degradation. Larger value of ϕ shows smaller extent of link capacity degradation. Route 1 consists of Link 1 and Link 2; Route 2 consists of Link 1, Link 3, and Link 5; and Route 3 consists of Link 4 and Link 5.
We conduct the sensitivity analysis based on a benchmark, where we assume θ = 0.95 (i.e., the risk-averse travelers), γ l = 5 , γ c = 5 (Route 1 and Route 3 can successfully achieve their target for travel cost, Route 2 cannot), α 1 = 3 , α 2 = 2 , and β b = β s = 1 (i.e., in the imperfect complementarity relationship). We solve the system of Equation (12) to obtain the target-oriented utilities, and obtain ξ 2 = 2 / 11 , ξ 3 = 3 / 11 , ξ 12 = 8 / 11 , and ξ 23 = 5 / 11 .

4.1.1. Sensitivity Analysis via Changing γ c

In this test, we compare three cases by changing the value of γ c , which are Case 1 ( γ c = 4 ), Case 2 ( γ c = 5 ), and Case 3 ( γ c = 6 ), and other parameter values follow the settings in the benchmark. Target for travel cost can be achieved only on Route 1 in Case 1, can be achieved on Route 1 and 3 in Case 2, and can be achieved on all routes in Case 3. That is, we aim to show the impact of target for travel cost. The equilibrium route utility is 0.6994 in Case 1, is 0.6997 in Case 2, and is 0.97 in Case 3, and the equilibrium route flows are shown in Figure 2a. From this figure, we see that, as the value of target for travel cost increases, i.e., this target can be achieved on more and more routes, and flows will shift from the route where this target cannot be achieved to the routes where this target can be achieved. Meanwhile, the larger target value of travel cost, the higher utility value of the equilibrium route, which is also due to the target achievement on more routes. However, the extent of this increase depends on the specific value of the target γ c , e.g., in our testing, the equilibrium route utility almost remains the same in Case 1 and Case 2.
We further discuss the changes of TAPt and TAPl by changing the value of γ c , which is shown in Figure 2b,c, respectively. From these figures, we find that the value of TAPt is about 0.35, when the target set for travel cost can be successfully reached. That is, although travelers are risk-averse, namely the on-time arrival probability is 0.95, they will become risk-seeking to achieve three targets simultaneously. We also see that TAPl on this kind of route is also not very large, around 0.6, which means that travelers bear the risk of violating the allowable delay of the company to achieve three targets simultaneously. Meanwhile, we find that the value of TAPt and TAPl is 0.95, when the target set for travel cost cannot be reached, i.e., travelers’ on-time arrival probability, and almost 1, respectively. Nevertheless, route flows on these routes are not very large, as shown in Figure 2a. By properly increasing the value of β b and β s , we can also run this testing in the situation that involves the perfect complementarity relationship, and obtain the similar trends as aforementioned.

4.1.2. Sensitivity Analysis via Changing θ

In this test, we show the impact of travelers’ on-time arrival probability by changing its value from 0.5 (i.e., the risk-neutral travelers) to 0.95, and other parameter values follow the settings in the benchmark. The value of equilibrium route utility always increases as the value of θ grows, and the minimum and maximum values are 0.45 and 0.7, respectively, i.e., the emergence of risk-aversion increases travelers’ equilibrium route utility. Next, we focus on the change of equilibrium route flows, which are shown in Figure 3a. From this figure, we find that the increase in θ will result in flows transforming from route 1 and 3 to route 2 initially, and then flows will shift back when the value of θ approaches 0.95. The reason for this is that a larger value of θ means larger budgeted travel time, i.e., larger target γ t , and thus travelers can obtain certain utility, which even holds up for routes where the target set for travel cost cannot be reached. Therefore, travelers will incline to select route 2. Nevertheless, when the value of θ becomes large enough, e.g., around 0.95, TAPl on route 2 becomes 1, i.e., the target set for LAP of route is large enough, and cannot increase any more, as shown in Figure 3c, and thus more travelers will choose route 1 and route 3 where the target set for travel cost of route could be reached for a larger utility.
We further discuss the changes of TAPt and TAPl by changing the value of θ , which is shown in Figure 3b,c, respectively. We can see that, with the increase in θ , TAPt and TAPl both increase. For some routes, where target for travel cost could be reached successfully, TAPt is less than travelers’ on-time arrival probability, and TAPl is also not very large, i.e., below 0.8, while for some routes where the target set for travel cost cannot be reached, TAPt is always travelers’ on-time arrival probability, and TAPl is close to 1. That is, in order to achieve the three targets simultaneously, travelers will become less risk-averse or even risk-seeking, and will bear the risk of violating the allowable delay of the company. By properly increasing the value of βb and β s , we can also run this testing in the situation that involves the perfect complementarity relationship, and obtain the similar trends as aforementioned.

4.1.3. Sensitivity Analysis via Changing γ l

In this test, we show the impact of travelers’ target for LAP by changing its value from 0 to 10, and other parameter values follow the settings in the benchmark. The value of equilibrium route utility keeps increasing, and the minimum and maximum values are 0.69 and 0.70, respectively, i.e., increase in the value of γ l brings little utility to the travelers. Next, we focus on the change of equilibrium route flows, which are shown in Figure 4a. From this figure, we find that the increase in γ l will result in flows transforming from route 2 to route 1 and 3. As shown in Figure 4b,c, for route 2, as the value of γ l increases, the value of TAPt is always 0.95, and the value of TAPI quickly converges to 1, while TAPl on route 1 and route 3 become larger and larger, i.e., route utility on route 1 and route 3 become larger and larger. Therefore, travelers will incline to select route 1 and 3, and TAPt on these routes will decrease.
We further discuss the changes of TAPt and TAPl by changing the value of γ l , which is shown in Figure 4b,c, respectively. From these figures, we find that the value of TAPt of the routes where the target set for travel cost cannot be reached, is always 0.95, i.e., travelers’ on-time arrival probability, and TAPl on this kind of route quickly becomes 1. Meanwhile, we see that on some routes where the target set for travel cost can be successfully reached, TAPl will become larger and larger as the value of γ l increases. In other words, to achieve the three targets simultaneously, travelers tend to transform into less risk-averse or even risk-seeking, and an increase in the value of γ l will decrease travelers’ risk of violating the allowable delay of the company, i.e., travelers will arrive within the allowable delay of the company with a high probability. By properly increasing the value of β b and β s , we can also run this testing in the situation that involves the perfect complementarity relationship, and obtain the similar trends as aforementioned.

4.1.4. Sensitivity Analysis via Changing α 1 and α 2

In this test, we show the impact of changing the value of α1 and α 2 , i.e., changing the utility ratios or the relative importance between different target achievements, and other parameter values follow the settings in the benchmark. The value of α 1 is changed from 2 to 4, and the value of α 2 is changed from 1 to 3. Given the specific values of α 1 and α 2 , we can calculate the value of all target-oriented utilities after addressing equation system (12).
The change of equilibrium route utility is shown in Figure 5, where its minimum value is 0.58, and its maximum value is 0.76. From this figure, we find that the value of equilibrium route utility always decreases with the value of α 1 grows, while the value of equilibrium route utility always increases with the value of α 2 grows. Figure 5 shows the impact of different target achievement on the equilibrium route utility. For example, recalling the definition of utility ratios, increase in the value of α 1 shows that achieving the targets for the travel cost and travel time on route deserves more attention, in contrast to the achievement of target set for LAP on route, while increase in the value of α 2 shows that achieving the targets for the LAP and travel time deserves more attention, in contrast to the achievement of target set for travel cost on route. We also see that the utility change is moderate as the value of α 1 increases, and this change is relatively large as the value of α 2 increases.
The change of equilibrium route flows is presented in Figure 6a–c. From these figures, we can find that the increase in α 1 will result in flows transforming from route 2 to route 1 and 3, as the achievement of target set for travel cost on route deserves more attention as aforementioned, while the increase in α 2 will result in flows transforming from route 1 and 3 to route 2, as the achievement of target set for travel cost on route becomes rarely noticed as aforementioned, recalling that the target set for travel cost fails to be achieved on route 2.
We further discuss the changes of TAPt and TAPl by changing the value of α 1 and α 2 , which is shown in Figure 7a–c and Figure 8a–c. From these figures, we see that as the value of α 1 increases, the value of TAPt and TAPl on routes where the target set for travel cost can be reached successfully, will decrease, while the value of TAPt and TAPl on these routes always decreases with the value of α 2 grows. We also see that on these routes, travelers will transform into less risk-averse and are very likely to arrive within the allowable delay of the company when α 2 is larger, and travelers will become risk-seeking and bear the risk of violating the allowable delay when α 2 is smaller. However, for some routes where the target set for travel cost fails to be reached, we see that TAPt is always 0.95, i.e., the probability that travelers will arrive on time, as well as the value of TAPl is almost 1. By properly increasing the value of β b and β s , we can also run this testing in the situation that involves the perfect complementarity relationship, and obtain the similar trends as aforementioned.

4.1.5. Sensitivity Analysis via Changing β b and β s

In this test, we show the impact of changing the value of complementarity ratios β b and β s , i.e., we compare the imperfect complementarity relationship with the perfect complementarity relationship, and other parameter values follow the settings in the benchmark. The values of β b is changed from 1.5 + δ to 3, and the value of β s is changed from 1 to 2, where δ is a positive with small value. It can be verified that these settings of β b and β s satisfy the aforementioned relationship. The change of equilibrium route utility is shown in Figure 9, where its minimum value is 0.23, and its maximum value is 0.47. From this figure, and combining the equilibrium route utility shown in the above tests, we see that the value of equilibrium route utility on routes always decreases as the value of β b grows, and the value of equilibrium route utility on routes also decreases as the value of β s grows, but this decrease is smaller. Figure 9 shows the impact of different target interactions on the equilibrium route utility. For example, the interaction among three targets, i.e., changing the value of βb, can bring fundamental changes to the equilibrium route utility, while the interaction between two targets, i.e., changing the value of β s , can only bring relatively smaller changes to the equilibrium route utility.
The change of equilibrium route flows is presented in Figure 10a–c. From these figures, we can find that the increase in β b will result in flows transform from route 2 to route 1 and 3, as achieving the three targets simultaneously brings more utility to the travelers, i.e., travelers are willing to realize these three targets at the same time, while the increase in βs will result in flows transform from route 1 and 3 to route 2, as interaction between two targets can weaken the impact of interaction among three targets, or weaken travelers’ desire to achieve the three targets at the same time as shown in Equation (7).
We further discuss the changes of TAPt and TAPl by changing the value of β b and β s , which is shown in Figure 11a–c and Figure 12a–c, respectively. From these figures, we see that the value of TAPt and TAPl on routes where the target set for travel cost can be reached successfully will decrease, while the value of β b grows, i.e., travelers tend to transform into more risk-seeking and bear the risk of violating the allowable delay of the company to achieve the three targets simultaneously, while the value of TAPt and TAPl on these routes always increases as the value of β s grows, i.e., interaction between two targets can weaken the impact of interaction among three targets. However, we see that TAPt or TAPI on the routes where the target set for travel cost fails to be reached is always 0.95, i.e., travelers’ on-time arrival probability, and almost 1, respectively.

4.2. Test on the Nguyen and Dupuis’s Traffic Network

Finally, we apply our proposed model on a Nguyen and Dupuis’s traffic network [48] shown in Figure 13, which can be seen as a more general network structure with several O D pairs. This traffic network consists of 13 nodes and 19 links. In this section, four OD pairs are taken into account, including 1 2 , 1 3 , 4 2 , and 4 3 , and the deterministic demands for all the O D pairs are all 800. Other characteristics of the network are shown in Table 2, where FFTT denotes free-flow travel time, and ϕ also reflects the link capacity degradation. Travelers’ on-time arrival probability is 0.95, the target for route LAP is 5, and the target set for travel cost of route is 31.
We test two situations, namely in the imperfect complementarity relationship, where βb = βs = 1, and in the perfect complementarity relationship, where β b = 2 and β s = 2 , α 1 = 3 , and α 2 = 2 . The convergence of our model is shown in Figure 14, where convergence of the first situation is shown in the left figure, and the convergence of the second situation is shown in the right figure. Specially, in order to better illustrate, we only display the first 500 iterations. As shown in these figures, we can find that the applicability and the performance of our proposed model, when it is embedded in a more general traffic network. As aforementioned, similar testing results can be received by properly changing the values of the relevant parameters, and we omit them for brevity.

4.3. Insights of the Sensitivity Analysis

The main insights of the sensitivity analysis into route choice decision of travelers in the stochastic tolled network are summarized as follows.
(1)
For some routes, where the target set for travel cost can be successfully reached (depending on travelers’ acceptable expense), target achievement probability for travel time is smaller than travelers’ on-time arrival probability, which means travelers become less risk-averse or even risk-seeking in order to achieve the three targets simultaneously. Meanwhile, on these kinds of routes, target achievement probability for late arrival penalty is also not very large, unless the allowable delay of the company is large; the relative significance between the achievements of the target set for travel cost and travel time is relatively large, or the extent of target interaction among three targets is small. That is, travelers usually bear the risk of violating the allowable delay of the company in order to achieve the three targets simultaneously, unless the aforementioned situations happen.
(2)
For some routes, where the target set for travel cost fails to be reached (depending on travelers’ acceptable expense), the probability of achieving the target for travel time is travelers’ on-time arrival probability, and the probability of achieving the target for late arrival penalty is almost 1, which are always valid in all the testing.
(3)
The equilibrium route flow distribution changes when the values of the three targets, the relative importance between different target achievements, and the target interaction change. Usually, travelers prefer the routes where the target set for travel cost can be successfully reached, which can be impacted by their risk-averse attitude, and the allowable delay of the company, but our testing shows that these impacts are not significant. If achieving the target for travel cost takes on greater importance, flows will shift from the routes where this target cannot be achieved to the routes where this target can be achieved in general, but the target achievements of travel time and late arrival penalty have impact on this flow shift. Meanwhile, the impact of the former one is larger. When the extent of the interaction among three targets become stronger, more travelers will choose the routes where the three targets can be achieved simultaneously, but stronger extent of the interaction between two targets could weaken this impact.
(4)
The equilibrium route utility changes when the values of the three targets, the relative importance between different target achievements, and the target interaction change. From our testing, travelers’ larger acceptable expense and stronger extent of risk-averse attitude leads to the larger equilibrium route utility, but the allowable delay of the company has almost no impact on this. When the relative significance between the achievements of the target set for travel cost and travel time becomes larger, the equilibrium route utility is larger, but the target interaction, among three targets or between two targets, can decrease the equilibrium route utility.

5. Conclusions

In this paper, we studied the target-oriented multi-attribute route choice decision of travelers in the stochastic tolled traffic network, considering the impact of three attributes, which are (stochastic) travel time, (stochastic) late arrival penalty, and (deterministic) travel cost. In particular, we proposed a new route choice model, termed as a target-oriented multi-attribute travel utility model, for this analysis based on our recent methodology. In this new model, each attribute is assigned a target by travelers, and travelers’ objective is to maximize their travel utility that is determined by the achieved targets. Moreover, the interaction between targets is interpreted as complementarity relationship between them, which can further affect their travel utility. Additionally, based on the travel utility model, target-oriented multi-attribute user equilibrium model is proposed, which is formulated as a variational inequality problem and solved with the method of successive average. Target for travel time is determined via travelers’ on-time arrival probability, while targets for late arrival penalty and travel cost are given exogenously. Lastly, we apply the proposed model to the Braess and Nguyen–Dupuis traffic networks, and conduct sensitivity analysis of the parameters, including these three targets and the target inter-action between them.
There are several directions that merit the further study. (1) In this paper, tolls on the traffic network are given, and we will study the optimal congestion tolling on the basis of our proposed model at a later stage. (2) We will also adopt more behaviorally consistent methods to determine the value of the target set for late arrival penalty, travel cost, and travel time at a later stage.

Author Contributions

Conceptualization, X.Z.; methodology, X.Z.; formal analysis, X.Z., Z.G. and J.M.; writing—original draft preparation, X.Z.; writing—review and editing, Y.Z. and X.J.; funding acquisition, Y.Z. supervision, Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (71974104).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Special thanks go to the anonymous reviewers for their suggestions which improve the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In the work of [40], link capacity is assumed to follow a uniform distribution. Upper bound of the capacity is the designed capacity c ¯ a and lower bound of the capacity is the product of fraction ϕ a and the designed capacity c ¯ a . Therefore, the distribution of link capacity is written as
C a U ( ϕ a c ¯ a , c ¯ a )
In Equation (A1), the link capacity degradation is reflected by ϕ a . The smaller the value of ϕa is, the less reliable the link is. If link capacity follows the uniform distribution constructed as above, route travel time will follow normal distribution [40], i.e.,
T p r s N ( E ( T p r s ) , σ ( T p r s ) )
where its mean and standard deviation are written as
E ( T p r s ) = a ( δ p a r s ( t a 0 + β t a 0 x a n 1 ϕ a 1 n c ¯ a n ( 1 ϕ a ) ( 1 n ) ) )
and
σ T p r s = α δ p a r s · β 2 t a 0 2 x a 2 n 1 ϕ a 1 2 n c ¯ a 2 n 1 ϕ a ( 1 2 n ) 1 ϕ a 1 n c ¯ a n 1 ϕ a ( 1 n ) 2

References

  1. Umar, M.; Ji, X.; Kirikkaleli, D.; Xu, Q. COP21 Roadmap: Do Innovation, Financial Development, and Transportation Infrastructure Matter for Environmental Sustainability in China? J. Environ. Manag. 2020, 271, 111026. [Google Scholar] [CrossRef]
  2. Wardrop, J.G. Road paper. Some theoretical aspects of road traffic research. Proc. Inst. Civ. Eng. 1952, 1, 325–362. [Google Scholar] [CrossRef]
  3. Xu, H.; Lou, Y.; Yin, Y.; Zhou, J. A prospect-based user equilibrium model with endogenous reference points and its application in congestion pricing. Transp. Res. Part B-Methodol. 2011, 45, 311–328. [Google Scholar] [CrossRef]
  4. Chen, A.; Zhou, Z. The α-Reliable mean-excess traffic equilibrium model with stochastic travel times. Transp. Res. Part B-Methodol. 2010, 44, 493–513. [Google Scholar] [CrossRef]
  5. Ji, X.; Ban, X.; Zhang, J.; Ran, B. Moment-based travel time reliability assessment with Lasserre’s relaxation. Transp. B-Transp. Dyn. 2019, 7, 401–422. [Google Scholar] [CrossRef]
  6. Khight, F.H. Risk, Uncertainty, and Profit; Houghton Mifflin: Boston, MA, USA, 1921. [Google Scholar]
  7. Ji, X.; Ban, X.; Zhang, J.; Ran, B. Subjective-utility travel time budget modeling in the stochastic traffic network assignment. J. Intell. Transport. Syst. 2017, 21, 439–451. [Google Scholar] [CrossRef]
  8. Ji, X.; Ban, X.; Li, M.; Zhang, J.; Ran, B. Non-expected route choice model under risk on stochastic traffic networks. Netw. Spat. Econ. 2017, 17, 777–807. [Google Scholar] [CrossRef]
  9. Li, X.; Natarajan, K.; Teo, C.-P.; Zheng, Z. Distributionally robust mixed integer linear programs: Persistency models with applications. Eur. J. Oper. Res. 2014, 233, 459–473. [Google Scholar] [CrossRef]
  10. Zhang, Y.; Song, S.; Shen, Z.-J.M.; Wu, C. Robust shortest path problem with distributional uncertainty. IEEE Trans. Intell. Transp. Syst. 2018, 19, 1080–1090. [Google Scholar] [CrossRef]
  11. Wang, Z.; You, K.; Song, S.; Zhang, Y. Wasserstein distributionally robust shortest path problem. Eur. J. Oper. Res. 2020, 284, 31–43. [Google Scholar] [CrossRef] [Green Version]
  12. Lo, H.K.; Luo, X.W.; Siu, B.W.Y. Degradable transport network: Travel time budget of travelers with heterogeneous risk aversion. Transp. Res. Part B-Methodol. 2006, 40, 792–806. [Google Scholar] [CrossRef]
  13. Abdel-Aty, M.; Kitamura, R.; Jovanis, P. Investigating effect of travel time variability on route choice using repeated measurement stated preference. Data. Transp. Res. Rec. 1995, 1493, 39–45. [Google Scholar]
  14. Brownstone, D.; Ghosh, A.; Golob, T.F.; Kazimi, C.; Van Amelsfort, D. Drivers’ willingness-to-pay to reduce travel time: Evidence from the San Diego I-15 Congestion Pricing Project. Transp. Res. Part A-Policy Pract. 2003, 37, 373–387. [Google Scholar] [CrossRef] [Green Version]
  15. Ji, X.; Ao, X. Travelers’ bi-attribute decision making on the risky mode choice with flow-dependent salience theory. Sustainability 2021, 13, 3901. [Google Scholar] [CrossRef]
  16. Xu, Q.; Ji, X. User equilibrium analysis considering travelers’ context-dependent route choice behavior on the risky traffic network. Sustainability 2020, 12, 6706. [Google Scholar] [CrossRef]
  17. Yang, H.; Huang, H.-J. The multi-class, multi-criteria traffic network equilibrium and systems optimum problem. Transp. Res. Part B-Methodol. 2004, 38, 1–15. [Google Scholar] [CrossRef]
  18. Watling, D. User equilibrium traffic network assignment with stochastic travel times and late arrival penalty. Eur. J. Oper. Res. 2006, 175, 1539–1556. [Google Scholar] [CrossRef] [Green Version]
  19. Wang, J.Y.T.; Ehrgott, M.; Chen, A. A bi-objective user equilibrium model of travel time reliability in a road network. Transp. Res. Part B-Methodol. 2014, 66, 4–15. [Google Scholar] [CrossRef] [Green Version]
  20. Ehrgott, M.; Wang, J.Y.T.; Watling, D.P. On multi-objective stochastic user equilibrium. Transp. Res. Part B-Methodol. 2015, 81, 704–717. [Google Scholar] [CrossRef] [Green Version]
  21. Sun, C.; Cheng, L.; Zhu, S.; Han, F.; Chu, Z. Multi-criteria user equilibrium model considering travel time, travel time reliability and distance. Transport. Res. Part D-Transport. Environ. 2019, 66, 3–12. [Google Scholar] [CrossRef]
  22. Kahneman, D.; Tversky, A. Prospect theory: An analysis of decision under risk. Econometrica 1979, 47, 263–291. [Google Scholar] [CrossRef] [Green Version]
  23. Tversky, A.; Kahneman, D. Advances in prospect theory: Cumulative representation of uncertainty. J. Risk Uncertain. 1992, 5, 297–323. [Google Scholar] [CrossRef]
  24. Kemel, E.; Paraschiv, C. Prospect theory for joint time and money consequences in risk and ambiguity. Transp. Res. Part B-Methodol. 2013, 56, 81–95. [Google Scholar] [CrossRef]
  25. Avineri, E. The effect of reference point on stochastic network equilibrium. Transp. Sci. 2006, 40, 409–420. [Google Scholar] [CrossRef] [Green Version]
  26. Schwanen, T.; Ettema, D. Coping with unreliable transportation when collecting children: Examining parents’ behavior with cumulative prospect theory. Transp. Res. Part A-Policy Pract. 2009, 43, 511–525. [Google Scholar] [CrossRef]
  27. Zhang, C.; Liu, T.-L.; Huang, H.-J.; Chen, J. A cumulative prospect theory approach to commuters’ day-to-day route-choice modeling with friends’ travel information. Transp. Res. Part C-Emerg. Technol. 2018, 86, 527–548. [Google Scholar] [CrossRef]
  28. Ghader, S.; Darzi, A.; Zhang, L. Modeling effects of travel time reliability on mode choice using cumulative prospect theory. Transp. Res. Part C-Emerg. Technol. 2019, 108, 245–254. [Google Scholar] [CrossRef]
  29. Hu, L.; Dong, J.; Lin, Z. Modeling charging behavior of battery electric vehicle drivers: A cumulative prospect theory based approach. Transp. Res. Part C-Emerg. Technol. 2019, 102, 474–489. [Google Scholar] [CrossRef] [Green Version]
  30. Ji, X.; Chu, Y. A target-oriented bi-attribute user equilibrium model with travelers’ perception errors on the tolled traffic network. Transp. Res. Part E-Logist. Transp. Rev. 2020, 144, 102150. [Google Scholar] [CrossRef]
  31. Ji, X. User-equilibrium-based investigation on the usage of tolled routes: A target-oriented perspective. Transp. Res. Part E-Logist. Transp. Rev. 2021. under review. [Google Scholar]
  32. Gao, S.; Frejinger, E.; Ben-Akiva, M. Adaptive route choices in risky traffic networks: A prospect theory approach. Transp. Res. Part C-Emerg. Technol. 2010, 18, 727–740. [Google Scholar] [CrossRef]
  33. Boyles, S.D.; Waller, S.T. Optimal information location for adaptive routing. Netw. Spat. Econ. 2011, 11, 233–254. [Google Scholar] [CrossRef]
  34. Gao, S.; Huang, H. Real-time traveler information for optimal adaptive routing in stochastic time-dependent networks. Transp. Res. Part C-Emerg. Technol. 2012, 21, 196–213. [Google Scholar] [CrossRef]
  35. Miller-Hooks, E. Adaptive least-expected time paths in stochastic, time-varying transportation and data networks. Networks 2001, 37, 35–52. [Google Scholar] [CrossRef]
  36. Unnikrishnan, A.; Waller, S.T. User equilibrium with recourse. Netw. Spat. Econ. 2009, 9, 575–593. [Google Scholar] [CrossRef]
  37. Unnikrishnan, A.; Lin, D.-Y. User equilibrium with recourse: Continuous network design problem. Comput-Aided Civ. Inf. 2012, 27, 512–524. [Google Scholar] [CrossRef]
  38. Wijayaratna, K.P.; Dixit, V.V.; Denant-Boemont, L.; Waller, S.T. An experimental study of the online information paradox: Does en-route information improve road network performance? PLoS ONE 2017, 12, e0184191. [Google Scholar] [CrossRef] [Green Version]
  39. Wijayaratna, K.P.; Dixit, V.V. Impact of information on risk attitudes: Implications on valuation of reliability and information. J. Choice Model. 2016, 20, 16–34. [Google Scholar] [CrossRef]
  40. Lo, H.K.; Tung, Y.-K. Network with degradable links: Capacity analysis and design. Transp. Res. Part B-Methodol. 2003, 37, 345–363. [Google Scholar] [CrossRef]
  41. Shao, H.; Lam, W.H.K.; Tam, M.L. A reliability-based stochastic traffic assignment model for network with multiple user classes under uncertainty in demand. Netw. Spat. Econ. 2006, 6, 173–204. [Google Scholar] [CrossRef]
  42. Ji, X.; Cheng, H.; Ao, X.; Zhou, X. Modeling travel time reliability and unreliability on the stochastic traffic network with target-oriented analysis. Ann. Oper. Res. 2021. under review. [Google Scholar]
  43. Isac, G.; Bulavsky, V.A.; Kalashnikov, V.V. Complementarity, Equilibrium, Efficiency and Economics; Springer: Boston, MA, USA, 2013. [Google Scholar]
  44. Nagurney, A. Network Economics: A Variational Inequality Approach; Springer: Boston, MA, USA, 2013. [Google Scholar]
  45. Facchinei, F.; Pang, J.-S. Finite-Dimensional Variational Inequalities and Complementarity Problems; Springer Series in Operations Research and Financial Engineering; Springer: Boston, MA, USA, 2003. [Google Scholar]
  46. Nikolova, E.; Stier-Moses, N. A mean-risk model for the traffic assignment problem with stochastic travel times. Oper. Res. 2014, 62, 366–382. [Google Scholar] [CrossRef] [Green Version]
  47. Sheffi, Y. Urban Transportation Networks: Equilibrium Analysis With Mathematical Programming Methods; Prentice-Hall: Englewood Cliffs, NJ, USA, 1984. [Google Scholar]
  48. Nguyen, S.; Dupuis, C. An efficient method for computing traffic equilibria in networks with asymmetric transportation costs. Transp. Sci. 1984, 18, 185–202. [Google Scholar] [CrossRef]
Figure 1. Braess traffic network.
Figure 1. Braess traffic network.
Sustainability 13 09992 g001
Figure 2. Testing results of changing γ c : (a) Equilibrium route utility; (b) TAPt; and (c) TAPl.
Figure 2. Testing results of changing γ c : (a) Equilibrium route utility; (b) TAPt; and (c) TAPl.
Sustainability 13 09992 g002
Figure 3. Testing results of changing θ: (a) Equilibrium route utility; (b) TAPt; and (c) TAPl.
Figure 3. Testing results of changing θ: (a) Equilibrium route utility; (b) TAPt; and (c) TAPl.
Sustainability 13 09992 g003
Figure 4. Testing results of changing γl: (a) Equilibrium route utility; (b) TAPt; and (c) TAPl.
Figure 4. Testing results of changing γl: (a) Equilibrium route utility; (b) TAPt; and (c) TAPl.
Sustainability 13 09992 g004
Figure 5. Testing results of equilibrium route utility with regard to changing α 1 and α 2 .
Figure 5. Testing results of equilibrium route utility with regard to changing α 1 and α 2 .
Sustainability 13 09992 g005
Figure 6. Testing results of equilibrium flow with regard to changing α 1 and α 2 : (a) route 1; (b) route 2; and (c) route 3.
Figure 6. Testing results of equilibrium flow with regard to changing α 1 and α 2 : (a) route 1; (b) route 2; and (c) route 3.
Sustainability 13 09992 g006
Figure 7. Testing results of TAPt with regard to changing α 1 and α 2 : (a) route 1; (b) route 2; and (c) route 3.
Figure 7. Testing results of TAPt with regard to changing α 1 and α 2 : (a) route 1; (b) route 2; and (c) route 3.
Sustainability 13 09992 g007
Figure 8. Testing results of TAPI with regard to changing α 1 and α 2 : (a) route 1; (b) route 2; and (c) route 3.
Figure 8. Testing results of TAPI with regard to changing α 1 and α 2 : (a) route 1; (b) route 2; and (c) route 3.
Sustainability 13 09992 g008
Figure 9. Testing results of equilibrium route utility with regard to changing β b and β s .
Figure 9. Testing results of equilibrium route utility with regard to changing β b and β s .
Sustainability 13 09992 g009
Figure 10. Testing results of equilibrium flow with regard to changing βb and β s : (a) route 1; (b) route 2; and (c) route 3.
Figure 10. Testing results of equilibrium flow with regard to changing βb and β s : (a) route 1; (b) route 2; and (c) route 3.
Sustainability 13 09992 g010
Figure 11. Testing results of TAPt with regard to changing β b and β s : (a) route 1; (b) route 2; and (c) route 3.
Figure 11. Testing results of TAPt with regard to changing β b and β s : (a) route 1; (b) route 2; and (c) route 3.
Sustainability 13 09992 g011
Figure 12. Testing results of TAPI with regard to changing β b and β s : (a) route 1; (b) route 2; and (c) route 3.
Figure 12. Testing results of TAPI with regard to changing β b and β s : (a) route 1; (b) route 2; and (c) route 3.
Sustainability 13 09992 g012
Figure 13. Nguyen and Dupuis’s traffic network.
Figure 13. Nguyen and Dupuis’s traffic network.
Sustainability 13 09992 g013
Figure 14. Convergence performance: (a) imperfect complementarity relationship; and (b) perfect complementarity relationship.
Figure 14. Convergence performance: (a) imperfect complementarity relationship; and (b) perfect complementarity relationship.
Sustainability 13 09992 g014
Table 1. Testing data on Braess traffic network.
Table 1. Testing data on Braess traffic network.
Link NumberFree-Flow Travel TimeCapacityTollϕ
1560020.8
21240020.7
3740020.9
41040030.7
5860020.8
Table 2. Characteristics of the Nguyen and Dupuis’s traffic network.
Table 2. Characteristics of the Nguyen and Dupuis’s traffic network.
Link NumberFree-Flow Travel TimeCapacityTollϕ
156003.50.8
21240080.9
3740050.8
41040060.8
5860060.7
6650040.8
71050070.6
81040080.8
9860050.8
10650040.9
111240090.9
121160080.9
13660040.7
14850050.8
15950060.8
161040060.9
1712400100.8
18860050.9
191050070.7
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, X.; Guo, Z.; Ma, J.; Zhong, Y.; Ji, X. Target-Oriented User Equilibrium Considering Travel Time, Late Arrival Penalty, and Travel Cost on the Stochastic Tolled Traffic Network. Sustainability 2021, 13, 9992. https://0-doi-org.brum.beds.ac.uk/10.3390/su13179992

AMA Style

Zang X, Guo Z, Ma J, Zhong Y, Ji X. Target-Oriented User Equilibrium Considering Travel Time, Late Arrival Penalty, and Travel Cost on the Stochastic Tolled Traffic Network. Sustainability. 2021; 13(17):9992. https://0-doi-org.brum.beds.ac.uk/10.3390/su13179992

Chicago/Turabian Style

Zang, Xinming, Zhenqi Guo, Jingai Ma, Yongguang Zhong, and Xiangfeng Ji. 2021. "Target-Oriented User Equilibrium Considering Travel Time, Late Arrival Penalty, and Travel Cost on the Stochastic Tolled Traffic Network" Sustainability 13, no. 17: 9992. https://0-doi-org.brum.beds.ac.uk/10.3390/su13179992

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