Urban Transportation Planning Trip Distribution Modelling Murtaza Haider

Trip Distribution Models • Problem Definition, Terminology • Growth Factor Model

• The Proportional Flow Model

• The Singly-constrained Gravity Model

• Bi-Proportional Updating (Fratar Method)

Urban Transportation Modelling System Population & Employment Forecasts

Trip Generation

Trip Distribution

Mode Split

Trip Assignment

Link & O-D Flows, Times, Costs, etc.

Transportation Network & Service Attributes

Forecasting Future Distributions • Trip Origins

• Trip Destinations 2004

D1

D2

D3

D4

D5

2020

D1

D2

D3

D4

D5

Oi

2004 O1

2020 O1

O2

O2

O3

O3

O4 O5

O4 O5

1

1 T11

2 T12

Dj 3 T13

2

T21

T22

T23

T24

T25

O2

3

T31

T32

T33

T34

T35

O3

4

T41 T51

T42 T52

T43 T53

T44 T54

T45 T55

O4 O5

D1

D2

D3

D4

D5

5

4 T14

5 T15

O1

Definitions Dest j

Org i

Tij

Oi

Dj

T

Definitions N = Number of zones in the urban area i j

= Subscript, used to denote origin zones

= Subscript, used to denote destination zones

Oi = Number of trips originating in zone i Dj = Number of trips destined for zone j

Tij = Number of trips (“flow”) from origin zone i to destination zone j

Tijm = Number of trips from i to j using mode m

Definitions, Cont’d

Oi

1

1 T11

2 T12

Dj 3 T13

2

T21

T22

T23

T24

T25

O2

3

T31

T32

T33

T34

T35

O3

4

T41 T51

T42 T52

T43 T53

T44 T54

T45 T55

O4 O5

D1

D2

D3

D4

D5

5

4 T14

5 T15

O1

O1= T11+ T12+ T13+ T14+ T15 (row vector) D1= [T11+ T21+ T31+ T41+ T51]’(column vector)

=

Objective in Trip Distribution Modeling • Need to to estimate the “most likely” origindestination (O-D) trip (or flow) matrix, given known origins and destinations for each zone (plus any additional relevant, known information)

O-D Table Example Zone 1 2 3 4 5 Base Dj Djf

Bal.Factor Sum

1 2 24.6 1 24.1 6.6 18.2 4.2 28.5 1.3 9.1 0.7 104.5 13.8 114

16

3 1 3.8 13.2 2.5 1.2 21.7

4 2.7 3.3 4.3 10.7 3.4 24.4

22

30

5 Base Oi 0.9 30.2 1.2 39 1.1 41 3.4 46.4 4.9 19.3 11.5 175.9

Oif

16

Balanced Trip matrix Matrix Bal.Factor Sum 27.9 1.2 1.0 3.5 1.3 0.00266 35.00 26.6 7.7 3.9 4.2 1.7 0.00206 44.00 20.4 4.9 13.6 5.5 1.6 0.00199 46.00 29.7 1.4 2.4 12.8 4.7 0.00168 51.00 9.4 0.8 1.1 4.0 6.7 0.00386 22.00 0.10695 0.80250 0.50923 0.46544 1.00000 114.00 16.00 22.00 30.00 16.00

35 44 46 51 22

198

Growth Factor Model • When a constant growth factor is applied to all OD pairs. • Let’s assume a 20% increase in the future year 1 2 3 4 5

1 24.6 24.1 18.2 28.5 9.1 104.5

2 1 6.6 4.2 1.3 0.7 13.8

3 1 3.8 13.2 2.5 1.2 21.7

4 2.7 3.3 4.3 10.7 3.4 24.4

5 0.9 1.2 1.1 3.4 4.9 11.5

30.2 39 41 46.4 19.3 175.9

1 2 3 4 5

1 29.5 28.9 21.8 34.2 10.9 125.4

2 1.2 7.9 5.0 1.6 0.8 16.6

3 1.2 4.6 15.8 3.0 1.4 26.0

4 3.2 4.0 5.2 12.8 4.1 29.3

5 1.1 1.4 1.3 4.1 5.9 13.8

36.2 46.8 49.2 55.7 23.2 211.1

Definitions Continued • T = Total trips = Σi Oi = Σj Dj

• Logical constraints which any feasible trip matrix must satisfy: – Σj Tij = Oi

– Σi Tij = Dj

i=1,…,N j=1,…,N

[1]

[2]

Proportional Flow Model • It is rather reasonable to assume that the flow between i and j is proportional to the total flow out of i and into j: – Tij ≅ Oi i=1,…,N – Tij ≅ Dj j=1,…,N – Tij = kOiDj

[3] [4] [5]

– k = 1/T

[6]

• where k is a constant of proportionality

• To solve for k, note that to be feasible [5] must satisfy constraints [1] & [2] ⇒ substitute [5] into [1] and solve for k; this yields: • Tij = OiDj/T

[7]

Σj Tij = Oi i=1,…,N Σj kOiDj = Oi Σj kDj = 1 k=1/ ΣjDj = 1/T

[1]

Proportional Flow Example Base Case Scenarios Zone 1 2 3 4 5 Total

1 24.6 24.1 18.2 28.5 9.1 104.5

2 1 6.6 4.2 1.3 0.7 13.8

3 1 3.8 13.2 2.5 1.2 21.7

4 2.7 3.3 4.3 10.7 3.4 24.4

5 Total 0.9 30.2 1.2 39 1.1 41 3.4 46.4 4.9 19.3 11.5 175.9

OiDj / T Tij1

Future Dj

Tij2 20.152 25.333 26.485 29.364 12.667 114

Tij3 2.828 3.556 3.717 4.121 1.778 16

Tij4 3.889 4.889 5.111 5.667 2.444 22

Future Oi

Tij5 5.303 6.667 6.970 7.727 3.333 30

2.828 3.556 3.717 4.121 1.778 16

35 44 46 51 22 198

Singly-Constrained Gravity Model • Proportional flow model does not generally represent travel well. It is insensitive to trip length, congestion, etc. – Add a third assumption:

Tij ≅ fij = f(1/tij) i,j=1,…,N [8] tij = travel time i to j fij = “impedance function”; ∂fij/ ∂ tij < 0

– Assumptions [3], [4] & [8] imply: Tij = kOiDjfij

[9]

k = 1/∑j’Dj’fij’

[10]

– As before, substitute [9] into [1] to solve for k – This yields:

Tij = OiDjfij / ∑ j’Dj’fij’

[11]

Gravity Model Cont’d • Eqn [11] does not, however, satisfy constraint [2], and so is not feasible. To ensure that both [1] and [2] are satisfied simultaneously, we must solve iteratively. • To do this, replace the Dj terms in [11] with “modified attraction” terms D*jx, which are defined for the xth iteration as follows

Gravity Model Balancing Procedure x=0

x for iteration

D*j1 = Dj x=x+1 Tijx = OiD*jxfij / Σj’D*j’xfij’ A

B

Balancing Procedure Cont’d A Rj = Dj / Σ iTijx |Rj - 1| < ε for all j?

Yes

No D*jx+1 = D*jxRj

B

Stop

Impedance Functions • Common impedance functions include: – fij = exp(βtij)

– fij = tijb

β<0

b<0

[12.1]

[12.2]

– β,b = parameters which must be estimated from observed data

“Entropy” Formulation • The gravity model balancing procedure can be rewritten more compactly in the following form • Define “balancing factors” as: – Ai = 1/ ∑ j’Bj’Dj’fij’

– Bj = 1/ ∑ i’Ai’Oi’fi’j – Tij = AiOiBjDjfij

[13.1]

[13.2]

[14]

• The iterative process begins with setting all Bj = 1 in eqn. [13.1]

“Entropy” Formulation Continued • It can be shown that: – D*j = BjDj

– and, hence, that Eqns [13] & [14] are equivalent to the gravity model balancing procedure

– It can also be shown that Eqns [13] & [14] are the “least-biased”, “most likely” estimates of the trip matrix which can be achieved, given the available data

Biproportional Updating • Gravity models often contain considerable errors

• This is not surprising, given that one is using a fairly simplistic model to represent very complex travel patterns • Instead of converting an observed trip matrix into a “synthetic” model, which is then used to predict the future, one can “update” the base directly

Synthetic Modelling Approach

Base O-D Matrix

Synthetic Model Tij = Model Calibration

Forecast Year Trip Ends

Oi

Dj

Forecast Year O-D Matrix

Base Data Updating Approach

Base O-D Matrix Forecast Year O-D Matrix Forecast Year Trip Ends

Oi

Dj

Base Data Updating Approach • Idea is to proportionally adjust the observed, base year trip matrix until it matches the forecast year row and column sums. An iterative procedure is required to balance rows and columns Let: – Oik = ∑ jTijk for the kth iteration

– Tij0 = base year O-D flow

– Oinew = forecast year row sum

– Similar definitions for Djk and Djnew

Biproportional Balancing Algorithm k=0 k=k+1 Tijk = Tijk-1(Oinew/Oik-1) No k=k+1 Tijk = Tijk-1(Djnew/Djk-1)

Convergence? Yes Stop

Balancing Iterations 1 Tij

1 2

2 12 34

22 11

Dj

46

33

Djnew

64

47

First Iteration Tij Dj Djnew/Dj Second Iteration Tij Dj Djnew/Dj

15.9 49.9 65.7

29.1 16.1 45.3

0.973

1.039

15.2 49.1 64.3

29.8 16.9 46.7

0.996

1.006

15.1 48.9 64.0

29.9 17.1 47.0

0.999

1.001

Third Iteration Tij Dj Djnew/Dj Fourth Iteration Tij Dj Djnew/Dj

15.1 48.9 64.0

29.9 17.1 47.0

1.000

1.000

Oi Oinew 34 45 45 66

Oi 45.0 66.0 Dj

Oi 45.0 66.0

Oi 45.0 66.0

Oi 45.0 66.0

15.5 48.5 64.0

30.2 16.8 47.0

15.2 48.8 64.0

30.0 17.0 47.0

15.1 48.9 64.0

29.9 17.1 47.0

15.1 48.9 64.0

29.9 17.1 47.0

Oinew/Oi 1.32 1.47

Oi 45.70 65.30

Oinew/Oi 0.98 1.01

Oi 45.11 65.89

Oinew/Oi 1.00 1.00

Oi 45.02 65.98

Oinew/Oi 1.00 1.00

Oi 45.00 66.00

Oinew/Oi 1.00 1.00

Gravity Model Parameter Estimation • If fij = exp(βtij) then one can show that the “most likely” value of β is the one which satisfies the equation: – ∑i ∑j tijTij / T = tavg

[15]

– where tavg is the observed average travel time and the left-hand-side of [15] is the average travel time predicted by the model – Eqn [15] can either be solved using trial & error or by applying the Newton-Raphson root-finding method

Goodness-of-Fit Statistics • Tij = Observed trips; T*ij = predicted trips

• T = Total trips; n = No of zones; T0 = T/n2 R2 =

∑i ∑j (Tij - T*ij)2

1- -------------------∑i ∑j (Tij - T0)2

χ2 = ∑i ∑j {(Tij - T*ij)2/T*ij}

MABSERR = {∑i ∑j |Tij - T*ij|}/n2

Normalized φ = ∑i ∑j (Tij/T)|log(Tij/T*ij)|

[16] [17] [18]

[19]

Model Evaluation Methods • Compare observed and predicted trip length frequency distributions (TFLD’s)

• Examine O-D residuals (perhaps on a super-zone basis) • Examine predicted versus observed screen-line counts

Observed versus Predicted Trip Lengths

30

obs pred

Trips

25

20

15

10

5

0 0

5

10

15

20

Travel Time

25

30

35

40

Model Estimation versus Calibration • Estimation is the process of finding model parameter values which cause the model to “best fit” observed data, according to some statistical procedure (regression, maximum likelihood, etc.)

• Calibration occurs post-estimation It involves ad hoc adjustments to model parameters which “force” the model to “better” fit observed data

Gravity Model Calibration: K-factors • The most common calibration procedure used for gravity models is to introduce K-factors into the equation: Tij = AiOiBjDjKijfij

Ai = 1/ ∑ j’Bj’Dj’Kij’fij’ Bj = 1/ ∑ i’Ai’Oi’Ki’jfi’j

[20]

[21.1] [21.2]

K-factors, cont’d • The Kij are chosen so as to reduce differences between observed and predicted: – screenline flows – TLFD’s and/or – key O-D pairs

K-factors, cont’d • The “rationale” for their use is to try to capture systematic differences in spatial flows not explained by travel time (or other terms in the gravity model)

– e.g., residential zones with high concentrations of white-collar workers are more likely to generate work trips to employment zones with high concentrations of white collar jobs than zones with predominantly blue collar jobs

• Use K-factors to link “white-collar” origin and destination zones together

Trip Distribution Models

Problem Definition, Terminology. • Growth Factor Model. • The Proportional Flow Model. • The Singly-constrained Gravity Model. • Bi-Proportional Updating ...

440KB Sizes 2 Downloads 242 Views

Recommend Documents

Distribution Forecasting in Nonlinear Models with ...
Nov 12, 2013 - A simulation study and an application to forecasting the distribution ... and Finance (Rotterdam, May 2013), in particular Dick van Dijk, for useful comments ...... below December 2008 forecasts”, and the Royal Bank of Scotland ...

Distribution Forecasting in Nonlinear Models with ...
Nov 12, 2013 - it lends itself well to estimation using a Gibbs sampler with data augmentation. ...... IEEE Transactions on Pattern Analysis and Machine Intelligence, 6:721–741, 1984. ... Business and Economic Statistics, 20:69–87, 2002.

Asymptotic distribution theory for break point estimators in models ...
Feb 10, 2010 - illustrated via an application to the New Keynesian Phillips curve. ... in the development of statistical methods for detecting structural instability.1.

The projection of species distribution models and the ...
... USDA Forest Service,. Southern Research Station, Asheville, NC 28804-3454, USA ... such novel conditions is not only prone to error (Heikkinen et al. 2006 ...

Sports Trip Driver and Field Trip Driver.pdf
Loading… Page 1. Whoops! There was a problem loading more pages. Main menu. Displaying Sports Trip Driver and Field Trip Driver.pdf.

magnificent - TOBU JAPAN TRIP
and yuba (bean curd skin) with them from their corners of Japan and those authentic flavors have been passed down unchanged through the generations.

Senior Trip 2018
Jun 1, 2018 - ... in your permission slips to Mrs. Sprague in Admin 1 or Ms. Newcomb in Admin 2 by. Thursday, May 24th (before Memorial Day weekend)​.

Field Trip Permission Form
Page 1. Sample Form.

magnificent - TOBU JAPAN TRIP
veloped in a fine mist giving it an air of mystery. In addition to Lake Chuzenji, the wondrous waters of Nikko can also be found at Kirikomi Lake, Saino Lake, and ...

Manual - Trip Advisor.pdf
Si puedes definir. el odio o el amor. Amigo que desilusión. No todo es blanco,. O negro: es gris. Todo depende del matiz,. Busca y aprende a distinguir. La luna puede calentar. Whoops! There was a problem loading this page. Retrying... Manual - Trip

Field Trip Request.pdf
(3) Transportation Department will forward approved form to the. Superintedent's Office. Page 1 of 1. Field Trip Request.pdf. Field Trip Request.pdf. Open. Extract.

Road Trip Instructions.pdf
Loading… Whoops! There was a problem loading more pages. Whoops! There was a problem previewing this document. Retrying... Download. Connect more apps... Try one of the apps below to open or edit this item. Road Trip Instructions.pdf. Road Trip Ins

Field-trip-permission.pdf
Page 1 of 1. Page 1 of 1. Field-trip-permission.pdf. Field-trip-permission.pdf. Open. Extract. Open with. Sign In. Main menu. Displaying Field-trip-permission.pdf.

Field Trip Forms.pdf
Loading… Whoops! There was a problem loading more pages. Whoops! There was a problem previewing this document. Retrying... Download. Connect more apps... Try one of the apps below to open or edit this item. Field Trip Forms.pdf. Field Trip Forms.pd

Choir Trip Combined.pdf
12:30 PM: Arrive at the Riverwalk – The downtown loop of the San Antonio River provides. unusual ways to see a portion of the city. Lunch on your own. 2:30 PM: Meet at THE ALAMO THEATER – Watch the IMAX “Alamo, the Price of Freedom.” (3pm sho

Don trip starlito
Revit 2016 pdf. Don trip starlito - Download.Don trip starlito.Good to great ... Mila milan balls.Iceage:continental drift hindi.Amy brookeand friends.Ar l.Xfacefox.

Field Trip Request.pdf
Loading… Whoops! There was a problem loading more pages. Whoops! There was a problem previewing this document. Retrying... Download. Connect more apps... Try one of the apps below to open or edit this item. Field Trip Request.pdf. Field Trip Reques

euro trip sub.pdf
There was a problem previewing this document. Retrying... Download. Connect more apps... Try one of the apps below to open or edit this item. euro trip sub.pdf.

Field Trip Form.pdf
Page 1 of 1. We understand the arrangement for taking your group to. We believe the necessary precautions and plans for the children's care have been made. We feel that reasonable vigilance in the care and supervision of the children during the. trip

Bathroom Trip Tallies.pdf
Page 1 of 1. Bathroom Trip Tallies.pdf. Bathroom Trip Tallies.pdf. Open. Extract. Open with. Sign In. Main menu. Displaying Bathroom Trip Tallies.pdf.

Senior Trip Policy.pdf
There was a problem previewing this document. Retrying... Download. Connect more apps... Try one of the apps below to open or edit this item. Senior Trip ...

Choir Trip Combined.pdf
Sign in. Loading… Whoops! There was a problem loading more pages. Whoops! There was a problem previewing this document. Retrying... Download. Connect ...

Field Trip Record.pdf
Loading… Whoops! There was a problem loading more pages. Whoops! There was a problem previewing this document. Retrying... Download. Connect more apps... Try one of the apps below to open or edit this item. Field Trip Record.pdf. Field Trip Record.