#### Abstract

This paper aims to discuss the trip mode choice problem by using cumulative prospect theory (CPT) rather than utility maximization from the network uncertainty perspective and evaluates the effect of the integrated service mode on taxi network equilibrium. The integrated service mode means taxis either are actively moving through traffic zones to pick up customers (cruising mode) or are queued at the center of a zone waiting for customers (dispatch mode). Based on this, CPT models are adopted to analyze the choice of customers’ trip mode. The travel time uncertainty of the network and the applicability of CPT are considered first, and the Nested Logit model was used to complete the trip mode split problem. Further, several relevant relationships including supply-demand equilibrium, network conditions, taxi behavior, and customer behavior perspectives were analyzed with respect to the integrated mode. Moreover, a network equilibrium model was established and its algorithm was designed. Finally, this paper presented a numerical example and discussed the taxi network equilibrium’s characteristic after introducing the integrated service mode.

#### 1. Introduction

Public transport via taxi has increasingly become an important service to modernity because of its personalized and flexible service, comfort, and available speed among other things. Despite its growing popularity, problems have emerged in the practical operation of the taxi market. For instance, it is well known that automobiles used in the taxi service are responsible for consuming a considerable amount of energy. This energy consumption is magnified under the conditions where there is a higher vacancy ratio in the taxi market. The pollution produced as a by-product of this industry is therefore much worse and seriously threatens energy sustainability as well as posing a serious risk to the environment. Based on this, two types of measures are needed in order to alleviate the negative impacts of such an important transportation service. To address the environmental impact of taxi trip, both the alternatively fueled vehicles and cleaner fuels have been proposed to provide taxi service [1]. Management policies which effectively advance the operational efficiency needed to be explored have also been proposed as a means of curbing the concerns associated with the industry. This investigation seeks to elucidate ways to improve the level of service and the operational situation through improvement to the management policies.

Different taxi policies have been heavily researched since the 1970s with current policies in the taxi field primarily concentrated on market regulations (i.e., entry restrictions and price control). Douglas, who first studied the aggregated model as a precursor in the realm of taxi networks, introduced the realm to economists that have proceeded to investigate the characteristics of different taxi markets, enabling them to understand the equilibrium between supply and demand in regulated markets [2–7]. Despite all the research, it is known and accepted that the taxi market is not the idealized market for conventional economic analyses since the spatial features of the road network are not considered [8]. From this standpoint, Yang and Wong et al. devised a series of models considering fluctuations in supply and demand over variations in time and space of Hong Kong for 1998–2011. In a previous attempt, Yang et al. sought to characterize taxi movements in a road network for a given origin-destination demand pattern, discussing the relationship between customer demand and taxi utilization by analyzing the impact of congestion, various markets, multiple user classes, stochastic travel time, and price formulation among other things [9–14]. Yang et al. [15] analyzed bilateral searching and meeting relationships between taxis and customers and studied fare control in competitive taxi markets. Further, by considering Pareto efficiency, Yang et al. [16] were able to evaluate quality of service and taxi utilization.

Despite all of the research over the years, models discussed typically account for the cruising service mode when determining regulations and the equilibrium problem; limited attention has been paid to other modes. In one of the first studies discussing multiple modes, Arnott [17] suggested that the dispatch mode was more effective in small cities and the cruising mode was more effective in large cities. von Massow and Canbolat [18] developed a model used to simulate the dispatch mode. One of the first to suggest this, Maria Salanora et al. [19], asserted that the integrated taxi service mode (e.g., cruising and dispatch modes) would be better in big cities.

Recently, in terms of prospect theory which assumes that lotteries are evaluated in a two-step process: an initial phase of editing and a subsequent phase of evaluation, a lot of achievements have been also published [20]. van de Kaa [21, 22] investigated the application of prospect theory and expected utility theory (EUT) where it is assumed that decision maker’s attitude toward risk can be rationalized by an expected utility function in trip choice behavior and concluded that, by extending the prospect theory further, it would be suitable to describe traveler behavior. Supporting this assumption, Li and Hensher [23] demonstrated that prospect theory was more suitable for researching travelers decision by analyzing the application of prospect theory in physics, behavior economics, and transportation fields and finally summarized the application constraints of prospect theory in transportation field.

By constructing and analyzing a supply-demand model, the impact factor on the market is evaluated as well as the taxis. The aim was to explore the essential characteristics of the taxi market and provide effective, valuable management policies based on the findings.

In this model, the expected utility theory was replaced with cumulative prospect theory to describe the split in customer trip mode. Cumulative prospect theory is regarded as a bounded rationality model reflective of the uncertainty and complication inherent to real road networks, expressing a decision maker’s attitude toward risk. Based on this, we consider the coexistence of cruising and dispatch modes to generate the model used for this paper. The equilibrium model describes the customer-search behavior of vacant-taxi drivers and the interactions between cruising and dispatch modes. Finally, parameters, such as customer-waiting time and taxi-searching time, were identified and evaluated based on variation in traffic zones and fleet sizes.

The rest of this paper is structured as follows: Section 2 is an analysis of the applicable characteristics of cumulative prospect theory and establishes the trip mode model for travelers based on CPT. Section 3 describes the relationship between cruising mode and dispatch modes presenting a supply-demand equilibrium relationship with customer demand variation in taxi market. Section 4 promotes a joint taxi network model of trip distribution and assignment. Section 5 explains the detailed steps of the algorithm developed to establish the supply-demand equilibrium mode. Section 6 describes the network conditions and parameters utilizing a numerical example and discusses the characteristic of equilibrium state with respect to the results. Finally, the last section makes a conclusion and proposes future extensions for the current taxi model.

#### 2. The Trip Mode Choice Model Based on CPT

It was assumed that there are two modes of taxi service available (i.e., cruising and dispatch) as well as other transit modes (bus, rail system, etc.) for customers on a given network. In terms of cruising mode, it means that taxi cruises on roadway to pick up next customers. In terms of dispatch mode, its implementation effect is greatly dependent on the installation rate of in-vehicle navigation system. However, there is a smaller popularizing rate of intelligent system in most current taxi markets. Thus, this paper adopts a simple dispatch policy which means that taxis were required to queue at the center of zone or the gravity point of zone to search for next customers by means of the achievements of von Massow and Canbolat [18]. Based on this assumption, the trip mode choice process may be broken down as follows.

From Figure 1, it can be seen that subsets of alternatives whose properties are similar are grouped in hierarchies or nests. The multinominal logit model (MNL) is used to estimate the probability of trip mode choice for the alternatives of the lower nest, including both cruising and dispatching service modes. And furthermore, at the higher nest, the MNL model consisting of composite taxi mode and other transit modes can be also estimated. Thus, the Nested Logit model is suitable to be used to study the trip mode choice problem described by this paper.

##### 2.1. The Descriptions of Travel Cost

The total travel time from zone to zone of customers involves the travel time on roadway and the waiting time which may be expressed as

In the preceding equations is the total time of cruising mode, is the total time of dispatch mode, and is the total time of other transit modes, while and are the travel time of taxis and is the travel time of other forms of transit from zone to zone . is the waiting time of customers that choose cruising mode in zone and is the waiting time of customers choosing cruising mode in zone , and refers to the waiting time of customers choosing other modes of transit in zone .

With defined as trip cost the variable indicates the cost brought by the perception time of mode while denotes the cost of the expected travel time relevant to a reference point and can be expressed as follows:

In the model, refers to the gain obtained by trip mode if , while is the loss obtained by trip mode if .

A set reference point makes distinguishing between gain and loss of a trip with respect to CPT possible for prospective customers. Customers typically perceive gain or loss by comparing trip cost with the reference point cost. In this context, it was assumed that a reference point may be related to the actual travel time experienced, typically about 30 minutes, and was thus used for this study [20].

##### 2.2. Cumulative Prospect Theory

CPT refers to the application of cumulative function to gains and losses. An uncertain prospect is a function that assigns to each state of a corresponding outcome. The prospect values are required to be arranged in increasing order to calculate the cumulative function. The prospect can be regarded as a sequence of pairs (), which yields if occurs, and is a partition of . The variable indicates the positive outcomes of and does the negative outcomes. For and , the relevant expressions are as follows:where is the projected gains and is the projected losses. Further, if is given by a probability distribution , can be regarded as an uncertain prospect . Here, can be expressed as

With the assumption of continuous traffic flow on a given roadway the distribution function can be introduced into this model as shown in

When combined with (8), (7) can be expressed as follows:

can be also defined by

From (8), (11) can also be expressed as

Thus cumulative prospect theory may be defined as follows:where and are the functions for gains and losses, respectively [24], and are defined by

From (4) and (5), it is known that is defined as the value function with the following formulation:where is the gain or loss; the variable is defined as a gain if and a loss if .

In these formulations, the parameters are , , , and [24].

Based on this, the cumulative prospect theory of trip modes can be obtained so that a customer may choose the trip mode based on maximal cumulative prospect value.

##### 2.3. The Relationship between Trip Mode Choices

Including cruising mode, dispatch mode, and other transit modes within the network the O-D pair (, , and ) is used to indicate the cumulative prospect values, respectively.

By adopting the Nested Logit model, the probability that a customer originating in zone chooses cruising mode to travel to zone is given as follows:

When the parameter is nonnegative, it may be used to reflect the degree of uncertainty in service from customers’ perspective.

Meanwhile, the probability of customers choosing dispatch mode is given as follows:

Therefore, the demand of customers choosing cruising service can be defined as follows:where is the customer demand for cruising mode and is the total customer demand for taxis in a network from zone to zone .

Similarly, the customer demand for dispatch service may be defined aswhere is the customer demand for dispatch mode. Based on this, the total customer demand for taxi service is expressed by the following formula:where is the total customer demand for taxi service from zone to zone .

#### 3. Characteristics of Taxi Movement in a Network

##### 3.1. Basic Assumptions of the Model

The basic assumptions of the model established in this paper are described as follows:(1)Assume that all occupied taxis follow the shortest path choice principle.(2)Once a taxi is occupied in zone , a driver will choose the shortest route to go to zone .(3)Given an O-D pair, (), for example, after a taxi has delivered customers to zone , the taxi driver either stays in the same zone or moves to other zones in search of the next customer in an attempt at minimizing search time.(4)Based on the context, all taxis in a network are divided into one of three types: occupied taxis from zone to zone , vacant taxis moving from zone to zone , or taxis searching for customers in zone .(5)The vacant taxi driver in zone can choose between cruising mode and dispatch mode to provide service for customers.

In other words, taxis on roadways can be divided into two types (occupied or vacant) and beginning to search for customers after arriving at zone ; those taxis can decide to make service mode choice.

##### 3.2. Introduction of Network Model Variables

Suppose that a road network with a taxi market is in a state of supply-demand equilibrium, defined as , where is the set of nodes and is the set of links. and are the set of origination and destination zones, respectively. is expressed as the travel time on link and is expressed as link flow including taxis and sources of traffic flow. Here, suppose that a taxi is only occupied by a customer.

Assume that is the travel time on route between O-D pair and is expressed as follows:

Further, if route between O-D pair uses link and 0 otherwise; is included in the set which is the set of routes between O-D pair .

Based on this context, is considered as the travel time via the shortest route from zone to zone :

##### 3.3. The Choice Relationship of Taxi Service Modes

In accordance with basic assumptions above, the probability that a vacant taxi arriving at zone chooses cruising mode to search for the next customer is given as follows: where is the probability of taxis choosing cruising mode, while is a nonnegative parameter, is the search time of taxi drivers that chose the cruising mode, and is the search time of taxi drivers that chose dispatch mode.

Similarly, the probability of vacant taxis choosing dispatch mode to search for the next customer is given as follows:where is the probability of taxis choosing dispatch mode.

##### 3.4. Taxi Service to Time Relationship

There is a given demand matrix OD and the 1 h interval is studied. This paper assumes that the total number of taxis in a network including occupied and vacant ones is and it is a constant.

The total number of occupied taxis is given as , with defined as the number of occupied taxis in the studied 1-h interval from zone to zone ; “ot” refers to occupied taxis.

Similarly, the total number of vacant taxis can be expressed as follows:where is the number of vacant taxis in the studied 1-h interval from zone to zone ; “vt” refers to vacant taxis. In the above relationship, is the number of vacant taxis on roadway; is the number of vacant taxis choosing cruising mode to search for customer in zone ; is the number of vacant taxis choosing dispatch mode in zone .

Therefore, the total number of taxis in a network can be expressed as

##### 3.5. Behavior Model of Taxi Drivers

The expected search time of drivers in one zone is assumed to be distributed in accordance with the Gumbel density function. Here is the total searching time of vacant taxis in zone and expressed as follows:The probability that a vacant taxi originating in zone eventually meets a customer in zone is given bywhere is the probability that a taxi meets the next customer in zone after having delivered a customer to zone , is a nonnegative parameter and its explanation is referred to by Wong et al. [13].

##### 3.6. The Relationship between Customer Wait Time and Taxi Search Time

According to Yang et al. [15], the function relationship between customer and taxi wait time can be described by applying the Cobb-Douglas model describing customer-taxi meeting function which is given bywhere is the average waiting time of customers choosing cruising mode in zone and is the average waiting time of customers choosing dispatch mode in zone .

##### 3.7. The Supply-Demand Equilibrium Relationship

When in an equilibrium state, the number of vacant taxis in a network should meet all customers’ demands in their respective origination zones. That is to say, taxi service is available to every customer. Therefore, the supply-demand equilibrium relationship in a network is given as follows:

In terms of cruising mode, its supply-demand equilibrium relationship is given by

In terms of dispatch mode, its supply-demand equilibrium relationship is expressed by

#### 4. A Mathematical Programming Model

##### 4.1. The Model Descriptions

A joint taxi network model of trip distribution and assignment can be given as follows:where is normal traffic flow in the studied 1-h interval from zone to zone and is the traffic flow on route in the studied 1-h interval.

In terms of the established supply-demand equilibrium model, the following Lagrangian function can be formed:

After differentiating the Lagrangian function and applying Kuhn-Tucker conditions, the following relationship can be obtained in an equilibrium state:

##### 4.2. The Deduction Process of the Model

Firstly, (36) can be rewritten as the following form of gravity model:where and .

By using and , it can be expressed as follows:

by first assuming and then alternatively solving (38) until converged.

Secondly, according to , (36) is also changed as follows:

Then, by combining (36) with (39), the following relationship is established:

By comparing (40) with (27), can now be interpreted as the taxi searching time in zone . After making a set of deductions, can be given as follows:

Here, the variable can be calculated by using

The theoretical process of (41) and (42) can be derived in detail [13]. By substituting the solution to in (42), the taxi searching time in zone can be then obtained.

Furthermore, by combining (28) with (41), the following form is also established as follows:

##### 4.3. The Existence of Equilibrium

In the taxi supply relationships, by combining (42) and (43), it can be seen that the taxi-searching time, , is a continuous function of customer demand vector . Furthermore, based on the assumption of the meeting function, (29) and (30), it can be discerned that customer wait time, , varies continuously with taxi search time, , in zone . Thus, it is concluded that the customer wait time is a continuous function of ; continuous mapping relationship can be expressed as . In the customer demand relationships, it can be seen that customer demand is a continuous function of the cumulative prospect vector . Additionally, based on the relationship between cumulative prospect and trip cost it may be deduced that cumulative prospect varies continuously with customer wait time . Such a continuous mapping relationship is also regarded as . Thus, the relationship is . In addition, the solution set of equilibrium model variables is a compact and convex set and the relationship exists. By applying Brouwer’s fixed-point theorem, it is easily concluded that has at least one fixed point in . Therefore the existence of the equilibrium solution is guaranteed [25].

#### 5. The Solution Algorithm

The iterative algorithm of taxi equilibrium problem of a network considering the trip mode choice model based on CPT is as follows.

*Step 1. *Initialize customer waiting time as follows.

Set an initial set of customer waiting times and , , then let .

*Step 2. *Updated customer demand value is as follows.

Substitute and into the cumulative prospect function and update , , and . Then and are updated in accordance with (16) and (17).

*Step 3. *Updated taxi search time is as follows. (1)Utilize the equilibrium model for taxis in a network.

By applying the iterative balancing algorithm, , and can be obtained. Substituting a set of variable values involving , , , , and , and so forth into (43) a formula involving and is established. The formula is defined by(2)Update the supply-demand equilibrium relationship.

Substitute and into (32) or substitute and into (33). Then update the appropriate variables. From this condition, another formulation is expressed as (45) and it includes and .(3)Combining (44) with (45) makes it possible to solve for the values of and .

*Step 4. *Update customer wait time as follows.

Substituting and into the relationship between taxi search time and customer wait time, (29) and (30), respectively, makes it permissible to update and .

*Step 5. *Check the following.

The standard to stop iteration is defined bywhere is a condition of iterative convergence.

If (46) is achieved, then stop. Otherwise, let and go to Step 2.

#### 6. A Numerical Example Analysis

A simple numerical example is given in order to check the convergence of the algorithm and present the resulting analysis of the equilibrium model.

##### 6.1. Introduction of a Numerical Example

For the model a six-zone network is utilized with nodes representing traffic zones, arcs representing relationship of adjacent nodes, and arrows indicating road direction (Figure 2). Further, each node represents a potential origin or a potential destination for customers, as well as a searching location for vacant taxis.

The travel time function for the links is defined as follows [26]:

Values for the parameters are obtained from Table 1.

A matrix of the total customer demands for taxis and other transit modes from zone to zone is supposed in Table 2.

The parameters used in the example were = 15 yuan/h, , , , , min. Additionally, the 30-minute time used in this example was regarded as a reference point for different trip modes in the CPT function. When making the sensitivity analysis, this paper set 1000, 1100, 1200, 1300, and 1400 as the values of , respectively.

##### 6.2. Result Analysis

Figure 3 shows the relationship of the taxi-searching times of dispatch mode versus the taxi fleet size in six zones. It can be seen that the taxi-searching times in zones are all positive and increase linearly with an increase in the taxi fleet size. This may be due to the fact that an increase in the taxi fleet size only enhances the availability of those taxis choosing dispatching service mode in taxi market but cannot decrease the searching time of dispatching service and achieve the improvement of dispatching service quality. That is to say, the changing trend of graphic in Figure 3 indicates that the demand for dispatch mode is less than its supply in taxi market. Specifically, taking a fixed taxi fleet size, for example, the taxi-searching times of dispatch mode in different zones are also different as shown in Figure 3. This is mainly because both the number of customer which chooses the dispatch mode and the one of taxi which provides the dispatching service are different in six zones. Figure 3 shows that, in terms of the taxi-searching time of dispatch mode, zone 4 is maximal and zone 6 is minimal. In addition, it can be seen from Figure 3 that the maximal gap of taxi-searching time in different zones reaches at least 0.7 h considering the taxi fleet size. It can be shown that the characteristic of a traffic zone should be discussed in terms of service and operation of taxi market.

Figure 4 presents the relationship of the taxi-searching times of cruising mode versus the taxi fleet size in six zones. These taxi-searching times in zones are all positive and increase nonlinearly with an increase in the taxi fleet size. It is attributed to the fact that an increase in the taxi fleet size leads to more greater supply of cruising mode and more worse service quality. Specifically, taking a fixed taxi fleet size, for example, the taxi-searching times of cruising mode in different zones are also different as depicted in Figure 4. It can be shown that, in terms of the taxi-searching time of cruising mode, zone 4 is maximal and zone 3 is minimal. The maximal gap of taxi-searching time in different zones reaches at least 0.6 h considering the taxi fleet size. In addition, combining Figure 3 with Figure 4, it can be seen that the taxi-searching times of cruising mode in zones are greater than the ones of dispatch mode, which reflects that the service quality of dispatch mode is better than the one of cruising mode in taxi market.

Figure 5 depicts the relationship between the probability of different modes chosen and an increase in the taxi fleet size in six zones. Figure 5 shows that, in terms of these two service modes, their values have no great changes with an increase in the taxi fleet size. And furthermore, compared with the probabilities of dispatch mode chosen in six zones, the ones of cruising mode chosen are relatively less, which is mainly attributed to the policy of dispatch mode and its parameters adopted by this paper.

**(a) The probability of cruising mode chosen**

**(b) The probability of dispatch mode chosen**

Figure 6 depicts the relationship of the total taxi-searching time versus the taxi fleet size in six zones. Figure 6 indicates that the total taxi-searching times in zones are all positive and increase linearly with an increase in the taxi fleet size. This is due to the fact that the increased taxis mainly search for customers instead of providing the occupied service. That is to say, an increase in the taxi fleet size can enhance taxi’s availability in taxi market but cannot indicate the improvement of service quality. The complicated relationship between taxi supply and customer demand is totally determined by the market and network conditions. This conclusion is consistent with the achievement from Yang et al. [15]. And furthermore, taking a fixed taxi fleet size, for example, total taxi-searching times in different zones are also different as shown in Figure 6. This is mainly because customer demands and network conditions (e.g., the travel time to traverse OD pair) are different in six zones. Figure 6 shows that, in terms of the total taxi-searching time, zone 4 is maximal, zone 6 is minimal, and zone 3 is nearly equal to zone 6. And moreover, Figure 6 indicates that the maximal gap of total taxi-searching time in different zones reaches at least thirty minutes considering the taxi fleet size. This means the characteristic of a traffic zone directly influences taxi’s availability and further the level of service quality in taxi market.

Figure 7 gives the relationship of the customer-waiting time in cruising mode versus the taxi fleet size in different zones. Figure 7 shows that the customer-waiting times in six zones are all positive and decrease nonlinearly with an increase in the taxi fleet size. This indicates that the marginal effect of advancing the taxi fleet size in improving service quality is diminishing. In terms of a given taxi fleet size, the customer-waiting times of cruising mode are different in six zones as shown in Figure 7. Figure 7 shows that, in terms of the customer-waiting time, zone 3 is maximal, zone 4 is minimal, and zone 4 is nearly equal to zone 5. And, moreover, Figure 7 indicates that the maximal gap of the customer-waiting time in different zones reaches at least three minutes. And, furthermore, this gap in six zones is decreasing with an increase in taxi fleet size, which means that the distribution and utilization of taxis in six zones will become similar with an increase in the taxi fleet size.

Figure 8 presents the relationship between the customer-waiting time of dispatch mode and an increase in the taxi fleet size in six zones. Figure 8 shows that the customer-waiting times in six zones are all positive and decrease nonlinearly with an increase in the taxi fleet size. In terms of a single taxi fleet size, the customer-waiting time of dispatch mode in zone 6 is maximal and the customer-waiting time of zone 4 is minimal. The maximal gap of this index in different zones reaches about six minutes. The tendency that the gap of customer-waiting time in dispatch mode in six zones is decreasing with an increase in the taxi number is then reflected.

#### 7. Conclusions

This paper considered the effect of the integrated service mode including cruising and dispatch modes and the travel time uncertainty on taxi network equilibrium. To facilitate this aim the cumulative prospect theory was adopted to study the trip mode choice problem. Based on the elasticity of customer demand, a mathematical programming model was established to describe the supply-demand equilibrium relationship of taxi service on a network with respect to the integrated service mode. From this an iterative algorithm was designed to present the optimal solutions of cruising mode’s parameters, dispatch mode’s parameters, and taxis network’s parameters in an equilibrium state. Finally, a numerical example was used to present the achievement of supply-demand equilibrium state and determine the outcomes of key parameters.

The results provided interesting insights and valuable recommendations for improving taxi market operation. Specifically, using taxi network equilibrium, the taxi-searching and customer-waiting time in six zones were obtained with an increase in the taxi fleet size. We found that the taxi-searching times of both dispatch mode and cruising mode in zones were positive and increased with an increase in the taxi fleet size, while the customer-waiting times of both dispatch mode and cruising mode in zones were also positive but decreased with an increase in the taxi fleet size. Moreover, the taxi-searching time of cruising mode was greater than the one of dispatch mode, while the customer-waiting time of cruising mode was also longer than the one of dispatch mode in six zones. Further, the algorithm indicated that the probability of dispatch mode being chosen by a customer was greater than cruising mode in the six zones considered. Finally, based on the above results from different service modes, the total taxi-searching times of integrated taxi mode were also obtained. Thus, from the practical point of view, these results should enable policymakers to know current operation situations in taxi market and to further adjust and formulate the reasonable operation policies in terms of integrated service mode.

However, to make this mathematical model more suitable for reality, several relationships in the model still need to be further improved. Firstly, the dispatch policy used in the model is a little simple so it cannot offer further insights into the interactions of cruising and dispatch modes. Therefore, applying a complicated dispatch policy will be more meaningful. Secondly, in the application of CPT, the reference points in those designed zones are only considered as constants instead of relevant functions, which need a suitable method to replace the current situation. And furthermore, an uncertainty of customers choosing trip mode should be reflected by adopting a continuous probability distribution function instead of the discrete function, which will be attempted to make further improvement.

#### Conflict of Interests

None of the authors has a financial relationship with a commercial entity that has an interest in the subject of this paper.

#### Acknowledgments

This research is supported by National Natural Science Foundation of China (Project no. 71390522). The authors would like to express her sincere thanks to the anonymous reviewers for their helpful comments and valuable suggestions on this paper.