Advances in reproducible simulation of parcel deliveries based on open data

Sebastian Hörl

6 December 2021

HEADS-UP at IRT SystemX

Classic transport planning

• Zones
• Flows
• Peak hours
• User groups

Aggregated

Agent-based models

0:00 - 8:00

08:30 - 17:00

17:30 - 0:00

0:00 - 9:00

10:00 - 17:30

17:45 - 21:00

22:00 - 0:00

• Discrete locations
• Individual travelers
• Individual behaviour
• Whole day analysis

Disaggregated

Icons on this and following slides: https://fontawesome.com

MATSim

• Flexible, extensible and well-tested open-source transport simulation framework
   
• Used by many research groups and companies all over the world
   
• Extensions for parking behaviour, signal control, location choice, freight, ...

matsim-org/matsim-libs

MATSim

Synthetic demand

MATSim

Mobility simulation

Synthetic demand

MATSim

Decision-making

10:00 - 17:30

17:45 - 21:00

22:00 - 0:00

Mobility simulation

Synthetic demand

MATSim

Decision-making

Mobility simulation

Synthetic demand

MATSim

Decision-making

Mobility simulation

Analysis

Synthetic demand

MATSim @ IRT SystemX

• On-demand mobility simulation
  
   

• First/last mile mobility simulation

Passenger transport

MATSim @ IRT SystemX

• On-demand mobility simulation
  
   

• First/last mile mobility simulation

Passenger transport

Freight transport

• First/last mile delivery
  
  
  
  
  
  
  
  
  
   
• VRP Solvers
  • Development of electric VRP for JSprit

LEAD Project

• H2020 Project from 2020 to 2023
   
• Six living labs with different innovative logistics concepts
  • Lyon, The Hague, Madrid, Budapest, Porto, Oslo
  • One partner for implementation and one for research each
     
• Development of a generic modeling library for last-mile logistics scenario simulation and analysis

Lyon Living Lab

• Peninsula Confluence between Saône and Rhône
   
• Interesting use case as there are limited access points
   
• Implementation of an urban consolidation center (UCC) to collect the flow of goods and organize last-mile distribution
  • using cargo-bikes
  • using electric robots
  • and others
     
• Analysis and modeling through
  • Flow estimation through cameras
  • Simulation of future scenarios
     
• Focus on parcel deliveries due to data availability

Modeling methodology

• Main question: What is the impact on traffic and population of implementing an Urban Consolidation Center in Confluence?
   
• Focus on B2C parcel deliveries due to data availability

Modeling methodology

Agent-based simulation of Lyon

• Main question: What is the impact on traffic and population of implementing an Urban Consolidation Center in Confluence?
   
• Focus on B2C parcel deliveries due to data availability

Modeling methodology

Agent-based simulation of Lyon

Demand in parcel deliveries

• Main question: What is the impact on traffic and population of implementing an Urban Consolidation Center in Confluence?
   
• Focus on B2C parcel deliveries due to data availability

Modeling methodology

Agent-based simulation of Lyon

Demand in parcel deliveries

Distribution by operators

• Main question: What is the impact on traffic and population of implementing an Urban Consolidation Center in Confluence?
   
• Focus on B2C parcel deliveries due to data availability

Modeling methodology

Agent-based simulation of Lyon

Demand in parcel deliveries

Distribution by operators

KPI Calculation

• Main question: What is the impact on traffic and population of implementing an Urban Consolidation Center in Confluence?
   
• Focus on B2C parcel deliveries due to data availability

Methodology: Synthetic travel demand

Synthetic travel demand has been generated for Lyon in order to perform agent-based mobility simulation of all residents' movements.

Methodology: Synthetic travel demand

+

=

Synthetic travel demand has been generated for Lyon in order to perform agent-based mobility simulation of all residents' movements.

Methodology: Synthetic travel demand

=

Paper published July 2021

+

=

Synthetic travel demand has been generated for Lyon in order to perform agent-based mobility simulation of all residents' movements.

Methodology: Synthetic travel demand

eqasim-org/eqasim-java

eqasim-org/ile-de-france

Methodology: Parcel demand

+

Based on sociodemographic attributes of the households, parcels are generated for the city on an average day.

Gardrat, M., 2019. Méthodologie d’enquête: le découplage de l’achat et de la récupération des marchandises par les ménages. LAET, Lyon, France.

Methodology: Parcel demand

+

Based on sociodemographic attributes of the households, parcels are generated for the city on an average day.

Gardrat, M., 2019. Méthodologie d’enquête: le découplage de l’achat et de la récupération des marchandises par les ménages. LAET, Lyon, France.

Methodology: Parcel demand

+

Based on sociodemographic attributes of the households, parcels are generated for the city on an average day.

Gardrat, M., 2019. Méthodologie d’enquête: le découplage de l’achat et de la récupération des marchandises par les ménages. LAET, Lyon, France.

Methodology: Parcel demand

Gardrat, M., 2019. Méthodologie d’enquête: le découplage de l’achat et de la récupération des marchandises par les ménages. LAET, Lyon, France.

+

Based on sociodemographic attributes of the households, parcels are generated for the city on an average day.

Gardrat, M., 2019. Méthodologie d’enquête: le découplage de l’achat et de la récupération des marchandises par les ménages. LAET, Lyon, France.

Methodology: Parcel demand

Gardrat, M., 2019. Méthodologie d’enquête: le découplage de l’achat et de la récupération des marchandises par les ménages. LAET, Lyon, France.

+

Based on sociodemographic attributes of the households, parcels are generated for the city on an average day.

• Based on synthetic population, find average number of purchases delivered to a household defined by socioprofesional class, age of the reference person and household size per day.

Methodology: Parcel demand

Gardrat, M., 2019. Méthodologie d’enquête: le découplage de l’achat et de la récupération des marchandises par les ménages. LAET, Lyon, France.

+

Based on sociodemographic attributes of the households, parcels are generated for the city on an average day.

• We now the average number of parcels, but we do not now the distribution of the number of parcels for a household on an average day.
   
• We know it must be non-negative, and we know the mean.
   
• Without additional data, we assume maximum entropy distribution, which is Exponential in this case.

Methodology: Parcel demand

Gardrat, M., 2019. Méthodologie d’enquête: le découplage de l’achat et de la récupération des marchandises par les ménages. LAET, Lyon, France.

+

Based on sociodemographic attributes of the households, parcels are generated for the city on an average day.

Methodology: Parcel demand

Gardrat, M., 2019. Méthodologie d’enquête: le découplage de l’achat et de la récupération des marchandises par les ménages. LAET, Lyon, France.

+

Based on sociodemographic attributes of the households, parcels are generated for the city on an average day.

Methodology: Parcel demand

Gardrat, M., 2019. Méthodologie d’enquête: le découplage de l’achat et de la récupération des marchandises par les ménages. LAET, Lyon, France.

+

Based on sociodemographic attributes of the households, parcels are generated for the city on an average day.

Methodology: Parcel demand

Gardrat, M., 2019. Méthodologie d’enquête: le découplage de l’achat et de la récupération des marchandises par les ménages. LAET, Lyon, France.

+

Based on sociodemographic attributes of the households, parcels are generated for the city on an average day.

First case study

Hörl, S. and J. Puchinger (2021) From synthetic population to parcel demand: Modeling pipeline and case study for last-mile deliveries in Lyon, Working paper.

Solve VRP-TW based on generated parcels and household presence

JSprit

First case study

Hörl, S. and J. Puchinger (2021) From synthetic population to parcel demand: Modeling pipeline and case study for last-mile deliveries in Lyon, Working paper.

Solve VRP-TW based on generated parcels and household presence

JSprit

Scaling up: City-wide baseline

Using shares of customer preference, parcels are assigned to operators and the closest distribution centers.

Using shares of customer preference, parcels are assigned to operators and the closest distribution centers.

Scaling up: City-wide baseline

Using shares of customer preference, parcels are assigned to operators and the closest distribution centers.

Scaling up: City-wide baseline

Using shares of customer preference, parcels are assigned to operators and the closest distribution centers.

Scaling up: City-wide baseline

Scaling up: City-wide baseline

• Allows us to cut out a better baseline scenario for Confluence

First main result: Estimate of the total distance driven for parcel deliveries in Lyon (by the most frequent operators).

Recap

Parcel optimization

Baseline movements

Cut-out of scenarios

Coming up: KPI Calculation

GHG/PM Emissions: Simulation based on Hbefa database with detailed driving conditions (cold start, current speed, ...)

(Alternative in LEAD: COPERT)

Coming up: KPI Calculation

GHG/PM Emissions: Simulation based on Hbefa database with detailed driving conditions (cold start, current speed, ...)

(Alternative in LEAD: COPERT)

Noise Emissions: Simulation simulates noise emissions and emissions to arrive at exposure values

Coming up: KPI Calculation

GHG/PM Emissions: Simulation based on Hbefa database with detailed driving conditions (cold start, current speed, ...)

(Alternative in LEAD: COPERT)

Noise Emissions: Simulation simulates noise emissions and emissions to arrive at exposure values

Aggregation

• Global emissions are summed up over the territory (CO2)
• Local emissions need to be looked at in terms of exposure (PM)

Coming up: KPI Calculation

GHG/PM Emissions: Simulation based on Hbefa database with detailed driving conditions (cold start, current speed, ...)

(Alternative in LEAD: COPERT)

Noise Emissions: Simulation simulates noise emissions and emissions to arrive at exposure values

Aggregation

• Global emissions are summed up over the territory (CO2)
• Local emissions need to be looked at in terms of exposure (PM)

• Noise is local by definition, impact is calculated in terms of exposure (of the population) and expected social cost (in EUR)

Coming up: KPI Calculation

Energy: Based on vehicle consumption as a factor on the distance driven)

• Sum over all vehicles and distance

Congestion: Comparison of travel time of the population with and without delivery traffic added

Optional to be developed: Explicit simulation of lane-blocking behaviour due to deliveries

• Sum over all delays for the population

Coming up: Scenario dashboard

Caution: Current test results to exemplify the methodology. Model needs to be improved in multiple stages (see slide after).

Scenarios

Policy

• Closing whole confluence for parcel delivery
• Closing Confluence except boulevard

Consolidation centers

• One vs. multiple centers
• Different locations of consolidation centers

Consolidation fleet

• Electric vans
• Cargo-bikes
• Delivery robots
• Mix

Demand

• 10% / 20% / 100% of parcels

Validation

Traffic

• Road counts (open data for Lyon)

Air quality

• Observatoire Air Rhone-Alpes

Deliveries

• Data will be provided by La Poste

Gaps

Synthetic travel demand

Parcel demand

Parcel distribution

Route optimization

KPI calculation

Validation

Sendcloud survey is only a proxy. Ideally: Obtain data on use of parcel delivery operators through survey.

Gaps

Synthetic travel demand

Parcel demand

Parcel distribution

Route optimization

KPI calculation

Validation

Delivery center and fleet specifications are based on best guess information. Do we have more detailed data? (In progress for La Poste)

Gaps

Synthetic travel demand

Parcel demand

Parcel distribution

Route optimization

KPI calculation

Validation

Use heuristic optimization algorithm from JSprit instead of best-response assignment (currently running optimization)

Gaps

Synthetic travel demand

Parcel demand

Parcel distribution

Route optimization

KPI calculation

Validation

Emissions not detailed for the moment
Detailed Hbefa data set has been ordered

Gaps

Synthetic travel demand

Parcel demand

Parcel distribution

Route optimization

KPI calculation

Validation

Some KPIs need extensions for the simulation (lane-blocking through deliveries)

Gaps

Synthetic travel demand

Parcel demand

Parcel distribution

Route optimization

KPI calculation

Validation

Validation data needs to be assessed and integrated

General Freight Traffic

Synthetic population

"Freight population"

Simulation

Analysis

General Freight Traffic (FRETURB)

Synthetic population

FRETURB

Simulation

Analysis

• Internship in collaboration with Future Cities Lab (November 2021 - February 2022)
   
• Generation of B2B movements based on aggregated published statistics
   
• Based on methodology developed at LAET, Lyon

General Freight Traffic (MASS-GT)

Synthetic population

MASS-GT

Simulation

Analysis

• Academic exchange with TU Delft
  (November 2021 - February 2022)
   
• Generation of B2B shipments and associated movements based on behavioural models of actors in the logistics system (producers, shippers, receivers)

General Freight Traffic (MASS-GT)

Synthetic population

MASS-GT

Simulation

Analysis

• Academic exchange with TU Delft
  (November 2021 - February 2022)
   
• Generation of B2B shipments and associated movements based on behavioural models of actors in the logistics system (producers, shippers, receivers)
   
• Potential simulations with feedback loop

MATSim + MASS-GT

• Simulation of Lyon (and Rhônes-Alpes) with influence of traffic, including freight traffic

• Simulation of shipper behaviour (B2C and B2B)

• For instance: What will be the effect of a zero-emission zone in the city on traffic?

Questions?

https://slides.com/sebastianhorl/headsup-6dec2021