Top Banner
I NTERDISCIPLINARY L IVELY A PPLICATIONS P ROJECT Travel Demand Forecasting
24

I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

May 03, 2018

Download

Documents

vocong
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

I N T E R D I S C I P L I N A R Y L I V E L Y A P P L I C A T I O N S P R O J E C T

Travel Demand Forecasting

Page 2: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

Project Intermath 1-23. ©Copyright 2002 by COMAP, Inc. All rights reserved. Permission to make digital or hard copies of part or all of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice. Abstracting with credit is permitted, but copyrights for components of this work owned by others than COMAP must be honored. To copy otherwise, to republish, to post on servers, or to redistribute to lists requires prior permission from COMAP.

Interdisciplinary Lively Application Project Title: Travel Demand Forecasting and Analysis for South Texas Authors: Lei Yu and Della D. Bell Department of Transportation Studies and Department of Mathematics, Texas Southern University, 3100 Cleburne Avenue, Houston, Texas 77004 Editor: David C. Arney Mathematics Classification: Algebra, Calculus, Statistics Disciplinary Classifications: Transportation Planning Prerequisite Skills:

1. Solving Linear Equations 2. Solving Quadratic Equations 3. Regression analysis concept 4. Elementary Optimization

Physical Concepts Examined:

1. Trip Attractions 2. Trip Distributions 3. Trip Modal Choice 4. Trip Route Choice

Materials:

1. Problem Statement; Student 2. Sample Solution; Instructor

Page 3: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

2 Project Intermath

3. Notes for the Instructor Computing Requirements:

1. Microsoft Excel or Similar Tools; 2. Ability to Use Calculator

Contents

1. Problem Statement 2. Sample Solution 3. Instructor Notes

1. Problem Statement Introduction Travel demand forecasting is the most important phase in the urban transportation planning process. The purpose of travel demand forecasting is to predict the travel demands on the roads in order to estimate the likely transportation consequences of transportation alternatives (including a do-nothing alternative) that are being considered for implementation. Usually, travel demand forecasting is performed using a 4-step sequential model described as follows:

Step 1: Trip Generation: Should I make a trip? Step 2: Trip Distribution: Where should I go? Step 3: Modal Choice: What mode of transportation should I use? Step 4: Trip Assignment: Which route in the network should I take?

Step 1 determines how many trips will be produced by each residential area and how many trips will be attracted by each commercial site. Step 2 determines where each trip generated in Step 1 will go. For example, for each residential area, there may exist multiple choices for the shopping purpose. The percentage of the total trips that will be attracted by each shopping site must be determined. Step 3 determines what mode of transportation each trip will use. The choices of transportation may include bus, light-rail (used only in the city), subway, private automobile, taxi, bicycle, walk and so on. Finally, Step 4 determines which route each trip will use.

Page 4: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

Travel Demand Forecasting 3

By the end of the travel demand forecasting process, the traffic volumes (the number of vehicles per unit time) on the roads will be produced, which provide useful information about the congestion on the streets. Thereafter, transportation planners can select the best transportation projects by reviewing the resulting levels of congestion from a series of transportation alternatives. The 4-step travel demand forecasting process is usually carried out using various mathematical models. PART 1: Step 1 - Trip Generation

General Information The objective of a trip-generation model is to forecast the number of trips that will begin from or end in each travel zone within the region for a typical day of the target year. The most widely used trip generation models are regression models which are expressed by the following two linear equations:

rri XaXaXaaP ++++= ...22110 (1) rri XaXaXaaA ++++= ...22110 (2)

Pi represents the total number of trips produced by a residential zone i and Ai

represents the total number of trips attracted by commercial activities in zone i. X1 to Xr are a series of independent variables, which are usually derived from the urban travel surveys. For example, X1 may represent the total zonal population, X2 may represent the average household income, and X3 may represent the average auto ownership of each zone. a0 to ar are constant values/coefficients of the independent variables, which are usually derived through the regression analysis.

It should be noted that the above regression equations for trip productions and attractions can be non-linear in many cases. However, since the regression analysis for non-linear equations is often performed by converting the non-linear form to a linear form first, the understanding of solving the linear equations is essential.

Example: Trip Attractions in South Texas A survey of the commercial activities was conducted for five zones in South Texas. Each zone represents a portion of the land that is used as the basic unit in an urban planning process. The data were collected based on three types of employment: manufacturing, retail and services, and others. The resulted zonal employment of three different commercial types and their respective trip attractions are listed in

Page 5: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

4 Project Intermath

the following table. The first column in Table 1 represents the Zone number. The second column is the total number of people employed in manufacturing related companies for each zone, while the third column is the total number of people employed in the retail and services sector. The fourth column is the total number of people employed in all companies other than manufacturing and retail and services, while the fifth column represents the total number of people employed for all companies. The sixth and also the last column lists the surveyed number of trips that are attracted by each zone. Variables X1, X2, X3, and X are used to express the employment types, while the variable Y is used to express the trip attractions.

Zonal Employment

Trip

Attraction Zone Manuf. Ret&Ser Others Total

X1 X2 X3 X Y 1 6820 2547 115 9482 9428 2 111 1899 0 2010 2192 3 228 87 259 574 330 4 0 127 0 127 153 5 2729 813 294 3836 3948

Table 1: The number of people employed in each sector and the total trip attractions This example intends to let students excise the concept of developing trip attraction equations through the regression analysis technique. In practice, the developed attraction equations are used to predict the actual trip attractions given the projected employment numbers for a future year. Trip productions should be predicted in the same way. To simplify the problem, this example is only designed to practice the regression analysis process, but not to predict the future trip attractions. In the Part 2: Step 2 of the next section, trip productions and attractions are assumed to have been predicted beforehand.

Requirement 1: Determine a single linear regression equation between dependent variable Y and each of independent variables X, X1, X2, and X3 individually. Requirement 2 (Optional): Determine a multiple linear regression equation between dependent variable Y and independent variables X1, X2, and X3. (Use of Microsoft Excel Data Analysis Tool is suggested.)

Page 6: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

Travel Demand Forecasting 5

Requirement 3: Select the equations from Requirements 1 and 2 that are acceptable for use in trip generation analysis. (At a minimum, a minimum value of R-squared, such as 0.65, should be set as the criteria and compared with the value of R-squared for each equation. Optimally, an F-test and/or t-test should be conducted.)

PART 2: Step 2 - Trip Distribution

General Information The objective of a trip distribution model is to determine the total number of trips between all pairs of zones i and j, where i is the trip-producing zone and j is the trip-attracting zone of the pair. The rationale of trip distribution is as follows: all trip-attracting zones j in the region are in competition with each other to attract trips produced by each zone i. Everything else being equal, more trips will be attracted by zones that have higher levels of “attractiveness.”

A widely used trip distribution model is the gravity model, which gets its name from the fact that it is conceptually based on Newton’s law of gravitation. Newton’s law of gravitation states that the force of attraction between two bodies is directly proportional to the product of the masses of the two bodies and inversely proportional to the square of the distance between them. It can be described by the following equation:

221

rMMkF = (3)

The variation of the Newton’s law of gravitation to the trip distribution takes the following form:

cij

jiij W

APkT = (4)

Tij is the total number of trips between zones i and j, Pi and Aj are the total trip productions for zone i and attractions for zone j, which can be derived from the trip generation step. Aj may also be the relative attractiveness of Zone j. For a given starting zone, k is a constant value. Wijc is the travel impedance, which can be defined as either the travel time or the travel distance. Rewriting Equation (4) in a form that is expressed by the trip productions and a probability factor results in the following Equation (5).

Page 7: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

6 Project Intermath

cij

ij

xixx

ijjiij W

FFA

FAPT 1 where =

=

∑ (5)

Fij is called the travel-time factor (or friction factor.) The term contained by the parentheses is a probability factor, which represents the proportion that trips produced by Zone i are attracted by Zone j. This term is usually expressed by pij. Usually, the value of exponent c in Equation (5) must be determined by calibration before the gravity model is applied in the trip distribution process.

Example: Trip Distribution in South Texas In a small town in South Texas, the land uses are divided into four zones as shown in Figure 1. Zone 1 is a residential zone, which does not have any commercial sites. Zone 2 and Zone 4 are commercial zones, which do not include any residential land uses. Zone 3 is a combined residential and commercial zone, which has both residential and commercial land uses. Zones 2, 3 and 4 compete in attracting trips that are produced by the households in Zones 1 and 3.

Based on a trip-generation analysis, the total trip productions and the relative zonal attractiveness for the target-year 2005 are generated and illustrated by Table 2. The inter-zonal and intra-zonal travel times are found to be the Table 3. A calibration of the gravity model has also found that c in the gravity Equation (5) equals 2.0.

Figure 1: A 4-zone network

Zone 1 Zone 2

Zone 3 Zone 4

Residential Zone Commercial Zone

Residential/Commercial Zone Commercial Zone

Page 8: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

Travel Demand Forecasting 7

Table 2: Trip productions Pi and attractiveness Aj for the 4-zone network

i\j 1 2 3 4 1 5 10 10 20 2 10 5 10 10 3 10 10 5 10 4 20 10 10 5

Table 3: Travel time Wij between Zone i and Zone j

Requirement 4: Use the gravity model Equation (5) to estimate the interchange trips between zones.

Requirement 5: Estimate the total trip attractions for Zones 2, 3 and 4.

PART 3: Step 3 - Modal Choice

General Information In a typical travel situation, trip-makers can select between several travel modes. These include driving, riding with someone else, taking the bus, walking, riding a motorcycle and so forth. A modal choice (sometimes called modal split) is concerned with the trip-maker’s behavior regarding the selection of travel mode.

The most widely used modal choice model is the multinomial Logit Model. The multinomial logit model calculates the proportion of travelers that will select a specific mode K according to the following relationship:

=

x

U

U

x

K

eeKp )( (6)

UK is a utility function, which measures the degree of satisfaction that people derive from their choices. The magnitude of UK depends on the characteristics of each choice, which are called attributes, and on the characteristics of the individual making that choice, which are called socioeconomic status. The utility function is typically expressed as the linear weighted sum of the independent variables or their transformation as the follow:

Zone Productions Attractiveness 1 20000 0 2 0 3 3 15000 2 4 0 6

Page 9: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

8 Project Intermath

rRKK XaXaXaaU ++++= ...2211 (7)

aK is the calibrated mode-specific constant for the mode K and Xi is the ith attribute weighted by the model parameter ai. The actual attributes may include such variables as the access and egress time, waiting time, line-haul time (actual time spent in the vehicle,) out-of-pocket cost, level of service and convenience.

Example: Modal Choice in South Texas In the same town as used in PART 2, people who go from Zone 1 to Zone 4 have two choices of transportation: private automobile and local bus, which is operated by a local public transportation authority. A calibration process has resulted in the following utility function for the auto and bus:

4321 001.0025.002.005.0 XXXXaU KK −−−−=

where X1 = waiting time, in minutes X2 = line-haul time, in minutes X3 = access time, in minutes X4 = out-of-pocket cost, in cents

The interchange trips between Zone 1 and Zone 4 for the target-year 2005 has been forecasted in the PART 2. The target-year service attributes of the two competing modes have been estimated as Table 4:

Attribute X1 X2 X3 X4 Automobile 0 15 5 80 Local Bus 15 35 10 50 Table 4: Target-year service attributes for automobile and local bus

Requirement 6: Assuming that the calibrated mode-specific constants aK are –0.14 for the automobile and –0.60 for the bus, use the logit model Equation (6) to estimate the target-year market shares of the automobile and bus and the resulting fare revenue of the bus system.

Requirement 7: A proposal is submitted to build a subway transit system (ST) between Zone 1 and Zone 4 for this town. A study has projected that the service attributes of the proposed subway system will be:

X1 (ST) = 5 X2 (ST) = 30 X3 (ST) = 10 X4 (ST) = 100

Page 10: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

Travel Demand Forecasting 9

Based on professional experience, the mode-specific constant for the subway system is –0.45. Estimate the new market shares of three modes that will result from implementing the subway system proposal and the effect on the revenues of the public transportation authority, which operates the bus and subway systems. Requirement 8: What other factors need to be considered to effectively decide if a subway should be built? Explain.

PART 4: Step 4 - Trip Assignment

General Information The last step in the sequential travel demand forecasting process is concerned with the trip-makers’ choice of path between pairs of zones by travel mode and with the resulting vehicular flows on the multimodal transportation network. The trip assignment problem can be described as follows: given traffic demand departure rates for each Origin Destination (OD) pair in the network, which are usually derived from the Steps 1-3 of the travel demand forecasting process, determine the likely path choices and path use probabilities for each OD pair in the network, and predict the resulting traffic flows and link travel times on each of the individual links that make up the network.

Most of the methods for solving the trip assignment problem are based on two trip assignment principles first presented by Wardrop. These principles are:

Wardrop’s 1st principle: each driver may strive to choose the path that will minimize his own travel time through the network. Wardrop’s 2nd principle: each driver may choose the path that will minimize the total network travel time of all drivers.

When the network trips are such assigned that no driver can reduce his own travel time by switching from his current route to an alternative, the traffic flows are said to be in a user-equilibrium state. This state coincides with Wardrop’s 1st principle. Similarly, when the network trips are such assigned that no driver can reduce the total network travel time by switching from his current route to an alternative, the traffic flows are said to be in a system equilibrium state. This state coincides with Wardrop’s 2nd principle.

The trip assignment problems can be formulated as non-linear mathematical programs based on the above two Wardrop’s assignment principles. The

Page 11: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

10 Project Intermath

following will use a simple two-node and two-link network to describe the trip assignment models based on the Wardrop’s principles.

It is assumed that there are only two nodes and two links in the network as shown below. Node 1 is an origin where the trips generate and Node 2 is a destination where the trips sink. There are two routes available for driving from Node 1 to Node 2 through either Link 1 or Link 2. The link travel time on either Link 1 or Link 2 is a function of the traffic flow on the link. In other words, link travel time will vary with the change of the number of vehicles on the link. If it is assumed that r12 vehicles per hour go from Node 1 to Node 2, the objective of the trip assignment process is to determine how much traffic of the r12 will use Link 1 and how much will use Link 2.

Figure 2: Illustration of a simple 2–node 2-link network

The link travel time functions are expressed as follows:

)( 11 Xftt = (8) )( 22 Xftt = (9)

where X1 and X2 are the total number of vehicles per hour on Link 1 and Link 2, and tt1 and tt2 are the travel times on Link 1 and Link 2. To perform the trip assignment for this simple network, two formulations can be made based on the Wardrop’s two principles.

Based on Wardrop’s 1st Assignment Principle - Based on the Wardrop’s 1st assignment principle, the travel times on Link 1 and Link 2 must be equal after all trips are assigned. Otherwise, some vehicles on the higher travel time link will automatically switch to the lower travel time link, as each individual driver seeks

1 2

Link 1

Link 2

Page 12: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

Travel Demand Forecasting 11

to minimize his/her own travel times. Based on this logic and the assumption that there are a total of r12 vehicles per hour going from Node 1 to Node 2, the following equations hold:

)()( 21 XfXf = (10) 1221 rXX =+ (11)

Equations (10) and (11) contain only two variables. Therefore, X1 and X2 can often be solved for using these two simultaneous equations if f(X1) and f(X2) are linear equations. Accordingly, the link travel times can also be calculated using Equation (8) or (9). Based on Wardrop’s 2nd Assignment Principle - Based on the Wardrop’s 2nd assignment principle, the total travel times of Link 1 and Link 2 should be minimum. Otherwise, some vehicles will switch the link in order to reduce the total travel times. The total travel times for Link 1 and Link 2 are expressed as follows:

)()( 2211 XfXXfXTT += (12)

Therefore, the following optimization formulation can be derived:

1221

2211

to subject)()( min

rXXXfXXfX

=++

(13)

The solution of Equation (13) will result in the values of X1 and X2. Thereafter, the individual Link travel times as well as the total link travel times can be calculated based on Equations (8), (9) and (12).

Example: Route Choice in South Texas In the same example network as the assignment in PART 2, it is assumed that the automobile users driving from Zone 1 to Zone 4 have two available routes to use, similar to what has been illustrated in Figure 2. During the morning peak hour, the travel times using the two routes are functions of the vehicles using each route, as indicated by Equations (8) and (9). It is assumed that travel time functions for Route 1 and Route 2 are expressed by the following equations:

minutes 02.010 11 Xtt += minutes 01.015 22 Xtt +=

Page 13: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

12 Project Intermath

where tt1 and tt2 are the actual travel times on Route 1 and Route 2 and X1 and X2 are the number of vehicles using Route 1 and Route 2 per hour. It is estimated from the trip distribution and modal split process that during the morning peak hour, about 2000 vehicles per hour will leave Zone 1 heading for Zone 4.

Requirement 9: Based on the Wardrop’s first user equilibrium principle and Equations (10) and (11), formulate the trip assignment problem. Solve the formulation and estimate the number of vehicles that will use Route 1 and the number of vehicles that will use Route 2. Calculate the resulting travel times on each route and the total vehicle travel times between Zone 1 and Zone 4.

Requirement 10: Based on the Wardrop’s second system optimal trip assignment principle and Equation (13), formulate the trip assignment problem. Solve the formulation and estimate the number of vehicles that will use Route 1 and the number of vehicles that will use Route 2. Calculate the resulting travel times on each route and the total vehicle travel times between Zone 1 and Zone 4.

Requirement 11: Compare the results from Question 1 and Question 2. Check if these results are consistent with the Wardrop’s trip assignment principles. Comment on your findings.

Page 14: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

Travel Demand Forecasting 13

2. Sample Solution PART 1: Step 1 – Trip Generation Trip Attractions in South Texas Requirement 1: Using the Excel Tool – Data Analysis – Regression, the linear regressions for Y-X, Y-X1, Y-X2, and Y-X3 result in the following output Table 5 to Table 8. Requirement 2: Using the Excel Tool – Data Analysis – Regression, the linear regression for Y-X1, X2 and X3 results in the following output Table 9 and Table 10. Requirement 3: Assuming that the value of R-square higher than 0.65 is acceptable for the analysis, the following equations are selected for the trip attraction analysis:

998124.0 001321.1163903.0 2 =+= RXY 954894.0 266307.19511.705 2

1 =+= RXY 69364.0 896011.222623.40 2

2 =+= RXY 998761.0 858613.0105403.1969715.0321931 2

321 =+++−= RXXXY If F-test and t-test are conducted to further examine the validity of each equation, the following equations are found acceptable statistically:

998124.0 001321.1163903.0 2 =+= RXY 954894.0 266307.19511.705 2

1 =+= RXY 99817.0 013371.1001821.17625.119 2

21 =++= RXXY

Page 15: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

14 Project Intermath

SUMMARY OUTPUT

Regression Statistics

Multiple R 0.999061 R Square 0.998124 Adjusted R Square 0.997498 Standard Error 190.2617 Observations 5

ANOVA

df SS MS F Significance F Regression 1 57775542 57775542 1596.03 3.45E-05 Residual 3 108598.6 36199.53 Total 4 57884141

Coefficients Standard

Error t Stat P-value Lower 95% Upper 95% Lower

95.0% Upper 95.0%

Intercept 0.163903 117.0305 0.001401 0.99897 -372.28 372.6074 -372.28 372.6074Total 1.001321 0.025064 39.95035 3.45E-05 0.921556 1.081087 0.921556 1.081087

Table 5: Simple linear regression result for trip attraction and total employment

SUMMARY OUTPUT

Regression Statistics Multiple R 0.977187 R Square 0.954894 Adjusted R Square 0.939859 Standard Error 932.9017 Observations 5

ANOVA

df SS MS F Significance F Regression 1 55273224 55273224 63.51013 0.004122 Residual 3 2610917 870305.6 Total 4 57884141

Coefficients Standard

Error t Stat P-value Lower 95% Upper 95% Lower

95.0% Upper 95.0%

Intercept 705.9511 522.3078 1.3516 0.269385 -956.267 2368.169 -956.267 2368.169Manuf. 1.266307 0.158898 7.969325 0.004122 0.760623 1.771991 0.760623 1.771991

Table 6: Simple linear regression result for trip attraction and manufacturing employment

Page 16: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

Travel Demand Forecasting 15

SUMMARY OUTPUT

Regression Statistics

Multiple R 0.832851 R Square 0.69364 Adjusted R Square 0.59152 Standard Error 2431.28 Observations 5

ANOVA

df SS MS F Significance F Regression 1 40150774 40150774 6.792412 0.079945 Residual 3 17733367 5911122 Total 4 57884141

Coefficients Standard

Error t Stat P-value Lower 95% Upper 95% Lower

95.0% Upper 95.0%

Intercept 40.22623 1631.45 0.024657 0.981877 -5151.78 5232.234 -5151.78 5232.234Ret&Ser 2.896011 1.11119 2.606226 0.079945 -0.64029 6.432316 -0.64029 6.432316

Table 7: Simple linear regression result for trip attraction and retail and service employment

Table 8: Simple linear regression result for trip attraction and others employment

SUMMARY OUTPUT

Regression Statistics Multiple R 0.087814 R Square 0.007711 Adjusted R Square -0.32305 Standard Error 4375.606 Observations 5

ANOVA

df SS MS F Significance F Regression 1 446360.4 446360.4 0.023314 0.888336 Residual 3 57437780 19145927 Total 4 57884141

Coefficients Standard

Error t Stat P-value Lower 95% Upper 95% Lower

95.0% Upper 95.0%

Intercept 2889.568 2870.336 1.006701 0.388241 -6245.13 12024.27 -6245.13 12024.27Others 2.399937 15.71792 0.152688 0.888336 -47.6216 52.42143 -47.6216 52.42143

Page 17: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

16 Project Intermath

Table 9: Multiple linear regression result for trip attraction and manufacturing, retail and service, and

others employment SUMMARY OUTPUT

Regression Statistics

Multiple R 0.999084 R Square 0.99817 Adjusted R Square 0.996339 Standard Error 230.1669 Observations 5

ANOVA

df SS MS F Significance F Regression 2 57778187 28889094 545.3159 0.00183 Residual 2 105953.6 52976.8 Total 4 57884141

Coefficients Standard

Error t Stat P-value Lower 95% Upper 95% Lower

95.0% Upper 95.0%

Intercept 119.7625 154.5094 0.775115 0.519369 -545.0381 784.563 -545.038 784.563Manuf. 1.001821 0.054921 18.24112 0.002992 0.765515 1.238128 0.765515 1.238128Ret&Ser 1.013371 0.14737 6.876347 0.020501 0.379286 1.647455 0.379286 1.647455

Table 10: Multiple linear regression result for trip attraction and manufacturing and retail and service

SUMMARY OUTPUT

Regression Statistics Multiple R 0.99938 R Square 0.998761 Adjusted R Square 0.995043 Standard Error 267.8365 Observations 5

ANOVA

df SS MS F Significance F Regression 3 57812404 19270801 268.6336 0.044814 Residual 1 71736.39 71736.39 Total 4 57884141

Coefficients Standard

Error t Stat P-value Lower 95% Upper 95% Lower

95.0% Upper 95.0%

Intercept -32.1931 284.1409 -0.1133 0.928177 -3642.53 3578.144 -3642.53 3578.144Manuf. 0.969715 0.079029 12.27035 0.051768 -0.03444 1.97387 -0.03444 1.97387Ret&Ser 1.105403 0.217178 5.089861 0.123503 -1.65409 3.864894 -1.65409 3.864894Others 0.858613 1.243212 0.690641 0.615216 -14.9378 16.65505 -14.9378 16.65505

Page 18: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

Travel Demand Forecasting 17

PART 2: Step 2 – Trip Distribution Trip Distributions in South Texas Requirement 4: For the origin i = 1, the trip production P1 = 20000. Based on the Equation (5) and Table 2, the friction factor F1j can be calculated as follows:

0025.020

11 01.010

11

01.010

11 04.0511

2214

142213

13

2212

122211

11

======

======

WF

WF

WF

WF

The interchange trips between Zone 1 and other zones can then be calculated based on the Equation (5). To simplify the calculation process, the following table is used to perform the calculation.

j Aj F1j AjF1j P1j T1j 1 0 0.04 0 0 0 2 3 0.01 0.03 0.462 9240 3 2 0.01 0.02 0.308 6160 4 6 0.0025 0.015 0.230 4600 ∑ = 0.065 1.000 For the origin i = 3, the trip production P3 = 15000. A similar calculation table to the above can be used to calculate the interchange trips between Zone 3 and other zones.

j Aj F3j AjF3j P3j T3j 1 0 0.01 0 0 0 2 3 0.01 0.03 0.176 2640 3 2 0.04 0.08 0.471 7065 4 6 0.01 0.06 0.353 5295 ∑ = 0.17 1.000 For Origins 2 and 4, since there are no trip productions, the interchange trips between these two zones and other zones become zero. Therefore, the final interchange trips can be summarized into the following table:

i\j 1 2 3 4 1 0 9240 6160 4600 2 0 0 0 0 3 0 2640 7065 5295 4 0 0 0 0

Page 19: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

18 Project Intermath

Requirement 5: If a row in the above interchange trip table is summed, it should be equal to the total trip productions for that zone. On the other hand, if a column in the table is summed, the total trip attractions for that can be derived. Therefore, the total trip attractions for Zones 2, 3 and 4 are calculated as follows:

98955295460013225706561601188026409240

*4

*3

*2

=+=

=+=

=+=

AAA

PART 3: Step 3 – Modal Choice Modal Choice in South Texas Requirement 6: Based on the trip distribution analysis in PART 2, the interchange trips between Zone 1 and Zone 4 are 4600. Given the utility function and the attribute values by Table 3, the utilities for automobile and bus can be calculated as

35.2)50(001.0)10(025.0)35(02.0)15(05.060.0)(095.1)80(001.0)5(025.0)15(02.0)0(05.014.0)(

−=−−−−−=−=−−−−−=

BusUAutoU

Therefore, the proportions of trips that will select the automobile and bus can be estimated as follows using the Equation (6):

%22)(

%78)(

35.2095.1

35.2

35.2095.1

095.1

=+

=

=+

=

−−

−−

eeeBusp

eeeAutop

The market shares of the automobile and bus are calculated as

trips/day 1012)4600)(22.0()(trips/day 3588)4600)(78.0()(

14

14

====

BusTAutoT

The fare revenue of the bus system = (1012 trips/day)($0.50/trip) = $506/day Requirement 7: By introducing the subway system, the utilities for automobile and bus will keep the same as those calculated in the above Question 1, which are U(Auto) = -

Page 20: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

Travel Demand Forecasting 19

1.095 and U(Bus) = -2.35. The utility of the proposed subway system is calculated as follows:

65.1)100(001.0)10(025.0)30(02.0)5(05.045.0)( −=−−−−−=SubwayU Therefore, the proportions of trips that will use each mode are estimated as:

%31)(

%15)(

%54)(

65.135.2095.1

65.1

65.135.2095.1

35.2

65.135.2095.1

095.1

=++

=

=++

=

=++

=

−−−

−−−

−−−

eeeeSubwayp

eeeeBusp

eeeeAutop

The new market shares after the introduction of the subway system are calculated as follows:

trips/day 1426)4600)(31.0()(trips/day 690)4600)(15.0()(

trips/day 2484)4600)(54.0()(

14

14

14

======

SubwayTBusTAutoT

The new fare revenue of both bus and subway systems = (690 bus trips/day)($0.50/bus trip) + (1426 subway trips/day)($1.00/subway trip) = $1,771/day. The revenue increase of the public transportation authority, which operates both the local bus and subway systems, due to the introduction of the subway system is $1,265 ($1,771 - $506). It should be noted that the increase of the fare revenue only does not justify the construction of the subway. In reality, a comprehensive cost-benefit analysis should be conducted before the decision is made. Requirement 8: Free thoughts.

PART 4: Step 4 – Trip Assignment Route Choice in South Texas Question 1: Based on the Wardrop’s first assignment principle, Equations (10) and (11), and the travel time functions for Route 1 and Route 2, the following linear simultaneous equations can be established:

Page 21: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

20 Project Intermath

200001.01502.010

21

21

=++=+

XXXX

The above equations can be solved as X1 = 833.33 and X2 = 1166.67, which means that 833.33 vehicles will actually use Route 1 and 1166.67 vehicles will use Route 2. The travel times are

minutes-vehicle 340,53)67.26)(67.1166()67.26)(33.833( times travel Totalminutes 67.26)67.1166(01.015 2 Route on time Travel

minutes67.26)33.833(02.010 1 Route on time Travel

=+==+=

=+=

Question 2: Based on the Wardrop’s second trip assignment principle, Equation (13), and the travel time functions, the following optimization formulation can be established:

2000 to subject)01.015()02.010( min

21

2211

=++++=

XXXXXXTT

where TT represents the total system travel times for all vehicles going from Zone 1 to Zone 4. The constraint can be re-written as X2 = 2000 – X1. Substituting this into the objective function results in

( )700005503.0

)2000(01.015)2000()02.010(

121

1111

+−=

−+−++=

XXXXXXTT

By setting the derivative of the above equation to zero, the optimal value for X1 can be solved.

12507502000 Therefore750

045)2)(03.0(

2

1

11

=−==

=−=

XX

XdXdTT

Therefore, based on the Wardrop’s second assignment principle, 750 vehicles will use Route 1 and 1250 vehicles will use Route 2. The travel times are calculated as follows:

Page 22: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

Travel Demand Forecasting 21

minutes-vehicle 125,53)5.27)(1250()25)(750( times travel Totalminutes 5.27)1250(01.015 2 Route on time Travel

minutes25)750(02.010 1 Route on time Travel

=+==+==+=

Question 3: Based on Wardrop’s first assignment principle, all the used routes should have equal travel times. In Question 1, both Routes 1 and 2 have travel time of 26.67 minutes, which indicates that this result conforms to the Wardrop’s first assignment principle. In Question 2, two routes have different travel times, which means that the result does not conform to the Wardrop’s first assignment principle. Based on the Wardrop’s second assignment principle, the total travel times should be minimum. Question 2 resulted in the total travel times of 53,125 vehicle-minutes, which are less than the total travel times of 53,340 in Question 1. This means that although Wardrop’s second assignment principle does not result in equal travel times for all routes, it does result in minimum total travel times.

Page 23: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

22 Project Intermath

3. Notes for the Instructor This project is designed to exercise some basic concepts in regression analysis, algebra, calculus, and optimization theory using a realistic example of transportation planning process. Understanding the travel demand forecasting process is critical in this project, while the mathematical calculation part for each assignment is not very difficult. The project is suitable for either individual or group work. With the understanding of the travel demand forecasting process concept by the students, the instructor can expand the assignment in each part to let students practice. The assignment in each part can also be used independently for exercise of mathematical topics. The requirements for PART 1 through PART 4 can also be used for a complete project for students. Some additional suggestions are as follows: 1. In PART 1, similar problem can be also designed for trip productions in which the

zonal population, household income and the household automobile ownership become the independent variables while the trip production becomes the dependent variable. If a single independent variable is used, students can be asked to do the simple linear regression analysis manually. Non-linear regression equations can be considered to let students practice how to transform a non-linear equation to a linear equation.

2. In PART 2, variations of the requirement can be designed. More zones can be added to the analysis areas. Residential zones and commercial zones can be redesigned. The attractiveness can be assigned different values to test how the changes will affect the trip distribution results. The travel time table can use different values for each cell. The exponential power c in the gravity equation can be assigned a different value other than 2 to examine how different values of c will change people’s choices of commercial zones.

3. In PART 3, students can be asked to change the attribute values for automobile and local bus, for example the out-of-pocket cost or waiting time, to examine how the changes will affect the choice of mode. The coefficients of attributes in the utility function can also be redesigned so that a separate problem can be created. In addition to the introduction of the subway transit system, more modes of travel can be considered, such as motorcycle, walking, ride-sharing and so on. If more modes are added to the problem, values of various attributes for each added mode should be assumed.

4. In PART 4, three-route and four-route problems can be designed, which will make the mathematical equations more complicated. In designing additional route(s) to

Page 24: I NTERDISCIPLINARY L IVELY A PPLICATIONS P … example, for each residential ... be attracted by zones that have higher levels of “attractiveness.” ... Based on a trip-generation

Travel Demand Forecasting 23

the network, new travel time function(s) should be assumed. Essentially, different forms of the travel time functions reflect different conditions of the roads. For example, some roads have one lane and some roads have two lanes in each direction. Some roads may have a lot of parked vehicles on the roadside, which significantly impede the moving of the vehicles, while some roads may have uphill and downhill slopes or short radius on curves. The travel time functions can also be made non-linear, which will make the optimization more complicated.