ART1
ART1
ART1
THÈSE
Présentée par
Walid BEHIRI
Pour obtenir le grade de
Docteur de l’université Paris-Est
Spécialité : Informatique
Composition du jury :
V.
VI.
Remerciements
Je tiens à exprimer ma profonde gratitude à tous ceux qui m’ont porté aide, assistance et
encouragements dans l'élaboration de cette recherche.
Mes remerciements s’adressent à Monsieur Mohamed Akil, pour son engagement en tant que
directeur de ma thèse, pour sa disponibilité, ses conseils et son soutien. Je tiens à exprimer mes
plus vifs remerciements à Mme Sana Berraf-Belmokhtar, qui en encadrant cette thèse a été au
plus près de mes questionnements, mais également de mes difficultés. Sa compétence, sa
rigueur intellectuelle et son esprit critique m’ont motivé et guidé tout au long de mon parcours
de recherche. Merci aussi pour les nombreuses relectures, les corrections, les encouragements
et l’intégration progressive dans ce monde qui m’était inconnu.
Que Messieurs Alexandre Dolgui, Professeur à l’IMT Atlantique, Chengbin Chu, Professeur à
l’ESIEE Paris, Eric Ballot, Professeur à Mines ParisTech et Joaquin Rodriguez, Directeur de
recherche à l’IFSTTAR de Lille, me permettent de leur exprimer ma gratitude pour avoir
accepté d’examiner ce manuscrit et pour leur participation au jury.
Je remercie également tous les doctorants et post-doctorants de l’ESIEE Paris, même si je les
ai trop peu croisés pendant la thèse, pour les échanges amicaux que nous avons pu avoir
ensemble. Je souhaite courage et réussite dans les plus brefs délais pour ceux qui n’ont pas
encore terminé.
Un remerciement particulier à toute ma famille. Mon père, mon frère et ma sœur ; mes premiers
reviewers et relecteurs. Merci de m’avoir toujours soutenu et encouragé, et surtout d’avoir cru
en moi et supporté tout au long de mon parcours.
Table des matières
Introduction Générale ................................................................................................................. 9
I. Chapitre I : Logistique et transport de marchandises en ville .......................................... 14
1. Introduction ................................................................................................................... 15
2. Logistique et transport de marchandises : définitions et concepts ................................ 15
3. Transport de marchandises en ville ............................................................................... 19
3.1. Enjeux du TMV ..................................................................................................... 20
3.2. Etat des lieux et pratiques actuelles ....................................................................... 21
3.3. Cadre légal et règlementaire .................................................................................. 24
3.4. Complexité du transport de marchandises en ville ................................................ 25
3.5. Evolution et perspectives futures du TMV ............................................................ 27
4. Exemple du TMV en Ile-de-France (Paris intra-muros et sa périphérie) ...................... 31
4.1. Infrastructures logistiques en Ile-de-France........................................................... 32
4.2. La gestion du transport de marchandises en région parisienne .............................. 34
5. Conclusion ..................................................................................................................... 36
II. Chapitre II : Transport de marchandises en ville par rail – mixité fret / voyageurs ......... 39
1. Introduction ................................................................................................................... 40
2. Pratiques et état des lieux : fret urbain par rail et mixité fret / voyageurs..................... 40
2.1. L’espace souterrain pour le transport de marchandises ......................................... 41
2.2. Mixité fret / voyageurs et utilisation du rail urbain : alternatives durables dans la
littérature .......................................................................................................................... 43
2.3. Mixité fret / voyageurs : expérimentations et solutions en cours d’exploitation ... 52
3. Conception d’un modèle de transport mixte fret / voyageurs dans le réseau ferroviaire
urbain .................................................................................................................................... 54
3.1. Composantes du réseau ferroviaire urbain de transport de voyageurs .................. 56
3.2. Identification des cas de mixité possibles .............................................................. 59
4. Problématiques décisionnelles liées à la mixité fret / voyageurs dans le réseau
ferroviaire urbain .................................................................................................................. 65
4.1. Conception de la solution de transport étudiée ...................................................... 65
4.2. Identification et description des problématiques décisionnelles ............................ 68
5. Conclusion ..................................................................................................................... 74
III. Chapitre III : Etat de l’art des méthodes de modélisation des problèmes de transport . 77
1. Introduction ................................................................................................................... 78
2
2. Problèmes classiques de la littérature liés à la logistique urbaine ................................. 78
2.1. Problèmes classiques liés au transport ................................................................... 78
2.2. Problème d’affectation généralisée ........................................................................ 80
3. Méthodes de résolution de la RO .................................................................................. 81
3.1. Méthodes exactes ................................................................................................... 82
3.2. Méthodes approchées ............................................................................................. 87
4. Problèmes de replanification ......................................................................................... 90
4.1. Le principe d’horizon glissant et résolution d’un problème de replanification ..... 91
4.2. La replanification dans la gestion de la chaine logistique ..................................... 92
4.3. La replanification dans le transport ferroviaire ...................................................... 94
5. Simulation à évènements discrets – ARENA ................................................................ 96
5.1. Généralités sur la simulation .................................................................................. 96
5.2. Démarche de réalisation d’un modèle de simulation ............................................. 97
5.3. Simulation avec ARENA ....................................................................................... 99
6. Couplage simulation / optimisation............................................................................. 102
7. Conclusion ................................................................................................................... 105
IV. Chapitre IV : Modélisation et simulation du FRTSP .................................................. 107
1. Introduction ................................................................................................................. 108
2. Hypothèses et description du FRTSP .......................................................................... 108
2.1. Hypothèses du FRTSP ......................................................................................... 109
2.2. Formalisation du FRTSP ...................................................................................... 111
3. Une modélisation par simulation du FRTSP ............................................................... 114
3.1. Génération des commandes ................................................................................. 114
3.2. Génération des trains ............................................................................................ 116
3.3. La ligne ferroviaire............................................................................................... 118
3.4. Implémentation des processus décisionnels sur ARENA .................................... 119
3.5. Modèle conceptuel du modèle de simulation ....................................................... 123
4. Conclusion ................................................................................................................... 126
V. Chapitre V : Modélisation mathématique et optimisation du FRTSP – cas prédictif et
replanification......................................................................................................................... 127
1. Introduction ................................................................................................................. 128
2. Un exemple numérique ............................................................................................... 128
3. Complexité du problème FRTSP ................................................................................ 131
4. Formulation d’un PLVM du FRTSP dans le cas prédictif .......................................... 132
3
5. L’adaptation des colonies de fourmis pour le FRTSP................................................. 136
5.1. La variante AS des ACO ...................................................................................... 136
5.2. Différentes variantes d’ACO ............................................................................... 139
6. Une approche par horizon glissant pour optimiser la replanification du FRTSP ........ 144
6.1. Un PLVM pour la minimisation des changements lors de la replanification du
FRTSP ............................................................................................................................ 147
6.2. L’approche de replanification par horizon glissant.............................................. 148
7. Couplage simulation / optimisation............................................................................. 149
7.1. Cas prédictif ......................................................................................................... 149
7.2. Cas de replanification .......................................................................................... 150
8. Conclusion ................................................................................................................... 151
VI. Chapitre VI : Expérimentations et analyse des résultats ............................................. 153
1. Introduction ................................................................................................................. 154
2. Génération des instances ............................................................................................. 154
2.1. Paramétrage du modèle de simulation ................................................................. 154
2.2. Le cas prédictif ..................................................................................................... 155
2.3. Le cas de la replanification .................................................................................. 155
3. Analyse des résultats du couplage simulation / optimisation ...................................... 156
3.1. PLVM dans le cas prédictif.................................................................................. 156
3.2. Les heuristiques ................................................................................................... 160
3.3. Comparaison et évaluation du gap des heuristiques VS PLVM .......................... 161
3.4. La replanification ................................................................................................. 165
4. Les colonies de fourmis............................................................................................... 168
4.1. Introduction à la méthode Taguchi ...................................................................... 168
4.2. Paramétrage des variantes de colonies de fourmis............................................... 171
4.3. Résultats des différentes variantes de colonies de fourmis .................................. 178
5. Conclusion ................................................................................................................... 181
Conclusion générale et perspectives....................................................................................... 183
Bibliographie .......................................................................................................................... 187
Liste des publications ............................................................................................................. 201
Résumé / Abstract .................................................................................................................. 203
4
Liste des Figures
Figure I.1 : Schématisation de la chaine logistique .................................................................. 16
Figure I.2 : schématisation de la relation entre les acteurs principaux de la chaine logistique à
travers des dépôts et plateforme de correspondance (Islam, et al., 2013) ................................ 16
Figure I.3 : Segmentation du transport de marchandises en France (Pimor & Fender, 2008) . 17
Figure I.4 : Evolution des parts modales de transport en France ............................................. 17
Figure I.5 : Les acteurs du TMV et leurs intérêts (Courivault, 2004) ...................................... 22
Figure I.6 : Répartition des livraisons / enlèvement au cours d’une journée ordinaire ............ 23
Figure I.7 : Répartition des activités commerciales en ville et mouvements générés
(Courivault, 2004) .................................................................................................................... 23
Figure I.8 : Les différentes formes de transport de marchandises en ville ............................... 26
Figure I.9 : Mode de livraison directe ...................................................................................... 27
Figure I.10 : Mode de livraison par tournée - optimisée par transporteur (Interface Transport,
2006)......................................................................................................................................... 27
Figure I.11 : Mode de livraison par tournée - optimisée par destinataire (Interface Transport,
2006)......................................................................................................................................... 28
Figure I.12 : Typologie des ELU en fonction de la couverture spatiale de la ville (Boudouin,
2006)......................................................................................................................................... 28
Figure I.13 : Evolution du trafic des véhicules de TMV à Paris et sa périphérie – moyenne
journalière................................................................................................................................. 32
Figure I.14 : Infrastructures et sites logistiques dans la région Ile-de-France (Paris Region,
2017)......................................................................................................................................... 33
Figure II.1 : Modèles 3D de la solution CargoCap : a- mise à l'echelle de la solution par
rapport à un metro, b- système de chargement / déchargement automatique, c- insertion de la
solution dans l'espace souterrain (CargoCap, 2002) ................................................................ 42
Figure II.2 : Modélisations de Cargo Sous Terrain : a- principe de fonctionnement, b- coupes
latérale et horizontale, c- modèle 3D des véhicules en développement (CST, 2016) .............. 43
Figure II.3 : Train de passagers modifié pour le transport de marchandises (Nuzzolo, et al.,
2008)......................................................................................................................................... 45
Figure II.4 : Photos de différentes expérimentations de fret urbain par rail ............................ 54
Figure II.5 : Réseau ferroviaire en Ile-de-France ..................................................................... 55
Figure II.6 : Les trois formes de stations ferroviaires en milieu urbain ................................... 56
Figure II.7 : Illustration des implantations possibles de voies ferrées ..................................... 57
Figure II.8 : Exemple de matériels roulant de chacun des modes de transport en région Ile-de-
France ....................................................................................................................................... 58
Figure II.9 : Les niveaux de mutualisation des différentes ressources du système de transport
ferroviaire pour une mixité fret / voyageurs ............................................................................. 60
Figure II.10 : Illustration du cas de mixité totale ..................................................................... 62
Figure II.11 : Schématisation d'un réseau de transport ferroviaire en milieu urbain ............... 67
Figure II.12 : Illustration du processus de transport de marchandises à travers le réseau
ferroviaire de la région Ile-de-France ....................................................................................... 68
5
Figure II.13 : Identification et décomposition des problématiques décisionnelles pour la
mixité des flux fret / voyageurs ................................................................................................ 69
Figure II.14 : Exemple d'une zone de stockage temporaire ..................................................... 72
Figure II.15 : Exemple de rangement de colis à l'intérieur du train (cas de la 2D) .................. 73
Figure III.1 : Schématisation du problème d'affectation généralisée dans le cas de
marchandises à livrer ................................................................................................................ 81
Figure III.2 : Méthodes de résolution pour l'optimisation combinatoire.................................. 82
Figure III.3 : Exemple d'une résolution graphique d'un PL à 5 contraintes ............................. 83
Figure III.4 : Illustration du principe d'optimalité de Bellman ................................................ 85
Figure III.5 : Illustration d'un graphe pour la recherche du chemin le plus court entre x et y . 86
Figure III.6 : Illustration de l'approche par horizon glissant .................................................... 92
Figure III.7 : Schématisation du processus de simulation (Pritsker, 1986) .............................. 97
Figure III.8 : Réalisation d’un modèle de simulation ............................................................... 98
Figure III.9 : Détail de la démarche pour la réalisation d’un modèle de simulation ................ 99
Figure III.10 : Evènements VBA lors de l’exécution de la simulation .................................. 102
Figure III.11 : Illustration de la démarche de couplage – optimisation pour simulation
(Borodin, et al., 2017) ............................................................................................................ 103
Figure III.12 : Illustration de la démarche de couplage – simulation pour optimisation
(Borodin, et al., 2017) ............................................................................................................ 103
Figure III.13 : Illustration de la démarche de couplage – simulation et optimisation (Borodin,
et al., 2017) ............................................................................................................................. 105
Figure IV.1 : Illustration des files d’attente dans le FRTSP .................................................. 112
Figure IV.2 : Exemple de conteneurs pour le transport de marchandises .............................. 113
Figure IV.3 : Schéma du modèle de simulation du FRTSP en utilisant ARENA .................. 115
Figure IV.4 : Modélisation de la création du plan de commandes périodique ....................... 115
Figure IV.5 : Modélisation de la génération de commandes en temps réel ........................... 116
Figure IV.6 : Modélisation du processus de génération des trains ......................................... 117
Figure IV.7 : Modélisation de la création des trains .............................................................. 117
Figure IV.8 : Modélisation des stations et des voies ferrées .................................................. 118
Figure IV.9 : Illustration de la construction d’une solution avec la 4ème heuristique ............. 121
Figure IV.10 : Modèle conceptuel du modèle de simulation ................................................. 125
Figure V.1 : Une illustration d’une instance numérique du FRTSP ...................................... 129
Figure V.2 : Réduction du problème FRTSP à un problème d’affectation généralisée ......... 132
Figure V.3 : Une étape durant le processus de construction d’une solution par une fourmi . 137
Figure V.4 : La replanification du FRTSP par horizon glissant ............................................. 149
Figure V.5 : Le schéma de couplage dans le cas prédictif ..................................................... 150
Figure V.6 : Le schéma de couplage dans le cas de la replanification ................................... 151
Figure VI.1 : Influence du nombre de trains sur les temps de résolution du PLVM ............. 159
Figure VI.2 : Influence du nombre de trains sur les temps d’attente ..................................... 159
Figure VI.3 : Comparaison des temps d’attente des heuristiques avec ceux du PLVM ........ 160
Figure VI.4 : Calculs des temps d’attente moyens ................................................................. 162
Figure VI.5 : Boites à moustaches des temps d’attente des commandes dans les stations –
(gauche) comparaison entre les différentes approches, (droite) heuristique 4 ....................... 163
Figure VI.6 : illustration du calcul de « UB » ........................................................................ 164
6
Figure VI.7 : Moyennes des temps d’attente des commandes – Colonies de fourmis VS
PLVM ..................................................................................................................................... 179
Figure VI.8 : Etude bi-critère des performances des colonies de fourmis ............................. 180
Figure VI.9 : Comparaison de convergence des différentes variantes ................................... 181
7
Liste des Tableaux
Tableau I.1 : Répartition des transporteurs selon le PTAC des véhicules utilisés (Patier, 2004)
.................................................................................................................................................. 22
Tableau II.1 : Comparaison entre les caractéristiques de transport de voyageurs /
marchandises (Delaitre, 2014).................................................................................................. 41
Tableau V.1 : Plan de transport des commandes – 1ère solution ............................................ 130
Tableau V.2 : Plan de transport des commandes – 2ème solution ........................................... 131
Tableau V.3 : Notations utilisées dans le PLVM ................................................................... 133
Tableau V.4 : Notations additionnelles pour le cas de replanification ................................... 147
Tableau VI.1 : Synthèse des résultats du PLVM dans le cas prédictif ................................... 157
Tableau VI.2 : Résultats de la résolution d’instances (150 commandes, 30 trains, 10 stations)
................................................................................................................................................ 157
Tableau VI.3 : Nombre d’instances infaisables ..................................................................... 158
Tableau VI.4 : Résultats des heuristiques VS solutions optimales ........................................ 161
Tableau VI.5 : Résultats de la simulation sur les temps d’attente des commandes ............... 161
Tableau VI.6 : Synthèse des cas d’infaisabilité ...................................................................... 163
Tableau VI.7 : comparaison des taux d’utilisation des stations ............................................. 165
Tableau VI.8 : Synthèse des résultats de la seconde itération de la replanification ............... 166
Tableau VI.9 : Résultats de la 1ère itération – plan initial....................................................... 166
Tableau VI.10 : Liste des nouvelles commandes ................................................................... 167
Tableau VI.11 : Résultats de 2ème itération – 1ère replanification ........................................... 167
Tableau VI.12 : Illustration d’une table de Taguchi 𝑳𝟖(𝟐𝟕)................................................. 170
Tableau VI.13 : Niveaux des paramètres pour les variantes AS et MMAS ........................... 171
Tableau VI.14 : Niveaux des paramètres pour la variante ACS ............................................ 172
Tableau VI.15 : Table de Taguchi 𝑳𝟏𝟔(𝟒𝟓).......................................................................... 173
Tableau VI.16 : Table de Taguchi 𝑳𝟐𝟓(𝟓𝟔).......................................................................... 173
Tableau VI.17 : Résultats des essais pour la variante AS ...................................................... 174
Tableau VI.18 : Ratio S/N pour la variante AS ...................................................................... 174
Tableau VI.19 : Résultats des essais pour la variante MMAS ............................................... 175
Tableau VI.20 : Ratio S/N pour la variante MMAS .............................................................. 175
Tableau VI.21 : Résultats des essais pour la variante ACS ................................................... 175
Tableau VI.22 : Ratio S/N pour la variante ACS ................................................................... 176
Tableau VI.23 : Résultats des essais pour la variante MMACS ............................................ 176
Tableau VI.24 : Ratio S/N pour la variante MMACS ............................................................ 177
Tableau VI.25 : Résultats des variantes colonies de fourmis VS solutions optimales ........... 179
8
Introduction générale
Introduction Générale
9
Introduction générale
Cette situation n’est pas propre à une région particulière, mais est largement répandue en
milieu urbain. En Europe, on attend le triplement des véhicules-kilomètres en milieu urbain à
l’horizon 2050, en particulier, pour le transport de marchandises (Van Audenhove, De Jongh,
& Durance, 2015). Ceci s’explique d’une part, par la densité de la population urbaine qui
représente 73% avec une tendance croissante (Nations United, 2015). D’autre part, les
changements du mode de consommation favorisent le e-commerce (Comi & Nuzzolo, 2016).
10
Introduction générale
Nos travaux sont structurés en 6 chapitres, que nous décrivons comme suit :
11
Introduction générale
12
Introduction générale
Développement de la solution de
transport mixte
Discussions avec des experts Revue de littérature Observations sur le terrain
Méthodologie de la recherche
13
Chapitre I
14
CHAPITRE I : Logistique et transport de marchandises en ville
1. Introduction
Au moment où la circulation de l’information et le transfert des données se sont affranchies
des distances physiques grâce à l’évolution des technologies de télécommunications, la
circulation des biens physiques se voit à son tour confrontée au même challenge. En effet,
aujourd’hui les produits ne sont plus consommés localement, mais traversent parfois des
milliers de kilomètres avant d’arriver chez le consommateur final. Pour répondre à cette
nouvelle réalité, le transport de marchandises a dû s’adapter et se développer, ce qui l’a fait
évoluer et a fait émerger de nouvelles problématiques. Le transport de marchandises fait
partie d’une discipline beaucoup plus complète qui est la logistique. Ainsi, la section 2 de ce
chapitre est dédiée à la définition des concepts principaux de la logistique, ainsi qu’à
l’exploration des spécificités du transport de marchandises et de ses enjeux de manière
générale.
L’idée qui a fait naitre ce travail de recherche, est la mutualisation du réseau de transport
ferroviaire de voyageurs en Ile-de-France, pour absorber une partie du transport urbain de
marchandises. Pour comprendre les motivations d’une telle proposition et le contexte d’une
possible mise en œuvre, dans la section 4 de ce chapitre, nous allons nous focaliser sur l’état
des lieux en matière de transport de marchandises dans cette région. Mais avant d’aborder ces
différents points, nous allons définir et rappeler quelques concepts clés.
Le Supply Chain Council définit la logistique comme étant : « la suite des étapes de
production et distribution d’un produit depuis les fournisseurs des fournisseurs du producteur
jusqu’aux clients de ses clients ». Aussi, la logistique peut se résumer en une équation (Islam,
et al., 2013) :
Une schématisation de la logistique est donnée par la Figure I.1. La logistique est composée
de 5 éléments clés (Islam, et al., 2013) : transport, entreposage, inventaire, emballage et
traitement de l’information. En général, le transport est la principale composante de la plupart
des services logistiques. Les aspects clés de la gestion des transports comprennent les modes
15
CHAPITRE I : Logistique et transport de marchandises en ville
de transport (tels que les voies par la route, ferroviaires, fluviales, aériennes, par pipelines,
multimodales ou intermodales), les infrastructures de transport, les conditions géographiques,
le type de livraison (par ex. express, normal, longue distance, dernier kilomètre…), la
planification du chargement, ainsi que l'itinéraire et la planification du transport.
Figure I.2 : schématisation de la relation entre les acteurs principaux de la chaine logistique à travers des dépôts et
plateforme de correspondance (Islam, et al., 2013)
Concernant les modes utilisés pour le transport de marchandises, en France, quatre modes se
partagent cette activité : la route, le ferroviaire, le fluviale, et par pipeline. La segmentation du
16
CHAPITRE I : Logistique et transport de marchandises en ville
marché du transport suivant la granularité des marchandises à transporter est donnée par la
Figure I.3. On remarque que le mode de transport ferroviaire assure, principalement, le
transport de marchandises de masse conditionnées et les marchandises en lots complets.
L’évolution de la part de chacun des modes de transport cités précédemment, sur les 30
dernières années est donnée par la Figure I.4 (données de (Insee, 2016)). On remarque que la
part du mode routier a augmenté tout au long de la période jusqu’en 2011 au détriment des
trois autres modes. Depuis 2010, on remarque le retour vers la croissance du mode de
transport ferroviaire, avec une moyenne annuelle de croissance de 3%.
Figure I.3 : Segmentation du transport de marchandises en France (Pimor & Fender, 2008)
100,00%
90,00%
80,00%
70,00%
60,00%
Pipelines
Transport ferroviaire
40,00%
Transport routier
30,00%
20,00%
10,00%
0,00%
1985
1986
1987
1988
1989
1990
1991
1992
1993
1994
1995
1996
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
2008
2009
2010
2011
2012
2013
2014
2015
17
CHAPITRE I : Logistique et transport de marchandises en ville
que les modes de consommation ont évolué (ex : produits personnalisés et juste à temps, e-
commerce…). Ainsi, le mode de transport routier est le plus adapté pour le transport de
produits en petite quantité (i.e. petits lots ou colis). En contrepartie, le mode de transport
ferroviaire qui est plus adapté au transport des produits de masse ou de lots importants, a subi
un recul de sa part dans le transport de marchandises.
Le transport routier domine les autres modes pour diverses raisons. L’association européenne
des véhicules de logistique mentionne quelques-uns de ses avantages comme suit :
Cependant, ce mode présente plusieurs inconvénients qui ont stoppé sa croissance depuis
2011, dont les principales :
- Les restrictions légales (au niveau Européen) et législatives (au niveau de chaque
pays) sur les temps de conduite, les tonnages et dimensions, les horaires d’accès aux
villes, etc.
- La fluctuation des coûts de transport en fonction des marchandises à transporter et de
la quantité. En effet, il n’est pas très adapté pour le transport de grandes quantités sur
de longues distances, où les prix deviennent moins compétitifs.
Quant au mode de transport ferroviaire, plusieurs de ses inconvénients ont fait obstacle à son
développement, selon l’association européenne des véhicules de logistique on peut citer :
Toutefois, ce mode de transport présente plusieurs avantages, qui font qu’il renoue avec la
croissance depuis 2010. Parmi ces avantages, on peut citer :
18
CHAPITRE I : Logistique et transport de marchandises en ville
Le transport intermodal nécessite d’avoir une vision globale du processus de transport, dans le
but d’identifier le moment et le lieu idéal pour procéder à une rupture de charge. La question
relative à la sélection du meilleur mode de transport, apporte une réponse plus complète, avec
des possibilités de faire du door-to-door, notamment avec l’utilisation de véhicules
électriques, pour assurer le dernier kilomètre par exemple. Le transport intermodal est défini
par (Macharis & Bontekoning, 2004) comme la combinaison d'au moins deux modes de
transport dans une seule chaîne de transport, sans changement de conteneur pour les
marchandises, avec la plupart des itinéraires parcourus par voies ferroviaires, par voies
navigables ou océaniques et avec les parcours initiaux et finaux les plus courts possibles par la
route.
Ces notions restent valables à l’échelle de la ville. Cependant, les spécificités des zones
urbaines nécessitent d’apporter quelques adaptations et de définir certaines notions
complémentaires.
19
CHAPITRE I : Logistique et transport de marchandises en ville
hétérogènes, qui sont impliqués dans le TMV de manière assez peu efficace, en raison
d’objectifs divergents ou bien d’un manque de coordination. Aujourd’hui, le TMV gagne de
l’importance aux deux niveaux :
- Economique : il doit absorber une activité économique qui connait une importante
croissance en milieu urbain.
- Social : il représente une part qui va continuer à croitre, de l’activité professionnelle
des citadins.
Le TMV est considéré comme étant nécessaire pour la ville, malgré les désagréments
engendrés en termes de réduction de la qualité de vie pour les citadins (ex : divers types de
pollutions, densification du trafic routier, etc.). Un autre effet néfaste est lié à la charge de
travail générée pour les collectivités (ex : gestion des flux supplémentaires de véhicules en
ville et des incidents qui y sont liés, adaptation des infrastructures de la ville aux spécificités
des véhicules lourds, gestion de l’impact environnemental, etc.).
Plus généralement, l’étude des flux de marchandises en ville reste très difficile et peu d’études
de terrain existent. En particulier, en raison des spécificités du mode routier (principal mode
de TMV) qui ne laisse entrevoir que peu de visibilité sur ses pratiques (particulièrement en ce
qui concerne : les types de marchandises transportées, la quantité par véhicules, le nombre de
livraisons pour compte propres, etc.), ainsi, les données disponibles sont peu précises
(Augereau, 2009).
Dans cette section, nous allons nous intéresser aux enjeux du TMV, aux pratiques actuelles et
cadre légal et juridique qui le régissent, ainsi qu’aux propositions d’évolution de ce dernier,
pour prendre en charge sa complexité et réduire ses nuisances.
20
CHAPITRE I : Logistique et transport de marchandises en ville
21
CHAPITRE I : Logistique et transport de marchandises en ville
des transporteurs selon le PTAC (Poids Total Autorisé en Charge) des véhicules
utilisés, ainsi que les usagers. Nous remarquons que la portion majoritaire du TMV
« 53% » est assurée par des véhicules légers (PTAC < 3,5 t), dont 78% pour compte
propre. Ceci s’explique par l’adéquation de ce type de véhicule avec les spécificités de
la ville.
- Une troisième caractéristique concerne un phénomène très répandu du TMV, à savoir
le stationnement illicite pour procéder aux livraisons / enlèvements (ex : en double file
ou sur les trottoirs). Certes, les temps d’arrêt sont faibles (60% des livraisons durent
moins de 5 minutes), cependant, en centre-ville la part des stationnements illicite est
d’en moyenne 60% et peut dépasser, dans certain cas, les 80%. Cette pratique
engendre des situations dangereuses, aussi bien pour les citadins, que pour les
transporteurs.
Tableau I.1 : Répartition des transporteurs selon le PTAC des véhicules utilisés (Patier, 2004)
22
CHAPITRE I : Logistique et transport de marchandises en ville
Figure I.6 : Répartition des livraisons / enlèvement au cours d’une journée ordinaire
Figure I.7 : Répartition des activités commerciales en ville et mouvements générés (Courivault, 2004)
- Dernier élément de la situation actuelle (et non des moindres), c’est la génération des
nuisances sous ses différentes formes. Comme nous l’avons évoqué au point
précèdent, l’impact environnemental du TMV est très important, cependant, d’autres
types de nuisances sont à déplorer :
Les nuisances sonores : à faible vitesse, les moteurs diesels des véhicules de
transport de marchandises génèrent beaucoup de bruit. De la même manière, à
23
CHAPITRE I : Logistique et transport de marchandises en ville
grande vitesse, c’est les roulements sur la chaussée qui en sont la cause.
D’autres bruits sont à constater, tel que le groupe frigorifique (des camions qui
en sont équipés), ouverture et fermeture des portes, l’utilisation du hayon
élévateur, etc.
Les mauvaises odeurs des gaz d’échappement.
L’accentuation de la congestion des rues.
La dégradation de la chaussée et des trottoirs lors du stationnement des
véhicules lourds.
Impact négatif sur le paysage urbain.
24
CHAPITRE I : Logistique et transport de marchandises en ville
comme cette activité fait intervenir trois composantes (les transporteurs, les
véhicules de transport et les espaces de livraisons), focalisons-nous sur ce que
mentionne la loi à leur sujet. D’abord, en ce qui concerne les professionnels qui
exercent dans cette activité, à travers l’utilisation de véhicules légers (<3,5 t),
ce n’est que depuis 2004 que le statut de chauffeur-livreur de véhicules
utilitaires légers existe. Ensuite, les véhicules utilitaires légers, qui sont
principalement utilisés pour les livraisons en ville, n’ont aucune définition
fonctionnelle. A la différence des poids lourds, leurs utilisations ne sont pas
encadrées par la loi ainsi, leur usage réel n’est pas officiellement connu. Enfin,
les espaces réservés à l’arrêt des véhicules effectuant un chargement ou un
déchargement de marchandises, n’ont été définis par la loi qu’en 2000.
Cependant, il manque de la précision en ce qui concerne la définition des temps
d’arrêt pour procéder aux livraisons.
- Au niveau local, la règlementation consiste en des textes adoptés par les instances
locales, à mettre en œuvre dans la localité correspondante. La circulation et le
stationnement en ville relèvent des prérogatives de la Mairie. Ainsi, c’est à cette
dernière de définir les différentes règlementations à appliquer à ce niveau. Que cela
soit pour les gabarits des véhicules, les horaires de livraison et les véhicules autorisés
pour chaque horaire, ou la durée de stationnement, chaque localité fixe ses propres
règles. Ce qui conduit sur le terrain, à des situations peu pratiques pour les
agglomérations. A titre d’exemple, en Ile-de-France, en 2004 il y avait 120
réglementations différentes (heures possibles d’entrée dans les différents quartiers et
gabarit de véhicules autorisés, itinéraires possibles en fonction des sens de circulation,
durée et lieux de stationnement pour livrer). Depuis, des guides d’harmonisation de la
règlementation ont été mis en place, mais un flou persiste au niveau des transporteurs
et logisticiens. Aussi, il faut noter que les réglementations locales diffèrent
énormément d’une agglomération à une autre, ce qui engendre pour des transporteurs
qui interviennent dans des agglomérations différentes, un surplus de complexité
technique, organisationnel et administrative.
25
CHAPITRE I : Logistique et transport de marchandises en ville
Les achats dans les commerces de détail, qui se situent dans les centres villes
(magasins de proximité).
- Les flux de gestion urbaine : ils représentent 10% des mouvements liés au transport de
marchandises. Ils correspondent à deux types de flux :
- Mode de livraison directe (Figure I.9) : il s’agit de transporter les marchandises d’une
commande, à partir d’un seul lieu d’expédition, à destination d’un seul lieu de
livraison.
26
CHAPITRE I : Logistique et transport de marchandises en ville
- Mode de livraison par tournée : il s’agit d’effectuer un certain parcours en ville, pour
livrer plusieurs établissements. Ce mode de livraison présente deux variantes, en
fonction de l’organisation du transporteur et des infrastructures logistiques de la ville.
Figure I.10 : Mode de livraison par tournée - optimisée par transporteur (Interface Transport, 2006)
27
CHAPITRE I : Logistique et transport de marchandises en ville
(Dufour, et al., 2007) présente une synthèse de certaines perspectives intéressantes, qu’on
présente à travers les points suivants.
Figure I.11 : Mode de livraison par tournée - optimisée par destinataire (Interface Transport, 2006)
Figure I.12 : Typologie des ELU en fonction de la couverture spatiale de la ville (Boudouin, 2006)
A l’échelle de la ville (ou partie de la ville, dans le cas des agglomérations très étendues),
c’est les CDU qui sont proposés. En effet, ils représentent l’interface logistique entre la ville
et l’extérieur et permettent une gestion des flux entrant et sortant homogènes, en groupant et
28
CHAPITRE I : Logistique et transport de marchandises en ville
dégroupant les marchandises de manière à optimiser leur transport. Selon (Browne, et al.,
2005), trois types de CDU sont à considérer :
- CDU privés / semi-privés : ils sont mis en place par des opérateurs ou chargeurs pour
leur propre compte, dans le but d’optimiser économiquement leurs activités.
- CDU mutualisés : ils sont mis en place par les autorités publiques, dans le but
d’améliorer la qualité de service de TMV et de minimiser ses impacts négatifs sur la
ville.
- CDU spécifiques : ils sont déployés de manière temporaire, pour assurer
l’approvisionnement d’activités temporaires en ville, tels que les chantiers de
construction. Il peut y avoir aussi des CDU de ce type mis en place de manière
permanente, dans le but de servir de plateforme de transbordement pour les aéroports
et ports fluviaux ou maritimes.
A l’échelle Européenne, plusieurs expériences de CDU ont été menées. On compte plus d’une
dizaine de CDU en cours d’exploitation en France, ainsi qu’en Italie, un autre à San Sebastian
en Espagne, etc. Cependant, plusieurs autres CDU ont été suspendus / arrêtés en raison de la
non rentabilité de ces derniers. Une étude de ces espaces logistiques a été réalisée dans
(Gonzalez-Feliu, et al., 2013), où trois facteurs de succès ont été identifiés :
- Lyon : limitation du temps de stationnement à 10 mn pour tous les usagers, ce qui crée
éventuellement, plus de places de stationnement pour les livraisons.
- Saint-Denis de la réunion : contrôle des zones dédiées aux livraisons par des agents,
dans le but d’assurer le respect de la règlementation.
- Amsterdam : utilisation d’un bateau par DHL comme arrière base pour les livraisons,
et livraison du dernier kilomètre par vélo, qui ne requiert pas de zone de
stationnement.
29
CHAPITRE I : Logistique et transport de marchandises en ville
- Barcelone : les voies latérales des rues avec une forte activité commerciale et une forte
circulation, sont réservées aux livraisons.
- Copenhague, Aalborg et Aarhus : mise en place d’une réglementation qui permet des
arrêts étendus pour les entreprises qui respectent certaines conditions, qui contribuent
à l’amélioration de la qualité de vie en ville.
Toutes ces expériences et plusieurs autres, dont certaines ont réussi et d’autres ont échoué,
montrent que les solutions à déployer dans les différentes villes seront forcément hétérogènes,
puisqu’ils devront s’adapter aux spécificités des lieux de leur mise en œuvre. Aussi, elles
montrent qu’il n’y a pas encore de consensus sur des solutions qui satisfont les besoins et les
contraintes des transporteurs, tout en réduisant l’impact sur la qualité de vie en ville.
30
CHAPITRE I : Logistique et transport de marchandises en ville
De plus, l’évolution du mode de consommation (achats plus fréquents toute l’année avec des
pics saisonniers) et l’impact qu’engendre la croissance continue du e-commerce, présagent
une explosion du nombre de kilomètres parcourus, pour transporter des colis de plus en plus
petits.
Les statistiques de la mairie de Paris (Mairie de Paris, 2009) montrent que 90% du TMV est
effectué par la route, représentant 20% de tout le trafic urbain. Ces 20% sont à l’origine d’un
tiers des émissions de polluants (tout type de polluants confondus), avec jusqu’à 50% des
émissions de particules, 50% des émissions de NOx et SO2, ainsi que 25% de CO2.
- Type de motorisation des véhicules routiers utilisés pour le TMV : 87% de diesel, 11%
d’essence et seulement 1% de véhicules propres. A noter qu’à l’horizon 2020, tous les
véhicules à motorisation diesel, seront interdits de circulation dans Paris intra-muros.
- Age des véhicules : seulement 40% des véhicules ont été mis en circulation après
2010.
- Nombre de zones de livraison à Paris : à la fin 2015, on comptait environs 9 000 zones
de livraison, dont 7 000 zones dites partagées (i.e. qui peuvent être utilisées par tous
les véhicules, sur des périodes définies de la journée et de la semaine) et 2 000 dites
sanctuarisées (i.e. qui sont exclusivement réservées à la livraison). Ce nombre peut
être considéré comme important, mais rapporté au nombre d’établissements (plus de
500 000 à Paris), il devient très faible (1 zone de livraison pour 55 établissements et
seulement 1 zone sanctuarisée pour 250 établissements).
- Evolution du trafic des véhicules utilitaires et des poids lourds dans Paris intra-muros
et sur le corridor périphérique : comme le montre la Figure I.13, sur les 15 dernières
31
CHAPITRE I : Logistique et transport de marchandises en ville
Figure I.13 : Evolution du trafic des véhicules de TMV à Paris et sa périphérie – moyenne journalière
32
CHAPITRE I : Logistique et transport de marchandises en ville
- 500 km de voies navigables, avec 70 ports urbains sur une surface de 1 000 ha. Ce qui
permet à cette région d’être la plus couverte en France par ce mode de transport et en
2ème position, au niveau européen.
- Réseau ferroviaire dédié au fret composé de 111 stations et connecté à 640
installations fret.
- 1 million m2 d’entrepôts de stockage, avec 6 plateformes multimodales et 6 terminaux
de containers.
La Figure I.14 montre la distribution des différentes infrastructures dédiées au fret dans toute
la région. Nous remarquons que la majeure partie de ces infrastructures se situe au nord et au
sud de la région parisienne, à plusieurs dizaines de kilomètres de la ville de Paris. A noter que
les stations ferroviaires de fret et les ports situés dans les principaux centres logistiques sont
principalement utilisés pour interconnecter la région au reste de la France (ou de l’Europe), et
non pas à desservir la région elle-même.
Figure I.14 : Infrastructures et sites logistiques dans la région Ile-de-France (Paris Region, 2017)
Il faut noter que la situation actuelle des sites logistiques dans cette région, est le fruit d’une
tendance d’externalisation des activités logistiques de la ville de Paris vers la banlieue, qui
s’est poursuivie durant plusieurs décennies (Dablanc & Andriankaja, 2011). Cette tendance a
été portée par la recherche permanente des acteurs économiques, de la réduction des coûts
logistiques. En effet, en banlieue il est possible d’avoir des surfaces plus importantes à des
coûts bien plus bas qu’au centre de Paris. Cependant, cette optimisation économique a eu un
33
CHAPITRE I : Logistique et transport de marchandises en ville
impact non négligeable sur les coûts environnementaux de cette activité (i.e. dégradation de la
qualité environnementale), ainsi que sur la croissance du flux de véhicules utilitaires en ville
et dans la périphérie.
4.2.1. La concertation
Il s’agit de réunir tous les acteurs du TMV au sein du même groupe de travail, dans le but
d’identifier les attentes et besoins de chaque acteur et de trouver les meilleurs compromis
entre ces derniers. Les acteurs principaux concernés par les concertations étaient : les conseils
généraux d’Ile-de-France, la ville de Paris, la préfecture de police, les représentants
professionnels du transport et les associations professionnelles. Les habitants de la ville ont
eux aussi été consultés.
Les résultats des concertations ont été regroupés sous forme de 7 orientations, dont la ville de
Paris devait tenir compte lors de l’établissement du règlement parisien des marchandises. Ces
orientations concernent :
34
CHAPITRE I : Logistique et transport de marchandises en ville
6) La mise en place d’une charte de bonnes pratiques entre la ville de Paris et les acteurs
du secteur des marchandises.
7) L’introduction du principe environnemental dans le règlement marchandises en
vigueur.
Ces différentes étapes ont donné lieu en 2007, à une nouvelle règlementation sur le TMV dans
la ville de Paris. Cependant, il y a toujours un certain nombre de divergence et un manque
d’harmonisation dans les règlementations des différentes localités de la région Ile-de-France,
malgré les efforts récents. De plus, bien que l’implantation et la gestion des zones de livraison
aient représenté un chapitre important, lors de l’établissement de la réglementation
« marchandises », le manque de contrôle sur le terrain de l’utilisation de ces zones, génère
toujours un nombre important de perturbations (tels que le stationnement illicite et le
dépassement du temps d’arrêt règlementaire).
35
CHAPITRE I : Logistique et transport de marchandises en ville
Plusieurs expérimentations ont abouti à des solutions innovantes à l’échelle de la région Ile-
de-France, mais aussi, ont permis l’implantation de ces mêmes solutions dans d’autres villes.
A titre d’exemple :
5. Conclusion
Dans ce chapitre, nous avons défini les différents concepts liés au transport de marchandises,
tout en situant cette activité dans la chaîne logistique. Nous avons étudié la situation actuelle
et les pratiques en vigueur, en matière de transport de marchandises en général, ce qui nous
permet d’en tirer les conclusions suivantes :
36
CHAPITRE I : Logistique et transport de marchandises en ville
Aussi, nous avons montré que la problématique de TMV est plus complexe que celle du
transport de marchandises en général. Cette complexité est principalement liée aux :
L’étude des pratiques actuelles en matière de TMV, quasi-dominé par la route, nous a permis
de mettre en évidence les multiples carences de ce mode de transport et ses limites, au vu des
modèles de développement durables, souhaités pour les villes. En effet, l’explosion de la
demande a eu des répercussions très négatives pour la ville, qu’on peut résumer comme suit :
D’où, la nécessité de faire évoluer le TMV. En premier, au niveau des modes de transport
utilisés, à travers la promotion des modes de transport alternatifs (ferroviaire, fluviale,
triporteurs et deux roues). Deuxièmement, au niveau organisationnel, pour réduire le nombre
de kilomètres assurés par le mode routier. L’une des idées étant le développement de
l’intermodalité.
37
CHAPITRE I : Logistique et transport de marchandises en ville
L’étude de la région Ile-de-France, nous a permis de montrer que les phénomènes observés en
général, dans les pratiques de TMV, sont accentués dans cette région. Cette accentuation est
due à l’importance de la densité de population dans cette région, ainsi qu’à la densité de
l’activité économique.
38
Chapitre II
Transport de marchandises en ville par
rail – mixité fret / voyageurs
39
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
1. Introduction
Dans le premier chapitre, nous avons montré que la pratique actuelle du TMV, quasi-
exclusivement assurée par le mode de transport routier, n’est pas à la hauteur des challenges
actuels pour accompagner le développement des villes. Plus encore, la saturation du réseau
routier urbain (en particulier, aux heures de pointes) et les différentes nuisances de ce mode,
poussent les autorités à durcir la règlementation et à encourager le report vers d’autres modes
de transport, plus durables. Ainsi et pour remédier à ce problème, plusieurs alternatives
modales et organisationnelles sont proposées.
Parmi les alternatives proposées, nous nous intéressons dans ce travail, à l’utilisation du
réseau ferroviaire urbain (initialement dédié au transport de voyageurs), pour le transport des
marchandises d’un point à un autre, en ville et ses environs. L’objectif de cette alternative,
n’est pas de se substituer au mode de transport routier, mais de transférer une partie des flux
de marchandises en ville, vers le mode ferroviaire. Ceci, permettra de réduire la pression sur
la route, ainsi que de rentabiliser la capacité non utilisée des trains de passagers et les
infrastructures y afférentes (en particulier, durant les heures creuses).
La première partie de ce chapitre, sera dédiée à une revue de littérature sur le transport mixte
fret / voyageurs, ainsi qu’à l’étude de différentes expérimentations et solutions déployées sur
le terrain (passées ou en cours d’exploitation). Dans la seconde partie, nous allons définir
l’activité de transport mixte fret / voyageurs, puis, nous allons identifier et étudier les
différents cas de mixité fret / voyageurs possibles, sur le réseau de transport ferroviaire
urbain. Enfin, dans la dernière partie, nous allons identifier et définir les problématiques
décisionnelles principales, liées à la mise en œuvre d’une telle solution de transport.
2. Pratiques et état des lieux : fret urbain par rail et mixité fret /
voyageurs
Le transport de marchandises par rail en ville existe déjà, mais sous une forme différente de
celle étudiée dans ce présent travail, et à une échelle assez réduite. En effet, il est mis en
œuvre sporadiquement et dans certaines circonstances exceptionnelles, mais aussi, à travers
certaines expériences, dont le but est d’évaluer son potentiel. D’un autre côté, le déploiement
d’une telle solution de transport en milieu urbain et périurbain, figure parmi les alternatives
suggérées par certains chercheurs.
Dans cette section, nous allons commencer par présenter les premières alternatives pour le fret
urbain, conçues sur la base du modèle de transport de voyageurs par rail, principalement, en
utilisant l’espace souterrain de la ville et ses environs. Puis, nous étudierons les quelques
40
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Voyageurs Marchandises
Voyagent de manière active Transportées de manière passive
Embarquent, débarquent et effectuent les Doivent être chargées, déchargées et
correspondances sans assistance transférées
Traitent l’information et l’utilisent sans L’information doit être traitée par des
assistance gestionnaires logistiques
Choisissent entre les différents moyens de
Les gestionnaires logistiques doivent choisir
transport sans assistance (parfois de manière
un mode de transport de manière rationnelle
irrationnelle)
L’idée d’utiliser l’espace souterrain urbain, pour un transport exclusivement dédié aux
marchandises (inspirée du modèle de transport de voyageurs), n’est pas récente. En 1991, une
équipe de chercheurs Japonais a proposé de développer un tel système, pour remédier à
plusieurs problèmes environnementaux, liés à la croissance continue de la fréquentation en
ville des véhicules de marchandises (Kashima, et al., 1993). Au Pays-Bas, des études de
faisabilité pour une telle solution de transport, ont été menées à partir de 1995, à travers le
projet « OLS-ASH ». Ces études ont été suivies par le déploiement d’expérimentations en
laboratoire, avec le développement de solutions technologiques avec des véhicules autonomes
« UFT » (Pielage, 2001). Les conclusions auxquelles ont abouti ces études et
expérimentations au sujet de cette solution de transport, peuvent être énoncées comme suit :
41
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Une autre étude expérimentale « CargoCap » (CargoCap, 2002), propose d’utiliser l’espace
souterrain en ville et ses environs (jusqu’à 150 km de distance), pour déployer un système de
transport automatisé par pipeline. Le système développé doit fonctionner 24 heures par jour,
toute la semaine (i.e. sans aucune interruption). Le caractère exclusif au transport de
marchandises, associé à l’utilisation d’une solution technologique automatisée, confère à cette
solution de transport, une grande fiabilité et la rend insensible aux influences des autres
modes de transport (tel que : routes saturées et circulation). Le véhicule de transport autonome
est dimensionné de sorte à pouvoir charger deux europalettes avec des dimensions :
800*1200*1050 mm (à savoir, l’europalette est très utilisé dans le transport de marchandises
et a des dimensions standardisées). Ainsi, chaque commande doit être contenue dans des
emballages respectant les dimensions de l’europalette, et est affecté à autant de véhicules
autonomes que nécessaire, transportant chacun deux palettes. La Figure II.1 montre les
éléments composant ce système et son implantation.
Figure II.1 : Modèles 3D de la solution CargoCap : a- mise à l'echelle de la solution par rapport à un metro, b-
système de chargement / déchargement automatique, c- insertion de la solution dans l'espace souterrain (CargoCap,
2002)
Selon (Delmastro, et al., 2016), bien que les études précédentes montrent un fort potentiel, à
l’utilisation d’une solution de transport souterrain dédiée uniquement aux marchandises, les
auteurs estiment qu’elle reste irréaliste, en raison du coût élevé de l’investissement nécessaire
et de la non compatibilité d’une partie des marchandises (particulièrement les marchandises
en vrac), avec le transport souterrain.
Cependant, un pas vers la mise en œuvre d’une telle solution a été franchi en Suisse. En effet,
une solution de transport souterrain « Cargo Sous Terrain » (CST) est en cours de
développement pour une mise en exploitation à l’horizon 2030 (CST, 2016). En effet, l’étude
de faisabilité finalisée en 2015, a montré que cette solution de transport est faisable, aussi bien
42
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Figure II.2 : Modélisations de Cargo Sous Terrain : a- principe de fonctionnement, b- coupes latérale et horizontale,
c- modèle 3D des véhicules en développement (CST, 2016)
43
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
la capacité résiduelle dans les bus pour voyageurs, pour transporter les marchandises. Quant
au déchargement, il est effectué dans les arrêts de bus habituels de la ligne, à partir desquelles,
un système de distribution capillaire est déployé, pour l’acheminement des livraisons vers leur
destination finale. Les principales conclusions auxquelles est parvenue cette étude sont :
44
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
- Les gains relatifs aux trois conclusions précédentes ont été évalués à 70€/t (i.e. le
surcoût indirect de l’utilisation du mode routier par rapport au mode ferroviaire).
- L’alternative ferroviaire présente un bénéfice global, relativement au mode routier,
équivalent à 10€/t.
Figure II.3 : Train de passagers modifié pour le transport de marchandises (Nuzzolo, et al., 2008)
Dans (Browne, et al., 2014), partant du constat de la prédominance du fret routier, en milieu
urbain, et de la nécessité de développer des alternatives plus durables, tel que le ferroviaire,
les auteurs ont adopté une démarche en deux phases :
- Maximiser son utilisation pour les marchandises en vrac, tel que : les matériaux de
construction, qui se déplacent vers les villes, et les déchets domestiques et industriels,
qui se déplacent vers l'extérieur des villes.
- Développer des solutions logistiques « vertes », pour le transport des biens de
consommation dans les zones urbaines, ceci en combinant l'utilisation du rail pour le
transport à partir des centres de distribution, vers des emplacements dans la zone
urbaine, et de là, l'utilisation de véhicules routiers « propres » (ex : véhicules
électriques), pour l’acheminement des marchandises jusqu’à leur destination finale.
- Dans les cas où le rail ne peut être utilisé directement, utiliser les infrastructures
ferroviaires existantes, ou de nouvelles (telles que les gares), comme des centres pour
la distribution des marchandises en ville.
45
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
- Transfert modal vers le fret ferroviaire limité. Ceci revient à des coûts excessifs
nécessaires à l’extension du réseau.
- Espace et infrastructures limités en ville pour le déploiement d'une telle solution (dans
le cas de nouvelles installations).
- Coûts élevés de l'infrastructure ferroviaire et des systèmes connexes.
- Flexibilité physique limitée.
- Croissance du besoin de transport de voyageurs par rail en ville, ce qui limite la
croissance de la part dédiée au fret.
Pour autant, ces barrières ne sont pas insurmontables, mais requiert des études
supplémentaires et surtout, une volonté des transporteurs et des autorités, pour investir dans le
développement du fret ferroviaire.
Dans (De Langhe, 2014), l’auteur réalise un travail d’investigation similaire au précèdent, en
identifiant aussi bien les facteurs de réussite d’un fret urbain par rail, que les obstacles à son
développement. Il ajoute quelques analyses pertinentes, telles que :
- Les caractéristiques de la ville pour une implantation d’une solution de fret ferroviaire,
sont un paramètre crucial, qui requiert une étude approfondie. En particulier, il faut
adapter le projet à la taille de la ville, la densité de population, la structure économique
de la ville, les infrastructures de transport existantes, etc.
- Lorsqu’on parle de transport de marchandises par rail en milieu urbain, le type de
marchandises à transporter peut varier d’une ville à l’autre et d’un projet à l’autre :
marchandises en vrac, livraison pour particuliers, livraison pour la grande distribution,
déchets, etc.
Dans (Diziain, et al., 2014), les auteurs ont étudié le transport urbain multimodal de
marchandises, en France et au Japon. Ils ont évalué le potentiel du report modal du mode
routier vers le ferroviaire et le fluvial. Ils ont conclu qu’aussi bien le mode ferroviaire, que le
fluvial, sont trop coûteux à mettre en œuvre, et qu’ils ne peuvent être rentables, que si
plusieurs conditions sont satisfaites (ex : réseau routier congestionné, existence
d’infrastructures multimodales, disponibilité d’un terminal ferroviaire en milieu urbain, etc.).
Les études qu’on vient de présenter concernent l’utilisation des trains de fret (au sens large),
pour le transport de marchandises en milieu urbain. Principalement, ce transport concerne
l’introduction des marchandises à partir de l’extérieur de la ville, ainsi que le transport des
produits et déchets de la ville vers l’extérieur (et non pas les livraisons en ville, d’un point à
un autre). Les études suivantes concernent, quant à elles, les livraisons de marchandises en
ville, moyennant un moyen de transport ferroviaire urbain (ex : métro ou tramway).
Dans (Motraghi & Marinov, 2012) et (Dampier & Marinov, 2015), les auteurs ont étudié la
faisabilité d’effectuer le transport de marchandises, en utilisant le métro à Newcastle upon
Tyne. Les deux études ont utilisé la modélisation pour évaluer le potentiel d’intégration du
flux de marchandises. La première, à travers un modèle à évènements discret, où l’intégration
de métros supplémentaires sur la ligne de l’aéroport a été évaluée, suivant plusieurs scénarios.
La seconde, à travers un logiciel de transport spécifique, qui permet d’évaluer l’impact
46
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
monétaire des accidents de la route et, ainsi, évaluer les bénéfices du report modal vers le
métro. Ces deux études ont conclu à la viabilité d’un tel système de transport de
marchandises. En effet, au-delà des avantages environnementaux et réduction des accidents et
perturbations sur la voirie urbaine, ce système permet :
Cependant, il est suggéré d’effectuer des études supplémentaires avec plus de données avant
de passer à une éventuelle implémentation.
Ces études ont été complétées par deux autres, (Brice, et al., 2015) et (Reece & Marinov,
2015), à travers un cas d’étude. Ce cas d’étude traite l’éventuelle utilisation du système
précèdent, pour le transfert de bagages à destination de l’aéroport de Newcastle, à partir du
centre-ville. Ces deux études ont montré la faisabilité technique de cette solution. Cependant,
elles ont mis en évidence l’important surcoût généré, comparativement au service actuel de
transfert de bagages par la route.
Dans (Gonzalez-Feliu, 2014), l’auteur développe une démarche sur la base de la méthode
d’analyse coût-bénéfice, pour évaluer la pertinence du déploiement d’une solution de
transport de marchandises, en utilisant le réseau ferroviaire urbain. Il a appliqué cette
démarche sur un cas d’étude, relatif à une ligne de tramway dans la ville de Paris. Les
conclusions de cette étude, présentent le bénéfice potentiel d’une telle solution de transport, et
pose les conditions nécessaires pour réaliser ce bénéfice. Ces conditions étant une forte
demande et l’accompagnement des autorités.
47
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
mise en œuvre d'un cas de mixité. Les modes de transport considérés étant : métro,
tramway et bus.
Les possibilités de mutualisation sont multiples. Elles ont été classées en trois niveaux
de mise en œuvre (ces niveaux permettent de mesurer l'étendu financier et temporel
quant à l'éventuelle mise en œuvre d'une telle solution) :
En utilisant l'existant tel qu'il est : mise en œuvre rapide, avec un coût
minimum. Cependant, cela réduit le champ des possibilités, étant donné que
l’activité de transport des marchandises, doit se conformer aux contraintes déjà
présentes, dans un système conçu exclusivement pour le flux passagers.
En utilisant l'existant moyennant quelques modifications : mise en œuvre plus
ou moins rapide, le temps d’effectuer un certain nombre de travaux et
d’aménagement, pour s’adapter aux contraintes du flux de marchandises. Ce
cas intermédiaire présente un équilibre intéressant en termes de : temps de mise
en œuvre, coûts d’implémentation et contraintes liées aux voyageurs.
En innovant : mise en œuvre lente et coûts élevés. Cependant, la solution
développée, considérant le flux des marchandises dès la phase de conception,
permettra de générer un maximum de profit, grâce à l’étendue des possibilités,
qui sera laissé au transport de marchandises.
Toutefois, une mise en œuvre d’une telle mixité, sur le réseau de transport de
voyageurs, doit faire face à plusieurs contraintes, qui ont été classées en quatre
catégories comme suit :
48
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Bien que contrastées, les conclusions de cette étude sont prometteuses, et peuvent être
résumées comme suit :
Cette étude s’est basée sur l’analyse : 1- du réseau de la RATP, 2- des types de
produits pouvant prétendre au transport sur ce réseau, 3- leur conditionnement, 4- les
49
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
La définition des critères pertinents, pour savoir si un espace RATP est apte à
accueillir des marchandises. Ceci à travers la définition d’une grille d’analyse
multicritère applicable à tous les espaces RATP (ce qui est appelé « Offre »).
La grille multicritère est composée de :
1- Généralités sur les espaces (stations et gares), telles que la localisation, les
numéros des lignes les desservant, la complexité de l’espace, etc.
2- Accessibilité de l’espace, telle que les dénivelés, le nombre de parcours, la
présence d’escaliers, etc.
3- La surface des espaces et leurs dimensions, telle que pour les quais,
couloirs, entrées / sorties, salles diverses, etc.
4- Les équipements de l’opérateur, tels que les guichets, les comptoirs
d’information, etc.
5- Les paramètres relatifs aux flux de voyageurs, tels que le nombre total
annuel d’entrants, le nombre journalier en hiver, l’heure la plus chargée du
matin, le nombre de voyageurs correspondants à l’heure la plus chargée du
matin, l’heure la plus chargée du soir, le nombre de voyageurs
correspondant à l’heure la plus chargée du soir, etc.
Cette étude a montré qu’il y avait un potentiel d’offre logistique de la part de la RATP,
en exploitant ses différentes ressources, pour répondre à une demande de service
logistique concrète. Toutefois, un ensemble de conditions doivent être préalablement
satisfaites (telles que celles énoncées par l’étude 1). Plusieurs services logistiques ont
été étudiés et à l’heure des conclusions de cette étude, seul le service de déploiement
de consignes des colis pour les particuliers, a été retenu. Il est spécifié que les autres
services requièrent des études complémentaires.
- Etude 3 : Transport collectif urbain : surcoûts consécutifs à l’ajout logistique dans un
projet (PREDIT GO n°6, 2012)
50
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
En raison des données stratégiques traitées par cette étude, il n’y a pas de version
publique de cette dernière. Cependant, sa fiche récapitulative, nous permet de saisir
ses grandes orientations.
D’abord, cette étude s’est articulée autour de quatre points majeurs, à savoir :
Pour les besoins de l’étude, il a été considéré la mixité fret / voyageurs dans un métro
automatique (type ligne 1 ou 14 du métro parisien). Aussi, les marchandises sont
conteneurisées dans des unités logistiques standardisées. Ainsi, trois modes
d'exploitation ont été proposées et étudiées (3 scénarios) :
Une évaluation du potentiel commercial, à l’image de celui réalisé par l’étude 2, mais
pour des zones urbaines plus étendues, a été réalisée. Elle a permis d’identifier les
secteurs de la grande distribution et de la messagerie, comme porteurs d’un grand
potentiel, en termes d’utilisation de la solution logistique étudiée. En effet, ces
secteurs répondent positivement à trois groupes de conditions, identifiées
préalablement comme nécessaires à toute mise en œuvre, d’une solution logistique
basée sur la mixité fret / voyageurs :
51
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Enfin, l’étude a été conclue par une application des concepts développés, à une partie
d’une ligne du futur « Grand Paris Express ». Une étude de la demande potentielle sur
le tracé de la ligne a été réalisée, puis, un placement des stations de marchandises a été
proposé, ainsi que des aménagements spécifiques sur la ligne. Cependant, l’évaluation
économique du surcoût d’intégration de la composante « marchandises », à un projet
de transport de voyageurs, n’a pas abouti. En effet, il est conclu que plusieurs outils
doivent être développés (tels que des outils mathématiques), dans le but d’effectuer
correctement l’évaluation des coûts de déplacement des marchandises et les surcoûts
d’intégration du nouveau service. Ceci requiert plus de ressources techniques et
temporelles, dont cette étude ne disposait pas.
Nous retenons de cette étude, la faisabilité technique dans différentes situations, ainsi
que le potentiel commercial et les secteurs cibles les plus prometteurs.
52
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
53
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Cette solution permet la collecte de 14 000 tonnes de déchets par an, avec 11 métros
dédiés, couvrant 359 stations et totalisant 567 arrêts.
Dans un premier temps, cela nous permet de proposer une définition du transport mixte fret /
voyageurs, sur le réseau ferroviaire de transport de voyageurs en milieu urbain, comme suit :
- Quand on parle de milieu urbain, on entend la ville à l’intérieur de ses limites, ainsi
que sa banlieue et les agglomérations limitrophes. A titre d’exemple, la région Ile-de-
France, constituée de la commune de Paris au centre et de plusieurs agglomérations
limitrophes.
- Quand on parle de réseau ferroviaire en milieu urbain, on entend les différents modes
de transport sur rail. En effet, dans une ville, on peut trouver le métro, le tramway,
ainsi que le train de banlieue qui traverse la ville et dessert les agglomérations
limitrophes.
54
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Réseau de métro*:
• 220 km
• 302 stations
Réseau de Tramway*:
• 104 km
• 181 stations
Réseau de RER**:
• 587 km
• 257 gares
Ce réseau est appelé à s’étendre dans les prochaines années, dans le cadre du projet « Grand
Paris Express ». Il s’agit de 4 nouvelles lignes de métro automatisées, sur une distance de 200
km et comportant 68 nouvelles stations. Ces lignes formeront une boucle reliant les
différentes régions d’Ile-de-France autour de la ville de Paris, avec des interconnexions aux
anciennes lignes. Ce projet présente une opportunité supplémentaire, d’intégration de
nouveaux services dès la réalisation, telle que la mixité fret / voyageurs.
55
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Dans le point suivant, nous allons identifier les différentes composantes d’un réseau
ferroviaire en milieu urbain.
- Les stations (ou gares) : la taille des stations varie selon le mode utilisé. Ainsi, les
stations de tramway sont les plus petites, puis celles dédiées au métro, enfin, celles
pour le train de banlieue. Trois formes de stations sont à recenser, comme le montre la
Figure II.6, à savoir :
Station aérienne
Station sous-terraine
Les stations souterraines : comme le métro est généralement déployé dans des
villes existantes très denses, ce type de stations est le plus fréquent, puisqu’il
n’interfère pas avec l’existant. Cependant, c’est le type de station dont le coût
de réalisation est le plus élevé. En ce qui concerne le train de banlieue, lorsqu’il
56
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
traverse les centres urbains (telle que la ville de Paris), les stations sont de ce
type. Enfin, on trouve très rarement des stations de ce type, tout au long des
réseaux de tramway (en raison de l’incompatibilité du coût de ce genre
d’infrastructure avec le moyen de transport qu’est le tramway).
- Les voies ferrées : ce qui donne son nom à ces trois modes de transport, c’est leur
spécificité à rouler sur une voie ferrée dédiée (dites en site propre, parce que
exclusivement dédiée au mode de transport en question). Même schéma que pour les
stations, on compte trois situations possibles (illustrées par la Figure II.7) :
Voies aériennes : ce type de voies est principalement utilisé pour les lignes de
métro (tel que la ligne 6 à Paris, une partie du réseau de New York et de
Chicago…). Plus rarement, on le trouve ponctuellement dans les deux autres
modes, pour passer un obstacle (tel qu’un fleuve ou une grande rue).
Voies au niveau de la chaussée : le réseau de tramway est quasi-exclusivement
constitué de ce type de voies. Les moins coûteuses à mettre en œuvre, elles ne
nécessitent aucune infrastructure supplémentaire (hormis les voies elles-
mêmes). Dans plusieurs villes, où la voierie n’est pas assez large dans certaines
parties de la ville, le tramway partage son chemin avec les voitures. Pour le
métro, ce type de voies est très peu utilisé. Quant au train de banlieue, en
général les voies entres les agglomérations voisines, se trouvent au niveau de la
chaussée. Cependant, ces dernières se trouvent derrière une clôture.
Voies aériennes
Voies sous-terraines
Figure II.7 : Illustration des implantations possibles de voies ferrées
Voies souterraines : c’est le type de voies qui coûte le plus cher, en raison de la
nécessité de creuser un tunnel. Ces voies sont principalement déployées dans
les centres urbains (très) denses pour les réseaux de métro, ainsi que, pour
permettre aux trains de banlieue de traverser la ville. Pour le tramway, on
57
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
- Le matériel roulant ferroviaire (trains) : il s’agit de tous les véhicules qui roulent sur
une voie ferrée. Dans le cas des trois modes de transport précédents, on compte trois
types de matériels roulants (d’ailleurs, qui partagent le nom du mode en question).
Tous ces véhicules roulent à l’énergie électrique. La Figure II.8 montre un exemple de
chacun de ces trois modes comme suit :
Tramway : les véhicules du tramway sont les plus petits parmi les trois modes
et ceux qui ont la plus faible vitesse de circulation. Les différentes
caractéristiques techniques de ces véhicules, varient d’une ville à une autre et
selon la génération du véhicule. Toutefois, pour apprécier ces caractéristiques,
un exemple à partir d’une ligne de tramway en région Ile-de-France, comme
suit : longueur : 30 m, largeur : 2,3 m, hauteur : 3,36 m, capacité maximum :
180 passagers et vitesse variant de 18 à 25 km/h.
Métro : les véhicules de métro sont assez variables, d’une ville à l’autre, voire
d’une ligne à l’autre. On peut trouver des véhicules automatisés, semi-
automatisés ou avec conducteur. Leur longueur varie en nombre de wagons.
Puis, on trouve des métros sur pneumatiques ou bien, plus classiques, avec des
roues en aciers. A titre d’exemple, l’un des métros parisiens est caractérisé
par : longueur : 90 m, largeur : 2,5 m, hauteur : 3,5 m, hauteur intérieur : 2,14
m, hauteur porte : 1,9 m, largeur porte : 1,65 m, capacité maximum : 976
passagers et vitesse maximum : 80 km/h.
Métro
RER
Tramway
Figure II.8 : Exemple de matériels roulant de chacun des modes de transport en région Ile-de-France
58
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Train de banlieue : les trains de banlieue sont assez grands, avec de grandes
capacités de charge. Ils varient suivant les destinations desservies. Ainsi, on
trouve des trains à double étage, des trains courts et des trains longs. Aussi,
suivant les lignes empruntées, la vitesse peut varier énormément. A titre
d’exemple, l’un des trains utilisé pour la ligne A du RER (en Ile-de-France) :
longueur : 112 m, largeur : 2,9 m, hauteur : 4,32 m, largeur porte : 2 m,
capacité maximum : 2 600 passagers et vitesse maximum : 120 km/h.
Pour un même mode de transport, ces éléments peuvent varier d’une ligne à une autre et d’une
ville à une autre. Ceci implique leur considération dans le cadre d’études techniques de
faisabilité, contextualisée au cas spécifique étudié. Dans le cadre de ce travail, l’objectif est de
proposer une démarche méthodologique, pour l’intégration du flux de marchandises, à un
mode de transport ferroviaire en milieu urbain, dans la perspective de proposer un nouveau
service de transport de marchandises. Ceci implique de se limiter aux éléments communs de
tels modes de transport.
Comme mentionné précédemment, la mixité fret / voyageurs signifie le partage d’au moins
l’une des trois composantes précédentes. Une quatrième composante est à prendre en
considération lors de l’étude de la mixité, il s’agit du moment auquel les marchandises sont
transportées (i.e. les heures d’exploitation du service de fret). Nous avons étudié les
différentes possibilités et avons identifié plusieurs cas de mixité possibles (en comptant les
différentes combinaisons, il y a 19 cas de mixité possibles). La Figure II.9 montre ces
différents cas.
59
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Transport de
marchandises sur le Partage de toutes les
réseau de transport de stations
Circulation des deux flux Partage des voies
voyageurs (une ligne)
durant les heures de ferrées de toute la
transport de voyageurs
Stations partagées
service ligne
Horaires d’exploitation partagés Voies ferrées partagées Partage de quelques
Moyens de transport disjoints: stations
• Moyen dédié fret
Stations partagées (Partiel)
• Moyen dédié voyageurs
Trains dédiés Partage de (quelques)
stations + stations
dédiées
Stations dédiées
Stations partagées
Figure II.9 : Les niveaux de mutualisation des différentes ressources du système de transport ferroviaire pour une mixité fret / voyageurs
60
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Ce cas de mixité est le plus complexe ayant été identifié, en raison du nombre de contraintes
liées au service de transport de voyageurs que le service de transport de marchandises doit
respecter. On retrouve aussi bien des contraintes techniques, qu’organisationnelles et
juridiques, telles que :
A savoir, une limitation juridique très forte existe actuellement (comme cité précédemment),
qui concerne l’impossibilité pour le transporteur public de voyageurs d’effectuer le service de
transport de marchandises et de faire circuler les deux flux, fret et voyageurs simultanément.
Cependant, la première étude du PREDIT avait conclu que bien que cette limitation soit très
forte, il est possible de l’abroger, si toutefois, d’autres conditions sont réunies.
Pour ce cas de mixité, nous proposons d’avoir des trains avec des wagons dédiés aux
marchandises, séparés des wagons dédiés aux voyageurs (i.e. une séparation physique entre
les deux). Le(s) wagon(s) dédiés aux marchandises est (seront) placé(s) à l’arrière du train).
Aussi, une petite zone de stockage clôturée sur le quai de la station, pour y mettre les
marchandises à transporter temporairement. La Figure II.10 montre une modélisation du
système proposé, avec un wagon dédié aux marchandises (où les sièges et équipements pour
les voyageurs sont supprimés) et les autres wagons standard dédiés aux voyageurs. La Figure
II.10 montre aussi un petit espace pour 6 colis sur le quai.
- Toutes les stations de voyageurs seront équipées et utilisées pour assurer le service de
transport de marchandises.
- Quelques stations seulement, choisies de manière stratégique (ce dernier cas est le plus
probable).
61
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Une 17ème possibilité de mixité est donnée par le cas où des espaces liés à l’activité de
transport de marchandises, sont implantées dans les stations ferroviaires.
Partage des horaires d’exploitation : il s’agit de la circulation des trains dédiés aux
voyageurs et des trains dédiés aux marchandises, au même moment. Nous avons identifié 4
cas de partage des stations, pour chacun 2 possibilités de partage des voies ferrées :
62
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
63
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
1- Pas de limite de temps lors des arrêts dans les stations, lié aux trains de
voyageurs.
2- Pas de limite dans la vitesse de circulation, lié à la synchronisation avec les
trains de voyageurs.
3- Pas de sécurité supplémentaire, pour considérer les risques d’accidents avec
les voyageurs sur les quais.
64
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
65
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
- Chaque ligne du réseau est physiquement indépendante des autres lignes (i.e. chaque
ligne dispose de ses propres voies ferrées, ses propres trains (métros, tramways ou
trains de banlieue), ainsi que ses propres quais dans les stations.
- Les lignes des différents modes de transport (i.e. les trois cités précédemment) sont
interconnectés entre eux de manière indirect, à travers des stations (ou gares)
d’interconnexions, avec rupture de charge (i.e. dans la station d’interconnexion, il est
nécessaire de quitter le moyen de transport de la première ligne, pour rejoindre celui
de la seconde, en traversant un certain parcours entre les quais des deux lignes, à
travers la station).
- En utilisant une seule ligne du réseau : les marchandises sont chargées dans un train de
la ligne, transportées entre les deux stations (départ-arrivée) et déchargées dans leur
station d’arrivée.
- En utilisant deux lignes (ou plus) du réseau, même mode de transport : les
marchandises sont chargées dans un train de la ligne, transportées entre la station
départ et la station de correspondance, déchargées et transférées à l’aide d’un moyen
de manutention vers le quai du train de la seconde ligne de correspondance, puis
chargées dans un train de cette seconde ligne, sachant que ce train, est du même mode
que le premier : métro, tramway ou train de banlieue. Enfin, elles sont transportées
entre les deux stations (correspondance-destination finale) et déchargées dans leur
« station arrivée ».
- En utilisant deux lignes (ou plus) du réseau, plusieurs modes de transport : les
marchandises sont chargées dans un train de la ligne, transportées entre la station de
départ et la station de correspondance, déchargées et transférées à l’aide d’un moyen
de manutention vers le quai du train de la seconde ligne de correspondance, puis
chargées dans un train de cette seconde ligne (sachant que, ce train est d’un autre
mode que le premier). Enfin, elles sont transportées entre les deux stations
(correspondance-destination finale) et déchargées dans leur station d’arrivée.
Nous remarquons que les deux derniers schémas, suivent un même processus dans le cadre du
transport de marchandises, à partir de la station de départ, jusqu’à sa livraison vers sa
destination finale.
66
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
voie retour
voie aller
Mode de transport 2 :
ligne 1
voie retour
voie retour
voie aller
Mode de transport 1 :
ligne 1
Mode de transport 2 :
voie aller
ligne 2
voie aller
Mode de transport 1 :
ligne 2
Ainsi, dans ce travail, nous allons nous limiter au cas composé d’une voie ferrée, tout au long
de laquelle se trouvent plusieurs stations, où il est possible de procéder au chargement /
déchargement des marchandises. Cette ligne est desservie par plusieurs trains de manière
régulière.
67
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Station de départ
Une partie du réseau ferroviaire de la région
Ile-de-France (RATP)
Processus de transport de
marchandises à partir de
“Creil” jusqu’à “Haussmann
Saint-Lazare” (≈52km):
1. Chargement des
marchandises dans la
station de départ dans le
RER D.
2. Transport jusqu’à la
station intermédiaire
“Gare du Nord” (≈50km).
3. Déchargement des
marchandises et leur
transfert à l’aide d’un
moyen de manutention
jusqu’au quai du RER E et Station
leur chargement. intermédiaire
4. Transport jusqu’à la
station d’ arrivée
“Haussmann Saint-Lazare”
(≈2km).
Station
d’arrivée
Figure II.12 : Illustration du processus de transport de marchandises à travers le réseau ferroviaire de la région Ile-
de-France
Ce schéma nous montre les différentes problématiques posées par la mise en exploitation
d’une telle solution de transport. En effet, nous avons identifié 9 problématiques aux trois
niveaux décisionnels (i.e. stratégique, tactique et opérationnel). A savoir, cette liste n’est pas
exhaustive, il peut être nécessaire de résoudre des problématiques supplémentaires, lors de
l’étude de cas particuliers, voire, pour certains autres cas de mixité identifiés dans la Figure
II.9.
68
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
& Problème de bin packing 2D/3D dans les zones de stockage temporaires [Niveau Opérationnel]
manutention
Station d’arrivée
Zone de
stockage Train (métro, tramway ou RER)
manutention
temporaire
Quai
manutention Zone de
Ligne d’interconnexion « rail » stockage
Quai
temporaire
Station de départ
manutention
Figure II.13 : Identification et décomposition des problématiques décisionnelles pour la mixité des flux fret / voyageurs
69
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Que cela soit lors de la réalisation d’une nouvelle station (comme dans le cadre du projet
« Grand Paris Express »), ou dans le cadre du réaménagement d’une station existante (encore
plus intéressant, au vu de l’intérêt fondamental de l’idée de mixité, qui est de pouvoir utiliser
l’existant, dans le but d’intégrer le nouveau service et de générer un bénéfice supplémentaire),
cette problématique est cruciale. La solution optimale est à privilégier, en raison de l’impact
sur le long terme et du manque à gagner correspondant à un éventuel mauvais
dimensionnement. Deux cas se présentent :
70
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
La modification de la configuration des trains, est une tâche qui nécessite un certain temps de
mise en œuvre, ainsi que la mobilisation de ressources humaines et matérielles. D’autre part,
elle est indispensable, pour maximiser la rentabilité dégagée par l’exploitation du service de
transport de marchandises. A ce titre, les trains peuvent être dimensionnés sur la base d’une
demande moyenne sur une période, ou bien, avoir des configurations différentes selon le jour
(jours de semaine ou weekend), ou selon les moments de la journée (matin, après-midi, soir).
A noter que le service de transport de marchandises est proposé durant les heures creuses de
la journée.
Le dimensionnement des trains nécessite aussi, la disponibilité de plusieurs données, dont les
principales étant les prévisions des deux flux fret / voyageurs et les temps d’attente.
- Problème de bin packing 2D/3D dans les zones de stockage temporaire (numéros
4 et 5 sur la Figure II.13)
Ce problème sera fixé (entre 2D et 3D) une fois la solution technique déterminée. En
effet, il s’agit d’abord de déterminer si les contenants des marchandises s’empilent,
puis, s’ils ont la même taille. Quelle que soit la solution adoptée, la problématique est
d’optimiser le rangement des marchandises dans la zone de stockage temporaire, avec
pour objectifs :
71
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Un exemple peut être illustré par la Figure II.14. Supposons que les prochaines
marchandises devant être transportées sont celles en rouge. Dans ce cas, pour y avoir
accès, il est nécessaire de sortir le coli en bleu puis les deux colis rouge, enfin remettre
le coli bleu à l’intérieur de la zone. Si les colis rouges étaient placés à l’avant, on
aurait gagné deux opérations de manutention.
72
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Nous avons évalué les différentes possibilités, dans le cas d’un seul niveau de
chargement des colis (i.e. pas d’empilement, soit le cas de la 2D). Ainsi, nous avons
identifié trois possibilités, comme le montre la Figure II.15.
Pour la première possibilité, le colis bleu est déchargé sans avoir à manipuler
les autres colis.
Pour la seconde possibilité, le déchargement du colis bleu requiert la
manipulation du colis qui le bloque.
Pour la troisième possibilité, les trois colis en orange peuvent être déchargés
sans manipuler les autres colis. Cependant, l’espace pour la manutention étant
réduit, l’opération doit se faire avec précaution.
73
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
5. Conclusion
Dans ce chapitre, nous avons étudié la mixité fret / voyageurs, sur le réseau ferroviaire de
transport de voyageurs en milieu urbain, en considérant le cas de la région Ile-de-France. Pour
cela, nous nous sommes basé sur les quelques références dans la littérature, qui ont traité des
problématiques du type fret ferroviaire en milieu urbain. Nous citons quelques
expérimentations et exploitations (passées ou en cours) de ce type de solutions de fret urbain,
que nous jugeons proche de notre cas de mixité.
Sous certaines conditions, la faisabilité technico-économique d’une solution de fret urbain par
voie ferroviaire, a été démontrée par plusieurs chercheurs. Cette faisabilité a été expérimentée
sur le terrain, et on compte aujourd’hui quelques exemples en cours d’exploitation. Les
motivations ayant encouragé le développement d’une telle solution, peuvent être décrites à
travers les points suivants :
74
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
Partant de ce constat, dans ce travail, nous proposons une démarche méthodologique, pour
l’intégration du flux de marchandises, à un mode de transport ferroviaire en milieu urbain,
dans la perspective de proposer un nouveau service de transport de marchandises. Le point de
départ étant l’identification de tous les cas possibles de mixité. Ainsi, nous avons identifié 19
cas de mixité, dont deux cas de mixité totale (i.e. toutes les composantes du moyen de
transport ferroviaire sont communes aux deux services, transport de voyageurs et de
marchandises).
On retiendra le cas de mixité le plus complexe, puis nous proposons une approche de
décomposition des problématiques décisionnelles, qui s’y posent. Ainsi, nous avons identifié
9 problématiques interdépendantes aux 3 niveaux décisionnels : 2 au niveau stratégique, 1 au
niveau tactique et 6 au niveau opérationnel. Nous avons montré que la résolution des
problématiques aux niveaux stratégique et tactique, nécessite la disponibilité d’un certain
nombre de données, issues de la future exploitation de cette solution de transport. En effet,
l’étude de la dynamique du système et de la circulation des flux est indispensable, à
l’évaluation de l’offre de transport potentielle. La vitesse à laquelle les marchandises sont
transportées à partir des stations, confrontée à la demande de transport (prévisionnelle),
permet d’évaluer les besoins en espace dans les zones de stockage temporaire. Aussi, elle est
indispensable au dimensionnement de la capacité des trains, en termes de transport de
marchandises.
Quant aux problématiques de bin packing, elles jouent un rôle important dans l’évaluation de
la faisabilité économique de la solution de transport étudiée. En effet, elles permettent
d’évaluer au plus juste les capacités du système, à travers la meilleure utilisation possible des
75
CHAPITRE II : Transport de marchandises en ville par rail – mixité fret / voyageurs
76
Chapitre III
Etat de l’art des méthodes de modélisation
des problèmes de transport
77
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
1. Introduction
Les activités de transport de personnes et de marchandises posent plusieurs problématiques à
tous les niveaux décisionnels (i.e. stratégique, tactique et opérationnel). Que cela soit pour la
construction de réseaux logistiques, leur gestion et leur organisation, l’optimisation de la
circulation des flux (personnes et marchandises), l’optimisation de la circulation des véhicules
(planification des tournées et optimisation des parcours), etc. La majeure partie des
problématiques relatives à la chaîne logistique et au transport, sont des problèmes
d’optimisation combinatoire et sont classées dans la catégorie des problèmes NP-difficile
(Gupta & Könemann, 2011).
Une grande variété de méthodes de résolution exactes et approchées ont été proposées pour
résoudre de tels problèmes. Les solveurs sont souvent en mesure de résoudre les petites et
moyennes instances de manière optimale. Cependant, des modèles riches ou des instances de
grande taille, ne peuvent être résolus de manière optimale, même par les solveurs les plus
performants, dans des délais acceptables. Ainsi, il est nécessaire de développer des
algorithmes pour une résolution approchée (Alumur , et al., 2015).
Ce chapitre est structuré en quatre sections. Il dresse l’état de l’art des méthodes de
modélisation et de résolution, des problèmes proches du FRTSP.
78
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
dispose d’une quantité limitée, 2- chaque client a une demande maximum et 3- chaque
paire « source-destination » a un coût propre, proportionnel à la quantité transportée.
L’objectif dans ce type de problème est de minimiser les coûts de transport/livraison.
Aussi, des plateformes de transbordement peuvent être introduites, pour réduire les
coûts de livraison, à travers la mutualisation du transport sur une certaine portion
commune des parcours.
- Le problème de localisation (Melo, et al., 2009) : ce problème se pose suite aux
problèmes issus de la classe précédente et dans le cas où les installations de production
(sources) n’existent pas encore (donc il s’agit de les réaliser). Le problème décisionnel
en question, est de définir les localisations de leur implantation, dans le but de
minimiser les futurs coûts de transport ou les futures distances à parcourir pour
procéder aux livraisons.
- Le problème de tournée de véhicules (Dantzig & Ramser, 1959) : cette classe de
problèmes classiques en optimisation combinatoire, consiste à déterminer les tournées
optimales de plusieurs véhicules, dans le but de livrer plusieurs clients (à noter que
lorsqu’il n’y a qu’un seul véhicule, ce problème se réduit au problème du voyageur de
commerce). Plusieurs variantes de problèmes de tournées de véhicules existent. A titre
d’exemple : 1- avec fenêtre de temps, où il s’agit de livrer chaque client, sur un
intervalle de temps réduit, 2- avec collecte et livraison, où il s’agit pour les véhicules
de transport, de procéder à des collectes de marchandises, en plus des livraisons, 3-
contraintes de capacité, où les véhicules ont une capacité limité en termes de poids /
volume…
- Le problème d’affectation (Schrijver, 2005) : c’est l’un des premiers problèmes
d’optimisation combinatoire. Il a été étudié par Gaspard Monge dès 1784. Ce
problème est le plus souvent assimilé au problème de transport décrit précédemment.
La description classique de ce problème consiste à déterminer l’affectation de tâches à
des agents, sachant que chaque affectation tâche-agent a un coût, l’objectif étant de
minimiser le coût total des affectations. Appliqué au transport, il s’agira de déterminer
pour chaque commande (produits ou marchandises), le véhicule qui va la transporter.
Une variante plus complexe de ce type de problème existe, il s’agit du problème
d’affectation généralisée, qui introduit des contraintes supplémentaires. Ce problème
sera décrit plus en détail dans la suite de ce chapitre.
- Le problème de bin packing (Johnson & Garey, 1985) : ce problème consiste à
déterminer le meilleur placement d’objets, à l’intérieur de plusieurs contenants, dont
l’objectif est la minimisation du nombre de contenants pour tous les objets. Appliqué
au transport, ce problème se traduit par l’optimisation du chargement de véhicules de
transport (ex : camions) ou containers (à transporter dans des bateaux ou trains), dans
le but de minimiser le nombre de véhicules à utiliser pour transporter tous le produits.
Des contraintes sur le volume du contenant et le poids maximum supporté par les
véhicules sont à considérer.
Au-delà des autres problèmes d’optimisation combinatoire, qui peuvent être rencontrés en
transport, on peut relever dans certains cas, l’interaction de plusieurs de ces problèmes, pour
traduire un problème réel. Par exemple : un problème de tournées de véhicules avec
contraintes de capacité, associé à un problème de bin packing pour optimiser le chargement
79
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
Ce problème peut être formulé par un modèle mathématique de type PLNE. Les notations
utilisées se présentent comme suit :
min 𝑓(𝑥) = ∑𝑀 𝑁
𝑚=1 ∑𝑛=1 𝑐𝑛𝑚 𝑥𝑛𝑚 (3.1)
s.c.
∑𝑁
𝑛=1 𝑏𝑛𝑚 𝑥𝑛𝑚 ≤ 𝑎𝑚 ∀𝑚 ∈ [1, 𝑀] (3.2)
∑𝑀
𝑚=1 𝑥𝑛𝑚 = 1 ∀𝑛 ∈ [1, 𝑁] (3.3)
80
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
L’objectif (3.1) étant la minimisation du coût total de l’affectation des tâches aux machines,
sous l’ensemble des contraintes (3.2) de capacité des machines et l’ensemble des contraintes
(3.3), qui permettent d’assurer l’affectation de chaque tâche à une seule machine. L’ensemble
des contraintes (3.4) représente les contraintes d'intégrité des variables de décision.
? ?
? ?
Figure III.1 : Schématisation du problème d'affectation généralisée dans le cas de marchandises à livrer
Il a été démontré que ce problème est NP-difficile dans (Fisher, et al., 1986). Toutefois,
plusieurs algorithmes de résolution exacte de ce problème, pour des instances de petite et
moyenne taille, ont été proposés par plusieurs auteurs dans : (Ross & Soland, 1975), (Fisher,
et al., 1986), (Guignard & Rosenwein, 1989), (Martello & Toth, 1990), etc. Dans
(Savelsbergh, 1997), l’auteur a proposé un algorithme qui fournit la solution optimale à un
problème composé de 20 machines et 200 tâches en moins de 1 160 secondes dans le pire des
cas. D’autre part, plusieurs heuristiques et métaheuristiques ont été proposées pour résoudre
des instances de taille plus importante. Dans (Amini & Racer, 1994), les auteurs ont proposé
une heuristique à recherche en profondeur variable, inspiré des travaux développés dans (Lin
& Kernighan, 1973) pour le problème du voyageur de commerce. Dans (Cattrysse & Van
Wassenhove, 1992), les auteurs ont formulé le problème sous forme d’un problème de
partitionnement et ont proposé une heuristique basée sur les techniques de génération de
colonnes. Dans (Chu & Beasley, 1997), un algorithme génétique est proposé pour améliorer la
qualité et la faisabilité des solutions explorées. Dans (Laguna, et al., 1995) et (Dıaz &
Fernández, 2001) des algorithmes de recherche tabous ont été proposé, pour la résolution du
problème de base et de certaines de ses extensions. Enfin, une autre classe de métaheuristique
a été utilisée par plusieurs auteurs, il s’agit des colonies de fourmis proposées, entre autres,
dans (Stützel, 1998), (Stützel & Hoos, 2000) et (Ramalhinho Lourenço & Serra, 2000).
3. Méthodes de résolution de la RO
En majorité, les outils utilisés par les chercheurs dans le cadre des études relatives aux
problématiques de transport, correspondent à des méthodes de résolution des problèmes
d’optimisation combinatoire. Ces méthodes de résolution sont classées en deux grandes
81
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
catégories, comme le montre la Figure III.2. D’un côté, nous avons les méthodes exactes qui
permettent de trouver et prouver l’optimalité. Cependant, pour la plupart des problèmes, les
temps de calcul augmentent exponentiellement en fonction de leur taille. Ainsi, d’autre part,
nous avons les méthodes approchées qui permettent de trouver des solutions faisables,
généralement de bonne qualité, avec des temps de calculs raisonnables.
Méthodes de
résolution
Exactes Approchées
- Programmation linéaire
Dans (Chvatal, 1983), une étude détaillée de la programmation linéaire, à partir de ces
origines, et à travers la présentation de différentes méthodes de résolution. Ainsi, la
programmation linéaire est une discipline mathématique relativement récente, qui doit
son émergence à la conception de la « méthode du simplexe » par G.B. Dantzig, en
1947, pour la résolution d’un problème de planification de l’US Air Force (qui a été
formulé sous la forme d’un programme linéaire).
Un problème de programmation linéaire est un problème de maximisation (ou
minimisation) d’une fonction linéaire, sujet à un nombre fini de contraintes linéaires.
En général, on utilise l’indice j pour les n variables « de décisions » et l’indice i pour
les m contraintes, ce qui nous permet d’écrire un PL sous la forme suivante :
𝑛
max ∑ 𝑐𝑗 𝑥𝑗 (3.5)
𝑗=1
82
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
s.c.
𝑛
∑ 𝑎𝑖𝑗 𝑥𝑗 ≤ 𝑏𝑖 𝑖 = (1,2, … , 𝑚) (3.6)
𝑗=1
𝑥𝑗 ≥ 0 𝑗 = (1,2, … , 𝑛) (3.7)
Cette écriture du PL est appelée la forme standard, avec (3.5) la fonction « objectif »
et (3.6) l’ensemble des m contraintes.
Une première méthode de résolution intuitive d’un PL, est la résolution graphique
(Nemhauser & Wolsey, 1988). Cette méthode ne peut être appliquée que si le
problème ne dépasse pas les 3 variables de décision. Chaque problème a un certain
nombre de solutions qui vont satisfaire toutes ces contraintes. On les appelle
l’ensemble des solutions réalisables S, cet ensemble peut être vide. Les autres
solutions, qui satisfont une partie des contraintes ou aucune d’elle, sont appelées des
solutions irréalisables. Ainsi, la représentation graphique des différentes contraintes,
permet de définir l’ensemble S (comme le montre la Figure III.3). La représentation de
l’objectif sous la forme d’une droite qu’on fait translater dans le sens de l’optimisation
(maximisation ou minimisation), nous permet l’obtention de la solution optimale (cette
solution peut être unique ou multiple). L’espace S est aussi appelé enveloppe convexe.
Aussi, la solution optimale se trouve forcément sur l’un de ses points extrêmes.
83
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
max ∑ 𝑐𝑗 𝑥𝑗 + ∑ 𝑑𝑘 𝑦𝑘 (3.8)
𝑗=1 𝑘=1
s.c.
𝑛 𝑝
∑ 𝑎𝑖𝑗 𝑥𝑗 + ∑ 𝑒𝑖𝑘 𝑦𝑘 ≤ 𝑏𝑖 𝑖 = (1,2, … , 𝑚) (3.9)
𝑗=1 𝑘=1
𝑥𝑗 ∈ ℤ+ 𝑗 = (1,2, … , 𝑛), 𝑦𝑘 ∈ ℝ+ 𝑘 = (1,2, … , 𝑝) (3.10)
Comme mentionné dans le préambule de cette sous-section, la résolution du PL obtenu
par relâchement de la contrainte d’intégrité, ne coïncide que très rarement avec la
solution optimale du PLNE. En effet, à la différence du PL, en général la solution
optimale du PLNE se trouve à l’intérieur de l’enveloppe convexe (Figure III.3) et non
pas sur sa frontière. De là, la résolution d’un PLNE requiert l’application d’approches
d’énumération et / ou de décomposition (Sun & Quilliot, 1993). Parmi ces approches
de résolutions : la procédure par séparation et évaluation (Lawler & Wood, 1966) et
(Narendra & Fukunaga, 1977), la programmation dynamique (Bellman, 1956), les
méthodes de coupes (Padberg & Rinaldi, 1991), les méthodes de génération de
colonnes (Vanderbeck & Wolsey, 1996), etc. Aussi, le problème formulé par un
PLNE, peut être formulé par la programmation par contraintes « PPC » (Mackworth,
1977).
- Procédure par séparation et évaluation « B&B »
Cette approche de résolution trouve ses racines dans (Land & Doig, 1960). Par la
suite, ses principes ont été développés par plusieurs auteurs tels que dans : (Bertier &
Roy, 1964), (Dakin, 1965) et (Roy, et al., 1965). Les algorithmes de type B&B
consistent à énumérer de manière implicite les solutions du PLNE, pour en trouver la
meilleure. Cependant, cette démarche implique des temps de calcul trop importants. A
titre d’exemple, l’exploration de toutes les solutions d’un PLNE avec 60 variables de
décision en (0,1), en utilisant un ordinateur qui met 10-9 secondes pour explorer
chaque solution, il faudrait 30 ans pour évaluer toutes les solutions possibles (Minoux,
1983). Ainsi, les schémas se basant sur l’énumération, procèdent de manière non
exhaustive.
La procédure de B&B a été décrite en détail par plusieurs auteurs tels que dans :
(Minoux, 1983), (Sakarovitch, 1984), (Nemhauser & Wolsey, 1988) et (Dell'Amico, et
al., 1997). Cette procédure se compose de deux étapes :
La séparation : cette étape consiste à diviser le problème en plusieurs sous-
problèmes. Chaque sous-problème a un ensemble de solutions réalisables, dont
l’union recompose l’ensemble des solutions du problème. Cette décomposition
se fait à travers la représentation de l’ensemble des solutions S sous la forme
d’une arborescence. Les nœuds de cet arbre correspondent aux sous-problèmes
(le nœud de niveau 0 est appelé la racine et correspond à l’ensemble S). Pour
construire un niveau de l’arborescence, il faut choisir une variable de manière
arbitraire et lui fixer une valeur. La suite est effectuée à travers la deuxième
étape de la procédure.
84
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
85
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
Figure III.5 : Illustration d'un graphe pour la recherche du chemin le plus court entre x et y
86
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
- Heuristiques
Une heuristique est une règle spécifique à un problème donné, qui permet de trouver
une, voire des solutions réalisables, souvent loin de l’optimalité, mais le temps de
réponse est très court, en contrepartie. Elles sont principalement utilisées dans les cas
où il n’y a pas d’autres alternatives, pour la résolution de problèmes avec des
contraintes de temps et d’espace.
- Méthodes de recherche locale
La recherche est initiée à partir d’une solution existante, suivie d’une recherche de
voisinage, pour améliorer la solution. Il y a trois types de recherche locale : 1- à
voisinage proche « LS », 2- à large voisinage « LNS » et 3- à voisinage variable
« VNS ». Ces méthodes sont décrites en détail dans (Siarry, 2014).
Ce qui est appelé voisinage, c’est l’ensemble des solutions pouvant être obtenues par
application d’une transformation (permutations ou mouvements) de la solution
courante. Parmi les procédures de transformation : 1- la recherche du meilleur
d’abord : qui consiste en la construction d’un arbre de recherche, où à chaque
itération, le voisin sélectionné est celui qui permet l’amélioration de la solution (le
premier à être trouvé). 2- la recherche du meilleur voisin : même chose que la
première procédure, sauf que le meilleur voisin choisi est sélectionné à partir de tous
les voisins de même niveau de l’arbre de recherche (i.e. tous les voisins sont testés et
celui qui améliore le plus la solution est sélectionné).
La solution obtenue à l’issue de l’application d’une méthode de recherche locale, est
appelée optimum local.
- Métaheuristiques
Les métaheuristiques permettent de trouver des solutions réalisables, de la meilleure
qualité possible, en un temps raisonnable. De plus, elles sont génériques, puisqu’elles
peuvent être adaptées à n’importe quel problème.
Créées au début des années 1980, elles ont connu un développement très important ces
dernières années. Elles sont conçues pour la résolution de problèmes d’optimisation
combinatoire difficiles, qui ne peuvent être résolus à travers des méthodes exactes ou
des heuristiques. Formellement, selon (Osman & Laporte, 1996), une métaheuristique
est un processus itératif qui pilote une heuristique subordonnée, combinée à diverses
techniques, dans le but d’exploiter et d’explorer au mieux l’espace de recherche. Des
stratégies d’apprentissage sont utilisées pour structurer l’information afin de trouver
des solutions proches de l’optimum (voire optimales dans certains cas).
87
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
88
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
89
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
phéromones. Tout cela, sans que les fourmis n’aient une vision globale du
trajet.
Appliqué à un problème d’optimisation, ce principe se présente comme un
mécanisme de construction d’une population de solutions à chaque
itération. En effet, chaque fourmi construit une solution complète et laisse
une trace de phéromone sur le chemin de construction de la solution, qui est
proportionnelle à la valeur de la fonction objectif (dans le cas d’un
problème de maximisation et inversement dans le cas d’une minimisation).
Lors de la construction de la solution, chaque fourmi sélectionne la
prochaine composante de la solution qu’elle construit, en fonction des
traces de phéromones précédentes (exploitation de l’expérience acquise),
ainsi qu’en fonction d’une heuristique (orientée ou aléatoire) qui permet
d’explorer d’autres solutions (diversification) et ne pas converger et rester
bloqué dans des optimums locaux.
Pour le problème d’affectation généralisée, on compte plusieurs tentatives
de résolution avec des algorithmes de type ACO, qui ont présenté de très
bonnes performances (comme cité dans la section 2.1 de ce chapitre).
Dans ce travail de recherche, nous allons adapter l’ACO pour résoudre le
FRTSP.
4. Problèmes de replanification
L’absorption des perturbations qui sont assez fréquentes, durant les périodes d’exploitation,
est l’un des challenges les plus importants pour les entreprises. Elle a un impact direct sur la
compétitivité et la performance de l’entreprise. Ainsi et pour faire face aux perturbations,
plusieurs outils d’aide à la décision existent, dont des solutions de replanification qui
améliorent sensiblement la gestion des perturbations (Cauvin, et al., 2009). Assez répondu en
gestion de production, la replanification trouve beaucoup d’applications, dans les différents
niveaux de la chaine logistique.
Dans ce travail de recherche, nous allons considérer les changements en cours d’exploitation,
en particulier celles concernant la demande, en utilisant une approche classique d’horizon
glissant, pour re-planifier le FRTSP. Cette approche est basée sur une formulation
mathématique de type programmation linéaire, avec des caractéristiques d’exécution
spécifiques et dont l’objectif est de minimiser les changements, par rapport aux décisions
prises lors de la prise de décision originale.
90
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
Dans toutes les activités liées à la production et à la logistique, les données concernant les
commandes des clients ne sont que rarement connues avec certitude. Lorsque c’est le cas, la
demande est connue sur un horizon de temps restreint (l’étendue de l’horizon varie en
fonction de l’activité et de l’industrie concernées). D’où, le premier paramètre à définir, c’est
l’horizon de planification « HP ». Pour comprendre le principe de cette approche dans le cas
général, on suppose que cet horizon est composé de « P » périodes.
De plus, même sur une période de temps donnée, aussi bien en production qu’en logistique
(plus globalement), les aléas et perturbations (telles que des défaillances humaines ou
matérielles, voire évènements exceptionnels externes) sont des éléments avec lesquels les
entreprises doivent composer. La préservation de leur performance est fortement liée à leur
capacité de prise en considération de ces derniers. C’est pourquoi chaque nombre donné de
périodes « ΔP », un réajustement des décisions prises au plan précédent, doit être effectué.
Pour cela, un nouveau paramètre est défini, c’est le pas de planification « r ». En théorie, la
valeur maximale que peut prendre r est HP , cependant et selon (Thierry, 2003), il est
préférable d’avoir une fréquence de planification plus importante, dans le but de prendre en
compte les aléas au plus tôt. Même si dans (Galasso, 2007), l’auteur est d’accord avec ce
raisonnement, il note toutefois que choisir un pas de planification long, entrainant un gel des
décisions lointaines, permet d’assurer la stabilité des plans précédents.
Enfin, même si une entreprise doit intégrer les conséquences de l’occurrence d’aléas de tous
les types, dans ses plans d’activités, il ne reste pas moins nécessaire d’assurer la stabilité des
systèmes sur un certain horizon, où les décisions prises lors du dernier plan, ne peuvent être
modifiées. D’où, la dernière notion à définir est l’horizon gelé « HG ». L’intérêt majeur de HG
est la limitation de la nervosité du système, à travers la non considération d’éventuelles
modifications tardives (Fontan, et al., 2001). A savoir que le complément de HG dans HP est
appelé l’horizon libre « HL » (i.e. les décisions prises pour une application dans HL peuvent
être remises en question par la nouvelle planification).
Les différentes notions définies, sont illustrées sur une échelle temporelle, par la Figure III.6.
Trois itérations sont schématisées, telles que :
- L’itération 1 représente le plan initial, qui est établi à travers la résolution d’un modèle
mathématique, dont l’objectif est d’obtenir le meilleur plan initial, au vu des données
disponibles.
- Les itérations 2 et 3 (et le suivantes) représentent les plans réajustés, en considérant les
changements durant la période 𝑟 écoulée et ayant un impact sur les périodes futures.
Ces plans sont obtenus à travers la résolution d’un modèle mathématique de
91
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
replanification, dont l’objectif est de minimiser les changements engendrés par la prise
en compte des perturbations.
r
Itération 1
HG
HP
²
r
Itération 2
HG
HP
r
Itération 3
HG
HP
92
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
identiques parallèles est considéré. Dans un premier temps, un plan de production initial est
établi, avec un objectif de réduction de la durée totale de production. Puis, durant la phase de
production, des perturbations liées aux pannes des machines sont considérées. Un modèle de
replanification est proposé, avec un objectif de réduction des changements par rapport au plan
initial. Pour cela, les dates d’achèvement de chaque tâche dans le premier plan, deviennent
des dates échues du modèle de replanification, ainsi, l’objectif de ce dernier devient la
réduction des retards « virtuels » de production.
Dans (Sabar, et al., 2009), les auteurs ont proposé une approche multi-agents, pour traiter un
problème de replanification de l’affectation du personnel dans un centre d’assemblage.
L’affectation s’effectue en fonction de la compétence des employés, de leur mobilité et de
leur préférence, avec un objectif de minimisation des coûts opérationnels et de l’insatisfaction
des employés. Les perturbations considérées pour la replanification de l’affectation sont celles
relatives à l’absentéisme des employés. Le principe adopté dans cet article est d’établir un
plan d’affectation initial, avant le début de la période de production. Puis, au début de la
production et après constat des personnes absentes de leur poste de travail, le modèle de
replanification n’est lancé que sur les activités devant être exécutées par les absents. Ces
activités doivent être affectées au personnel présent avec les mêmes contraintes et objectifs,
que dans le plan initial. Cette approche a montré une bonne performance, tant sur le plan de la
qualité de la solution proposée, que sur le temps de calcul nécessaire à l’élaboration du
nouveau plan.
Un autre problème de planification logistique a été étudié dans (Silva, et al., 2008). Ce
problème concerne la considération d’annulation de certaines commandes. L’article effectue
une comparaison entre deux méthodes de résolution : GA et ACO, dans la résolution de divers
problèmes d’optimisation (problème du voyageur de commerce, problème d’affectation
quadratique, problème de tournée de véhicules et problème de job-shop), dans un contexte
statique, puis dynamique (avec replanification dans le cas des perturbations). Cette étude a
montré que les deux métaheuristiques étaient performantes (i.e. elles fournissent des solutions
équivalentes), cependant, la métaheuristique GA fournit plus rapidement les solutions dans le
cas des problèmes statiques, à contrario, l’ACO est beaucoup plus rapide dans le cas
dynamique, puisqu’elle capitalise les calculs du cas statistique, grâce à la phéromone, pour
poursuivre la recherche et fournir un réajustement de la solution.
D’autre part, et dans un contexte où le challenge environnemental est au centre de toutes les
préoccupations, dans (Salido, et al., 2017), les auteurs proposent un modèle de replanification
dans un atelier de type job-shop. Mêmes préoccupations que celles présentées par les autres
travaux, à savoir, dans un environnement dynamique, les perturbations font partie du
quotidien de l’entreprise. Suite à l’occurrence de ces perturbations, les plans initiaux ne sont
plus optimaux et requiert des réajustements. Pour cela, il est proposé un modèle de
replanification implémenté à travers un algorithme mimétique, dont l’objectif est double. En
effet, il s’agit de préserver le temps de production total, tout en réduisant la consommation
énergétique sur la période impactée par les changements dus à la replanification.
93
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
Dans (Vieira, et al., 2003) et (Herrmann, 2006), une revue de littérature de plusieurs travaux
en replanification a été réalisée et une formalisation des concepts et définitions usités a été
proposée.
94
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
2012), les auteurs ont proposé un PL en variables mixtes pour minimiser les retards
des trains. Même chose dans (Rodriguez, 2007), où l’auteur a proposé deux modèles
PPC, pour planifier le passage des trains par des nœuds ferroviaires. L’objectif étant
de minimiser les retards liés aux conflits.
2- Replanification du matériel roulant (trains, wagons et locomotives) : la recherche au
niveau de la replanification des ressources utilisées par le transport ferroviaire, aussi
bien matérielles (infrastructure et matériel roulant) qu’humaines, n’est pas aussi dense
que dans le cas précédent (Nielsen, et al., 2012). Dans (Huisman, et al., 2005), les
auteurs ont fait une revue de littérature des travaux dans ce domaine.
Un premier modèle commercial a été élaboré pour le compte de la Société Nationale
des Chemins de Fer Français « SNCF », par SABRE Technology Solutions (Ben-
Khedher, et al., 1998), pour la planification aux niveaux stratégique et tactique, de
l’activité de transport ferroviaire à grande vitesse. Dans un premier temps, ce modèle
propose des plannings basés sur les prévisions de la compagnie. Dans un second
temps, un réajustement des plannings et du matériel roulant utilisé, est effectué en
fonction de la variation de la demande effective et des dernières prévisions.
Dans (Nielsen, et al., 2012), une approche par horizon glissant est proposée, pour la
replanification du matériel roulant, dans le cas d’une perturbation majeure, au niveau
du système ferroviaire. L’objectif du modèle de replanification développé, étant la
minimisation de la déviation dans l’inventaire de tous les types de matériels roulants et
dans toutes les stations de la ligne, par rapport à l’inventaire du plan d’activité initial.
3- Replanification du personnel : la replanification du personnel est un axe de recherche
assez récent, qui intervient plus, dans les cas de perturbations de type majeur (ou
interruption de service). Un premier travail de recherche (Walker, et al., 2005), s’est
intéressé à l’établissement des plannings de circulation des trains, en planifiant le
matériel roulant et les équipages des trains (organisés par équipe) simultanément, avec
un objectif de minimisation des temps morts des trains, durant les heures d’activité de
chaque équipage. Un modèle mathématique de type PLNE a été développé pour cela.
Puis, lors de l’occurrence de perturbations majeures, un second PLNE, construit sur la
base du premier, a été établi avec un objectif de minimisation des changements par
rapport au plan initial. Le problème est résolu par une approche de type génération de
colonnes.
Dans un autre article (Veelenturf, et al., 2012), la replanification du personnel se réduit
à celle des conducteurs de trains. La replanification à ce niveau, peut nécessiter une
replanification au niveau du matériel roulant, voire, au niveau du planning de
circulation des trains. En effet, s’il n’y a pas de conducteur pour conduire un train
planifié, il sera nécessaire de procéder à un réajustement. Cette étude se focalise sur
l’occurrence de perturbations majeures. Pour cela, le modèle de replanification se
présente sous la forme d’un PLNE, dont la fonction « objectif » à minimiser, se
compose de trois membres : la déviation par rapport au plan des conducteurs initial,
les pénalités liées à l’annulation de certaines tâches des conducteurs, ainsi que les
pénalités liées aux retards. Le problème est résolu par une approche de type génération
de colonnes, couplée à une heuristique lagrangienne. Pour réduire les temps de calcul,
seul un nombre limité de trains est considéré lors de chaque replanification.
95
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
- L’étude de l’intégration d’un nouveau service, dans un système existant, dont une
expérimentation sur le terrain serait trop coûteuse.
- Formaliser la dynamique du nouveau service et étudier les flux y afférents.
- L’étude de différents scénarii, traduisant plusieurs stratégies de gestion du nouveau
service.
- Le couplage avec des modèles d’optimisation, pour évaluer les performances de la
solution optimale au vu de la dynamique du système, en considérant d’autres aspects.
- Utiliser le modèle de simulation pour valider les modèles d’optimisation.
- Les résultats visuels sont un moyen de compréhension et de confiance, qui permettent
de faire évoluer le modèle et, par conséquent, le nouveau service proposé.
- Le modèle de simulation est un démonstrateur pour les décideurs des parties
prenantes.
96
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
l’étude des systèmes de grande taille composés de plusieurs éléments en interaction. Elle
permet ainsi, de répondre à certains problèmes, à chaque fois que l’expérimentation en
grandeur nature, se révèle impossible et/ou trop coûteuse (Carrie, 1988).
Simuler revient à créer un modèle informatique fidèle au système étudié de par son
comportement. A l’aide du modèle défini, on peut mettre au point de véritables plans
d’expériences pour faire des tests sur certaines décisions à prendre, certains paramétrages
modifiés (Zeigler, 1976). Selon l’article (Pritsker, et al., 1989), l’objectif de la simulation sur
ordinateur est de reproduire et anticiper, prédire le comportement dynamique d’un système
existant (ou futur), modélisé de façon virtuelle, donc sans aucun risque.
Le processus de simulation détaillé a été décrit dans (Pritsker, 1986), il se compose de dix
étapes, suivant un processus incrémental. Ce dernier peut être schématisé à travers la Figure
III.7.
Formulation du problème
Modélisation
Vérification du modèle
Validation du modèle
Exécution de la simulation
97
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
La démarche de réalisation d’un modèle de simulation est structurée en plusieurs étapes, dont
le nombre varie d’un auteur à un autre (Le Moigne, 1990), (Anu, 1997), (Vernadat, 1999), etc.
Nous proposons de regrouper les différentes étapes, en trois étapes principales :
Le système
- Objectifs de l’étude
Un modèle Des résultats
1- La première étape consiste à décrire précisément les entités du modèle et ses règles de
gestion. On décrit ensuite la dynamique du système. Dans le cas d’un système simple
avec un faible nombre d’entités, une représentation graphique peut s’avérer suffisante.
Dans le cas de systèmes plus complexes, on peut être amené à utiliser des outils de
formalisation comme les réseaux de Petri. Cette étape de formalisation est très
importante car elle évite de remettre en cause le modèle de simulation ultérieurement.
2- La deuxième étape concerne le codage. Elle s’appuie sur la description du système
réalisée précédemment.
3- La troisième étape concerne la vérification du programme afin de s’assurer qu’il fait
ce que l’on veut qu’il fasse. Il s’agit de vérifier le code et sa logique. La validation du
modèle est une étape très importante et délicate car, dans le cas d’un système
complexe, il n’existe pas de méthode permettant de s’assurer que le modèle représente
correctement le système de départ. On s’appuie alors sur :
Des techniques d’animation du modèle montrant les entités en déplacement
conformément aux règles de fonctionnement.
Des comparaisons avec des résultats mathématiques théoriques lorsque ces
derniers sont disponibles.
Des comparaisons avec les résultats fournis par des modèles conçus
indépendamment.
Enfin, dans le cas d’un modèle de système existant, il faut souvent recourir à
des spécialistes du domaine pour interpréter les résultats obtenus par
simulation, et vérifier qu’ils sont en conformité avec le comportement du
système réel.
Cette démarche peut être formalisée par la Figure III.9.
98
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
Hypothèses simplificatrices
Objectifs de l’étude FORMALISER
Formalisation
DEVELOPPER
Programme informatique
VALIDER
Résultats numériques
Résultats statistiques
- Les entités : ce sont les premiers éléments à identifier dans un système, pour pouvoir
les modéliser. Les entités sont les objets dynamiques de la simulation. Elles sont
créées et peuvent être dupliquées, regroupées ou disposées du système. Selon le
99
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
système à modéliser, les entités représentent des éléments réels ou fictifs. Elles
peuvent modéliser les pannes d’une machine donnée. Comme elles peuvent être
assimilées à des pièces mécaniques, des personnes, des encours de production, etc.
- Les attributs : ils servent à personnaliser les entités. Un attribut est généralement
commun entre toutes les entités, mais sa valeur change d’une entité à une autre. De ce
fait, les attributs sont considérés comme des variables locales pour chaque entité.
- Les variables (globales) : les variables contiennent des informations sur le système.
Nous distinguons deux types de variables globales :
Les variables prédéfinies (états des files d’attente, états des machines, etc.).
Les variables définies par l’utilisateur (variables spécifiques au modèle).
Contrairement aux attributs, les variables n’appartiennent pas à une entité spécifique,
mais à tout le système.
- Les ressources : une ressource est affectée à une entité, afin de la transformer, de lui
rajouter de la valeur ou juste pour la faire attendre. Une ressource peut être une
personne, une machine, etc. Elle peut contenir un ensemble de serveurs identiques en
parallèle.
- Les files d’attente : lorsqu’une entité doit attendre la libération d’une unité
appartenant à une ressource, elle est placée dans une file d’attente. Cette dernière se
caractérise par une capacité et une règle de priorité.
- Les accumulateurs des statistiques : ils servent à mesurer les différentes grandeurs
et indices de performance. Parmi eux, nous distinguons :
Le temps moyen de séjour des entités dans le système.
Le temps de séjour des entités dans une file d’attente.
Le nombre moyen d’entités dans une file d’attente.
Le taux d’occupation d’une ressource.
Les accumulateurs se basent sur le temps actualisé (TNOW), pour leur calcul.
- Les évènements : un évènement est une action de durée nulle, suite à laquelle, les
valeurs des attributs, des variables ou des accumulateurs sont changées. Nous citons
comme exemples d’évènements :
L’entrée de nouvelles entités dans le système.
La fin de traitement d’une tâche.
Dans ARENA, les évènements programmés pour une simulation sont stockés dans un
calendrier d’évènements. Pour chaque futur évènement, on enregistre dans le
calendrier des évènements :
Une identification de l’entité concernée.
Le temps d’exécution de l’évènement.
Le type d’évènement.
Les évènements dans le calendrier sont triés selon leurs dates d’occurrence, de
manière à ce que les prochains évènements soient les premiers dans la liste. Lorsque la
date d’un évènement arrive, son enregistrement est effacé de la liste du calendrier des
évènements, et les informations que contient son enregistrement sont utilisées pour
produire l’évènement voulu.
- L’horloge de la simulation : le temps actualisé de la simulation est contenu dans la
variable « horloge de la simulation ». Cette variable ne balaie pas toutes les valeurs du
100
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
temps. Elle interagit avec le calendrier des évènements et elle actualise sa valeur
uniquement lorsqu’un évènement se produit. Ceci pour éviter un comptage intégral en
vain, puisque généralement nous n’avons pas besoin des valeurs du temps entre
évènements. Cependant, cette manière de faire présente l’inconvénient de « sauter »
les conditions temporelles.
- Le générateur de nombres aléatoires « GNA » : l’efficacité du simulateur à donner
des résultats valides et proche de la réalité, est directement liée à sa capacité à générer
l’aléa. La génération de l’aléa par un programme informatique, permet d’obtenir des
nombres pseudo-aléatoires. Le principe consiste à générer des nombres, d’une
distribution continue et uniforme entre 0 et 1. Les nombres générés sont des
observations. Par la suite, et selon la loi choisie par l’utilisateur, les valeurs
correspondantes des variables aléatoires seront calculées par la fonction inverse de la
loi dans le cas continu, et par une décomposition en intervalles, des probabilités
d’occurrence de chaque valeur que peut prendre la variable aléatoire, dans le cas
discret.
Les dernières versions d’ARENA utilisent un nouveau code de génération de nombres
aléatoires, dit CMRG (combined multiple recursive generator). Il utilise un système
de flots et sous-flots, pour subdiviser un cycle de génération d’une longueur de
3.1 ∗ 1057 . Lorsque nous aurons à effectuer plusieurs réplications, ARENA bascule
automatiquement d’un sous-flot à l’autre, lors de chaque nouvelle réplication.
La programmation en VBA peut être utilisée dans ARENA de deux manières différentes :
- En insérant un bloc VBA dans le modèle logique. Ainsi, à chaque fois qu’une entité
passe par ce bloc, le code correspondant s’exécute.
- En l’ouvrant avec l’éditeur de VBA directement à partir d’ARENA, et en
implémentant le code au moment opportun (la Figure III.10 montre les moments
opportuns d’implémentation d’un code).
Les primitives de VBA sont reliées aux différents évènements apparaissant lors d’une
simulation, voire, avant le début de la simulation. Les évènements sont classés en trois
grandes familles :
101
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
1. RunBegin
2. ARENA teste et initialise le modèle
3. RunBeginSimulation
4. RunBeginReplication
5. ARENA exécute la réplication
OnKeyStroke
UserFunction
6. RunEndReplication
7. RunEndSimulation
8. ARENA termine la simulation
9. RunEnd
Figure III.10 : Evènements VBA lors de l’exécution de la simulation
102
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
Figure III.11 : Illustration de la démarche de couplage – optimisation pour simulation (Borodin, et al., 2017)
Figure III.12 : Illustration de la démarche de couplage – simulation pour optimisation (Borodin, et al., 2017)
103
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
Dans ce travail de recherche, le problème étudié met en interaction plusieurs types d’entités et
de files d’attentes (relatifs aux marchandises et aux trains). De plus, les flux relatifs aux deux
entités sont indépendants et, à travers des règles décisionnelles, ils sont associés sur des
portions différentes du parcours. Le problème étudié est assez complexe et les éventuelles
perturbations auxquelles il est soumis, accroissent cette complexité. Pour prendre en charge
ses particularités et sa complexité, nous allons utiliser le dernier type de couplage.
104
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
Figure III.13 : Illustration de la démarche de couplage – simulation et optimisation (Borodin, et al., 2017)
7. Conclusion
Dans ce chapitre, nous avons présenté les différents outils d’aide à la décision, qui seront
utilisés dans le cadre de ce travail de recherche. Nous avons commencé par relever à travers la
littérature, la complexité des problèmes liés à la logistique en général, et au transport en
particulier. L’optimisation combinatoire tient une place importante dans la prise en charge de
ces problématiques. Ainsi, nous avons présenté un panorama des méthodes de résolution
exactes et souligné certains de leurs inconvénients, tels que :
- Des temps de calcul trop importants, dans le cas des grandes instances.
- Une incapacité à prendre en compte certains aspects des problèmes réels.
Pour compléter les méthodes exactes, nous avons présenté les approches de résolutions
approchées, telles que les heuristiques et les métaheuristiques. En effet, ces dernières
permettent de fournir des solutions rapidement, avec pour inconvénient, l’incapacité d’évaluer
ou garantir, un seuil de qualité ou de performance des solutions obtenues.
L’activité de transport en général, et celle des marchandises en particulier, est caractérisée par
de fortes perturbations, sous différentes formes (défaillance d’un transporteur, commandes de
dernière minute, circulation importante, pannes techniques, etc.). Ceci requiert des outils
d’aide à la décision adaptés à cet environnement, dans le but de maintenir la performance des
décisions prises en cours d’exploitation. Dans ce chapitre, nous avons présenté l’approche par
horizon glissant pour la replanification en prenant compte des perturbations. Cette approche a
été utilisée dans le cadre de plusieurs problèmes en logistique, mais aussi dans le transport
ferroviaire. Nous avons montré que cette approche dispose de plusieurs paramètres, pour
ajuster sa performance et lui permettre de proposer une solution satisfaisante avec les
nouvelles données. Cependant, elle compte quelques inconvénients parmi lesquels :
105
CHAPITRE III : Etat de l’art des méthodes de modélisation des problèmes de transport
Enfin, le couplage simulation / optimisation permet de bénéficier des avantages des deux
approches. La simulation permet de :
Et les prises de décision étant déduites de la résolution du modèle d’optimisation, avec les
données du modèle de simulation.
Ainsi, dans la suite de ce travail, les outils décrits dans ce chapitre seront utilisés comme suit :
- La nouvelle solution de transport étudiée sera modélisée sur ARENA, dans le but de
comprendre son comportement et d’étudier sa dynamique.
- Les problèmes décisionnels étudiés seront modélisés par des PLVM. Ces derniers,
dans les cas prédictifs, seront résolus par des méthodes exactes (pour les instances de
taille modérée) et des méthodes approchées (pour les instances de grande taille).
- Les perturbations seront prises en charge par l’approche de replanification à horizon
glissant.
- Un couplage simulation / optimisation sera développé. Dans un premier temps, il
permettra de valider les PLVM, en faisant tourner la simulation avec les solutions
issues de leur résolution exacte. Dans un second temps, il permettra la prise en charge
des perturbations enregistrées en cours de simulation.
106
Chapitre IV
107
CHAPITRE IV : Modélisation et simulation du FRTSP
1. Introduction
Dans le présent chapitre, nous nous concentrerons sur la problématique qui est au cœur du
système de transport ferroviaire intégrant le flux de marchandises avec le flux voyageurs.
Nous nous plaçons dans le cas de mixité maximal des ressources du système, car si ce niveau
de mixité est évalué, comme une perspective techniquement et économiquement plausible,
alors tous les autres le seront également. Comme cela a été mis en évidence dans le second
chapitre dans la figure II.13, le FRTSP est central et critique par rapport aux autres problèmes
décisionnels qui se posent. Ainsi, la faisabilité de la solution de mixité est conditionnée par les
performances réalisées à ce niveau. Plus précisément, si les trains de passagers démontrent
leur capacité à absorber une quantité de marchandises suffisamment significative, alors la
rentabilité de cette solution sera prouvée. Par conséquent, il nous a semblé plus pertinent
d’aborder en premier le FRTSP, pour ensuite s’intéresser aux autres problématiques
identifiées et directement liées à celui-ci.
Cette solution de mixité n’existe pas encore, hormis quelques expérimentations à échelle
réduite, dans certaines villes. Pour comprendre l’impact de l’intégration du flux de
marchandises, sur le réseau ferroviaire de transport de voyageurs, il est nécessaire de passer à
une échelle supérieure. Cependant, une expérimentation à grande échelle, serait trop coûteuse
et nécessiterait la mobilisation d’importantes ressources. De plus, comme mentionné dans le
chapitre 2, l’opérateur de transport de voyageurs n’a pas encore le droit d’assurer ce service,
d’où, l’impossibilité de procéder à une telle expérimentation. Dans ce cas, la démarche
usuelle suggère l’utilisation de la simulation. Quelques arguments supplémentaires
soutiennent cette proposition :
108
CHAPITRE IV : Modélisation et simulation du FRTSP
- Montrer la faisabilité de tous les autres niveaux, à travers la séparation de chacune des
ressources partagées initialement, pour aller vers une ressource dédiée à chacun des
flux.
- L’implémentation de ce cas de mixité, ne requiert pas la mobilisation de trains
supplémentaires, exclusive pour le transport de marchandises.
- Le transport de marchandises ne requiert pas la mobilisation de conducteurs
supplémentaires.
- Sachant que les trains doivent circuler pour assurer l’activité de transport de
voyageurs, le transport de marchandises ne requiert pas une consommation
supplémentaire de ressources (à noter que le transport d’un poids supplémentaire,
implique une surconsommation énergétique. Toutefois, cette dernière reste
négligeable, comparativement à la mise en circulation d’un nouveau train).
109
CHAPITRE IV : Modélisation et simulation du FRTSP
certaines stations. Malgré le peu de place disponible en milieu urbain, il peut être
réalisable d’augmenter les capacités des stations dans le cas de certaines stations.
- Hypothèse 2 : 3ème problématique - dimensionnement des trains
L’attribution d’un espace dédié au transport de marchandises, à l’intérieur d’un train
de voyageurs, dépend de deux éléments :
Comme dans un système de transport, une des difficultés majeures est de minimiser
l’écart entre la capacité disponible (quasi constante) et la demande correspondant au
flux de fret, qui est intrinsèquement variable, la capacité dédiée au transport de
marchandises sera fixée pour tous les trains. Pour cela, nous allons nous baser sur les
estimations préliminaires, de l’opérateur de transport en région Ile-de-France (études
du PREDIT présentées au chapitre 2).
- Hypothèse 3 : 4ème et 5ème problématiques - problèmes de bin packing dans les
zones de stockage temporaire
Dans le chapitre 2, nous avons soulevé l’importance du rangement des marchandises,
dans les zones de stockage temporaire. Cette problématique a un impact direct sur le
temps de disponibilité des colis pour être chargés, en raison des temps de transfert et
d’accès à ces colis. Ces temps dépendent directement des moyens de manutention et
de la politique de placement des colis dans les zones de stockage. Une modélisation
plus fine pourrait nous amener à définir un temps de disponibilité par colis et non plus
par commande (qui peut être composée de plusieurs colis), afin de se rapprocher plus
de la réalité. Dans un premier temps, nous considérons que la commande est
disponible entièrement au même moment pour être chargée et que la définition du
temps de disponibilité est déterminé en amont du FRTSP.
- Hypothèse 4 : 6ème problématique - planification des trains
La planification des trains, peut être un levier d’ajustement de l’offre de transport à la
demande. Cependant, à ce stade, l’objectif est de montrer le potentiel d’absorption
d’un flux additionnel, avec le planning actuel de circulation des trains de passagers,
durant les heures creuses de la journée.
- Hypothèse 5 : 7ème problématique - problèmes de bin packing à l’intérieur des
trains
Il est indéniable que les temps de chargement et de déchargement des commandes,
sont variables et qu’ils dépendent du nombre de colis constituant les commandes et de
leur placement dans le train. Dans le présent travail, nous allons faire une estimation
d’un temps moyen de chargement / déchargement de chaque colis.
La détermination de ce temps pourra être plus précise une fois que les technologies à
adopter seront connues.
- Hypothèse 6 : 8ème problématique - livraison des marchandises dans la station de
départ
110
CHAPITRE IV : Modélisation et simulation du FRTSP
La mise en place d’un système de régulation des arrivées des marchandises, est
indispensable à une rationalisation de l’utilisation des différentes ressources
mobilisées.
A ce stade préliminaire, nous admettons qu’aucune action n’est entreprise, pour influer
sur les moments d’arrivée des commandes, aux stations de départ.
A noter que toutes ces hypothèses, ne remettent pas en question la possibilité de mettre en
œuvre une telle solution de transport. Toutefois, elles conduiront à une exploitation sous-
optimale de cette dernière. D’autre part, elles peuvent être retenues comme première étape
d’une mise en œuvre progressive. Tout en sachant, qu’une fois la dynamique de circulation
des flux de marchandises étudiée, à travers l’étude du FRTSP, ces problématiques pourront
être étudiées, au vu des dépendances citées précédemment.
111
CHAPITRE IV : Modélisation et simulation du FRTSP
- Une meilleure utilisation de l’espace dédié aux marchandises dans les trains,
accélérant la rentabilité de la nouvelle solution
- L’augmentation de la fluidité dans le système de transport permettra la prise en charge
de nouvelles demandes de transport plus rapidement.
Pour rappel, le FRTSP correspond au cas mono-ligne, ce qui correspond au cas particulier qui
doit être étudié en premier, avant de pouvoir envisager une généralisation, vers le cas multi-
ligne, voire multimodes. En pratique, ce cas devrait être celui ayant le plus de chances d’être
déployé, au vu des inconvénients liés au changement de ligne, voire changement de mode.
Une représentation graphique est proposée dans la Figure IV.1, pour mettre en évidence le
phénomène de files d’attente qui modélise la dynamique des évènements du problème. En
particulier, les trains constituent eux même une file d’attente, de même pour chaque station de
chargement, pour laquelle il y a une file d’attente des colis en attente d’être transportés. Les
caractéristiques de cette ligne ferroviaire et des différentes composantes du service de
transport étudié, peuvent être décrites comme suit :
t=T ( )
t=2 ( )
t=1 ( )
j=1 ( ) j=3 ( 5)
j=7 ( 1) j=5 ( 3)
j=9 ( 2) j=8 ( 2)
j=10 ( 2)
Station 1
Station 4
File d’attente des
File d’attente des
commandes
commandes
Facilité de manutention, aussi bien dans les stations, qu’à l’intérieur des trains.
112
CHAPITRE IV : Modélisation et simulation du FRTSP
Il est déjà utilisé dans l’approvisionnement des magasins en ville, ainsi que par
les services postaux (transport groupé de colis).
- Les demandes de transport transmises par les clients seront appelées « commandes »
dans la suite. Durant chaque période d’exploitation considérée, nous aurons 𝐽
commandes à transporter. Chaque commande sera caractérisée par :
L’ensemble des commandes est supposé fixé et connu à l’avance, pour l’établissement
du plan de transport prédictif. Lors de l’extension en considérant la replanification de
ce plan de transport, nous autoriserons la prise en compte d’une liste de commandes
pouvant subir des changements.
- L’ensemble des trains est connu également à l’avance, il y a 𝑇 trains. On considère
uniquement les trains de voyageurs, disposant d’un espace dédié aux marchandises.
Ainsi, les trains exclusivement dédiés au transport de voyageurs, qui pourront
continuer à circuler, sont exclus de cet ensemble de 𝑇 trains). Les caractéristiques des
trains et de leur parcours se présentent comme suit :
113
CHAPITRE IV : Modélisation et simulation du FRTSP
Le temps d’attente effectif dans chaque station pour chaque train, dépendra des temps
de chargement et de déchargement des colis dans la station en question. Néanmoins,
ce temps est contraint par une borne inférieure qui est wait min et une borne supérieure
qui est wait max .
- La présence des processus d’arrivé des trains, ainsi que des marchandises, impose
d’elle-même leur modélisation par un logiciel de simulation de flux, qui est prévu à cet
effet.
- Le couplage avec un modèle d’optimisation permet d’une part, la validation du modèle
mathématique et d’autre part, l’évaluation des performances de la solution optimale
dans un milieu incertain.
Dans la Figure IV.3, qui schématise le modèle de simulation de cette solution, on peut
distinguer les deux sous-systèmes correspondant aux composantes dynamiques qui sont les
trains et les commandes. De plus, le modèle est constitué des différentes stations composant la
ligne simulée, ces dernières sont reliées par des blocs temporels simulant les temps
nécessaires à la traversée des stations. Dans ce qui suit, nous détaillons chacune des parties du
modèle comme suit :
114
CHAPITRE IV : Modélisation et simulation du FRTSP
115
CHAPITRE IV : Modélisation et simulation du FRTSP
- Le bloc « VBA » : les données relatives à chaque commande, sont enregistrées dans
un tableau dans l’environnement VBA. Ce tableau sera utilisé pour procéder aux
différents calculs, nécessaires à la simulation du processus d’affectation des
commandes aux trains. Cette dernière opération sera détaillée dans les points suivants.
- Le bloc « Dispose » : une fois les données du plan de commandes enregistrées, les
entités modélisant ces commandes sont supprimées à travers le bloc « Dispose ».
116
CHAPITRE IV : Modélisation et simulation du FRTSP
Les différents blocs utilisés dans cette séquence permettent d’effectuer les opérations
suivantes :
- « Create : Génération des trains » : il s’agit de créer autant d’entités, que de trains
devant parcourir la ligne, durant la période d’exploitation. Les entités sont toutes
générées à l’instant 0.
- « Assign » : plusieurs attributs sont affectés pour chaque train, pour représenter ses
caractéristiques. Dans la Figure IV.7, un exemple illustre les caractéristiques
attribuées, telles que la date de mise en circulation du train (l’intervalle de temps avec
le train qui le précède), l’initialisation du chargement et un numéro d’identification.
- La séquence « Queue », « Seize », « Delay » et « Release », modélise la mise en
circulation effective des trains. Le bloc « Delay » temporise chaque train, le temps que
sa date de mise en circulation arrive.
- Les blocs « VBA » : le premier enregistre dans l’environnement VBA, les données
relatives à chaque train. Ces données sont utilisées dans les calculs nécessaires à la
simulation du processus d’affectation des commandes aux trains. Le second bloc,
permet d’ajuster les dates de mise en circulation des trains, en modifiant la valeur de la
variable fixant le délai d’attente dans le bloc « Delay ».
117
CHAPITRE IV : Modélisation et simulation du FRTSP
- Les stations : elles sont modélisées par une séquence de blocs « VBA » et « Process »,
comme le montre la Figure IV.8.
Bien que le bloc « Process » dispose d’une file d’attente intégrée, le modèle développé ne
l’utilisera pas. En effet, dans le but de reproduire le comportement réel dans le transport
ferroviaire, si un train rattrape son prédécesseur (i.e. il se trouve à une distance limite de ce
dernier), un mécanisme d’ajustement de la distance est activé. Dans notre cas, il s’agit
d’augmenter le temps de parcours « TTs,s+1 » de ce train, en utilisant le bloc « VBA » (les
tests de distance sont réalisés dans VBA).
118
CHAPITRE IV : Modélisation et simulation du FRTSP
Ces trois heuristiques, nous ont permis de valider le schéma organisationnel proposé, avec des
acteurs du transport ferroviaire, ainsi que d’avoir une référence sur les conséquences d’une
potentielle mise en exploitation.
Bien que simple à mettre en œuvre, il est indispensable que chaque heuristique, considère
toutes les contraintes d’exploitation. Ces contraintes consistent en :
Ces heuristiques sont implémentées en utilisant le langage VBA. Ainsi, lors du passage de
chaque train à chaque station (représentée par un bloc VBA suivi d’un bloc Process dans la
Figure IV.8), et selon l’heuristique appliquée, certaines commandes sont affectées au train.
Cette affectation s’effectue suivant la règle de priorité et en tenant compte des différentes
contraintes. Plus précisément, les calculs pour la sélection des commandes à affecter à un
train, sont effectués lors de son passage par le bloc « VBA ». L’affectation résultante, se
traduit en un temps d’arrêt nécessaire au chargement des commandes et l’occupation d’un
certain espace, à l’intérieur du train. Ces trois heuristiques se présentent comme suit :
119
CHAPITRE IV : Modélisation et simulation du FRTSP
tailles, devrait permettre la libération rapide de l’espace dans les stations, ainsi que
l’amélioration du taux d’utilisation de l’espace.
Le principe de cette règle est de hiérarchiser l’ordre d’affectation des commandes aux
trains, en ordre décroissant, selon le nombre de colis les composant. Toutefois, comme
pour le FIFO, s’il est impossible d’affecter la commande prioritaire, l’affectation des
commandes suivantes est évaluée.
- 3ème heuristique : priorité aux commandes de plus petite taille
Les capacités des trains et leurs temps d’attente dans les stations étant fortement
contraints, l’idée étant de privilégier les commandes de petites tailles. Cette règle
s’inspire de la charge la plus faible d’abord, ce qui devrait augmenter la fluidité dans
le système. En effet, il y a plus de chances de pouvoir charger les petites commandes,
car cela prend moins de place et moins de temps également.
Les commandes ayant le moins de colis sont dans ce cas prioritaires. Comme pour les
deux règles précédentes, un processus itératif permet de tester les candidats à
l’affectation dans le cas de l’échec de l’affectation de la commande prioritaire. Le
processus se termine une fois le train saturé, ou le temps d’attente du train dans la
station épuisé.
A noter que pour chacune de ces trois heuristiques, une règle de priorité supplémentaire
s’additionne aux spécificités de chacune d’entre elle. Il s’agit de la priorité relative à l’ordre
des stations desservies par les trains. En effet, chaque heuristique est appliquée au niveau de
chaque station, ce qui octroie une priorité aux commandes, selon leur station de départ.
3.4.2. 4ème heuristique : basée sur l’optimisation du problème d’un seul train
Cette heuristique procède par décomposition et résolution exacte des sous-problèmes, tel que
chaque sous problème est relatif au chargement optimal d’un seul train. Ceci, à travers un
processus d’évaluation de différentes possibilités d’affectation des commandes à chaque train,
pour maximiser le nombre de commandes pris en charge par le train correspondant. La
résolution de chaque sous-problème est basée sur un algorithme de type séparation et
évaluation. L’intérêt de cette décomposition est lié directement à la réduction de la complexité
du problème initial en raison de la considération d’un sous ensemble des commandes
disponibles, lors du passage du train en question. Ceci réduit sensiblement, le nombre de
possibilités à explorer. En effet, comme les dates de disponibilité des commandes, ainsi que
leur taille, sont disponibles à l’avance, cela permet de construite les sous-ensembles des
commandes disponibles pour chaque train au niveau de chaque station. La solution optimale
de chaque sous-problème permet de construire une solution approchée du problème.
Pour illustrer l’algorithme de résolution de cette 4ème heuristique, nous proposons l’exemple,
dont l’arbre des possibilités est schématisé par la Figure IV.9. Nous considérons le train t1 ,
pour lequel trois commandes seront disponibles lors de son passage par la ligne ferroviaire
(j1 , j2 et j3). Chaque nœud de l’arbre correspond à une instanciation de la variable binaire de
décision d’affectation de la commande. Le problème mono-train revient à déterminer les
commandes à transporter, en maximisant leur nombre, tel que les contraintes d’espace et de
temps sont respectées. Lorsqu’il y a plusieurs solutions équivalentes, celle qui minimise le
temps d’attente des commandes, avant leur chargement dans le train, est privilégiée. Dans le
120
CHAPITRE IV : Modélisation et simulation du FRTSP
Sur la Figure IV.9, nous illustrons un exemple d’un problème mono-train en présence de trois
commandes à transporter :
1 0
1 0 1 0
1 0 1 0 1 0 1 0
ème
Figure IV.9 : Illustration de la construction d’une solution avec la 4 heuristique
Le principe de fonctionnement de la 4ème heuristique est donné par l’Algorithme IV.1. Nous
introduisons de nouveaux paramètres comme suit :
- La liste des commandes pouvant être transportées par le train 𝑡 est notée 𝐽𝑡 .
- Date d’arrivée du train 𝑡 à 𝑑𝑒𝑝𝑗 est notée : 𝑅𝑗 .
- La borne supérieure de la date d’arrivée du train 𝑡 à 𝑑𝑒𝑝𝑗 , correspond à un temps
d’attente maximum du train 𝑡 dans toutes les stations qui précèdent 𝑑𝑒𝑝𝑗 : 𝑅𝑗𝑚𝑎𝑥 =
𝑑𝑒𝑝 −1
𝑙𝑡 + ∑𝑠=1𝑗 𝑡𝑡𝑠,𝑠+1 + (𝑑𝑒𝑝𝑗 − 1) ∗ 𝑤𝑎𝑖𝑡𝑚𝑎𝑥
- La borne supérieure courante de la fonction objectif est notée : 𝐹𝑏𝑒𝑠𝑡 .
121
CHAPITRE IV : Modélisation et simulation du FRTSP
} alors
7. Affecter 𝑗 à 𝑡 (𝑥𝑗𝑡 = 1)
8. 𝐹𝑖𝑛𝑡 𝐹𝑖𝑛𝑡 calculée au nœud précèdent + [𝑅𝑗 − 𝑟𝑗 ]
9. Sinon
10. 𝑗 n’est pas affectée à 𝑡: 𝑥𝑗𝑡 = 0 (couper la branche correspondant à 𝑥𝑗𝑡 = 1)
11. Fin Si
12. Si 𝐹𝑖𝑛𝑡 > 𝐹𝑏𝑒𝑠𝑡 alors
13. Couper la branche actuelle et aller à la branche suivante
14. Fin Si
15. Si (toutes les commandes de 𝐽𝑡 ont été évaluées) et (𝐹𝑖𝑛𝑡 < 𝐹𝑏𝑒𝑠𝑡 ) alors
16. 𝐹𝑏𝑒𝑠𝑡 𝐹𝑖𝑛𝑡
17. Fin Si
18. Si (toutes les commandes de 𝐽𝑡 ont été évaluées) et (il reste des branches à
explorer) et (∃ une branche avec un nombre de commandes pouvant être affectée à
𝑡 ≥ nombre de commandes affectées en considérant la meilleure solution actuelle)
alors
19. Retour arrière jusqu’à la dernière commande affectée à 𝑡 (dernier 𝑥𝑗𝑡 = 1)
20. Ne pas affecter cette commande à 𝑡 (mettre 𝑥𝑗𝑡 = 0).
21. Fin Si
22. Prochaine commande
23. Considérer l’affectation correspondant à 𝐹𝑏𝑒𝑠𝑡
122
CHAPITRE IV : Modélisation et simulation du FRTSP
Les trois sous-processus en haut à droite de la figure, sont utilisés pour trier les
commandes, en fonction de la règle de priorité choisie.
Quant au quatrième sous-processus, il est utilisé pour récupérer la solution de
l’exécution de la 4ème heuristique.
Enfin, le diagramme encadré n°2 est dédié à la 4ème heuristique et, il doit
remplacer dans le diagramme principal, la partie encadrée n°1, qui est dédiée
au trois autres heuristiques.
Ainsi, le processus global de transport de marchandises par les trains, à travers la ligne
ferroviaire, se résume comme suit :
123
CHAPITRE IV : Modélisation et simulation du FRTSP
dernier à partir du dépôt de départ. Ainsi, lors de son passage par chaque
station, il s’agira de :
Vérifier pour chaque commande disponible dans la station, si elle doit être
transportée par le train actuel. Si c’est le cas, la charger.
Mettre à jour le temps d’attente du train dans la station actuelle.
Mettre à jour la capacité actuelle du train.
124
CHAPITRE IV : Modélisation et simulation du FRTSP
Aller à la prochaine
Non
Dernière station Sous-processus - Heuristique 1
station de la ligne ?
Trier les
commandes en
Arrivée du train à la ordre décroissant
station selon leur taille
Réinitialisation du
temps d’attente :
wait=0 Sous-processus - Heuristique 2
Trier les
Commandes commandes en
dans le train ordre croissant
Temps d’attente du selon leur taille
train dans la
station actuelle =
wait Sous-processus - Heuristique 3
Non
Oui
Récupérer le
résultat de
Wait = waitmin l’exécution de
Considérer la Oui l’algorithme
prochaine
commande
Sous-processus - Heuristique 4
1
2
Station
Non Oui
d’arrivée de la Non wait>=waitmin
Non
commande
Toutes les
commandes
Oui évaluées Oui
Toutes les
Temps d’attente :
commandes
wait = wait + time * Q
Capacité du train : évaluées
ChgT = ChgT + Q
Non
Non
Non Temps d’attente :
Capacité du train : wait = wait + time * Q
Capacité du train :
ChgT = ChgT - Q ChgT = ChgT + Q
Non Non
Oui Non
Non
Test temps d’attente
dans la station d’arrivée
Temps d’attente :
Toutes les wait = wait + time * Q
commandes
évaluées Oui
Oui
Test temps d’attente
dans la station actuelle
Oui
Si la commande
Oui doit être transportée
Commandes par ce train
dans la station Test
disponibilité de
la capacité
Oui Considérer la
Sous- Considérer la prochaine
processus prochaine commande
(heuristiques) commande
125
CHAPITRE IV : Modélisation et simulation du FRTSP
4. Conclusion
L’intérêt porté à ce problème, découle de sa capacité à évaluer le potentiel de la solution de
transport proposée, en évaluant les volumes qui peuvent être transportés. Par ailleurs, l’étude
du FRTSP permet également de fournir, les données nécessaires à la résolution de plusieurs
autres problématiques décisionnelles, identifiées dans le second chapitre. Pour de nombreuses
raisons cités dans le présent chapitre, le FRTSP est jugé comme critique et prioritaire
notamment au vu de la considération du niveau de mixité totale. Dans ce chapitre, nous avons
proposé une formalisation et une première approche de modélisation et de résolution du
FRTSP, au vu des hypothèses considérées.
Le problème FRTSP est inspiré du projet du Grand Paris Express. Le problème revient à
définir un plan de transport des marchandises en utilisant une seule ligne ferroviaire pour
passagers. Il s’agit de déterminer le train et le moment de chargement de chaque marchandise,
en attente dans sa station de départ. Plusieurs contraintes sont considérées, telles que le
respect de l’espace disponible dans les trains, le temps d’arrêt des trains dans les stations et les
temps de disponibilité des marchandises aux stations. La fonction objectif permet de
minimiser les temps d’attente des marchandises aux stations.
Dans un second temps, nous avons présenté le modèle de simulation à évènements discrets
développé sur ARENA. Ce dernier est composé de trois parties interdépendantes, à savoir, la
génération des commandes, la génération des trains et la ligne ferroviaire. Les avantages de
développer un modèle de simulation sont multiples. Entre autres, il nous permet de
comprendre la dynamique de ce système de transport ferroviaire, en représentant les différents
processus de files d’attente présents. Par ailleurs, les règles décisionnelles pour l’affectation
des commandes aux trains, sont implémentées dans l’environnement VBA. Par la suite, nous
allons montrer comment coupler ce modèle de simulation, avec le modèle de programmation
linéaire en variables mixtes, qui permet de déterminer un plan optimal pour le FRTSP.
126
Chapitre V
Modélisation mathématique et
optimisation du FRTSP – cas prédictif et
replanification
127
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
1. Introduction
Après avoir défini le problème FRTSP et étudier son comportement dynamique grâce au
modèle de simulation à évènements discrets, nous allons nous intéresser dans le présent
chapitre à l’optimisation de son processus décisionnel. D’abord nous étudions la nature
combinatoire du FRTSP en montrant que la relaxation de plusieurs de ses contraintes le
ramène au problème d’affectation généralisée. Un exemple numérique est présenté afin
d’illustrer l’ensemble des contraintes et données du problème, plusieurs solutions sont
également proposées pour saisir les enjeux de ce problème. Ensuite, nous nous intéressons à
l’établissement du plan prédictif en formulant un PLVM pour réduire les temps d’attente des
marchandises. Pour compléter cette approche, nous adapterons un algorithme de colonies de
fourmis pour les instances de grande taille. Par la suite, nous proposons une approche par
horizon glissant pour la replanification du FRTSP afin de prendre en compte des changements
relatifs à la demande. Cette mise à jour du plan initial permet d’ajuster les décisions initiales
de manière à minimiser les écarts dus au recalcul. Enfin, nous décrivons les processus de
couplage simulation / optimisation, aussi bien dans le cas prédictif, que dans le cas de la
replanification.
2. Un exemple numérique
Avant de procéder à la formalisation mathématique du problème FRTSP, dans cette section,
nous proposons un exemple numérique sur une instance réduite afin d’appréhender sa
complexité.
Nous considérons une ligne ferroviaire composée de 4 stations de voyageurs, où il est possible
de procéder aux chargements / déchargements des marchandises. Cette ligne est parcourue par
3 trains mixtes, à intervalle de 10 minutes. Enfin, 10 commandes doivent être transportées sur
cette ligne, suivant des parcours différents (i.e. les stations de départs et d’arrivées varient
d’une commande à une autre). L’exemple est illustré par la Figure V.1.
128
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
t=3 ( )
Temps d’attente :
t=2 ( ) secondes
secondes
t=1 ( ) Temps de chargement d’1 colis :
secondes
File d’attente des trains
Nous proposons deux plans différents de transport des commandes, élaborés en appliquant des
règles décisionnelles aléatoires. Ceci nous permettra de constater l’impact de toutes les
contraintes sur l’activité de transport de marchandises, ainsi que l’importance du schéma
décisionnel adopté. Le Tableau V.1 et le Tableau V.2 présentent ces deux plans. Les notations
et formules suivantes sont utilisées :
Le temps d’attente moyen de chaque commande dans sa station de départ, avant son transport
par un train, est égal à : 7 min.
129
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
- La commande 1 n’a pas été chargée dans le train 1 à la station 1, parce que son
chargement aurait nécessité une durée d’arrêt totale égale à 70 secondes (> 𝑤𝑎𝑖𝑡𝑚𝑎𝑥 ).
- La commande 5 n’a pas été chargée dans le train 1 à la station 2, parce que son
chargement aurait nécessité une durée d’arrêt totale égale à 70 secondes (> 𝑤𝑎𝑖𝑡𝑚𝑎𝑥 ).
Aussi, la capacité du train a été atteinte (𝑐𝑎𝑝1 = 10).
- Les commandes 6 et 7 n’ont pas été chargée dans le train 1, parce que le déchargement
des commandes 3 et 4 a nécessité 60 secondes (= 𝑤𝑎𝑖𝑡𝑚𝑎𝑥 ).
- La commande 5 n’a pas été chargée dans le train 2, bien que le temps d’attente fût
suffisant dans la station 2. Cette décision, a permis le transport des commandes 6 et 7
par le train 2. En effet, le chargement de la commande 5 aurait requis 10 secondes
d’attente supplémentaires dans la station 4, nécessaires à son déchargement.
Tableau V.1 : Plan de transport des commandes – 1ère solution
Bien que le temps d’attente moyen des commandes transportées soit égal à : 5 min, la
commande 10 (non transportée), le ramène à l’infini. En effet, l’incapacité à prendre en
charge l’ensemble des commandes conduit à une infaisabilité. En pratique, sous condition
d’une demande plus importante que la capacité disponible, la formulation du FRTSP peut
évoluer pour déterminer le plan pouvant transporter le maximum de commandes possibles
(pour générer le maximum de bénéfice).
130
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
l’interaction entre les différentes contraintes et la multitude de possibilités, même pour une
instance de taille réduite. Ceci met en lumière la difficulté de résolution du FRTSP.
Tableau V.2 : Plan de transport des commandes – 2ème solution
131
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
En ce qui concerne les contraintes, nous retrouvons à travers (5.2) l’ensemble de contraintes
relatives à l’obligation d’affecter chaque commande à un seul train (i.e. chaque tâche à une
seule machine). Aussi, à travers (5.3) l’ensemble des contraintes relatives à la capacité des
trains (i.e. capacité des machines).
Comme mentionné dans le chapitre 3, le problème d’affectation est considéré comme étant
NP-difficile (Fisher, et al., 1986), ainsi, notre problème l’est aussi.
j=1 ( 2)
t=1 ( )
j=2 ( 2)
j=3 ( 2)
t=2 ( )
j=4 ( 2)
j=5 ( 2)
t=T ( )
j=J ( 2)
S=1 S=2
132
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
d’exploitation. Ainsi, le carnet des commandes sera connu à l’avance et aucun changement ne
sera considéré en cours d’exploitation. Aussi, le planning de circulation des trains est fixé et
respecté durant toute la période d’exploitation.
Les notations utilisées pour la définition des données considérées dans ce problème sont
résumées dans le Tableau V.3.
Tableau V.3 : Notations utilisées dans le PLVM
S Nombre de stations
𝑡𝑡𝑠,𝑠+1 Temps de trajet entre deux stations successives
𝑄𝑗 Nombre de colis standard composant la commande j
Capacité de chargement du train t en nombre de colis
𝐶𝑎𝑝𝑡
standard
𝑊𝑎𝑖𝑡𝑚𝑎𝑥 Temps d’attente maximum des trains dans chaque station
𝑊𝑎𝑖𝑡𝑚𝑖𝑛 Temps d’attente minimum des trains dans chaque station
time Temps de chargement / déchargement d’un colis standard
M Un grand nombre positif
1 si la commande j est présente à l’intérieur du train t dans la
𝑥𝑗𝑡𝑠
Variables de
station s, 0 sinon
décision
133
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
134
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
est considéré dans 𝐽′𝑠 qui est l’ensemble des commandes, dont les stations de départ
et d’arrivée sont la station s.
- Respect de la date de disponibilité des commandes
Les commandes ne sont pas toutes disponibles à l’instant 0. De ce fait, elles ne
peuvent être transportées par un train, qu’après qu’elles soient disponibles dans leur
station de départ. Ceci donne lieu à :
𝑙𝑡 + ∑𝑠−1
𝑠′ =1(𝐶𝑡𝑠′ + 𝑡𝑡𝑠′,𝑠′+1 ) ∀𝑗, 𝑡 𝑒𝑡 𝑠 = 𝑑𝑒𝑝𝑗 𝑠𝑖 𝑑𝑒𝑝𝑗 > 1
𝑥𝑗𝑡𝑠 ∗ 𝑟𝑗 ≤ { (5.12)
𝑙𝑡 ∀𝑗, 𝑡 𝑒𝑡 𝑠 = 𝑑𝑒𝑝𝑗 𝑠𝑖 𝑑𝑒𝑝𝑗 = 1
Dans cet ensemble de contraintes (5.12), la date d’arrivée du train t à la station de
départ 𝑑𝑒𝑝𝑗 de la commande j, est calculée sur la base de :
𝑙𝑡 + ∑𝑠−1
𝑠′ =1(𝐶𝑡𝑠′ + 𝑡𝑡𝑠′,𝑠′+1 ) − 𝑀(1 − 𝑥𝑗𝑡𝑠 ) ∀𝑗, 𝑡 𝑎𝑛𝑑 𝑠 = 𝑑𝑒𝑝𝑗 𝑠𝑖 𝑑𝑒𝑝𝑗 > 1
𝑅𝑗𝑡𝑠 ≥{ (5.13)
𝑙𝑡 − 𝑀(1 − 𝑥𝑗𝑡𝑠 ) ∀𝑗, 𝑡 𝑎𝑛𝑑 𝑠 = 𝑑𝑒𝑝𝑗 𝑠𝑖 𝑑𝑒𝑝𝑗 = 1
L’introduction de M dans l’ensemble des contraintes (5.13) est nécessaire à la
linéarisation de la forme initiale de cet ensemble de contraintes. En effet, la forme de
base de (5.9) est quadratique et s’écrit comme suit : 𝑅𝑗𝑡𝑠 ≥ 𝑥𝑗𝑡𝑠 ∗ (𝑙𝑡 + ∑𝑠−1
𝑠′ =1(𝐶𝑡𝑠′ +
𝑡𝑡𝑠′ ,𝑠′ +1 )).
- Domaines de définition des variables
𝑥𝑗𝑡𝑠 ∈ {0,1}, 𝐶𝑡𝑠 ≥ 0, 𝑅𝑗𝑡𝑠 ≥ 0 ∀𝑗, 𝑡, 𝑠 (5.14)
A ce stade de l’étude, nous avons défini l’objectif le plus important comme celui relatif à la
maximation du taux de rotation des marchandises, ce qui revient à la minimisation de leur
temps d’attente au niveau des stations de départ. La formalisation mathématique de cet
objectif se présente comme suit :
𝐽 𝑇
Dans (5.15), la somme des temps d’attente de toutes les commandes dans leur station de
départ est considérée.
135
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
Nous allons développer plusieurs variantes de colonies de fourmis, ce qui nous permettra
d’identifier la variante qui serait en mesure de considérer au mieux, les spécificités du
problème étudié.
Ce problème s’intéresse à la recherche du plus court chemin, reliant n villes. Chaque ville
devant être visitée exactement une seule fois. Ce problème est formalisé par un graphe
complet, tel que les sommets représentent les villes et les arrêtes modélisent les chemins les
reliant. L’algorithme AS va procéder de manière à ce qu’à chaque itération, chaque fourmi va
construire une solution complète (i.e. un chemin parcourant toutes les villes). Plusieurs
mécanismes sont mis en place pour reproduire la recherche du meilleur chemin par chaque
fourmi, en s’inspirant du phénomène biologique. Dans notre adaptation des AS au FRTSP,
nous nous sommes également inspirés de cette modélisation, en apportant les ajustements
nécessaires, pour prendre en compte les spécificités de notre problème.
136
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
?
Fourmi n° k
Figure V.3 : Une étape durant le processus de construction d’une solution par une fourmi
137
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
que celles qui concernent l’affectation de commandes à des trains, dont la date de
passage est avant la date de disponibilité de ces mêmes commandes).
𝑘 : l’identifiant d’une fourmi.
𝛼 : paramètre qui détermine l’importance relative des traces de phéromones, c’est-à-
dire l’intensification ou l’exploitation.
𝛽 : paramètre qui détermine l’importance relative de la règle heuristique, représentant
la diversification ou l’exploration.
Les deux paramètres 𝛼 et 𝛽 permettent d’équilibrer la balance intensification /
diversification, en définissant la part relative de chacun de ces deux mécanisme.
- A la fin de la construction d’une solution complète, par chaque fourmi à chaque
𝑘
itération, chacune dépose une quantité de phéromones ∆𝜏(𝑗,𝑡)(𝑗 ′ ,𝑡 ′ ) (𝑖𝑡), qui est calculée
comme suit :
𝑄
⁄𝐿𝑘 (𝑖𝑡) 𝑠𝑖 𝑙 ′ 𝑎𝑟𝑟ê𝑡𝑒 𝑟𝑒𝑙𝑖𝑎𝑛𝑡 (𝑗, 𝑡) 𝑒𝑡 (𝑗 ′ , 𝑡 ′ ) 𝑒𝑠𝑡 𝑢𝑡𝑖𝑙𝑖𝑠é𝑒 𝑝𝑎𝑟 𝑘
𝑘
∆𝜏(𝑗,𝑡)(𝑗 ′ ,𝑡 ′ ) (𝑖𝑡) ={ (5.18)
0 𝑠𝑖𝑛𝑜𝑛
𝑘
𝐿 (𝑖𝑡) : c’est le temps d’attente total de toutes les commandes dans leur station de
départ (i.e. la valeur de la fonction objectif, en considérant la solution construite par la
fourmi 𝑘). Comme c’est un problème de minimisation, nous considérons l’inverse de
cette valeur. Ainsi, plus le temps d’attente sera grand, moindre sera la quantité de
phéromone déposée.
𝑄 : un paramètre à fixer
- La mise à jour de la quantité de phéromones, présente sur chaque arrête du graphe des
possibilités, s’effectue en considérant deux éléments :
Ceci donne lieu à la formule de calcul des quantités de phéromones sur chaque arrête,
comme suit :
𝑘
𝜏(𝑗,𝑡)(𝑗′ ,𝑡 ′ ) (𝑖𝑡 + 1) = 𝜌. 𝜏(𝑗,𝑡)(𝑗′ ,𝑡 ′ ) (𝑖𝑡) + ∑𝑚
𝑘=1 ∆𝜏(𝑗,𝑡)(𝑗 ′ ,𝑡 ′ ) (𝑖𝑡) (5.20)
A noter que la trace de phéromones est initialisée de la même manière sur toutes les
arrêtes du graphe, avec 𝜏0 ≥ 0.
L’algorithme de la variante de base AS, appliqué au problème FRTSP est présenté par
l’Algorithme V.1.
138
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
Les modifications principales qu’apporte la variante MMAS à la variante AS, peuvent être
résumées à travers les points suivants :
139
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
- Seule une fourmi est autorisée à déposer une quantité de phéromones sur le parcours
de la solution qu’elle a construit. Bien que deux possibilités sont suggérées par la
littérature, à savoir :
Ces modifications sont accompagnées par de nouvelles formules, pour calculer les nouveaux
paramètres et réajuster d’autres, comme suit :
- Le calcul des bornes de l’intervalle des quantités de phéromones, s’effectue suivant les
formules suivantes (il s’agit d’une adaptation des formules proposées dans (Stützel &
Hoos, 2000), aux spécificités du problème FRTSP) :
𝜏𝑚𝑎𝑥 = (1⁄1 − 𝜌) (1⁄ 𝑏𝑒𝑠𝑡 ) (5.21)
𝐿 (𝑖𝑡)
(1 − 𝑇√𝑝𝑏𝑒𝑠𝑡 )
𝜏𝑚𝑖𝑛 = 𝜏𝑚𝑎𝑥 ( ⁄𝑇 ) (5.22)
( ⁄2 − 1) 𝑇√𝑃𝑏𝑒𝑠𝑡
𝐿𝑏𝑒𝑠𝑡 (𝑖𝑡) : la meilleure valeur de la fonction objectif à l’itération 𝑖𝑡, en considérant la
solution construite par la meilleure fourmi « 𝑏𝑒𝑠𝑡 ».
𝑃𝑏𝑒𝑠𝑡 : il s’agit de la probabilité qu’une fourmi ait fait le meilleur choix d’affectation
des commandes, pour chaque train. Dans (Stützel & Hoos, 2000), les auteurs
proposent de la fixer à 0,05.
Le calcul de ces bornes s’effectue à chaque fois qu’une meilleure solution est trouvée.
- La mise à jour de la trace de phéromones, à la fin de chaque itération (qui était donnée
par (5.20)), devient :
k
𝜌. 𝜏(𝑗,𝑡)(𝑗′ ,𝑡 ′ ) (𝑖𝑡) + ∆𝜏(𝑗,𝑡)(𝑗 ′ ,𝑡 ′ ) (𝑖𝑡) 𝑠𝑖 𝑘 = 𝑏𝑒𝑠𝑡
(𝑖𝑡
𝜏(𝑗,𝑡)(𝑗′ ,𝑡 ′ ) + 1) = { (5.23)
𝜌. 𝜏(𝑗,𝑡)(𝑗′ ,𝑡 ′ ) (𝑖𝑡) 𝑠𝑖𝑛𝑜𝑛
140
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
best
∆𝜏(𝑗,𝑡)(𝑗 ′ ,𝑡 ′ ) (𝑖𝑡) : représente la quantité de phéromones calculée, pour la fourmi ayant
k
construit la meilleure solution à l’itération 𝑖𝑡. A savoir, ∆𝜏(𝑗,𝑡)(𝑗 ′ ,𝑡 ′ ) (𝑖𝑡) = 0, ∀𝑘 ≠
𝑏𝑒𝑠𝑡 (i.e. pour toutes les arrêtes qui ne correspondent pas au chemin parcouru par la
fourmi « best », il n’y a que le processus d’évaporation qui est considéré).
Si la quantité de phéromones mise à jour sur une arrête est supérieure à 𝜏𝑚𝑎𝑥 , elle est
initialisée à 𝜏𝑚𝑎𝑥 . D’autre part, si elle est inférieure à 𝜏𝑚𝑖𝑛 , elle est initialisée à 𝜏𝑚𝑖𝑛 .
141
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
Les modifications principales qu’apporte la variante ACS à la variante AS, peuvent être
résumées à travers les points suivants :
Ces modifications requièrent la définition de nouvelles formules, dans le but de calculer les
nouveaux paramètres. Une description de ces formules pour le problème du voyageur de
commerce, est donnée dans (Dorigo & Gambardella, 1997) et (Dorigo, et al., 2006). Nous
proposons une adaptation aux spécificités du problème FRTSP, comme suit :
- Nous avons décrit pour AS, le processus qui pilote à chaque étape de la construction
de la solution, le déplacement d’une fourmi dans le graphe, d’un sommet (𝑗, 𝑡) à un
autre (𝑗 ′ , 𝑡′). Dans le cas de la variante ACS, ce processus est enrichi avec une règle
aléatoire. Ainsi, la sélection du prochain sommet (𝑗 ′ , 𝑡′) s’effectue selon la
formulation suivante :
𝛽
arg 𝑚𝑎𝑥(𝑗𝑚 ,𝑡𝑚) {𝜏(𝑗,𝑡)(𝑗𝑚 ,𝑡𝑚 ) [𝜂(𝑗,𝑡)(𝑗𝑚 ,𝑡𝑚) ] } 𝑠𝑖 𝑞 ≤ 𝑞0
(𝑗 ′
, 𝑡′) = { (5.24)
𝑠𝑢𝑖𝑣𝑎𝑛𝑡 𝑙𝑎 𝑓𝑜𝑟𝑚𝑢𝑙𝑒 (5.17) 𝑠𝑖𝑛𝑜𝑛
𝑞 : variable aléatoire choisi à l’intérieur de l’intervalle [0,1].
𝑞0 : paramètre de la variante ACS, à fixer lors de l’initialisation. A noter : 0 ≤ 𝑞0 ≤ 1.
La valeur prise par le paramètre 𝑞0 permet de faire pencher la balance intensification /
diversification, soit d’un côté soit de l’autre.
- La mise à jour locale de la trace de phéromones, s’effectue durant la construction de la
solution par une fourmi. Une fois une arrête sélectionnée, la quantité de phéromones
qui y est associée est mise à jour, suivant la formule suivante :
𝜏(𝑗,𝑡)(𝑗′ ,𝑡 ′ ) (𝑖𝑡) = φ. 𝜏(𝑗,𝑡)(𝑗′ ,𝑡 ′ ) (𝑖𝑡) + (1 − 𝜑). 𝜏0 (5.25)
𝜑 : le facteur de persistance local (i.e. (1 − 𝜑) représente le facteur d’évaporation). Il
est compris dans l’intervalle [0,1].
Cette mise à jour locale favorise la diversification, à travers la diminution de la
désirabilité de l’arrête choisie, pour les prochaines fourmis de la colonie, durant la
même itération.
- La mise à jour globale de la trace de phéromones, présente deux légères variations
dans la formule, par rapport à la variante MMAS. Ces variations se présentent comme
suit :
142
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
A noter que la convergence de la variante ACS a été prouvée théoriquement dans (Stützle &
Dorigo, 2002).
Algorithme V.3 : L’adaptation de l’algorithme ACS pour le FRTSP à partir de AS
- Seule une fourmi est autorisée à déposer une nouvelle quantité de phéromones, à la fin
de chaque itération. Ce qui implique pour les itérations suivantes, des solutions
potentiellement voisines de la meilleure solution précédente.
- Pour éviter la stagnation du processus de recherche, les deux variantes mettent en
œuvre un mécanisme régulant la variation de la quantité de phéromones, sur chaque
arrête. En effet, pour la variante MMAS, les limites sont explicites à travers
l’intervalle [𝜏𝑚𝑖𝑛 , 𝜏𝑚𝑎𝑥 ]. Alors que pour la variante ACS, elles sont implicites et
découlent de la forme des deux formules (5.25) et (5.26). En ce qui concerne la limite
inférieure, aussi bien la mise à jour locale, que globale, n’autorisent une quantité de
phéromones inferieure à 𝜏0 . En ce qui concerne la limite supérieure, la formule de
143
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
mise à jour globale nous permet de déduire sa borne supérieure, qui est égale à
Q
⁄𝐿𝑘 (𝑖𝑡) avec 𝑘 correspondant à la fourmi ayant trouvé la meilleure solution globale
Cependant, une différence majeure existe entre ces deux variantes. Il s’agit de la règle
décisionnelle qui dirige le processus de recherche. En particulier, pour la variante MMAS,
c’est la décision probabiliste de la variante de base qui est utilisée. Cette dernière présente un
équilibre entre exploration de nouvelles solutions (diversification) et exploitation de
l’expérience acquise par les fourmis (intensification). En ce qui concerne la variante ACS, la
règle proportionnelle pseudo-aléatoire utilisée, accentue le processus d’intensification, ce qui
conduit principalement, à l’exploration de solutions voisines.
La variante hybride que nous proposons, est dérivée d’une variante proposée dans
(Ramalhinho Lourenço & Serra, 2000), et développée pour le problème d’affectation
généralisée. Nous y avons apporté quelques ajustements, pour prendre en compte les
spécificités du problème FRTSP.
Cette variante propose de combiner une partie des caractéristiques de la variante MMAS et
une autre partie de la variante ACS, comme suit :
144
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
perturbations liées aux commandes. En ce qui concerne les perturbations liées aux trains, elles
seront traitées dans d’autres travaux.
Algorithme V.4 : L’adaptation de l’algorithme MMACS pour le FRTSP à partir de AS
Les perturbations relatives aux commandes peuvent être classées comme suit :
145
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
- Les commandes annulées : il s’agit des commandes planifiées pour être transportées,
mais qui devront être retirés des plans de transport. En fonction du modèle
économique choisi, la nécessité d’utiliser au mieux l’espace libéré sera plus ou moins
critique, pour la rentabilité de la solution.
- Les commandes réceptionnées en retard : il s’agit des commandes qui arrivent à leur
station de départ, à une date ultérieure à leur date de disponibilité prévisionnelle.
Le modèle mathématique développé dans le cas prédictif, fourni un plan initial de transport
des marchandises. Ce dernier reste optimal, si aucune perturbation n’est enregistrée. Ainsi, il
détermine le train 𝑡 qui va transporter chaque commande 𝑗 à partir de sa station de départ
𝑠 = 𝑑𝑒𝑝𝑗 , avec la valeur de la variable de décision 𝑥𝑗𝑡𝑠 . Aussi, il détermine le moment auquel
cette commande sera transportée, avec la valeur de 𝑅𝑗𝑡𝑠 . En pratique, il est très rare, de ne pas
avoir de changements entre le moment ou le plan initial est calculé et le moment ou les
décisions doivent être exécutées. L’objectif du modèle de replanification, est de réadapter le
plan de transport des commandes précédemment établi, au vu des changements opérés,
comme suit :
146
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
Paramètres
Date de disponibilité de la commande j (il peut s’agir
𝑛𝑟𝑗 d’une nouvelle commande, ou une commande retardée.
Dans le cas où il n’y a aucun changement 𝑛𝑟𝑗 = 𝑟𝑗 )
Variation maximale entre le premier plan de transport et le
𝑉𝑎𝑟𝑚𝑎𝑥
nouveau
Variable de décision
Nouveau moment auquel le train t charge la commande j à
𝑁𝑅𝑗𝑡𝑠
la station s (après replanification)
L’objectif de la replanification est la minimisation des variations avec le plan initial. Ainsi, la
fonction objectif se présente comme suit :
Dans (5.27), il s’agit de minimiser la somme des variations, entre le plan précédent et celui
résultant de la replanification. L’indice 𝑗 est utilisé pour l’ensemble 𝐽′ , correspondant à
l’ensemble des commandes non transportées et pouvant faire l’objet d’une replanification
(𝐽′ = (𝐽 ∖ 𝐽𝑡 ∖ 𝐽𝑐 ) ∪ 𝐽𝑛 ). L’indice 𝑡 est relatif aux trains de l’ensemble 𝑇 ′ , ne considérant que
ceux qui seront mis en circulation après l’horizon gelé (∀𝑡, 𝑙𝑡 ≥ 𝐻𝐺 ). A noter que pour les
nouvelles commandes (𝑗 ∈ 𝐽𝑛 ), 𝑅𝑗𝑡𝑠 n’existe pas encore. Pour préserver la performance du
système global, à travers la réduction du temps d’attente des commandes dans leur station de
départ avant leur transport, nous considérons : ∑𝑡∈𝑇 ′ 𝑅𝑗𝑡𝑠 = 𝑛𝑟𝑗 , ∀𝑗 ∈ 𝐽𝑛 𝑎𝑣𝑒𝑐 𝑠 = 𝑑𝑒𝑝𝑗 .
L’ensemble des contraintes du PLVM dans le cas prédictif (de (5.5) jusqu’à (5.14)), sont
considérées pour le modèle de replanification. Les modifications à y apporter se présentent
comme suit :
147
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
Dans le but de limiter les variations lors de la replanification, nous définissons les deux
ensembles de contraintes supplémentaires comme suit :
avec 𝐽′′ = 𝐽 ∖ 𝐽𝑡 ∖ 𝐽𝑐 \ 𝐽𝑟
A ces différents paramètres, s’ajoute le nombre d’itérations à considérer, durant une période
d’exploitation (qui peut être une journée, ou une semaine par exemple).
Lors de la première itération, c’est le modèle du cas prédictif qui est utilisé. La solution
proposée est ce qu’on qualifie de plan de transport initial. Toutes les itérations suivantes,
utilisent le modèle de replanification. A noter que le plan initial reste la référence pour les
commandes du plan prédictif. Pour les nouvelles commandes, leur première planification
devient la référence. Enfin, si aucune perturbation n’a été enregistrée durant le pas de
planification actuelle, le recalcul est inutile.
148
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
Considérer la
prochaine itération
Non
Mettre à jour les
données
Occurrence de d’affectation des
perturbations durant commandes aux
l’actuelle itération trains
Oui
149
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
Le couplage mis en œuvre est schématisé par la Figure V.5. On peut le résumer comme suit :
Simulation
Données externes :
Optimisation
• Liste des commandes
• Planning de circulation Modèle
des trains Problème
Initialisation
Données de conception du
modèle de simulation Approche de résolution
Réponse
Plan de transport optimal
Fin
Figure V.5 : Le schéma de couplage dans le cas prédictif
Le couplage mis en œuvre est schématisé par la Figure V.6. En plus de toutes les étapes du
couplage dans le cas prédictif, nous introduisons les étapes supplémentaires suivantes :
150
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
Simulation Optimisation
Données externes :
Modèle Initial
• Liste des commandes Problème initial
• Planning de circulation Instant 0
des trains
Initialisation Approche de résolution
Données de conception du
modèle de simulation
Modèle de replanification
Génération de perturbations:
• Nouvelles commandes Perturbations
• Modifications dans la liste Chaque pas de planification
initiale des commandes Approche de résolution
Fin
Figure V.6 : Le schéma de couplage dans le cas de la replanification
Lors de la circulation des trains à travers la ligne ferroviaire, c’est le plan de transport optimal
actuel, fourni par la dernière itération de résolution du modèle de replanification, qui est
considéré.
8. Conclusion
Dans ce chapitre nous avons développé la méthodologie de résolution du problème FRTSP,
dans les deux cas : prédictif et replanification. Nous avons commencé ce chapitre par l’étude
d’un exemple illustrant les données et contraintes du FRTSP afin d’appréhender sa
complexité. En particulier, à travers cet exemple, nous avons constaté :
151
CHAPITRE V : Modélisation mathématique et optimisation du FRTSP
- L’impact direct des contraintes de temps d’attente dans les stations, à travers le temps
nécessaire au déchargement des commandes, puis, celui nécessaire au chargement
d’autres commandes.
- L’impact de la limitation de l’espace à l’intérieur des trains, sur les possibilités de
transport de chaque train, en fonction des commandes chargées.
- La difficulté d’appréhender les conséquences des différentes possibilités
combinatoires, sur la performance des décisions de transport prises. Ce qui résulte de
la combinaison entre les deux contraintes : temps d’attente et capacité des trains,
associée à la combinaison entre les deux caractéristiques des commandes : parcours
spécifique pour chaque commande et tailles différentes.
Nous avons étudié la complexité du FRTSP en montrant que c’est une extension du problème
d’affectation généralisé. Ensuite, nous avons formalisé ce problème en utilisant un modèle
mathématique de type PLVM, considérant toutes les contraintes techniques et
organisationnelles. L’objectif de cette étude préliminaire est de montrer le bénéfice potentiel,
en quantifiant les volumes qu’il est possible de transporter. Ainsi, nous avons proposé de
minimiser le temps d’attente des commandes avant leur transport. Nous avons adapté les
colonies de fourmis au FRTSP, en raison des bonnes performances de cette métaheuristique,
pour la résolution du problème d’affectation généralisée.
De plus, le transport de marchandises étant sujet à divers types de perturbations, nous avons
proposé une approche par horizon glissant pour replanifier le FRTSP. Cette démarche permet
de considérer périodiquement, toutes les perturbations enregistrées en cours de la période
écoulée, pour générer un nouveau plan de transport optimal. Le modèle d’optimisation
mathématique de replanification, ayant pour objectif, la minimisation des variations relatives,
par rapport au plan de transport initial.
Enfin, nous avons décrit le protocole de couplage simulation / optimisation, dans le but :
152
Chapitre VI
153
CHAPITRE VI : Expérimentations et analyse des résultats
1. Introduction
Ce chapitre est dédié à la phase expérimentale afin de présenter les résultats obtenus des
différents modèles et approches proposées pour le FRTSP. Dans un premier temps, nous
abordons la génération aléatoire des instances du FRTSP, qui est justifiée par l’absence de
données réelles, notamment concernant la demande prévisionnelle pour le transport de
marchandises. Dans la section suivante, nous rapportons les résultats du PLVM, des quatre
heuristiques, ainsi que ceux de l’approche de replanification, qui sont couplés au modèle de
simulation dans le but d’évaluer les taux d’occupation des stations et des trains. Ces résultats
sont cruciaux pour le dimensionnement des espaces et de la capacité nécessaire à la mise en
place de ce nouveau service de transport.
Une section est complètement dédiée aux colonies de fourmis. D’abord, on y présente
l’utilisation de la méthode Taguchi pour définir les paramètres des différentes variantes
étudiées. Ensuite, on présente les résultats obtenus en comparant les performances des
différentes variantes en termes de gap, de temps de calcul et de vitesse de convergence.
2. Génération des instances
L'une des principales difficultés, lors de l'étude d'une nouvelle problématique, est l’absence de
données, telles que la prévision de la demande de transport de marchandises (vu que ce
service n’existe pas encore). Cependant, pour développer un modèle de simulation, des
instances numériques sont nécessaires avec idéalement des données proches des projections
de ce nouveau service. Dans un premier temps, nous proposons de générer des instances
aléatoirement, en se basant sur les données actuelles qui caractérisent le réseau ferroviaire de
transport pour passagers. Entre autres, la fréquence de circulation des trains et leur capacité,
correspondent aux données des lignes franciliennes. Pour les autres paramètres, les valeurs
aléatoires sont générées sur la base des projections obtenues suite aux études du PREDIT.
- Une ligne de transport ferroviaire est composée de 10 stations (𝑆 = 10). Les stations
sont indexées de 1 à 10, en ordre croissant. Ainsi, les trains parcourent la ligne en
traversant les stations de l’index 1, jusqu’à celle avec l’index 10.
- La distance entre deux stations consécutives, requiert un temps de parcours égal à 5
minutes (∀𝑠, 𝑡𝑡𝑠,𝑠+1 = 5).
- Les temps d’attente limites des trains dans chaque station sont données par :
𝑤𝑎𝑖𝑡𝑚𝑖𝑛 = 30 𝑠𝑒𝑐𝑜𝑛𝑑𝑒𝑠 et 𝑤𝑎𝑖𝑡𝑚𝑎𝑥 = 60 𝑠𝑒𝑐𝑜𝑛𝑑𝑒𝑠.
154
CHAPITRE VI : Expérimentations et analyse des résultats
- Les dates de disponibilité des commandes sont générées suivant une loi uniforme,
ainsi : ∀𝑗, 𝑟𝑗 est issue de la distribution 𝑈(0,240), avec 0 correspondant au début de la
période (soit 10h) et 240 (soit 14h), le dernier moment auquel il est possible de
réceptionner une commande pour son transport. Nous arrêtons la réception des
commandes, à 50 minutes de la mise en circulation du dernier train, pour leur
permettre d’être transportées.
- Les commandes sont composées d’un nombre de colis variant entre 1 et 5 colis
standard, générées suivant une loi uniforme, ainsi : ∀𝑗, 𝑄𝑗 est issue avec la loi 𝑈(1,5).
- Les stations de départ et d’arrivée de chaque commande, sont générées aléatoirement
entre la 1ère et la 10ème station. Sachant que la station 1 ne peut pas être une station
d’arrivée et la station 10 ne peut pas être une station de départ.
10 familles d’instances sont considérées, telles que chacune des familles correspond à un
nombre de commandes à transporter, allant de 10 à 100 (avec un pas de 10). Chaque famille
contient 25 instances.
155
CHAPITRE VI : Expérimentations et analyse des résultats
En ce qui concerne le nombre d’itérations simulées, nous allons considérer deux seulement.
L’objectif est d’expliciter la démarche de replanification et d’évaluer la performance du
PLVM de replanification.
Dans un premier temps, nous allons évaluer le modèle d’optimisation dans le cas prédictif.
Puis, nous comparons les résultats des règles décisionnelles implémentées sur ARENA.
Dans un second temps, nous analysons l’impact des différentes stratégies décisionnelles, sur
les différentes composantes du système de transport étudié. Enfin, nous présentons les
résultats de l’approche de replanification.
Pour évaluer les limites du PLVM, en ce qui concerne la considération d’un nombre de
commandes plus important, nous avons testé des instances avec 150 commandes (les autres
paramètres étant fixes). Le PLVM (composé de 49 800 variables et 50 850 contraintes) n’a
pas été en mesure de fournir la solution optimale après 30 minutes de calcul. La décision
d’arrêter le processus de résolution après 30 minutes, découle des spécificités de la solution de
transport étudié. En effet, l’établissement du plan de transport est une tâche qui doit être
renouvelée assez fréquemment (le plus probable, chaque début de journée). Ainsi, il n’est pas
possible d’avoir des temps de calcul très importants.
156
CHAPITRE VI : Expérimentations et analyse des résultats
Temps de
Nombre de Nombre de Nombre de
résolution
commandes variables contraintes
moyen
10 3 600 4 510 0,1
20 6 900 7 820 0,3
30 10 200 11 130 0,4
40 13 500 14 400 0,6
50 16 800 17 750 0,8
60 20 100 21 060 1,2
70 23 400 24 370 1,8
80 26 700 27 680 3,2
90 30 000 30 990 4,2
100 33 300 34 300 17,9
Inst. 1 2 3 4 5 6 7 8 9 10 11 12
Gap% x 3,86 4,79 8,89 4,74 2,89 8,01 9,58 2,25 5,73 x 9,11
13 14 15 16 17 18 19 20 21 22 23 24 25
11,01 13,37 2,91 12,19 5,11 3,17 6,09 2,93 11,18 x 4,32 2,74 8,55
A noter que sur les 25 instances, 3 sont infaisables (instances : 1, 11 et 22). L’étude en détail
de ces instances, montre que la configuration de transport adoptée, ne permet pas la prise en
charge de toutes les commandes. En effet, la date de disponibilité tardive de certaines
commandes, ne leur permet pas d’être transportées par les trains qui restent encore en
circulation, pour deux raisons :
L’un des avantages du PLVM est de vérifier la faisabilité de l’instance très rapidement
(quelques millisecondes). En pratique, cette information peut permettre aux managers, de
prendre des dispositions rapidement, pour ajuster les paramètres du système de transport, dans
le but de satisfaire la demande de transport. Les décisions à prendre peuvent consister à
augmenter le nombre de trains en circulation, ou bien augmenter l’espace dédié aux
marchandises, à l’intérieur des trains planifiés initialement.
Pour évaluer l’influence du nombre de trains, sur la performance du PLVM, nous avons fait
varier leur nombre entre 10 et 100, avec un pas de 10. D’abord, nous avons enregistré un
157
CHAPITRE VI : Expérimentations et analyse des résultats
certain nombre d’instances infaisables (Tableau VI.3). Cette infaisabilité s’explique par les
mêmes raisons citées précédemment.
Tableau VI.3 : Nombre d’instances infaisables
Nombre de trains
Nombre de
10 20 30 40 50 60 70 80 90 100
commandes
10 0 0 0 0 0 0 0 0 0 0
20 0 0 0 0 0 0 0 0 0 0
30 2 0 0 0 0 0 0 0 0 0
40 8 0 0 0 0 0 0 0 0 0
50 18 0 0 0 0 0 0 0 0 0
60 21 0 0 0 0 0 0 0 0 0
70 25 0 0 0 0 0 0 0 0 0
80 25 2 0 0 0 0 0 0 0 0
90 25 2 0 0 0 0 0 0 0 0
100 25 11 0 0 0 0 0 0 0 0
En ce qui concerne l’évolution des temps de résolution, la Figure VI.1 montre leur évolution
en fonction du nombre de trains et du nombre de commandes. Nous relevons que les temps de
résolution augmentent avec l’augmentation du nombre de commandes. Toutefois, pour un
même nombre de commandes, lorsque le nombre de trains est moindre, les temps de
résolution augmentent. En effet, pour les instances faisables, le temps de résolution décroit
sensiblement avec l’augmentation du nombre de trains jusqu’à 40 trains. Ce phénomène est
lié à la réduction du nombre de solutions faisables et à la difficulté de trouver la combinaison
optimale. Puis, à partir de 40 jusqu’à 100 trains, nous observons l’effet inverse, mais avec une
intensité moindre. A noter que le point d’inversion de la tendance se décale vers la droite,
avec l’augmentation du nombre de commandes à transporter.
D’autre part, nous avons enregistré la valeur moyenne de la fonction objectif pour chaque
famille. La Figure VI.2 montre son évolution, en fonction du nombre de trains en circulation.
Nous relevons qu’avec l’augmentation du nombre de trains, le temps d’attente des
commandes se réduit rapidement entre 10 et 50 trains. Ce comportement s’explique par le fait
que l’augmentation du nombre de trains, a un impact positif conséquent sur la performance du
système de transport, jusqu’à l’atteinte d’un point d’équilibre, où la mise en circulation de
trains supplémentaires, n’a qu’un impact réduit. Dans le cas avec 100 commandes à
transporter, ce point d’équilibre se trouve entre 50 et 60 trains. En effet, l’étude en détail des
résultats pour chaque instance avec 100 commandes, montre qu’à partir de 50 trains,
quasiment toutes les commandes sont transportées par le premier train qui les croise. Avec 60
trains, c’est toutes les commandes qui sont prises en charge par le premier train. Au-delà (i.e.
avec un nombre de trains supérieur à 60), la réduction du temps d’attente des commandes
résulte de la réduction de l’intervalle, entre chaque deux trains consécutifs. En effet, cela a
pour conséquence directe, la réduction pour certaines commandes, de leur temps d’attente de
quelques secondes.
158
CHAPITRE VI : Expérimentations et analyse des résultats
70,00
60,00
Nombre de
50,00 trains " T "
temps de résolution en secondes
10
20
40,00 30
Note : aucune instance faisable pour les combinaisons 40
(10 trains, 70-80-90-100 commandes) 50
30,00 60
70
80
20,00
90
100
10,00
-
10 20 30 40 50 60 70 80 90 100
Nombre de commandes " J "
Figure VI.1 : Influence du nombre de trains sur les temps de résolution du PLVM
140 000
Note : aucune instance faisable pour les combinaisons
(10 trains, 70-80-90-100 commandes)
120 000
Nombre de
Valeurs de la fonction objectif (secondes)
80 000 30
40
50
60 000 60
70
80
40 000
90
100
20 000
-
10 20 30 40 50 60 70 80 90 100
Nombre de commandes " J "
Ce résultat apporte une information cruciale, pour la résolution du problème n°6, qui consiste
à déterminer la fréquence optimale des trains en circulation (ou planning de circulation des
trains).
159
CHAPITRE VI : Expérimentations et analyse des résultats
Les histogrammes présentés par la Figure VI.3, montrent la corrélation entre le nombre de
commandes et le temps d’attente total des commandes. Aussi, nous observons la performance
de chaque heuristique, comparativement à la solution optimale. Le Tableau VI.4 donne le gap
entre les solutions fournies par chaque heuristique et les solutions optimales, pour chaque
famille d’instances.
1200
1100
Moyennes des temps d'attente des commandes (secondes)
1000
900
800
PLVM
700 Heuristique 1
Heuristique 2
600 Heuristique 3
Heuristique 4
500
400
300
200
10 20 30 40 50 60 70 80 90 100
Nombre de commandes " J "
Figure VI.3 : Comparaison des temps d’attente des heuristiques avec ceux du PLVM
Les résultats des différentes heuristiques sont assez éloignés de ceux obtenus grâce au PLVM.
De plus, la qualité de leurs solutions se dégrade, d’autant plus que le nombre de commandes à
transporter augmente. Si on définit des classes d’instances telle que le 1er groupe est celui des
20-40 commandes, le 2ème des 50 à 70 et le 3ème celui des 80 à 100, alors la courbe des écarts
est strictement croissante pour chaque heuristique. A noter que l’effort de décomposition
permet à la 4ème heuristique d’obtenir de meilleures solutions comparativement aux 3
premières. Aussi, nous avons relevé que l’écart avec la solution optimale atteint, dans le pire
des cas, 40% sur l’ensemble des instances comprenant celles avec 150 commandes
160
CHAPITRE VI : Expérimentations et analyse des résultats
Des analyses approfondies seront réalisées dans la prochaine section, pour étudier l’impact
respectif de chaque heuristique, sur le système de transport.
Tableau VI.4 : Résultats des heuristiques VS solutions optimales
Nombre de
Heur. 1 Heur. 2 Heur. 3 Heur. 4
commandes
20-40 27% 28% 28% 27%
50-70 31% 33% 30% 31%
80-100 49% 53% 46% 34%
3.3.1. Commandes
Pour caractériser le temps d’attentes des commandes dans le Tableau VI.5, nous rapportons :
- Le temps d’attente moyen des commandes (colonne « Moy. »).
- Le temps d’attente moyen des commandes dont la station de départ est l’une des 3
premières stations de la ligne (colonne « MS »).
- Le temps d’attente moyen des commandes dont la station de départ est l’une des 3
stations au milieu de la ligne (colonne « MM »).
- Le temps d’attente moyen des commandes dont la station de départ est l’une des 3
dernières stations de la ligne - hors dernière station - (colonne « ME »).
Tableau VI.5 : Résultats de la simulation sur les temps d’attente des commandes
La Figure VI.4 montre le mode de calcul des différentes moyennes de temps d’attente.
161
CHAPITRE VI : Expérimentations et analyse des résultats
s=1 s=2 s=3 s=4 s=5 s=6 s=7 s=8 s=9 s=10
Moy. Temps Moy. Temps Moy. Temps Moy. Temps Moy. Temps Moy. Temps Moy. Temps Moy. Temps Moy. Temps
de séjours de séjours de séjours de séjours de séjours de séjours de séjours de séjours de séjours
sta. 1 sta. 2 sta. 3 sta. 4 sta. 5 sta. 6 sta. 7 sta. 8 sta. 9
Moy =
Nous remarquons que l'optimisation du temps d’attente total des commandes, produit des
temps d’attente non équilibrés entre les différentes stations. En effet, le temps d’attente dans
les stations augmente à mesure que la ligne est parcourue. Ce phénomène pourrait s'expliquer,
en partie, par le fait que les trains sont vides au début et se remplissent au fur et à mesure
qu’ils parcourent la ligne. De plus, dans les premières stations, généralement, il n’y a pas de
déchargement, tandis qu'au milieu de la ligne, les deux processus de chargement et de
déchargement, peuvent s’accumuler. Une conséquence directe de cette superposition est la
réduction du temps résiduel pouvant potentiellement permettre des opérations de chargement.
Enfin, dans les dernières stations de la ligne, il y a plus de déchargement que de chargement,
ce qui implique que les commandes qui sont en attente de transport, à partir de ces stations,
mettent plus de temps à être prises en charge.
Une tentative d’équilibrage des temps d’attente des commandes dans les stations, implique la
dégradation du temps d’attente total. A noter qu’en pratique, la majeure partie des
commandes, aurait comme station de départ, les premières stations de la ligne. Ce qui aura
pour conséquence, le rééquilibrage de ces temps d’attente.
D’autre part, il est intéressant de noter la capacité de la 4 ème heuristique à reproduire le même
comportement que le PLVM. Les boîtes à moustaches de gauche de la Figure VI.5 montrent
que la variabilité du temps d’attente des commandes, entre les différentes stations de la ligne,
donnée par l'heuristique 4, est la plus maitrisée. En outre, elles montrent, pour les autres
heuristiques, que cette variabilité est répartie sur un plus large intervalle. Ce constat permet de
souligner l’intérêt de l’heuristique 4, basée sur la décomposition pour les instances difficiles
(voir le max des heuristiques 1, 2 et 3 dans les boîtes à moustaches à gauche de la Figure
VI.5). L'analyse des résultats détaillés dans le Tableau VI.5, montre que pour les trois
premières heuristiques, les temps d’attente des commandes augmentent fortement, lorsque
l'indice de la station augmente. Les boîtes à moustaches à droite de la Figure VI.5 présentent
la variation des temps d’attente des commandes, en fonction du nombre de commandes à
transporter, lorsque la 4ème heuristique est appliquée. Elles montrent que les temps d’attentes
des commandes et leur variabilité entre les stations, augmentent proportionnellement avec
l’augmentation du nombre de commandes à transporter.
162
CHAPITRE VI : Expérimentations et analyse des résultats
1800 1800
1600 1600
1400 1400
1200
1200
1000
1000
800
800
600
600
400
400
200
200
0
10 20 30 40 50 60 70 80 90 100
0
Nombre de commandes - Heuristique 4
PLVM Heur. 1 Heur. 2 Heur. 3 Heur. 4
Figure VI.5 : Boites à moustaches des temps d’attente des commandes dans les stations – (gauche) comparaison entre
les différentes approches, (droite) heuristique 4
Pour rappel, nous avons mentionné la capacité du PLVM, à évaluer rapidement la faisabilité
des instances. Ainsi, le PLVM a fourni la solution optimale, pour les 250 instances résolues,
ce qui signifie que toutes ces instances sont faisables. Toutefois et en ce qui concerne les 4
heuristiques, nous avons enregistré quelques cas, où ces dernières fournissent une solution
incomplète. Plus exactement, à la fin du processus de simulation, certaines commandes n’ont
pas pu être transportées, en raison des règles décisionnelles adoptées et leur impact sur
l’impossibilité technique de les charger en raison d’un manque d’espace ou de temps. Le
Tableau VI.6 fourni le nombre d’instances infaisables et le nombre de commandes non
transportées pour les 250 instances testées.
Là aussi, la 4ème heuristique est la plus intéressante, étant donné qu’un seul cas d’infaisabilité
a été enregistré. A noter que la 3ème heuristique qui permet le transport des commandes de
petite taille en priorité, se classe juste après. Cette dernière, en privilégiant les commandes de
petite taille, permet d’exploiter l’espace à l’intérieur des trains, de manière plus performante
que les 2 premières heuristiques.
Tableau VI.6 : Synthèse des cas d’infaisabilité
3.3.2. Trains
L’évaluation de la performance d’utilisation des trains, pour l’activité de transport de
marchandises, est effectuée à travers le calcul des taux d’utilisation de la capacité dédiée aux
marchandises. Ce taux est relatif à chaque train et, est obtenu par le rapport entre le nombre
total de colis effectivement transportés, et le nombre maximum de colis (noté 𝑈𝐵) pouvant
être transportés par ce dernier. Pour évaluer 𝑈𝐵, nous considérons le pire cas, où un nombre
maximum de colis, doivent être transportées entre chaque deux stations consécutives. A noter
que chaque colis chargé dans un train, doit être déchargé avant que le train ne quitte la ligne,
163
CHAPITRE VI : Expérimentations et analyse des résultats
sachant que les temps de chargement et de déchargement sont égaux. Le pire des cas en
termes de nombre maximum de colis transportés, va correspondre au cas dans lequel chaque
colis a pour station d’arrivée, la station qui suit celle de sa station de départ. En commençant
par un chargement au maximum de la capacité d’accueil et en supposant que tous les colis
sont tous disponibles, cela se ramène au cas illustré dans la Figure VI.6 (i.e. le processus de
chargement / déchargement se répète S⁄2 fois tout au long de la ligne). Ainsi, la formule (6.1)
nous permet de calculer sa valeur, comme suit :
𝑤𝑎𝑖𝑡𝑚𝑎𝑥⁄ 𝑆
𝑈𝐵 = 𝑚𝑖𝑛 (𝑐𝑎𝑝, 𝑡𝑖𝑚𝑒) ∗ ⁄2 (6.1)
Légende :
: Nombre de colis chargés
: Nombre de colis déchargés
: Temps d’attente du train dans la station
s=1 s=2 s=3 s=4 s=5 s=6 s=7 s=8 s=9 s=10
0 6 0 6 0 6 0 6 0 6
Pour les 25 instances de la famille de 100 commandes à transporter, le nombre de colis total
moyen à transporter est égal à 300. Sachant qu’il y a 30 trains en circulation, chacun disposant
d’une capacité maximum UB = 30, nous obtenons un taux d’utilisation moyen de la capacité
des trains égal à 33%. Ceci signifie qu’on aurait pu transporter trois fois plus de colis.
Toutefois et comme indiqué précédemment, selon l’algorithme adopté, plusieurs cas avec des
commandes non transportées, ont été relevés. Sans être en contradiction avec la valeur du taux
d’utilisation des trains, les dates de disponibilité des commandes, combinées à leur taille et
leur parcours (station de départ – station d’arrivée), justifient l’occurrence d’infaisabilités. En
particulier, dans certains cas, les commandes ne sont disponibles que vers la fin du service de
transport de marchandises. Ceci a pour conséquence, la réduction du nombre de possibilités
de transport de ces dernières, voire l’impossibilité de les transporter, alors que des trains en
début de service, ont circulé avec des taux de remplissage très faibles. La résolution de la
problématique décisionnelle n°8, que nous avons identifiée au chapitre 2 (figure II.13),
devrait permettre de réduire l’intensité de ce phénomène, en améliorant la synchronisation des
commandes avec celles des espaces libres dans les trains.
Le temps d’arrêt minimum des trains dans les stations, pour permettre aux voyageurs de
monter et descendre, est fixé à 30 secondes. Pour les instances avec 100 commandes, le temps
d’arrêt moyen des trains dans chaque station est égal à 37 secondes. Ceci implique une
augmentation du temps d’attente de 23%. De plus et pour la même famille, le temps de
parcours moyen des trains, de la première station jusqu’à la dernière, dans le cas où des
marchandises ont été transportées, s’établi à 50,5 minutes. A noter que s’il n’y avait pas
d’activité de transport de marchandises, ce temps s’établira à 49,5 minutes (en considérant un
164
CHAPITRE VI : Expérimentations et analyse des résultats
temps d’arrêt des trains égal à wait min ). Ainsi, la simulation a montré que l’augmentation du
temps de parcours moyen des trains, en raison de l’activité de transport de marchandises, était
de 2% seulement.
3.3.3. Stations
Le Tableau VI.7 présente les nombres moyens / maximums des commandes / colis, présents
simultanément dans les stations, en attente de transport par un train.
Tableau VI.7 : comparaison des taux d’utilisation des stations
L’analyse des résultats présentés dans le Tableau VI.7, nous permet de constater que le
nombre de commandes simultanément présents dans une station, augmente
proportionnellement avec le nombre de commandes total à transporter. Le même constat est
fait pour les colis. Toutefois, en relevant ces résultats en cours du processus de simulation,
nous avons observé que les pics ne durent pas très longtemps, en général quelques secondes
seulement, jusqu’à 5 minutes dans quelques rares cas.
A noter que les trois premières heuristiques sont équivalentes, alors que la quatrième, présente
des résultats légèrement meilleurs.
Cette évaluation de l’espace, notamment des pics atteints dans le temps, devrait apporter des
informations très pertinentes pour le dimensionnement des zones de stockage temporaire.
Ainsi, l’étude de la dynamique de circulation des flux de marchandises, nous permet d’évaluer
les besoins en espace de stockage dans chacune des stations de la ligne, en fonction de la
demande de transport.
3.4. La replanification
Dans le cas de la replanification, la première itération correspond au plan prédictif de
transport. Ainsi, les résultats sont similaires à ceux présentés dans la section 3.1. Le modèle
de replanification est exécuté, après l’écoulement de l’intervalle de recalcul (le pas de
replanification : r = 60 minutes).
Pour chaque famille d’instances, le Tableau VI.8 présente les moyennes des résultats
enregistrés, pour chaque famille d’instances. Nous relevons que le modèle de replanification
165
CHAPITRE VI : Expérimentations et analyse des résultats
présente les mêmes caractéristiques que le modèle du cas prédictif. Ainsi, avec un nombre de
commandes inférieur à planifier, pour chaque famille d’instances, les temps de résolution sont
moins importants.
Tableau VI.8 : Synthèse des résultats de la seconde itération de la replanification
- Itération 1
Le plan de transport initial est présenté par le Tableau VI.9. Pour chaque commande,
nous avons l’identifiant du train qui doit la transporter, ses dates de disponibilité et de
transport (exprimées en secondes), ainsi que son temps d’attente dans la station avant
son transport.
Tableau VI.9 : Résultats de la 1ère itération – plan initial
N° commande 1 2 3 4 5 6 7 8 9 10
N° train 6 9 18 13 20 12 5 23 9 20
Date dispo. 5359 7040 11625 8427 12014 8549 2320 12635 5604 13466
Date transp. 5359 7110 12510 8520 12746 8910 2730 13200 5790 13466
Temps séjour 0 70 885 93 732 361 410 565 186 0
N° commande 11 12 13 14 15 16 17 18 19 20
N° train 8 10 13 17 16 1 4 8 3 11
Date dispo. 4609 7761 8766 10899 8888 530 3145 7077 2129 7554
Date transp. 4917 7761 9180 11250 9330 1320 3145 7077 2190 7650
Temps séjour 308 0 414 351 442 790 0 0 61 96
N° commande 21 22 23 24 25 26 27 28 29 30
N° train 2 21 3 4 19 21 1 4 17 2
Date dispo. 497 12228 2443 3365 12227 13700 0 2017 12203 604
Date transp. 600 12330 2870 3805 12227 13700 1670 2130 12240 1260
Temps séjour 103 102 427 440 0 0 1670 113 37 656
166
CHAPITRE VI : Expérimentations et analyse des résultats
Les différents évènements que nous avons enregistrés, en cours de la première heure
d’exploitation (correspondant au pas de replanification, r = 60 minutes), se
présentent comme suit :
- Itération 2
Les résultats de la replanification sont présentés dans le Tableau VI.11. Les cellules
des commandes en fond rouge sont les commandes dont la date de disponibilité est
retardée, celles en fond bleu, sont les nouvelles commandes. Les données écrites en
rouge, sont celles qui ont subi un changement suite au processus de replanification.
Tableau VI.11 : Résultats de 2ème itération – 1ère replanification
N° commande 2 3 4 5 6 8 9 10 13 14
N° train (initial) 9 18 13 20 12 23 9 20 13 17
N° train (Repl.) 9 17 13 22 12 23 16 21 13 17
Date transp. initial 7110 12510 8520 12746 8910 13200 5790 13466 9180 11250
Date transp. Repl. 7110 11910 8520 13920 8910 13200 9990 14040 9180 11250
Variations - -600 - 1174 - - 4200 574 - -
15 18 20 22 25 26 29 31 32 33 34 35
16 8 11 21 19 21 17 - - - - -
16 8 11 21 19 22 18 21 17 18 15 20
9330 7077 7650 12330 12227 13700 12240 - - - - -
9330 7077 7650 12330 12227 14260 12840 13340 10260 10200 8400 13050
- - - - - 560 600 - - - - -
167
CHAPITRE VI : Expérimentations et analyse des résultats
Les différentes simulations réalisées, nous ont permis de constater la capacité du modèle de
replanification, à insérer, annuler et retarder les commandes, tout en réduisant leur impact par
rapport au plan initial. Aussi, nous avons relevé la réactivité du modèle de replanification
(implémenté sous la forme d’une application Java, faisant intervenir CPLEX pour la
résolution du PLVM, et à laquelle fait appel le modèle ARENA, en cours de simulation). En
effet, il suffit de quelques secondes au plus, pour fournir le nouveau plan de transport, qui est
mis en œuvre 30 minutes après son calcul (en raison de l’horizon gelé).
Les paramètres des colonies de fourmis permettent de moduler le comportement à adopter lors
de la construction des solutions. Suivant la variante utilisée et ses spécificités, il peut être plus
performant d’accentuer son comportement d’intensification ou de diversification, voire à
trouver un équilibre entre les deux. La structure du problème auquel est appliquée la
métaheuristique, a un impact direct sur le choix à faire.
Compte tenu du nombre de paramètres et des valeurs qu’ils peuvent prendre, il n’est pas
envisageable d'analyser l’ensemble des combinaisons possibles, à travers un plan factoriel
complet. Pour déterminer les meilleures valeurs des paramètres de chaque variante, nous
utilisons la méthode Taguchi.
168
CHAPITRE VI : Expérimentations et analyse des résultats
développée pour le secteur industriel, car elle permet de réduire les variabilités des processus
de fabrication qui sont soumis à différents aléas, à travers un paramétrage optimal de ces
derniers. Le principe étant d’identifier les différents facteurs (paramètres) contrôlables, de
définir leurs niveaux, c’est à dire les valeurs qu’ils peuvent prendre. Puis, il faut en déduire la
combinaison la plus robuste pour faire face aux aléas (i.e. un paramétrage qui présente une
sensibilité minimum face aux aléas). Depuis, le champ d’application de cette méthode s’est
élargie. Les différentes notions de cette méthode et ses différentes applications, sont
présentées dans (Taguchi, 1986) et (Roy, 2010).
Il n’est pas évident de développer ce type de tables, toutefois, Taguchi a développé un certain
nombre de tables, qui s’avèrent suffisantes. Ainsi, au besoin, il s’agira de sélectionner la table
qui correspond au mieux aux besoins de l’étude.
Telle que décrit dans (Roy, 2010), la sélection d’une table de Taguchi requiert au préalable,
une étude du processus concerné, à travers :
169
CHAPITRE VI : Expérimentations et analyse des résultats
D’où,
𝑑𝑀 = (2 − 1) + (2 − 1) + (2 − 1)(2 − 1) = 3
Ainsi, pour cet exemple, il faut au moins 3 essais.
- La recherche de la table de Taguchi qui assure au moins le nombre d’essais minimum,
dont le nombre de colonnes est supérieur ou égal au nombre de facteurs et
d’interactions entre facteurs, avec le nombre de niveaux de tous les facteurs.
Un exemple d’une table de Taguchi est donné dans le Tableau VI.12. Dans cet exemple, il
s’agit de trouver la combinaison optimale de 7 facteurs sans interactions, à 2 niveaux chacun.
Dans ce cas dM = 7. Donc le nombre minimum d’essais est égal à 7. Toutefois, comme il y a
une condition sur le fait que les tables de Taguchi doivent être orthogonales, la première table
qui satisfait toutes ces caractéristiques est celle qui compte 8 essais. A noter que le plan
factoriel complet aurait nécessité (27 = 128) essais.
Facteurs contrôlables
N° essai
A B C D E F G
1 1 1 1 1 1 1 1
2 1 1 1 2 2 2 2
3 1 2 2 1 1 2 2
4 1 2 2 2 2 1 1
5 2 1 2 1 2 1 2
6 2 1 2 2 1 2 1
7 2 2 1 1 2 2 1
8 2 2 1 2 1 1 2
La notation d’une table de Taguchi se présente sous la forme LT (nC ) (Sabre, 2007), avec :
Les résultats obtenus à partir des plans d’expériences de Taguchi, sont évalués à travers des
ratios signal/bruit « S/N ». L’avantage de cet indicateur de performance résulte du fait qu’il
considère simultanément, l’objectif recherché (signal), ainsi que la dispersion de cette valeur
(le bruit). Son calcul est donné par la formule suivante (dans le cas d’un critère à
minimiser) (Roy, 2010) :
1
𝑆/𝑁 = −10. log10 [𝑛 ∑𝑛𝑖=1 𝑦𝑖2 ] (6.4)
170
CHAPITRE VI : Expérimentations et analyse des résultats
Enfin, pour chaque paramètre, la moyenne des ratios S/N pour chaque niveau est calculée.
Puis, le ratio le plus élevé va déterminer le meilleur niveau. A noter que la différence entre les
S/N le plus élevé et le plus faible, pour chaque paramètre, nous permettra de définir « 𝑑𝑒𝑙𝑡𝑎 »
qui représente l’influence du paramètre sur la performance de la variante. Plus 𝑑𝑒𝑙𝑡𝑎 sera
faible, moins le paramètre correspondant aura de l’influence.
Pour les autres paramètres et selon les différentes possibilités suggérées dans la littérature,
pour chacune des variantes, nous proposons ce qui suit :
171
CHAPITRE VI : Expérimentations et analyse des résultats
A noter que dans notre cas, nous recherchons le paramétrage optimal d’algorithmes, ainsi, les
variations au niveau de chacun des paramètres (ou facteurs) définis sont faciles. D’autre part,
selon la littérature, il n’y a pas d’interactions entre les différents paramètres.
En ce qui concerne les variantes AS et MMAS, le modèle qui traduit la relation entre les
différents paramètres, s’écrit comme suit :
𝑦 =𝛽+𝜌+𝑄 (6.5)
Le degré de ce modèle est donné par :
𝑑𝐴𝑆/𝑀𝑀𝐴𝑆 = (4 − 1) + (4 − 1) + (4 − 1) = 9
Ainsi, il faut au moins 9 essais. La recherche de la table de Taguchi la plus adéquate, nous a
permis d’identifier la table L16 (45 ) qui est donnée par le Tableau VI.15. Cette dernière
permet d’étudier le cas avec au plus 5 facteurs à 4 niveaux. Dans notre cas, nous avons 3
facteurs seulement, donc les deux dernières colonnes de la table, ne seront pas considérées.
Cette table propose d’effectuer 16 essais pour trouver la combinaison optimale, alors que le
plan factoriel complet requiert (43 = 64) essais.
172
CHAPITRE VI : Expérimentations et analyse des résultats
En ce qui concerne les deux variantes ACS et MMACS, les modèles qui traduisent la relation
entre les différents paramètres, s’écrivent respectivement comme suit :
𝑦 = 𝛼 + 𝛽 + 𝜌 + 𝑄 + 𝑞0 + 𝜑 (6.6)
𝑦 = 𝛼 + 𝛽 + 𝜌 + 𝑄 + 𝑞0 (6.7)
Les degrés de ces deux modèles sont donnés par :
𝑑𝐴𝐶𝑆 = (5 − 1) + (5 − 1) + (5 − 1) + (5 − 1) + (5 − 1) + (5 − 1) = 24
𝑑𝐴𝑀𝑀𝐴𝐶𝑆 = (5 − 1) + (5 − 1) + (5 − 1) + (5 − 1) + (5 − 1) = 20
Ainsi, il faut au moins 24 essais, pour la variante ACS et 20 essais, pour la variante MMACS.
La recherche de la table de Taguchi la plus adéquate pour ces deux cas, nous a permis
d’identifier la table L25 (56 ) (elle est donnée par le Tableau VI.16). Elle permet la
considération au plus de 6 facteurs à 5 niveaux. Elle propose d’effectuer 25 essais, au lieu de
(56 = 15 625) essais pour la variante ACS et (55 = 3 125) essais pour la variante MMACS.
A noter que pour la variante MMACS, la dernière colonne n’est pas considérée.
Tableau VI.16 : Table de Taguchi 𝑳𝟐𝟓 (𝟓𝟔 )
173
CHAPITRE VI : Expérimentations et analyse des résultats
- Variante AS
Le Tableau VI.17 présente les résultats des essais de la variante AS. La colonne « Moy »
donne la moyenne des gap d’optimalité des différentes mesures de chaque essai.
Tableau VI.17 : Résultats des essais pour la variante AS
N° N°
𝛽 𝜌 𝑄 Moy S/N 𝛽 𝜌 𝑄 Moy S/N
essai essai
1 2 0,25 0,1 41,26% 7,690 9 4 0,25 10 15,33% 16,292
2 2 0,5 1 42,15% 7,504 10 4 0,5 100 14,79% 16,601
3 2 0,75 10 44,01% 7,129 11 4 0,75 0,1 15,27% 16,325
4 2 0,98 100 44,17% 7,097 12 4 0,98 1 15,48% 16,205
5 3 0,25 1 22,29% 13,037 13 5 0,25 100 12,43% 17,110
6 3 0,5 0,1 20,71% 13,677 14 5 0,5 10 12,18% 18,286
7 3 0,75 100 19,92% 14,015 15 5 0,75 1 12,50% 18,393
8 3 0,98 10 20,55% 13,742 16 5 0,98 0,1 12,48% 18,076
Le Tableau VI.18 présente la moyenne des ratios S/N pour chaque niveau de chaque
paramètre, ainsi que le delta. En vert, le niveau des paramètres optimaux.
Tableau VI.18 : Ratio S/N pour la variante AS
Niveau 𝛽 𝜌 𝑄
1 7,355 55,129 55,768
2 13,618 56,068 55,139
3 16,356 55,862 55,448
4 18,216 55,121 55,824
𝑑𝑒𝑙𝑡𝑎 10,861 0,947 0,685
𝛼 = 1 ; 𝛽 = 5 ; 𝜌 = 0,5 ; 𝑄 = 100
- Variante MMAS
174
CHAPITRE VI : Expérimentations et analyse des résultats
N° N°
𝛽 𝜌 𝑄 Moy S/N 𝛽 𝜌 𝑄 Moy S/N
essai essai
1 2 0,25 0,1 11,35% 18,898 9 4 0,25 10 7,24% 22,811
2 2 0,5 1 4,36% 27,367 10 4 0,5 100 7,81% 22,151
3 2 0,75 10 4,76% 26,452 11 4 0,75 0,1 9,91% 20,078
4 2 0,98 100 25,25% 11,956 12 4 0,98 1 12,65% 17,957
5 3 0,25 1 6,37% 23,913 13 5 0,25 100 7,79% 22,164
6 3 0,5 0,1 10,17% 19,850 14 5 0,5 10 8,20% 21,720
7 3 0,75 100 7,76% 22,202 15 5 0,75 1 8,65% 21,261
8 3 0,98 10 15,22% 16,354 16 5 0,98 0,1 11,74% 18,607
Le Tableau VI.20 présente la moyenne des ratios S/N pour chaque niveau de chaque
paramètre, ainsi que le delta. En vert, le niveau des paramètres optimaux.
Tableau VI.20 : Ratio S/N pour la variante MMAS
Niveau 𝛽 𝜌 𝑄
1 21,168 87,786 77,434
2 20,580 91,088 90,498
3 20,749 89,994 87,337
4 20,938 64,875 78,474
𝑑𝑒𝑙𝑡𝑎 0,588 26,213 13,064
𝛼 = 1 ; 𝛽 = 2 ; 𝜌 = 0,5 ; 𝑄 = 1
- Variante ACS
175
CHAPITRE VI : Expérimentations et analyse des résultats
Le Tableau VI.22 présente la moyenne des ratios S/N pour chaque niveau de chaque
paramètre, ainsi que le delta. En vert, le niveau des paramètres optimaux.
Tableau VI.22 : Ratio S/N pour la variante ACS
Niveau 𝛼 𝛽 𝜌 𝑄 𝑞0 𝜑
1 22,147 16,795 20,979 20,962 22,327 22,250
2 22,344 19,411 21,129 21,457 25,127 20,531
3 22,573 24,105 21,197 21,244 20,218 20,218
4 20,590 22,770 21,447 21,447 20,002 19,435
5 20,460 25,031 20,730 20,655 20,440 23,933
𝑑𝑒𝑙𝑡𝑎 2,113 8,236 0,717 0,802 5,125 4,498
- Variante MMACS
176
CHAPITRE VI : Expérimentations et analyse des résultats
Le Tableau VI.24 présente la moyenne des ratios S/N pour chaque niveau de chaque
paramètre, ainsi que le delta. En vert, le niveau des paramètres optimaux.
Tableau VI.24 : Ratio S/N pour la variante MMACS
Niveau 𝛼 𝛽 𝜌 𝑄 𝑞0
1 20,133 19,640 21,545 21,492 23,371
2 21,792 20,034 21,808 22,181 25,331
3 22,727 23,832 22,103 21,243 22,123
4 21,935 22,602 22,582 22,582 19,967
5 24,680 25,160 21,981 21,883 20,475
𝑑𝑒𝑙𝑡𝑎 4,547 5,52 1,037 1,401 5,364
Même chose que la variante ACS, nous notons que les influences de ρ et Q sont les plus
faibles.
177
CHAPITRE VI : Expérimentations et analyse des résultats
4.2.5. Validation
Le paramétrage optimal de chaque variante est le résultat d’une étude théorique, qui s’appuie
principalement, sur l’hypothèse d’addition des effets moyens. Ainsi, il est indispensable de
valider ce résultat théorique, en effectuant le même nombre de mesures pour chaque essai des
plans d’expériences de Taguchi. A l’issue de cette opération, il convient de vérifier que ce
paramétrage, permet effectivement de produire la meilleure solution. Dans certains cas, il peut
s’avérer nécessaire de réajuster certains paramètres.
Dans notre cas et pour chacune des quatre variantes, cette étape de validation a permis la
confirmation de la performance des paramètres sélectionnés. Dans la section suivante, nous
présenterons les résultats des quatre variantes de colonies de fourmis, avec les paramètres
fixés précédemment.
D’abord, nous relevons que la variante de base AS, est beaucoup moins performante que les
trois autres variantes. Aussi, jusqu’à 50-60 commandes à transporter, les quatre variantes sont
plus ou moins équivalentes. Cependant, la performance de la variante MMAS se dégrade,
avec l’augmentation du nombre de commandes. Pour les deux variantes ACS et MMACS,
leur performance est relativement stable, indépendamment du nombre de commandes.
Toutefois, pour les familles d’instances avec 90 et 100 commandes, nous avons relevé deux
instances difficiles à résoudre, qui dégradent la moyenne des écarts observés.
En ce qui concerne les temps de résolution, en secondes dans le Tableau VI.25, nous pouvons
observer qu’ils augmentent suivant une tendance linéaire, avec l’augmentation du nombre de
commandes. Toutefois, cette augmentation n’est pas liée uniquement au nombre de
commandes, mais aussi au nombre de fourmis de la colonie. Pour rappel, nous avons décidé
de lier le nombre de fourmis au nombre de commandes. D’autre part, nous remarquons pour
les familles d’instances avec 100 commandes, que la variante MMAS a les temps de
résolution les plus faibles, à contrario, la variante ACS a les temps de résolution les plus
importants. Pour une meilleure étude de dominance au vu des deux critères (gap et temps de
résolution), la Figure VI.8 confronte les temps de résolution et les résultats (gap). Nous
relevons que les deux variantes AS et ACS sont dominées par MMAS et MMACS. Par
178
CHAPITRE VI : Expérimentations et analyse des résultats
ailleurs, on peut départager ces deux dernières au vu du gap obtenu avec la variante MMACS
pour moins de 2 secondes de temps de calcul supplémentaires. Ainsi, on peut dire que la
variante MMACS présente de meilleures performances.
900
Moyennes des temps d'attente des commandes (secondes)
800
700
600
PLVM
AS
MMAS
500
ACS
MMACS
400
300
200
10 20 30 40 50 60 70 80 90 100
Nombre de commandes " J "
Figure VI.7 : Moyennes des temps d’attente des commandes – Colonies de fourmis VS PLVM
179
CHAPITRE VI : Expérimentations et analyse des résultats
12,00%
10,00%
8,00%
AS
Gap%
6,00%
MMAS
4,00% ACS
MMACS
2,00%
0,00%
0 2 4 6 8 10 12
Temps de résolution (secondes)
Comme mentionné précédemment, le principe des colonies de fourmis est d’exploiter les
résultats de chaque itération, pour guider la recherche des itérations suivantes. Ce principe est
mis en œuvre, à travers le mécanisme de dépôt de quantités de phéromones, proportionnelles
à la qualité de la solution produite par chaque fourmi (ou la meilleure fourmi, suivant la
variante). Les différentes variantes permettent une exploitation différente de ce mécanisme.
L’objectif étant d’accentuer son apport dans le processus de recherche de la solution optimale,
ou bien, de l’atténuer, pour une exploration plus large du champ des possibilités. La Figure
VI.9 montre la vitesse de convergence de chaque variante. Pour cela, nous avons considéré
une instance particulière avec 100 commandes.
180
CHAPITRE VI : Expérimentations et analyse des résultats
- La variante MMACS combine les avantages des deux variantes précédentes, à savoir :
elle permet une convergence assez rapide, et grâce au mécanisme de réinitialisation de
la trace de phéromones, elle arrive à quitter les optimaux locaux, pour chercher
d’autres solutions de meilleure qualité. En effet, le résultat final le confirme, puisqu’on
est à 0,09% de la solution optimale.
- Après 100 itérations, les variantes AS et ACS proposent des solutions équivalentes.
Cependant, 200 itérations plus tard, la variante ACS a réussi à améliorer la solution de
13%.
Gap :
17,51%
6,97%
4,58%
0,09%
A noter que d’autres instances ont été testées, avec des nombres de commandes différents. Le
comportement des différentes variantes a été le même que celui décrit ci-dessus. Enfin, nous
pouvons conclure que la variante MMACS, présente la meilleure performance (qualité de la
solution / temps de résolution), au vu des spécificités de la problématique étudiée.
5. Conclusion
Dans ce chapitre, nous avons étudié les résultats numériques des différentes approches
proposées, pour le FRTSP. Pour cela, nous avons généré un certain nombre d’instances
aléatoirement de différentes tailles.
Les résultats numériques ont montré que le PLVM, dans le cas prédictif, est en mesure de
résoudre rapidement des instances de taille modérée, avec un nombre de commandes allant
jusqu’à 100. En ce qui concerne les instances de grande taille, ce PLVM a montré ses limites
pour fournir la solution optimale, en un temps raisonnable pour un problème de niveau
opérationnel.
181
CHAPITRE VI : Expérimentations et analyse des résultats
D’autre part, les règles de priorité, qui sont souvent privilégiées en milieu industriel, pour la
facilité de leur mise en œuvre, présentent des résultats assez médiocres dans le cas du système
de transport étudié. De plus, la qualité des solutions qu’elles fournissent, se dégrade très
rapidement, avec l’augmentation du nombre de commandes. A noter que la décomposition du
FRTSP, en plusieurs sous-problèmes mono-train, apporte les meilleurs résultats. Mais même
ces derniers, ne sont pas à la hauteur des efforts de calcul.
Les résultats numériques de l’approche de replanification, ont confirmé que ce dernier est en
mesure de fournir rapidement, un plan de transport réajusté, qui considère les différentes
perturbations, en minimisant l’écart par rapport au calcul initial. A noter que le point de départ
de la replanification, étant le plan de transport initial, il est possible de générer ce dernier en
utilisant l’algorithme des colonies de fourmis, notamment pour les instances de grande taille.
Concernant les colonies de fourmis, nous avons adapté différentes variantes de ces dernières.
Nous avons fait appel à la méthode Taguchi qui garantit la qualité du paramétrage des
colonies de fourmis, en déterminant les meilleures combinaisons pour ces algorithmes.
L’étude expérimentale et l’analyse des résultats, mettent en évidence la dominance de la
variante MMACS au vu du gap obtenu et du temps de résolution.
182
Conclusion générale et perspectives
183
Conclusion générale et perspectives
Dans cette thèse, nous avons développé une méthodologie pour modéliser et optimiser, une
nouvelle solution de transport de marchandises en milieu urbain, basée sur la mutualisation du
réseau de transport ferroviaire urbain, pour le transport des marchandises et des voyageurs.
D’abord, nous avons montré la nécessité de déployer de nouvelles solutions de transport de
marchandises en milieu urbain. Ces dernières doivent être en mesure d’absorber une partie
des mouvements de marchandises, tout en réduisant l’impact de cette activité sur la qualité de
vie en ville (moins de pollution et moins de circulation routière). Dans ce contexte,
l’utilisation du réseau ferroviaire, actuellement exploité pour le transport de voyageurs, est
une idée qui est en mesure de répondre favorablement à ces attentes, d’une part, d’autre part,
elle permet d’améliorer le rendement des équipements et infrastructures de ce mode de
transport durant les heures creuses. L’état de l’art dans ce domaine permet de conforter cette
idée, puisque plusieurs études ont mis en évidence le potentiel et la faisabilité d’une telle
solution de mixité, pour intégrer le fret urbain au flux de passagers.
Dans le chapitre 2, nous avons proposé une identification et une classification des 19
différentes configurations de mixité, permettant de définir le niveau de difficulté de mise en
œuvre correspondant. Cette classification se base sur le degré de partage des ressources du
système de transport. Dans une perspective de déploiement de cette solution dans d’autres
villes, il serait pertinent de se doter d’une méthodologie qui permettrait d’identifier le cas de
mixité à adopter. Cette identification devrait s’appuyer sur les caractéristiques de la ville
d’implantation, le type de marchandises à transporter et particulièrement, les caractéristiques
du réseau ferroviaire.
Dans cette thèse, nous avons souhaité montrer la faisabilité de cette solution, c’est pourquoi
nous avons retenu le cas de mixité le plus critique, à savoir le degré maximal de partage des
ressources. Une seconde contribution a consisté en une approche par décomposition
temporelle pour lier l’ensemble des 9 problématiques décisionnelles (stratégique, tactique et
opérationnelle). Une seconde perspective serait de généraliser l’approche de décomposition, à
184
Conclusion générale et perspectives
l’ensemble des cas de mixité possibles, en identifiant les problématiques communes et celles
qui sont spécifiques pour chacun des cas de mixité.
Nous avons retenu le problème FRTSP, car si ce dernier ne serait pas capable d’absorber un
flux supplémentaire, qu’est celui des marchandises, la solution n’aurait aucune chance d’être
adoptée. De plus, si le modèle économique de la nouvelle offre de service venait à se préciser,
l’étude de ce problème permettrait à l’opérateur de transport d’effectuer une première
évaluation de ses gains potentiels. Ceci pourrait se faire de manière assez immédiate, pour
passer d’une quantification des volumes pouvant être transporté à une évaluation financière
du flux correspondant. Par ailleurs, cet outil de simulation pourrait également être adapté,
pour dimensionner la capacité nécessaire pour maximiser les bénéfices.
Dans le chapitre 4, nous avons formalisé le problème FRTSP, en définissant ses différentes
hypothèses. Ensuite, nous avons proposé un modèle de simulation à évènements discrets, qui
se justifie par l’étude du comportement dynamique de cette solution, non encore déployée à
l’échelle réelle. Ce modèle est enrichi grâce à l’implémentation de plusieurs règles
décisionnelles. La validation de ce modèle de simulation a été faite, en s’appuyant sur le
retour d’experts du transport ferroviaire en milieu urbain. Par la suite, nous avons formulé un
modèle mathématique de type PLVM qui a été validé lui-même, grâce au modèle de
simulation.
L’analyse de l’impact du plan de transport optimal, nous a permis de relever que l’utilisation
de l’espace dans les stations, bien que corrélée au nombre de commandes, est relativement
maîtrisée. Ceci est une conséquence directe de l’optimisation du taux de rotation, ce qui
permet d’accélérer la dynamique, en prenant en charge le plus rapidement possible les
commandes, dès leur réception. Toutefois, cette conclusion doit être modérée par l’occurrence
d’incidents majeurs, qui risquent de créer un tassement des commandes, si aucune action n’est
entreprise. D’où une troisième perspective de ce travail, la considération de perturbations sur
la ligne ferroviaire, telles que l’annulation de trains, l’interruption du service durant un laps de
temps assez important, l’impact de ces incidents sur l’accroissement temporaire de la
demande de transport de voyageurs, etc.
Cette thèse ouvre la voie à plusieurs autres perspectives très intéressantes. Par exemple, une
quatrième perspective consisterait en la généralisation du FRTSP au cas multi-ligne avec
rupture de charge (correspondance). La pertinence de cette perspective relève de sa capacité à
accroitre les taux de couverture du réseau de transport de fret. La généralisation au cas
multimodal s’inscrit également dans le même objectif de l’amélioration du taux de couverture,
185
Conclusion générale et perspectives
186
Bibliographie
Bibliographie
187
Bibliographie
Akkan, C., 2015. Improving schedule stability in single-machine rescheduling for new
operation insertion. Computers&OperationsResearch, Volume 64, pp. 198-209.
Almeder, C., Preusser, M. & Hartl, R., 2009. Simulation and optimization of supply chains:
alternative or complementary approaches?. OR spectrum, 31(1), pp. 95-119.
Altazin, E., Dauzère-Pérès, S., Ramond, F. & Tréfond, S., 2017. Rescheduling through stop-
skipping in dense railway systems. Transportation Research Part C: Emerging Technologies,
Volume 79, pp. 73-84.
Altiok, T. & Melamed, B., 2007. Simulation Modeling and Analysis with ARENA. Burlington:
Academic Press.
Alumur , S., Kara , B. & Melo , M., 2015. Location and Logistics. In: Laporte G., Nickel S.,
Saldanha da Gama F. (eds) Location Science. Springer, Cham.
Amaran, S., Sahinidis, N., Sharda, B. & Bury, S., 2016. Simulation optimization: a review of
algorithms and applications. Annals of Operations Research, Volume 240, pp. 351-380.
Amini, M. & Racer, M., 1994. A rigorous computational comparison of alternative solution
methods for the generalized assignment problem. Management Science, 40(7), pp. 868-890.
Anu, M., 1997. Introduction to modeling and simulation. Proceedings of the 1997 Winter
Simulation Conference.
Bakalem, M., 1996. Modélisation et simulation orientées objet des systèmes manufacturiers.
Thèse de doctorat en Électronique - Électrotechnique - Automatique, Université de Savoie,
France.
Baptiste, P., Le Pape, C. & Nuijten, W., 2003. Constraint-based scheduling applying
constraint programming to scheduling problems. London: Kluwer Academic Publisher.
Bellman, R., 1956. Dynamic programming and lagrange multipliers. Proceedings of the
National Academy of Sciences of the United States of America, 42(10), pp. 767-769.
Bellman, R. & Dreyfuce, S., 1962. Applied dynamic programming. Princeton: Princeton
university press.
Ben-Khedher, N., Kintanar, J., Queille, C. & Stripling, W., 1998. Schedule optimization at
SNCF: from conception to day of departure. INTERFACES, 28(1), pp. 6-23.
Bertier, P. & Roy, B., 1964. Une Procédure de résolution pour une classe de problèmes
pouvant avoir un caractère combinatoire. Cahiers du Centre d’Études de recherche
Opérationnelle, Volume 6, pp. 202-208.
188
Bibliographie
Bhattacharya, A., Kumar, S., Tiwari, M. & Talluri, S., 2014. An intermodal freight transport
system for optimal supply chain. Transportation Research Part C: Emerging Technologies,
Volume 38, pp. 73-84.
Bierlaire, M., 2015. Simulation and optimization: A short review. Transportation Research
Part C: Emerging Technologies, Volume 55, pp. 4-13.
Borodin, V., Bourtembourg, J., Hnaien, F. & Labadie, N., 2017. COTS software integration
for simulation optimization coupling: Case of ARENA and CPLEX products. Preprint
submitted to Mathematics and Computers in Simulation.
Boudouin, D., 2006. Guide de mise en oeuvre d’espaces logistiques urbains, s.l.: CRET-Log,
Action fédérative ELU, PREDIT.
Boudouin, D. & Morel, C., 2002. L'optimisation de la circulation des biens et services en
ville: logistique urbaine : Programme national "marchandises en ville". s.l.:Documentation
française.
Brice, D., Marinov, M. & Rüger, B., 2015. A newly designed baggage transfer system
implemented using event-based simulations. Urban Rail Transit, 1(4), pp. 194-214.
Browne, M., Allen, J., Woodburn, A. & Piotrowsk, M., 2014. The potential for non-road
modes to support environmentally friendly urban logistics. 1st International Conference
Green Cities 2014 – Green Logistics for Greener Cities, London, UK, pp. 29-36.
Browne, M., Sweet, M., Woodburn, A. & Allen, J., 2005. Urban freight consolidation : Final
Report, University of Westminster for the Department for Transport: Transport Studies
Group.
Cacchiani, V. et al., 2014. An overview of recovery models and algorithms for real-time
railway rescheduling. Transportation Research Part B: Methodological, Volume 63, pp. 15-
37.
Carosi, S., Gualandi, S., Malucelli, F. & Tresoldi, E., 2015. Delay management in public
transportation: service regularity issues and crew re-scheduling. Transportation Research
Procedia: 18th Euro Working Group on Transportation, Delft, The Netherlands, Volume 10,
pp. 483-492.
Carrie, A., 1988. Simulation of manufacturing systems. Chichester: John Wiley & Sons.
Carson, Y. & Maria, A., 1997. Simulation optimization: methods and applications.
Proceedings of the 1997 Winter Simulation Conference. Atlanta, Georgia: IEEE, pp. 118-126.
Cattrysse, D. & Van Wassenhove, L., 1992. A survey of algorithms for the generalized
assignment problem. European Journal of Operational Research, 10(3), pp. 260-272.
189
Bibliographie
Cauvin, A., Ferrarini, A. & Tranvouez, E., 2009. Disruption management in distributed
enterprises: A multiagent modelling and simulation of cooperative recovery behaviors.
International Journal of Production Economics, Volume 122, pp. 429-439.
Chand, S., Hsu, V. & Sethi, S., 2002. Forecast, solution and rolling-horizons in operations
management problems: A classified bibliography. Manufacturing and Service Operations
Management, 4(1), pp. 25-43.
Christofides, N., Mingozzi, A. & Toth, P., 1981. Exact algorithms for the vehicle routing
problem, based on spanning tree and shortest path relaxations. Mathematical Programming,
Volume 20, pp. 255-282.
Chu, P. & Beasley, J., 1997. A genetic algorithm for the generalised assignment problem.
Computers & Operations Researc, 24(1), pp. 17-23.
Chvatal, V., 1983. Linear programming. W. H. Freeman and Company ed. New York:
Macmillan .
Comi, A. & Nuzzolo, A., 2016. Exploring the Relationships between E-Shopping Attitudes
and Urban Freight Transport. Transportation Research Procedia, Volume 12, pp. 399-412.
Corman, F., D’Ariano, A., Pacciarelli, D. & Pranzo, M., 2012. Bi-objective conflict detection
and resolution in railway traffic management. Transportation Research Part C: Emerging
Technologies, 20(1), pp. 79-94.
CST, 2016. Cargo sous terrain: The goods transport system of the future. Press release 26th
January 2016. [Online] Available at: http://www.cargosousterrain.ch [Accessed 28 07 2017].
Dablanc, L., 2007. Goods transport in large European cities: difficult to organize, difficult to
modernize. Transportation Research Part A, 41(3), pp. 280-285.
Dablanc, L., 2007. Goods transport in large European cities: difficult to organize, difficult to
modernize. Transportation Research Part A: Policy and Practice, 41(3), pp. 280-285.
Dablanc, L., 2009. Freight transport for development toolkit: Urban freight. Department for
International Development, The World Bank.
190
Bibliographie
Dakin, R., 1965. A Tree Search Algorithm for Mixed Integer Programming Problems.
Computer Journal, 8(3), pp. 250-255.
Dampier, A. & Marinov, M., 2015. A study of the feasibility and potential implementation of
metro-based freight transportation in Newcastle upon Tyne. Urban Rail Transit, 3(1), pp.
164-182.
Dantzig, G., 1951. Activity analysis of production and allocation - Chapter XXI.
Maximization of a linear function of variables subject to linear inequalities. New York: John
Wiley & Sons.
Dantzig, G. & Ramser, J., 1959. The Truck Dispatching Problem. Management Science , 6(1),
pp. 80-91.
De Langhe, K., 2014. Analysing the role of rail in urban freight distribution. Next Generation
Supply Chains: Trends and Opportunities, epubli GmbH, Berlin.
Delaitre, L., 2014. Introduction à la logistique urbaine - Entre intérêt public et optimisation
industrielle. Logistic Systems & Software, SAVOYE.
Dell'Amico, M., Maffioli, F. & Martello, S., 1997. Annotated bibliographies in combinatorial
optimization. s.l.:John Wiley & Sons.
Delmastro, C., Lavagno, E. & Schranz, L., 2016. Underground urbanism: Master plans and
sectorial plans. Tunnelling and Underground Space Technology, Volume 55, pp. 103-111.
Denardo, E. & Mitten, L., 1967. Elements of sequential decision processes. Journal of
industrial engineering, Volume 18, pp. 106-112.
Dessouky, M. & Leachman, R., 1995. A simulation modeling methodology for analyzing
large complex rail network. Simulation, 65(2), pp. 131-142.
Dıaz, J. & Fernández, E., 2001. A tabu search heuristic for the generalized assignment
problem. European Journal of Operational Research, 132(1), pp. 22-38.
Diziain, D., Taniguchi, E. & Dablanc, L., 2014. Urban logistics by rail and waterways in
France and Japan. 8th International Conference on City Logistics, Procedia - Social and
Behavioral Sciences, Volume 125, pp. 159-170.
Dollevoet, T. & Huisman, D., 2014. Fast heuristics for delay management with passenger
rerouting. Public Transport, 6(1-2), pp. 67-84.
Dorigo, M., Birattari, M. & Stutzle, T., 2006. Ant colony optimization. IEEE computational
intelligence magazine, 1(4), pp. 28-39.
Dorigo, M. & Gambardella, L., 1997. nt colony system: a cooperative learning approach to
the traveling salesman problem. IEEE Transactions on evolutionary computation, 1(1), pp.
53-66.
191
Bibliographie
Dorigo, M., Maniezzo, V. & Colorni, A., 1991. Positive feedback as a search strategy.
Technical report 91-016, Dipartimento di Elettronica, Politecnico di Milano, Milan.
Dorigo, M. & Stützle, T., 2004. Ant Colony Optimization. A Bradford Book ed. Cambridge,
Massachusetts: The MIT Press.
Dufour, J., Patier, D. & Routhier, J., 2007. Du transport de marchandises en ville à la
logistique urbaine. Techniques de l'ingénieur, p. ag8210.
Eberhart, R. & Kennedy, J., 1995. A new optimizer using particle swarm theory. Proceedings
of the sixth international symposium on micro machine and human science, Volume 1, pp. 39-
43.
European Commission, 2011. White Paper - Roadmap to a single European transport area -
Towards a competitive and resource efficient transport system, COM/2011/0144 final, s.l.:
European Commission publishing.
Feo, T., Resende, M. & Smith, S., 1994. A greedy randomized adaptive search procedure for
maximum independent set. Operations Research, 42(5), pp. 860-878.
Fisher, M., Jaikumar, R. & Van Wassenhove, L., 1986. A multiplier adjustment method for
the generalized assignment problem. Management Science, 32(9), pp. 1095-1103.
Flodén, J., 2007. Modelling intermodal freight transport: the potential of combined transport
in Sweden, s.l.: PhD Thesis, Goteborg University, BAS Publisher.
Flodén, J. & Williamsson, J., 2016. Business models for sustainable biofuel transport: the
potential for intermodal transport. Journal of Cleaner Production, Volume 113, pp. 426-437.
Fontan, G., Mercé, C. & Lasserre, J., 2001. La planification des flux de production. In:
Performance Industrielle et Gestion des flux (P. Burlat, J.P. Campagne). s.l.:Hermers Science
Publications Traité IC2.
Fu, M., 1994. Optimization via simulation: a review. Annals of Operations Research, 53(1),
pp. 199-204.
Fu, M., 2002. Optimization for simulation: theory vs. practice. INFORMS Journal on
Computing, 14(3), pp. 192-215.
Galasso, F., 2007. Aide à la planification dans les chaînes logistiques en présence de demande
flexible. Thèse de Doctorat, Institut National Polytechnique de Toulouse.
192
Bibliographie
Glover , F., 1986. Future paths for integer programming and links to artificial intelligence.
Computers and Operations Research, 13(5), pp. 533-549.
Gonzalez-Feliu, J., 2014. Costs and benefits of railway urban logistics: a prospective social
cost benefit analysis. <halshs-01056135>.
Gonzalez-Feliu, J., Malhéné, N., Morganti, E. & Trentini, A., 2013. Développement des
espaces logistiques urbains. CDU et ELP dans l'europe du sud-ouest, s.l.: halshs-00862009.
Guignard, M. & Rosenwein, M., 1989. An improved dual-based algorithm to the knapsack
problem. Operations Research, 37(4), pp. 658-663.
Guo, Y., Huang, M., Wang, Q. & Leon, V., 2016. Single-machine rework rescheduling to
minimize maximum waiting-times with fixed sequence of jobs and ready times. Computers &
Industrial Engineering, Volume 91, pp. 262-273.
Gupta, A. & Könemann, J., 2011. Approximation algorithms for network design: a survey.
Surveys in Operations Research and Management Science, 16(1), pp. 3-20.
Haralick, R. & Elliott, G., 1980. Increasing tree search efficiency for constraint satisfaction
problems. Artificial intelligence, 14(3), pp. 263-313.
Holland, J., 1975. Adaptation in natural and artificial systems : an introductory analysis with
applications to biology, control, and artificial intelligence. Michigan : University of Michigan
Press.
Huisman, D., Kroon, L., Lentink, R. & Vromans, M., 2005. Operations Research in
passengerrailway transportation. Statistica Neerlandica, 59(4), pp. 467-497.
Insee, 2014. Comparateur de territoire : Région d’Île-de-France (11). [Online] Available at:
https://www.insee.fr/fr/statistiques/1405599?geo=REG-11 [Accessed 29 06 2017].
Insee, 2015. Institut national de la statistique et des études économiques : Transports routiers
de fret de proximité. [Online] Available at: http://www.insee.fr [Accessed 13 05 2015].
Islam, D. et al., 2013. Logistics and supply chain management. Research in Transportation
Economics, 41(1), pp. 3-16.
Jacobson, S. & Schruben, L., 1989. Techniques for simulation response optimization.
Operations Research Letters, 8(1), pp. 1-9.
193
Bibliographie
Johnson, D., 1990. Local optimization and the Traveling Salesman Problem. 17th
International Colloquium on Automata, Languages, and Programming, Volume 443, pp. 446-
461.
Johnson, D. & Garey, M., 1985. A 71/60 theorem for bin packing. Journal of complexity,
Volume 1, pp. 65-106.
Kirkpatrick, S., Vecchi, M. & Gelatt Jr, C., 1983. Optimization by simmulated annealing.
Science, 220(4598), pp. 671-680.
Kleijnen, J. & Wan, J., 2007. Optimization of simulated systems: OptQuest and alternatives.
Simulation Modelling Practice and Theory, 15(3), pp. 354-362.
Köhler, J., Whitmarsh, L., Nykvist, B. & Schilpero, M., 2009. A transitions model for
sustainable mobility. Ecological Economics, 68(12), pp. 2985-2995.
Laguna, M., Kelly, J., González-Velarde, J. & Glover, F., 1995. Tabu search for the
multilevel generalized assignment problem. European journal of operational research, 82(1),
pp. 176-189.
Land , A. & Doig, A., 1960. An automatic method of solving discrete programming problems.
Econometrica, 28(3), pp. 497-520.
Lawler, E. & Wood, D., 1966. Branch-and-bound methods : A survey. Operations research,
14(4), pp. 699-719.
Lindholm, M. & Behrends, S., 2012. Challenges in urban freight transport planning - a review
in the Baltic Sea Region. Journal of Transport Geography, Volume 22, pp. 129-136.
Lin, S. & Kernighan, B., 1973. An efficient heuristic algorithm for the traveling salesman
problem. Operations Research, Volume 21, pp. 498-516.
Macharis, C. & Bontekoning, Y., 2004. Opportunities for OR in intermodal freight transport
research: a review. European Journal of Operational Research, Volume 153, pp. 400-416.
Macharis, C., Caris, A., Jourquin, B. & Pekin, E., 2011. A decision support framework for
intermodal transport policy. European Transport Research Review, 3(4), pp. 167-178.
MacHaris, . C., Van Mierlo, J. & Van Den Bossche, P., 2007. Combining intermodal transport
with electric vehicles: towards more sustainable solutions.. Transportation Planning and
Technology, 30(2-3), pp. 311-323.
Mackworth, A., 1977. Consistency in networks of relations. Artificial Intelligence, 8(1), pp.
99-118.
194
Bibliographie
Mairie de Paris, 2009. Le transport de marchandises dans la ville de Paris, Paris: Direction
de la Voirie & des Déplacements, Agence de la Mobilité - Usagers & Marchandises.
Mairie de Paris, 2009. Le transport de marchandises dans la ville de Paris, s.l.: Direction de
la Voirie & des Déplacements, Agence de la Mobilité - Usagers & Marchandises.
Martello, S. & Toth, P., 1990. Knapsack problems: algorithms and computer
implementations. Chichester, England: John Wiley & Sons Ltd.
Martí, R., Laguna, M. & Glover, F., 2006. Principles of scatter search. European Journal of
Operational Research, 169(2), pp. 359-372.
Meketon, M., 1987. Optimization in simulation: a survey of recent results.. Proceedings of the
1987 Winter Simulation Conference. Piscataway, NJ: IEEE, pp. 58-67.
Melo, M., Nickel, S. & Saldanha-da-Gama, F., 2009. Facility location and supply chain
management – A review. European Journal of Operational Research, 196(2), pp. 401-412.
MONOPRIX, 2007. L’acheminement des marchandises par voie ferrée et véhicules roulant au
GNV vers les magasins monoprix et MONOP’ de Paris. Dossier de Presse Ferroviaire, July
04, MONOPRIX.
Montreuil, B., Meller , R. & Ballot, E., 2012. Physical Internet Foundations. Proceedings of
the 14th IFAC Symposium on Information Control Problems in Manufacturing Bucharest,
Romania.
Motraghi, A. & Marinov, M., 2012. Analysis of urban freight by rail using event based
simulation. Simulation Modelling Practice and Theory, Volume 25, pp. 73-89.
MTA, 2014. NYCT Trash Can Free Stations Pilot Update. Report Presentation, January,
Metropolitan Transportation Authority "MTA".
Narendra, P. & Fukunaga, K., 1977. A branch and bound algorithm for feature subset
selection. IEEE transactions on computers, 26(9), pp. 917-922.
Nations United, 2015. World Urbanization Prospects - The 2014 Revision, s.l.: United
Nations Department of Economic and Social Affairs/Population Division.
Nemhauser, G. & Wolsey, L., 1988. Integer and Combinatorial Optimization. s.l.:John Wiley
& Sons.
Nielsen, L., Kroon, L. & Maróti, G., 2012. A rolling horizon approach for disruption
management of railway rolling stock. European Journal of Operational Research, Volume
220, pp. 496-509.
195
Bibliographie
Nuzzolo, A., Crisalli, U. & Comi, A., 2008. Metropolitan freight distribution by railways. In:
R. G. T. Eiichi Taniguchi, ed. Innovations in City Logistics. New York: Nova Science
Publishers, pp. 351-368.
OECD, 2014. OECD Economic Surveys: European Union 2014, s.l.: OECD Publishing.
Offer, G., Howey, D., Contestabile, M. & Clague , R., 2010. Comparative analysis of battery
electric, hydrogen fuel cell and hybrid vehicles in a future sustainable road transport system.
Energy Policy, 38(1), pp. 24-29.
Padberg, M. & Rinaldi, G., 1991. A branch-and-cut algorithm for the resolution of large-scale
symmetric traveling salesman problems. SIAM review, 33(1), pp. 60-10.
Paris Region, 2017. Paris region: key figures 2017. IAU Ile-de-France.
Patier, D., 2004. La place du transport de marchandises en compte propre, Lyon: Laboratoire
d'économie des transport.
Pellegrini, P., Marlière, G. & Rodriguez, J., 2012. Real time railway traffic management
modeling track-circuits. In ATOMOS, 12th Workshop on Algorithmic Approaches for
Transportation Modelling, Optimization, and Systems, Volume 25, pp. 23-34.
Pielage, B., 2001. Underground Freight Transportation - A new development for automated
freight transportation systems in the Netherlands. IEEE Intelligent Transportation Systems
Conference Proceedings - Oakland (CA) USA.
Pimor, Y. & Fender, M., 2008. Logistique - 5ème édition - Production - Distribution -
Soutien. l'usine nouvelle ed. Paris: Dunod.
PREDIT GO n°3, 2012. Faisabilité de l’intégration logistique des espaces transport, Paris:
Rapport finale.
PREDIT GO n°6, 2012. Transport collectifs urbains : surcouts consécutifs à l’ajout logistique
dans un projet, Paris: Rapport intermédiaire.
Pritsker, A., 1986. Introduction to simulation and SLAM II. New York: Halsted Press.
Pritsker, A., Sigal, C. & Hammesfahr, R., 1989. SLAM II : Network Models for decision
Support. Englewood Cliffs: Prentice Hall.
196
Bibliographie
Quak, H., 2008. Sustainability of urban freight transport: Retail distribution and local
regulations in cities. PhD Thesis, Erasmus Universiteit Rotterdam..
Ramalhinho Lourenço, H. & Serra, D., 2000. Adaptative search heuristics for the generalized
assignment problem. Mathware & soft computing, Volume 7, pp. 1-15.
Reece, D. & Marinov, M., 2015. Modelling the implementation of a Baggage Transport
System in Newcastle upon Tyne for passengers using Mixed-Mode Travel. Transport
Problems, 10(4), pp. 149-155.
Ripert, C., 2004. Transport de marchandises dans Paris. Mairie de Paris, présentation
Pollutec.
Rodriguez, J., 2007. A constraint programming model for real-time train scheduling at
junctions. Transportation Research Part B: Methodological, 41(2), pp. 231-245.
Rossetti, M., 2015. Simulation Modeling and Arena, 2nd ed. New York: John Wiley & Sons,
Incorporated.
Ross, G. & Soland, P., 1975. A branch and bound based algorithm for the generalized
assignment problem. Mathematical Programming, Volume 8, pp. 91-103.
Routhier, J., Segalou, E. & Durand, S., 2001. Mesurer l'impact du transport de marchandises
en ville : le modèle de simulation FRETURB (V.1), s.l.: LET, MTL, ADEME.
Roy, B., Bertier, P. & Nghiem, P., 1965. Programmes linéaires en nombre entiers et
procédures SEP. Metra, 4(3).
Roy, R., 2010. A primer on the Taguchi method, 2nd edition. s.l.:Society of Manufacturing
Engineers.
Runge, D. & Becker, H., 2008. Urban Commercial Transport - Background Paper for the
Commission 4 Meeting on 22nd October 2008. World Association of the Major Metropolises
- Urban Mobility Management, Metropolis, Berlin.
Sabar, M., Montreuil, B. & Frayret, J., 2009. An Agent-based Algorithm for Personnel
Scheduling and Rescheduling in Assembly Centers. Proceedings of the 13th IFAC Symposium
on Information Control Problems in Manufacturing, Moscow, Russia.
Sabre, R., 2007. Plans d’expériences - Méthode de Taguchi. Techniques de l'ingénieur, F1006
V1.
Salido, M., Escamilla, J., Barber, F. & Giret, A., 2017. Rescheduling in job-shop problems for
sustainable manufacturing systems. Journal of Cleaner Production, Volume 162, pp. S121-
S132.
197
Bibliographie
Sato, K., Tamura, K. & Tomii, N., 2013. A MIP-based timetable rescheduling formulation
and algorithm minimizing further inconvenience to passengers. Journal of Rail Transport
Planning & Management, Volume 3, pp. 38-53.
Savelsbergh, M., 1997. A branch-and-price algorithm for the generalized assignment problem.
Operations research, 45(6), pp. 831-841.
Schrijver, A., 2005. On the History of Combinatorial Optimization (Till 1960). Handbooks in
Operations Research and Management Science, Volume 12, pp. 1-68.
Siarry, P., 2014. Métaheuristiques : recuits simulé, recherche avec tabous, recherche à
voisinage variables, méthodes GRASP, algorithmes évolutionnaires, fourmis artificielles,
essaims particulaires et autres méthodes d'optimisation. Paris: Editions Eyrolles.
Silva, C., Sousa, J. & Runkler, T., 2008. Rescheduling and optimization of logistic processes
using GA and ACO. Engineering Applications of Artificial Intelligence, 21(3), pp. 343-352.
Southworth, F. & Peterson, B., 2000. Intermodal and international freight network modeling.
Transportation Research Part C: Emerging Technologies, 8(1-6), pp. 147-166.
Stützel, T., 1998. Local search algorithms for combinatorial problems analysis, improvements
and new applications. PhD thesis, Departement of computer science, Darmstadt university of
technology, Germany.
Stützel, T. & Hoos, H., 2000. Max-Min ant system. Future generation computer systems,
Volume 16, pp. 889-914.
Stützle, T. & Dorigo, M., 2002. A short convergence proof for a class of ant colony
optimization algorithms. IEEE Transactions on evolutionary computation, 6(4), pp. 358-365.
Sun, X. & Quilliot, A., 1993. Résolution de la programmation linéaire en nombre entiers.
Clermont-Ferrand: Thèse de Doctorat à l'Université Blaise Pascal.
Taguchi, G., 1986. Introduction to quality engineering: designing quality into products and
processes. s.l.:Asian Productivity Organization.
Thierry, C., 2003. Gestion de chaînes logistiques, Modèles et mise en oeuvre pour l’aide à la
décision à moyen terme. Mémoire HDR, Université Toulouse 2, Le Mirail.
Toilier, F. et al., 2016. How can Urban Goods Movements be Surveyed in a Megacity? The
Case of the Paris Region. Transportation Research Procedia, Volume 12, pp. 570-583.
Törnquist Krasemann, J., 2012. Design of an effective algorithm for fast response to the re-
scheduling of railway traffic during disturbances. Transportation Research Part C: Emerging
Technologies, 20(1), pp. 62-78.
198
Bibliographie
Trentini, A., 2012. Proposition d’un système de transport urbain mixte. Application dans le
cadre de la ville moyenne de La Rochelle. Thèse de doctorat de l’école nationale supérieure
des mines de Paris.
Vamanan, M., Wang, Q., Batta, R. & Szczerba, R., 2004. Integration of COTS software
products ARENA & CPLEX for an inventory/logistics problem. Computers & Operations
Research, 31(4), pp. 533-547.
Van Audenhove, F., De Jongh, S. & Durance, M., 2015. Urban Logistics – How to unlock
value from last mile delivery for cities, transporters and retailers, s.l.: Arthur D. Little Future
of Urban Mobility Lab.
Van Duin, J., Tavasszy, L. & Quak, H., 2013. Towards E(lectric)- urban freight: first
promising steps in the electric vehicle revolution. European Transport \ Trasporti Europei,
Volume 54, p. Paper n° 9.
Van Wee, B., Maat, K. & De Bont, C., 2012. Improving Sustainability in Urban Areas:
Discussing the Potential for Transforming Conventional Car-based Travel into Electric
Mobility. European Planning Studies, 20(1), pp. 95-110.
Vanderbeck, F. & Wolsey, L., 1996. An exact algorithm for ip column generation. Operations
research letters, 19(4), pp. 151-159.
Veelenturf, L., Potthoff, D., Huisman, D. & Kroon, L., 2012. Railway crew rescheduling with
retiming. Transportation Research Part C: Emerging Technologies, 20(1), pp. 95-110.
Vieira, G., Herrmann, J. & Lin, E., 2003. Rescheduling manufacturing systems: A framework
of strategies, policies, and methods. Journal of Scheduling, Volume 6, pp. 39-62.
Visser, J., 2017. Automated freight transport. 1th Chinas Symposium on Underground Freight
Transportation - Shanghai.
Walker, C., Snowdon, J. & Ryan, D., 2005. Simultaneous disruption recovery of a train
timetable and crew roster in real time. Computers & Operations Research, 32(8), pp. 2077-
2094.
Wang, L. & Shi, L., 2013. Simulation optimization: a review on theory and applications. Acta
Automatica Sinica, 39(11), pp. 1957-1968.
Wang, Y. et al., 2016. Towards enhancing the last-mile delivery: An effective crowd-tasking
model with scalable solutions. Transportation Research Part E, Volume 93, pp. 279-293.
Wu, S., Storer, R. & Chang, P., 1992. A rescheduling procedure for manufacturing systems
under random disruptions. In: New directions for operations research in manufacturing (G.
Fandel, T. Gulledge, & A. Jone ). Berlin: Springer, pp. 292-308.
199
Bibliographie
Xie, M. & Li, X., 2012. Railway Timetable Rescheduling Based on the Feedback of. Procedia
- Social and Behavioral Sciences: 8th International Conference on Traffic and Transportation
Studies Changsha, China, Volume 43, pp. 781-789.
Yin, Y., Cheng, T. & Wang, D., 2016. Rescheduling on identical parallel machines with
machine disruptions to minimize total completion time. European Journal of Operational
Research, Volume 252, pp. 737-749.
Zeigler, B., 1976. Theory of Modeling and simulation. New York: John Wiley & Sons.
200
Liste des publications
201
Liste des publications
202
Résumé / Abstract
Résumé / Abstract
203
Résumé / Abstract
Titre : Une méthodologie pour modéliser et optimiser la mutualisation du transport ferroviaire urbain
de marchandises et de passagers
Malgré la prédominance actuelle du mode routier, pour le transport de marchandises en milieu urbain,
une alternative durable est nécessaire, au vu des enjeux environnementaux et sociétaux. Dans cette
thèse, nous proposons l’étude d’une des perspectives possibles, pour absorber une partie de ce flux de
marchandises toujours plus dense, en utilisant le réseau ferroviaire urbain, initialement dédié aux
voyageurs. Une méthodologie intégrant le fret dans ce dernier est proposée, avec comme première
étape, l'identification et la classification de tous les niveaux de mixité fret / voyageurs possibles. Le
niveau le plus contraint est retenu, car sa faisabilité induira celle des autres. Notre seconde
contribution est relative à une approche par décomposition du problème d’insertion du flux de fret en
plusieurs sous-problèmes interdépendants, selon les trois horizons temporels (long, moyen et court).
Dans le but d’évaluer la capacité du système global, à absorber un flux supplémentaire de nature
différente, le problème de détermination du meilleur plan de transport des marchandises est identifié
comme central et critique. La troisième contribution concerne la simulation du système de transport,
puis sa formalisation par un PL en variables mixtes, pour affecter chaque commande à un train, en
déterminant le moment auquel elle sera chargée et en minimisant les temps d’attente cumulés des
commandes. Plusieurs variantes de colonies de fourmis sont développées, pour la résolution
d’instances de grande taille. La quatrième contribution concerne le couplage du modèle de simulation,
qui permet l’évaluation des performances de cette nouvelle solution de transport, avec les différents
algorithmes optimisant le plan de transport. Enfin, nous proposons une approche de replanification par
horizon glissant, pour absorber les perturbations de la demande, en minimisant les changements du
plan de transport.
Title: A modeling methodology to introduce freight into urban passenger rail network
Keywords: urban freight, rail transport, optimization, discrete event simulation, simulation
optimization coupling.
Urban freight transport is almost exclusively carried out by truck. Beyond the drawbacks caused in the
city, this transport mode is nearly saturated. This study discusses an alternative way of transporting
freight by using urban rail infrastructure. The first contribution deals with the identification and
classification of all different sharing possibilities of mixing freight with passenger’s traffic using rail
network. The second contribution is the definition of global freight/passenger transport problem,
which is decomposed into several optimization interdependent sub-problems with different temporal
decision horizon. In order to show the capacity of the global system to absorb an additional flow with
different nature, the Freight Rail Transport Schedule Problem “FRTSP” is identified as the bottleneck
of transportation system and is formalized with MIP model. As third contribution, this problem
determines train and loading time for each demand to be assigned respecting several constraints while
minimizing total waiting time. The fourth contribution deals with a discrete event simulation
approach, which studies this alternative and validates several proposed decision algorithms. Finally,
the fifth contribution consists in a dynamic approach based on a rolling horizon, which is proposed in
order to update the initial plan. The updated plan allows to determine a new assignment regarding new
demand such as the modifications from the previous plan are minimized.
204