# Distributed Monitoring Based onP-Time Petri Nets and Chronicle Recognition of the Tunisian Railway Network

Authors: Anis M’halla; Dimitri Lefebvre; Mouhaned Gaied
Journal: Journal of Electrical and Computer Engineering
 (2020)

---

## Abstract

This paper falls under the problems of the monitoring of a Discrete Event System (DES) with time constraints. Among the various techniques used for online and distributed monitoring, we are interested in the chronicle recognition. Chronicles are temporal patterns that represent the system’s possible evolutions. The proposed models are based onP-time Petri nets that are suitable to represent with accuracy and modularity the Tunisian railway network. These models are scalable and may be used to represent a large variety of railway networks. Then, monitoring is based on the generation of chronicles that are suitable to detect and isolate traffic incidents in a distributed setting. Consequently, the proposed approach is tractable for large networks. Finally, to demonstrate the effectiveness and accuracy of the approach, an application to the case study of the Tunisian railway network is outlined.


---

## Body

## 1. Introduction

Transport systems need to be monitored online to avoid critical situations and temporal disturbances. These disturbances concern either the railway infrastructure or the traffic management and may lead to service disruptions due to weather, obstacles on the tracks, malice, social movements, and so on. In the systems under consideration, the processing times are interval-valued with parameters that depend on the operation to be performed. Any deviations from the specified interval will characterize a traffic disturbance. Consequently, the monitoring of the time intervals will be used as the main principle to detect and isolate the disturbances that affect the system. In this context, we propose the chronicles as a formalism for interpreting events to monitor temporal disruptions. A chronicle is composed of a set of events and a set of temporal constraints linking the pairs of events. The proposed approach is applied to the railway transport network in Tunisia. The main contribution of this paper is to use the combination of Petri net formalism and chronicles for the analysis of traffic disruptions. The overall goal is to maintain the stability and efficiency of the studied railway network.Over the past decade, the chronicle formalism has been developed and used by numerous authors, particularly for faults diagnosis [1–4]. It was also used by the Telecom operators for fault diagnosis in telecommunication networks [5] or for the diagnosis of web service failures [6, 7]. Lately, to deal with the time aspects inherent to the chronicles, the authors in [8] propose an automated translation of chronicles into a set of Labeled Time Petri Nets with Priorities (LTPNPr). The aim of this study was to improve the characterization of diagnosability for a discrete event system by taking into account the notion of finite durations between the occurrence of two events. The diagnosability study was done by comparing the chronicles associated with the normal behavior to the chronicle associated with faulty behaviors. The diagnosability analysis is based on LTPNPr in which priority relationships on transitions are used to complete the firing conditions. This approach is different from our contribution that uses P-time Petri nets for modelling issues and that excludes the notion of priority for transition firings. In [9], the authors proposed a diagnosis approach based on chronicles and modular temporal analysis. In fact, a failure was allocated to each set of chronicles, and each chronicle recognized the fault signature which was obtained from the state diagnosed with the finite state automata associated to each process module.To describe clearly all desired behaviors of an automated production system, an incremental learning approach, based on Causal Temporal Signatures (CTS) was proposed in [10]. The CTS was defined as “a subset of partially-ordered observable events that characterizes the system’s faulty behavior”. The formalism of CTS, based on chronicles, was used to diagnose faults in Discrete Event Systems (DES). The main benefits of this formalism were the readability and the expressiveness. In [11], the authors proposed a chronicle model allowing integrating prohibition constraints representing the absence of events in the model. Based on this formalism, criteria were proposed to characterize the chronicles and to compare them in order to evaluate the relevance of the chronicle database used for the diagnosis. The chronicles were modeled by temporal Petri nets, and the exclusivity test was solved by a model checking approach.The association of the Petri net formalism and chronicles for the monitoring and diagnosis is widely proposed in the literature [12–15]. In these works, the notions of causality, conflict, competition, unfolding, time branching formulated in Petri net setting lead to the modelling, identification, and discrimination of fault patterns. These approaches can be extended to more complex monitoring problems where the places not only represent faults but also process states. This context was also exploited in our approach: no specific assumption is made about what the chronicle represents, and it can be associated with normal or faulty behaviors.In transportation systems, a large number of observations are regularly collected and should be processed. Such a large amount of data cannot be treated online efficiently by a human operator. On the contrary, chronicles are suitable to analyze such data in order to recognize normal and faulty behaviors. The main advantage of this tool is the ability to describe and represent the possible evolutions of the transport systems and to recognize these behaviors in a flow of events.The contributions presented in this paper are devoted to distributed monitoring of rail transport networks and deal with the detection of traffic disturbance symptoms in these systems.P-time Petri nets are used for modelling aspects. Among the various techniques used for the distributed monitoring, we are interested in the chronicles recognition.The paper is organized as follows. Section two presents the railway transport system in Tunisia and gives some basic notions aboutP-time Petri nets. The third section presents the chronicle formalism and illustrates how chronicles are applied to the distributed monitoring. In Section four, the approach is applied to the Tunisian railway transport network. The main purpose of this study is to avoid time constraints violation and to guarantee the stability and safety of railway traffic. Finally, the last section concludes the paper. Some perspectives for further works are also proposed.


## 2. Modelling of the Tunisia Sahel Railway Network

### 2.1. Presentation of the Studied Transport Network

The railway network of the Sahel Tunisia, Figure1, ensures the transportation of passengers between the main cities of the Sahel from Sousse to Mahdia. It starts at the Sousse Bab Jdid station, crosses the industrial zone of Sousse, and bypasses the Sahline city until the Monastir station by the serving the airport, the tourist zone, and the University of Monastir; Figure 1. South of Monastir, the line continues to the Ksar Hellal city going through the industrial zone, and the main agglomerations located on this axis are Khniss, Ksibet Elmedouini, Bouhjar, and Sayada. Leaving Moknine, it serves the industrial zone of Teboulba and continues to Mahdia. The particularity of this railway network is the single line connecting the Bekalta station to Mahdia tourist zone station. Therefore, the segments are unidirectional in this area, while it has double tracks in the rest of the network.Figure 1
Sahel railway network.With an average of 40 minutes of frequency, the metro of Sahel ensures daily 44 journeys scheduled between 5 : 00 am until 10 : 00 pm and ensures the transport of more than 9 million passengers per year with an average of 27,000 passengers daily. Figure2 shows the travel times between the various stations deserved by the metro of the Sahel.Figure 2
Travel times between stations belonging to the Sahel railway transport network.


### 2.2. Model of the Studied Railway Network

#### 2.2.1.P-Time Petri Net

Several modelling approaches have been proposed for transport networks. Models may be graphical or analytical or both as the same time as Petri nets. In this work,P-time Petri Nets (P-TPNs) are used for modelling the transport network. In such models, travel durations are included within a minimum value and a maximum one.Definition 1 (see Reference [16]).
The formal definition of a P-TPN is given by a pair <R, IS>, where R is a marked Petri net:(1)IS:P⟶Q+×Q+∪+∞,pi⟶ISi=ai,biwith 0≤ai≤bi.
ISi defines the static interval of sojourn time of a mark in the place pi belonging to the set of places P (Q+ is the set of positive rational numbers). A token in the place pi fires an output transition when it has stayed in pi a duration at least ai and at most bi. After the duration bi, there is a death of mark in pi.


#### 2.2.2. Modelling

The Sahel railway network consists of 3 main stations (Sousse, Monastir, and Mahdia) and 28 other stations, Figure3. The travel from Sousse to Monastir is represented by the places p63 to p45, Figure 3. In the considered application, the white places represent the paths between two stations, the red places represent the stations, and the gray places represent resources introduced to avoid catching-up. A specific module for bidirectional segments is detailed in [17] and used in Figure 3.Figure 3
P-TPN model for the Sahel railway network.The sojourn durationqi in each place pi represents either the travelling time between the station or the parking times in stations. Such duration has an expected value (qie) which should be computed in order to satisfy the planned schedule (Table 1).Table 1
Sectors of the railway network between Sousse and Monastir.
Section 1: Sousse Med 5 sectorSection 2: Sousse Sud sectorSection 3: Sousse Z.Ind sectorSection 4: Sahline sectorSection 5: Sahline Sabkha sectorSection 6: Hôtels sectorSection 7: Les Hôtels sectorSection 8: L’Aeroport sectorSection 9: La faculté sectorSection 10: Monastir sectorTable 2
The values of the static intervals ISi and the expected sojourn time of the token in the place pi.
piISiqei1IS1 = [60, inf]992IS2 = [230, 250]2363IS3 = [60, 120]804IS4 = [172, 188]1775IS5 = [60, 120]806IS6 = [174, 186]1787IS7 = [60, 120]808IS8 = [173, 187]1779IS9 = [60, 120]8010IS10 = [168, 192]17612IS12 = [60, 120]8013IS13 = [170, 190]17615IS15 = [60, 120]8016IS16 = [227, 253]23517IS17 = [60, 120]8018IS18 = [111, 129]11719IS19 = [60, 120]8020IS20 = [168, 192]17621IS21 = [60, 120]8022IS22 = [112, 128]11723IS23 = [60, 120]8024IS24 = [110, 130]11625IS25 = [60, 120]8026IS26 = [113, 127]11727IS27 = [60, 120]8028IS28 = [107, 133]11529IS29 = [60, 120]8030IS30 = [109, 131]11631IS31 = [60, 120]8032IS32 = [110, 130]11633IS33 = [60, 120]8034IS34 = [172, 188]17735IS35 = [60, 120]8036IS36 = [165, 195]17537IS37 = [60, 120]8038IS38 = [106, 134]11539IS39 = [60, 120]8040IS40 = [110, 130]11641IS41 = [60, 120]8042IS42 = [110, 130]11643IS43 = [60, 120]8044IS44 = [832, 848]83745IS45 = [60, inf]7846IS46 = [233, 247]23747IS47 = [60, 120]8048IS48 = [351, 369]35749IS49 = [60, 120]8050IS50 = [108, 132]11651IS51 = [60, 120]8052IS52 = [166, 194]17553IS53 = [60, 120]8054IS54 = [230, 250]23655IS55 = [60, 120]8056IS56 = [173, 187]17757IS57 = [60, 120]8058IS58 = [171, 189]17759IS59 = [60, 120]8060IS60 = [170, 190]17661IS61 = [60, 120]8062IS62 = [113, 127]11763IS63 = [60, inf]7164IS64 = [110, 130]11665IS65 = [60, 120]8066IS66 = [168, 192]17667IS67 = [60, 120]8068IS68 = [166, 194]17569IS69 = [60, 120]8070IS70 = [167, 193]17571IS71 = [60, 120]8072IS72 = [233, 247]23773IS73 = [60, 120]8074IS74 = [176, 184]17875IS75 = [60, 120]8076IS76 = [110, 130]11677IS77 = [60, 120]8078IS78 = [349, 371]35679IS79 = [60, 120]8080IS80 = [228, 252]23681IS81 = [832, 848]83782IS82 = [60, 120]8083IS83 = [111, 129]11784IS84 = [60, 120]8085IS85 = [110, 130]11686IS86 = [60, 120]8087IS87 = [107, 133]11588IS88 = [60, 120]8089IS89 = [173, 187]17790IS90 = [60, 120]8091IS91 = [170, 190]17692IS92 = [60, 120]8093IS93 = [105, 135]11594IS94 = [60, 120]8095IS95 = [115, 125]11896IS96 = [60, 120]8097IS97 = [109, 131]11698IS98 = [60, 120]8099IS99 = [108, 132]116100IS100 = [60, 120]80101IS101 = [108, 132]116102IS102 = [60, 120]80103IS103 = [105, 135]115104IS104 = [60, 120]80105IS105 = [174, 186]178106IS106 = [60, 120]80107IS107 = [108, 132]116108IS108 = [60, 120]80109IS109 = [230, 250]236110IS110 = [171, 189]177111IS111 = [167, 193]175112IS112 = [175, 185]178113IS113 = [60, 120]80114IS114 = [172, 188]177115IS115 = [60, 120]80116IS116 = [170, 190]176117IS117 = [60, 120]80118IS118 = [233, 247]237The main parameters of the models are defined as follows:A time interval [Li, Hi] expressed with time unit TU is associated to each segment between two successive stations. Its lower bound Li indicates the minimum time required for the travel and the upper bound Hi fixes the maximum time not to be exceeded in order to avoid any temporal disturbance on railway traffic.A train can park in a principal station for at least one minute. The static intervals associated with the three main stations areISp1 = ISp45 = ISp63 = [60, +∞]; Figure 3.The sojourn time of a metro in any other station is estimated from one to two minutes: ISi = [60, 120]; Figure 3.At the beginning of each day, it is assumed that 4 trains start from the Mahdia station, and 2 others are stationed at the Sousse station; Figure3.The static intervals ISi and the effective sojourn time qie associated with the stations and with the segments between two successive stations are detailed in Table 2 [17]. The static intervals are defined based on the SNCFT traffic dataset.


## 2.1. Presentation of the Studied Transport Network

The railway network of the Sahel Tunisia, Figure1, ensures the transportation of passengers between the main cities of the Sahel from Sousse to Mahdia. It starts at the Sousse Bab Jdid station, crosses the industrial zone of Sousse, and bypasses the Sahline city until the Monastir station by the serving the airport, the tourist zone, and the University of Monastir; Figure 1. South of Monastir, the line continues to the Ksar Hellal city going through the industrial zone, and the main agglomerations located on this axis are Khniss, Ksibet Elmedouini, Bouhjar, and Sayada. Leaving Moknine, it serves the industrial zone of Teboulba and continues to Mahdia. The particularity of this railway network is the single line connecting the Bekalta station to Mahdia tourist zone station. Therefore, the segments are unidirectional in this area, while it has double tracks in the rest of the network.Figure 1
Sahel railway network.With an average of 40 minutes of frequency, the metro of Sahel ensures daily 44 journeys scheduled between 5 : 00 am until 10 : 00 pm and ensures the transport of more than 9 million passengers per year with an average of 27,000 passengers daily. Figure2 shows the travel times between the various stations deserved by the metro of the Sahel.Figure 2
Travel times between stations belonging to the Sahel railway transport network.


## 2.2. Model of the Studied Railway Network

### 2.2.1.P-Time Petri Net

Several modelling approaches have been proposed for transport networks. Models may be graphical or analytical or both as the same time as Petri nets. In this work,P-time Petri Nets (P-TPNs) are used for modelling the transport network. In such models, travel durations are included within a minimum value and a maximum one.Definition 1 (see Reference [16]).
The formal definition of a P-TPN is given by a pair <R, IS>, where R is a marked Petri net:(1)IS:P⟶Q+×Q+∪+∞,pi⟶ISi=ai,biwith 0≤ai≤bi.
ISi defines the static interval of sojourn time of a mark in the place pi belonging to the set of places P (Q+ is the set of positive rational numbers). A token in the place pi fires an output transition when it has stayed in pi a duration at least ai and at most bi. After the duration bi, there is a death of mark in pi.


### 2.2.2. Modelling

The Sahel railway network consists of 3 main stations (Sousse, Monastir, and Mahdia) and 28 other stations, Figure3. The travel from Sousse to Monastir is represented by the places p63 to p45, Figure 3. In the considered application, the white places represent the paths between two stations, the red places represent the stations, and the gray places represent resources introduced to avoid catching-up. A specific module for bidirectional segments is detailed in [17] and used in Figure 3.Figure 3
P-TPN model for the Sahel railway network.The sojourn durationqi in each place pi represents either the travelling time between the station or the parking times in stations. Such duration has an expected value (qie) which should be computed in order to satisfy the planned schedule (Table 1).Table 1
Sectors of the railway network between Sousse and Monastir.
Section 1: Sousse Med 5 sectorSection 2: Sousse Sud sectorSection 3: Sousse Z.Ind sectorSection 4: Sahline sectorSection 5: Sahline Sabkha sectorSection 6: Hôtels sectorSection 7: Les Hôtels sectorSection 8: L’Aeroport sectorSection 9: La faculté sectorSection 10: Monastir sectorTable 2
The values of the static intervals ISi and the expected sojourn time of the token in the place pi.
piISiqei1IS1 = [60, inf]992IS2 = [230, 250]2363IS3 = [60, 120]804IS4 = [172, 188]1775IS5 = [60, 120]806IS6 = [174, 186]1787IS7 = [60, 120]808IS8 = [173, 187]1779IS9 = [60, 120]8010IS10 = [168, 192]17612IS12 = [60, 120]8013IS13 = [170, 190]17615IS15 = [60, 120]8016IS16 = [227, 253]23517IS17 = [60, 120]8018IS18 = [111, 129]11719IS19 = [60, 120]8020IS20 = [168, 192]17621IS21 = [60, 120]8022IS22 = [112, 128]11723IS23 = [60, 120]8024IS24 = [110, 130]11625IS25 = [60, 120]8026IS26 = [113, 127]11727IS27 = [60, 120]8028IS28 = [107, 133]11529IS29 = [60, 120]8030IS30 = [109, 131]11631IS31 = [60, 120]8032IS32 = [110, 130]11633IS33 = [60, 120]8034IS34 = [172, 188]17735IS35 = [60, 120]8036IS36 = [165, 195]17537IS37 = [60, 120]8038IS38 = [106, 134]11539IS39 = [60, 120]8040IS40 = [110, 130]11641IS41 = [60, 120]8042IS42 = [110, 130]11643IS43 = [60, 120]8044IS44 = [832, 848]83745IS45 = [60, inf]7846IS46 = [233, 247]23747IS47 = [60, 120]8048IS48 = [351, 369]35749IS49 = [60, 120]8050IS50 = [108, 132]11651IS51 = [60, 120]8052IS52 = [166, 194]17553IS53 = [60, 120]8054IS54 = [230, 250]23655IS55 = [60, 120]8056IS56 = [173, 187]17757IS57 = [60, 120]8058IS58 = [171, 189]17759IS59 = [60, 120]8060IS60 = [170, 190]17661IS61 = [60, 120]8062IS62 = [113, 127]11763IS63 = [60, inf]7164IS64 = [110, 130]11665IS65 = [60, 120]8066IS66 = [168, 192]17667IS67 = [60, 120]8068IS68 = [166, 194]17569IS69 = [60, 120]8070IS70 = [167, 193]17571IS71 = [60, 120]8072IS72 = [233, 247]23773IS73 = [60, 120]8074IS74 = [176, 184]17875IS75 = [60, 120]8076IS76 = [110, 130]11677IS77 = [60, 120]8078IS78 = [349, 371]35679IS79 = [60, 120]8080IS80 = [228, 252]23681IS81 = [832, 848]83782IS82 = [60, 120]8083IS83 = [111, 129]11784IS84 = [60, 120]8085IS85 = [110, 130]11686IS86 = [60, 120]8087IS87 = [107, 133]11588IS88 = [60, 120]8089IS89 = [173, 187]17790IS90 = [60, 120]8091IS91 = [170, 190]17692IS92 = [60, 120]8093IS93 = [105, 135]11594IS94 = [60, 120]8095IS95 = [115, 125]11896IS96 = [60, 120]8097IS97 = [109, 131]11698IS98 = [60, 120]8099IS99 = [108, 132]116100IS100 = [60, 120]80101IS101 = [108, 132]116102IS102 = [60, 120]80103IS103 = [105, 135]115104IS104 = [60, 120]80105IS105 = [174, 186]178106IS106 = [60, 120]80107IS107 = [108, 132]116108IS108 = [60, 120]80109IS109 = [230, 250]236110IS110 = [171, 189]177111IS111 = [167, 193]175112IS112 = [175, 185]178113IS113 = [60, 120]80114IS114 = [172, 188]177115IS115 = [60, 120]80116IS116 = [170, 190]176117IS117 = [60, 120]80118IS118 = [233, 247]237The main parameters of the models are defined as follows:A time interval [Li, Hi] expressed with time unit TU is associated to each segment between two successive stations. Its lower bound Li indicates the minimum time required for the travel and the upper bound Hi fixes the maximum time not to be exceeded in order to avoid any temporal disturbance on railway traffic.A train can park in a principal station for at least one minute. The static intervals associated with the three main stations areISp1 = ISp45 = ISp63 = [60, +∞]; Figure 3.The sojourn time of a metro in any other station is estimated from one to two minutes: ISi = [60, 120]; Figure 3.At the beginning of each day, it is assumed that 4 trains start from the Mahdia station, and 2 others are stationed at the Sousse station; Figure3.The static intervals ISi and the effective sojourn time qie associated with the stations and with the segments between two successive stations are detailed in Table 2 [17]. The static intervals are defined based on the SNCFT traffic dataset.


## 2.2.1.P-Time Petri Net

Several modelling approaches have been proposed for transport networks. Models may be graphical or analytical or both as the same time as Petri nets. In this work,P-time Petri Nets (P-TPNs) are used for modelling the transport network. In such models, travel durations are included within a minimum value and a maximum one.Definition 1 (see Reference [16]).
The formal definition of a P-TPN is given by a pair <R, IS>, where R is a marked Petri net:(1)IS:P⟶Q+×Q+∪+∞,pi⟶ISi=ai,biwith 0≤ai≤bi.
ISi defines the static interval of sojourn time of a mark in the place pi belonging to the set of places P (Q+ is the set of positive rational numbers). A token in the place pi fires an output transition when it has stayed in pi a duration at least ai and at most bi. After the duration bi, there is a death of mark in pi.


## 2.2.2. Modelling

The Sahel railway network consists of 3 main stations (Sousse, Monastir, and Mahdia) and 28 other stations, Figure3. The travel from Sousse to Monastir is represented by the places p63 to p45, Figure 3. In the considered application, the white places represent the paths between two stations, the red places represent the stations, and the gray places represent resources introduced to avoid catching-up. A specific module for bidirectional segments is detailed in [17] and used in Figure 3.Figure 3
P-TPN model for the Sahel railway network.The sojourn durationqi in each place pi represents either the travelling time between the station or the parking times in stations. Such duration has an expected value (qie) which should be computed in order to satisfy the planned schedule (Table 1).Table 1
Sectors of the railway network between Sousse and Monastir.
Section 1: Sousse Med 5 sectorSection 2: Sousse Sud sectorSection 3: Sousse Z.Ind sectorSection 4: Sahline sectorSection 5: Sahline Sabkha sectorSection 6: Hôtels sectorSection 7: Les Hôtels sectorSection 8: L’Aeroport sectorSection 9: La faculté sectorSection 10: Monastir sectorTable 2
The values of the static intervals ISi and the expected sojourn time of the token in the place pi.
piISiqei1IS1 = [60, inf]992IS2 = [230, 250]2363IS3 = [60, 120]804IS4 = [172, 188]1775IS5 = [60, 120]806IS6 = [174, 186]1787IS7 = [60, 120]808IS8 = [173, 187]1779IS9 = [60, 120]8010IS10 = [168, 192]17612IS12 = [60, 120]8013IS13 = [170, 190]17615IS15 = [60, 120]8016IS16 = [227, 253]23517IS17 = [60, 120]8018IS18 = [111, 129]11719IS19 = [60, 120]8020IS20 = [168, 192]17621IS21 = [60, 120]8022IS22 = [112, 128]11723IS23 = [60, 120]8024IS24 = [110, 130]11625IS25 = [60, 120]8026IS26 = [113, 127]11727IS27 = [60, 120]8028IS28 = [107, 133]11529IS29 = [60, 120]8030IS30 = [109, 131]11631IS31 = [60, 120]8032IS32 = [110, 130]11633IS33 = [60, 120]8034IS34 = [172, 188]17735IS35 = [60, 120]8036IS36 = [165, 195]17537IS37 = [60, 120]8038IS38 = [106, 134]11539IS39 = [60, 120]8040IS40 = [110, 130]11641IS41 = [60, 120]8042IS42 = [110, 130]11643IS43 = [60, 120]8044IS44 = [832, 848]83745IS45 = [60, inf]7846IS46 = [233, 247]23747IS47 = [60, 120]8048IS48 = [351, 369]35749IS49 = [60, 120]8050IS50 = [108, 132]11651IS51 = [60, 120]8052IS52 = [166, 194]17553IS53 = [60, 120]8054IS54 = [230, 250]23655IS55 = [60, 120]8056IS56 = [173, 187]17757IS57 = [60, 120]8058IS58 = [171, 189]17759IS59 = [60, 120]8060IS60 = [170, 190]17661IS61 = [60, 120]8062IS62 = [113, 127]11763IS63 = [60, inf]7164IS64 = [110, 130]11665IS65 = [60, 120]8066IS66 = [168, 192]17667IS67 = [60, 120]8068IS68 = [166, 194]17569IS69 = [60, 120]8070IS70 = [167, 193]17571IS71 = [60, 120]8072IS72 = [233, 247]23773IS73 = [60, 120]8074IS74 = [176, 184]17875IS75 = [60, 120]8076IS76 = [110, 130]11677IS77 = [60, 120]8078IS78 = [349, 371]35679IS79 = [60, 120]8080IS80 = [228, 252]23681IS81 = [832, 848]83782IS82 = [60, 120]8083IS83 = [111, 129]11784IS84 = [60, 120]8085IS85 = [110, 130]11686IS86 = [60, 120]8087IS87 = [107, 133]11588IS88 = [60, 120]8089IS89 = [173, 187]17790IS90 = [60, 120]8091IS91 = [170, 190]17692IS92 = [60, 120]8093IS93 = [105, 135]11594IS94 = [60, 120]8095IS95 = [115, 125]11896IS96 = [60, 120]8097IS97 = [109, 131]11698IS98 = [60, 120]8099IS99 = [108, 132]116100IS100 = [60, 120]80101IS101 = [108, 132]116102IS102 = [60, 120]80103IS103 = [105, 135]115104IS104 = [60, 120]80105IS105 = [174, 186]178106IS106 = [60, 120]80107IS107 = [108, 132]116108IS108 = [60, 120]80109IS109 = [230, 250]236110IS110 = [171, 189]177111IS111 = [167, 193]175112IS112 = [175, 185]178113IS113 = [60, 120]80114IS114 = [172, 188]177115IS115 = [60, 120]80116IS116 = [170, 190]176117IS117 = [60, 120]80118IS118 = [233, 247]237The main parameters of the models are defined as follows:A time interval [Li, Hi] expressed with time unit TU is associated to each segment between two successive stations. Its lower bound Li indicates the minimum time required for the travel and the upper bound Hi fixes the maximum time not to be exceeded in order to avoid any temporal disturbance on railway traffic.A train can park in a principal station for at least one minute. The static intervals associated with the three main stations areISp1 = ISp45 = ISp63 = [60, +∞]; Figure 3.The sojourn time of a metro in any other station is estimated from one to two minutes: ISi = [60, 120]; Figure 3.At the beginning of each day, it is assumed that 4 trains start from the Mahdia station, and 2 others are stationed at the Sousse station; Figure3.The static intervals ISi and the effective sojourn time qie associated with the stations and with the segments between two successive stations are detailed in Table 2 [17]. The static intervals are defined based on the SNCFT traffic dataset.


## 3. Chronicles Formalism

Chronicles provide formalism for monitoring transport systems. A chronicle is composed of a set of events, a set of temporal constraints linking the pairs of events, and the diagnosis test that describes the recognized situation [18].


### 3.1. Basic Definition

Preliminary definitions, useful for the rest of this paper, are given in order to explain the distributed detection principles.An event is a stimulus to which the system can react by a state change [19]. An event can occur after a message has been sent by the process at the beginning or at the end of any operation.Definition 2 (see Reference [19]).
The occurrence date is the time corresponding to an event issued from the process. LetO be the occurrence function which associates to each event ei its occurrence date O(ei); then,(2)O:E⟶Q+,ei⟶Oei,where E is the set of events.Definition 3.
A constraint is a relationship expressed by the duration between event occurrences. Two types of constraints can be distinguished; Figure4:Local constraints link timed events in the same monitoring siteGlobal constraints link timed events in different monitoring sites
In our study, the monitoring system (site)Si manages a physical zone (i.e., a set of sensors and resources) in the system considered as a Discrete Event System. A monitoring site is composed of six functions: failure detection, diagnosis, prognosis, follow, resumption, and emergency. Our study is limited to the detection function, whose role consists in recognizing a deviation from the normal (expected) functioning, locally or in a distributed manner by communication and cooperation between different monitoring sites.
In many cases, the local information generated by each monitoring site is insufficient, and communications between sites become essential to compensate this deficiency. When these systems are connected, we have a distributed monitoring system. The ease of implementation, test, and maintenance and the reduction of the software complexity in distributed systems are furthered by a high degree of modularity and a low degree of coupling between modules [20]. This type of systems can easily be extended if the existing elements change or if new elements are added; the reconfiguration is, therefore, encouraged.
In a distributed monitoring architecture, the communication aspects are significant and specific problems must be taken into account, such as the clocks synchronization, the reconstitution of the exchanged messages order, and the communication delays. In this context, the main contribution of this paper is the evaluation of the influence of the communication delays between modules, on the temporal constraints verification. In particular, the proposed solution, inspired from the work of Boufaied, evaluates the impact of the uncertainties that concern communications delays, on the global constraints in railway transport networks.Figure 4
Global and local constraints.


### 3.2. Time Constraints Verification

The time constraints verification is based on some assumptions:The communication delays Ω belong to the interval [Ωmin, Ωmax]Each monitoring site has its own clock, and all clocks have the same time base. Therefore, an event is dated only in the time frame associated to the receiving site. The occurrence function represents its temporal coordinate in the associated time frame.The proposed technique is suitable for any system in which the operational constraints and the communication time uncertainties are of the same order, such as railway transport systems.The problem to be solved can be summarized as follows:LetCX,Y be a global constraint linking an event eX to eY; Figure 5. Whatever be the nature of the constraint, it is always possible to express it as an interval. Thus, CX,Y will be defined by(3)dX,Y≤OeX−OeY≤fX,Y.Figure 5
Operating delays between monitoring sites.Since eventseX and eY are received and dated by two different sites, it is not possible to directly evaluate the duration O(eX) − O(eY). Therefore, it is necessary that the site SY informs the site SX when it received the event eY (if we consider that only the site X is responsible for the time constraint monitoring). This is performed by sending the event ek; Figure 5. This communication between the two sites is generally affected by a communication delay Ω with Ω∈Ωmin,Ωmax. In these conditions, the following are known:The global constraint CX,YThe occurrences times of eventsex and ek on the site XThe communication delay Ω between the two sitesX and Y (O(eX) − O(eY) = O(eX) − O(ek) + Ω)The question is that is it possible to know if the constraint CY,X is satisfied (i.e., if the occurrence time O(eY) satisfies the CY,X constraint)? It is possible to reformulate the initial global constraint into a new constraint which is a local constraint, for which it is possible to evaluate the involved durations.(4)OeX−Oek+Oek−OeY=Φ+Ω.As Ωmin ≤ Ω ≤ Ωmax, we obtain(5)Φ+Ωmin≤OeX−OeY≤Φ+Ωmax,which leads to(6)OeX−Oek+Ωmin≤OeX−OeY≤OeX−Oek+ΩmaxwithΦ∈−∞,+∞.The verification of the interval constraint consists, by means of the measurable duration Φ, of looking for the durationsO(eX) − O(eY) that verify both:(7)dX,Y≤OeX−OeY≤fX,Y,Φ+Ωmin≤OeX−OeY≤Φ+Ωmax.A graphical representation of these two constraints is shown in Figure6. In the plane (O(eY) − O(ek), O(eY) − O(eA)), the inequalities (1) and (2) define two bands. The searched durations belong to the intersection of these two bands which defines a polygon noted PO. The position of PO depends not only on the duration O(eY) − O(eX) but also on the terminal positions of the duration Ω with respect to dY,X and fY,X. For any duration Ω, PO is defined by 4 points: X1, X2, X3, and X4 of respective coordinates X1 = dY,X − Ωmax; dY,X − Ωmin; fY,X − Ωmin; and fY,X − Ωmax.Figure 6
Representation of the admissible domain [19].In many usual diagnosis approaches, the notions of uncertainty, inaccuracy, and incompleteness of information are modeled through bounded intervals. Whether the observation belongs to the expected interval leads to a binary reasoning about the detected inconsistencies. This type of binary reasoning may be insufficient for monitoring and diagnosis, especially for the evaluation of marginal deviations that can occur during transient phases. To enrich the reasoning, fuzzy sets and the theory of possibilities were integrated and lead to the definition of new fuzzy models. In order to quantify the set of possible durations Φ, a graphical representation, inspired from [19], is proposed; Figure 7. Considering a bounded delay, the possibility, that a time constraint is satisfied, belongs to the interval [0, 1]. Consequently, the verification of time constraints can be represented by fuzzy membership functions. This result makes it possible to highlight zones of certainty for the detection function: if the constraints to be verified are associated with a normal behavior of the monitored system, a high possibility value ensures a normal behavior of the monitored system; Figure 7. On the contrary, a low possibility value implies the detection of an illegal behavior.Figure 7
Possibility function for an interval constraint withX1 ≥ 0, X2 ≥ 0, X3 ≥ 0 and X4 ≥ 0.


## 3.1. Basic Definition

Preliminary definitions, useful for the rest of this paper, are given in order to explain the distributed detection principles.An event is a stimulus to which the system can react by a state change [19]. An event can occur after a message has been sent by the process at the beginning or at the end of any operation.Definition 2 (see Reference [19]).
The occurrence date is the time corresponding to an event issued from the process. LetO be the occurrence function which associates to each event ei its occurrence date O(ei); then,(2)O:E⟶Q+,ei⟶Oei,where E is the set of events.Definition 3.
A constraint is a relationship expressed by the duration between event occurrences. Two types of constraints can be distinguished; Figure4:Local constraints link timed events in the same monitoring siteGlobal constraints link timed events in different monitoring sites
In our study, the monitoring system (site)Si manages a physical zone (i.e., a set of sensors and resources) in the system considered as a Discrete Event System. A monitoring site is composed of six functions: failure detection, diagnosis, prognosis, follow, resumption, and emergency. Our study is limited to the detection function, whose role consists in recognizing a deviation from the normal (expected) functioning, locally or in a distributed manner by communication and cooperation between different monitoring sites.
In many cases, the local information generated by each monitoring site is insufficient, and communications between sites become essential to compensate this deficiency. When these systems are connected, we have a distributed monitoring system. The ease of implementation, test, and maintenance and the reduction of the software complexity in distributed systems are furthered by a high degree of modularity and a low degree of coupling between modules [20]. This type of systems can easily be extended if the existing elements change or if new elements are added; the reconfiguration is, therefore, encouraged.
In a distributed monitoring architecture, the communication aspects are significant and specific problems must be taken into account, such as the clocks synchronization, the reconstitution of the exchanged messages order, and the communication delays. In this context, the main contribution of this paper is the evaluation of the influence of the communication delays between modules, on the temporal constraints verification. In particular, the proposed solution, inspired from the work of Boufaied, evaluates the impact of the uncertainties that concern communications delays, on the global constraints in railway transport networks.Figure 4
Global and local constraints.


## 3.2. Time Constraints Verification

The time constraints verification is based on some assumptions:The communication delays Ω belong to the interval [Ωmin, Ωmax]Each monitoring site has its own clock, and all clocks have the same time base. Therefore, an event is dated only in the time frame associated to the receiving site. The occurrence function represents its temporal coordinate in the associated time frame.The proposed technique is suitable for any system in which the operational constraints and the communication time uncertainties are of the same order, such as railway transport systems.The problem to be solved can be summarized as follows:LetCX,Y be a global constraint linking an event eX to eY; Figure 5. Whatever be the nature of the constraint, it is always possible to express it as an interval. Thus, CX,Y will be defined by(3)dX,Y≤OeX−OeY≤fX,Y.Figure 5
Operating delays between monitoring sites.Since eventseX and eY are received and dated by two different sites, it is not possible to directly evaluate the duration O(eX) − O(eY). Therefore, it is necessary that the site SY informs the site SX when it received the event eY (if we consider that only the site X is responsible for the time constraint monitoring). This is performed by sending the event ek; Figure 5. This communication between the two sites is generally affected by a communication delay Ω with Ω∈Ωmin,Ωmax. In these conditions, the following are known:The global constraint CX,YThe occurrences times of eventsex and ek on the site XThe communication delay Ω between the two sitesX and Y (O(eX) − O(eY) = O(eX) − O(ek) + Ω)The question is that is it possible to know if the constraint CY,X is satisfied (i.e., if the occurrence time O(eY) satisfies the CY,X constraint)? It is possible to reformulate the initial global constraint into a new constraint which is a local constraint, for which it is possible to evaluate the involved durations.(4)OeX−Oek+Oek−OeY=Φ+Ω.As Ωmin ≤ Ω ≤ Ωmax, we obtain(5)Φ+Ωmin≤OeX−OeY≤Φ+Ωmax,which leads to(6)OeX−Oek+Ωmin≤OeX−OeY≤OeX−Oek+ΩmaxwithΦ∈−∞,+∞.The verification of the interval constraint consists, by means of the measurable duration Φ, of looking for the durationsO(eX) − O(eY) that verify both:(7)dX,Y≤OeX−OeY≤fX,Y,Φ+Ωmin≤OeX−OeY≤Φ+Ωmax.A graphical representation of these two constraints is shown in Figure6. In the plane (O(eY) − O(ek), O(eY) − O(eA)), the inequalities (1) and (2) define two bands. The searched durations belong to the intersection of these two bands which defines a polygon noted PO. The position of PO depends not only on the duration O(eY) − O(eX) but also on the terminal positions of the duration Ω with respect to dY,X and fY,X. For any duration Ω, PO is defined by 4 points: X1, X2, X3, and X4 of respective coordinates X1 = dY,X − Ωmax; dY,X − Ωmin; fY,X − Ωmin; and fY,X − Ωmax.Figure 6
Representation of the admissible domain [19].In many usual diagnosis approaches, the notions of uncertainty, inaccuracy, and incompleteness of information are modeled through bounded intervals. Whether the observation belongs to the expected interval leads to a binary reasoning about the detected inconsistencies. This type of binary reasoning may be insufficient for monitoring and diagnosis, especially for the evaluation of marginal deviations that can occur during transient phases. To enrich the reasoning, fuzzy sets and the theory of possibilities were integrated and lead to the definition of new fuzzy models. In order to quantify the set of possible durations Φ, a graphical representation, inspired from [19], is proposed; Figure 7. Considering a bounded delay, the possibility, that a time constraint is satisfied, belongs to the interval [0, 1]. Consequently, the verification of time constraints can be represented by fuzzy membership functions. This result makes it possible to highlight zones of certainty for the detection function: if the constraints to be verified are associated with a normal behavior of the monitored system, a high possibility value ensures a normal behavior of the monitored system; Figure 7. On the contrary, a low possibility value implies the detection of an illegal behavior.Figure 7
Possibility function for an interval constraint withX1 ≥ 0, X2 ≥ 0, X3 ≥ 0 and X4 ≥ 0.


## 4. Distributed Monitoring of the Sahel Railway Networks

### 4.1. Principle

Centralized monitoring consists of associating a single diagnoser to the whole process model. The latter collects the various process information before making a final decision on the process operating state. Although powerful in terms of diagnosis, the centralized structure is difficult to use for large and distributed systems such as transport systems since the building of a global model generates, in most cases, combinatorial explosion problems.In distributed monitoring structures, the process is decomposed into several local models. To each model is associated a local diagnoser. Each diagnoser makes his decision based on local observation. In the case of global specifications, a protocol allows the communication between the different diagnosers in order to make a final decision. The distributed monitoring, based on chronicles, verifies that each observed event is consistent with the chronicles specified by the time constraints. If the chronicles represent an abnormal behavior model, chronicles recognition is used for monitoring purposes.


### 4.2. Distributed Monitoring Architecture

In the studied railway networks, every set of sensors providing useful information has its own monitoring sites. Then, these sites are connected in distributed architecture. According to Figure6, the chronicles are distributed into several subchronicles associated to the different monitoring sites Si. Each diagnoser monitors locally a subsystem and communicates with other diagnosers in order to get necessary information in order to take decisions relative to the distributed diagnosis; Figure 8. Therefore, the detection function of each site has to recognize an evolution of the monitored process by subchronicles which are defined as sets of patterns of observable events temporally constrained. The recognition is performed through the verification of the time constraints associated to the considered subchronicles. The distributed detection function monitors the system evolution through the recognition of chronicles; Figure 8.Figure 8
Distributed monitoring architecture based on chronicles.


### 4.3. Distributed Monitoring of the Railway Sousse-Monastir

In order to help the supervisor in charge of managing the studied railway networks (i.e., detecting traffic perturbations, alert traveler claims, and maintain stability and security of the networks), a monitoring task is needed. The overall aim of the proposed monitoring approach is to control that the railway traffic proceeded well to avoid undesirable situations. The distributed monitoring of the Sahel railway networks based on the chronicle method involves the following steps.


#### 4.3.1. Step 1: Subsystems Identification

The first task is to subdivide the railway network system into sectors (subsystems).The railway network between Sousse and Monastir has been split into 10 sectors, as shown in Figure9. Each sector is indicated by the dashed lines. These sectors areFigure 9
Travel time from Sousse to Monastir stations.


#### 4.3.2. Step 2: Sensor Recognition

This step consists of identifying the sensors needed to perform the system monitoring. In the proposed approach, each siteSi, Figure 9, monitors a specific zone of the railway from Sousse to Monastir, through the observation of a set of events generated by the process. In the railway network, the travelling times should be within two bounds. Any deviations (occurrence of a temporal disturbance) from the allowed lower (resp. upper) bounds will lead to a low service quality and can lead to disaster scenarios. Therefore, the proposed monitoring approach, based on chronicle formalism, uses the additional information provided by the knowledge of interval constraints and allows detecting disruptions when a constraint is violated. Figure 9 gives journey times between the railway stations and parking times of the metro in stations. These durations are represented by the static intervals “ISi”.


#### 4.3.3. Step 3: Distributed Detection Function

To each event is associated a diagnoserSi, Figure 9, and the detection function is distributed on the different sites:eS1: Beginning journey from the Sousse Bab Jdid station (Sensor S1).eS2: Metro arrival at Sousse Med.5 (Sensor S2).eS3: Departure from the subway station Sousse Med.5 (Sensor S3).eS4: Metro arrival at Sousse Sud Station (Sensor S4).eS5: Departure from the subway station Sousse Sud (Sensor S5).eS6: Metro arrival at Sousse Z.Ind Station (Sensor S6).eS7: Departure from the subway station Sousse Z.Ind (Sensor S7).eS8: Metro arrival at Sahline Ville station (Sensor S8).eS9: Departure from the subway station Sahline Ville (Sensor S9).eS10: Metro arrival at Sahline Sabkha station (Sensor S10).eS11: Departure from the subway station Sahline Sabkha (Sensor S11).eS12: Metro arrival at the Les Hôtels station (Sensor S12).eS13: Departure from the subway station Les Hôtels (Sensor S13).eS14: Metro arrival at the L’Aeroport station (Sensor S14).eS15: Departure from the subway station L’Aeroport (Sensor S15).eS16: Metro arrival at the La faculté station (Sensor S16).eS17: Departure from the subway station La faculté (Sensor S17).eS18: Metro arrival at the Monastir station (Sensor S18).


#### 4.3.4. Step 4: Scenario Generation

A scenario is defined for each possible deviation in each section bringing the system to an erroneous situation. Deviations can occur due to unavailability of transportation systems, temporal disruptions, or traffic management and disturbances.Let us detail the monitoring between the two events,eS18 (metro arrival at the Monastir station) and eS1 (departure from Sousse Bab Jdid); Figure 9. To monitor this duration, it is necessary to check the time constraint linking the occurrences of the two events eS18 and eS1. This timing constraint is a global one; therefore, the verification of this constraint can be performed through the measure of the travelling time between the station and parking time of a metro in stations. As long as these durations are included in the mentioned intervals, no disturbance is detected. Otherwise, a traffic disturbance is detected. As previously mentioned, the global constraint to compute is an interval constraint type defined as CS18,S1: dS18,S1 ≤ O(eS18) − O(eS1) ≤ fS18,S1, with(8)dS18,S1=∑i=4563Li,fS18,S1=∑i=4563Hi.According to the time intervals (Figure3), the minimum time (dS18,S1) of the travel between the Sousse and Monastir stations is 2315 s, whereas the maximum time (fS18,S1) is 3085 s.


#### 4.3.5. Tested Scenario: Delay Departure of the Metro from the Sousse Bab Jdid Station

In the studied railway network, the detection function monitors the system evolution through the verification of time constraints. Let us suppose that(9)Φes1⟶es18=∑i=018Φi,where(10)Φ1=71s,Φ2=118s,Φ3=83s,Φ4=178s,Φ5=88s,Φ6=177s,Φ7=78s,Φ8=177s,Φ9=85s,Φ10=236s,Φ11=83s,Φ12=175s,Φ13=81s,Φ14=116s,Φ15=80s,Φ16=357s,Φ17=80s,Φ18=237s,Φ19=78s.From an experimental point of view, at the end of each journey, the current real-time values are collected and stored in an SNCFT database. The only data, from which the identification is computed is, therefore, a sequence of vectors indicating the planned schedule and the actual train departure times at each station. In order to illustrate our approach, a small part of the real data has been extracted from Table3. This table shows the planned and measured times associated to Sousse-Mahdia direction. These measurements have been recorded at the month of June, 2018.Table 3
Planned and measured times: Sousse to Mahdia direction.
StationPlanned timeReal timeSousse Bab Jdid05 : 40 : 0005 : 41 : 22Sousse Mohamed V05 : 42 : 0005 : 43 : 43Sousse Sud05 : 45 : 0005 : 48 : 28Sousse zone industrielle05 : 48 : 0005 : 52 : 34Sahline Ville05 : 51 : 0005 : 55 : 33Sahline Sabkha05 : 54 : 0005 : 57 : 33Les Hôtels05 : 58 : 0006 : 01 : 01Aeroport06 : 00 : 0006 : 02 : 01La faculté06 : 06 : 0006 : 08 : 00Monastir06 : 10 : 0006 : 21 : 25La faculté 206 : 24 : 0006 : 24 : 24Monastir zone industrielle06 : 26 : 0006 : 28 : 24Frina06 : 28 : 0006 : 30 : 05Khnis Bembla06 : 30 : 0006 : 33 : 24Ksibet Bennane06 : 34 : 0006 : 37 : 23Bouhadjar06 : 37 : 0006 : 40 : 23Lamta06 : 39 : 0006 : 43 : 22Sayada06 : 41 : 0006 : 44 : 43Ksar Hellal Z.I.06 : 43 : 0006 : 47 : 03Ksar Hellal06 : 45 : 0006 : 48 : 03Moknine Gribaa06 : 47 : 0006 : 50 : 02Moknine06 : 50 : 0006 : 52 : 03Téboulba Z.I.06 : 53 : 0006 : 56 : 20Téboulba06 : 56 : 0006 : 59 : 21Bekalta07 : 00 : 0007 : 06 : 02Baghdadi07 : 10 : 0007 : 13 : 04Mahdia Z.T.07 : 15 : 0007 : 18 : 03Sidi Massaoud07 : 18 : 0007 : 20 : 34Borj Arif07 : 21 : 0007 : 22 : 33Ezzahra07 : 24 : 0007 : 26 : 54Mahdia07 : 30 : 0007 : 28 : 54Let us suppose a late departure of the metro from the Sousse Bab Jdid station (departure at 05 : 41 : 22; see Table3), observed by the monitor S3: Φes1⟶es3 = 620 s. This delay may involve an illegal behavior and can lead to a degraded service. In fact, according to Figure 10, we deduce that the possibility of railway traffic delay is of 0.39 (possibility of violation of the global constraint CS18,S1). This delay can affect the stability of the studied railway network: according to Table 3, the metro arrives in the Monastir station with a delay of 685 s: violation of the constraint CS18,S1 (planned arrival time 06 : 10 : 00/measured arrival time at 06 : 21 : 25). Consequently, the distributed monitoring, based on chronicle, allows an early detection of traffic disturbance, to avoid catastrophic scenarios and preserve stability and security of the studied railway networks.Figure 10
Possibility function considering Φes1⟶es3.


## 4.1. Principle

Centralized monitoring consists of associating a single diagnoser to the whole process model. The latter collects the various process information before making a final decision on the process operating state. Although powerful in terms of diagnosis, the centralized structure is difficult to use for large and distributed systems such as transport systems since the building of a global model generates, in most cases, combinatorial explosion problems.In distributed monitoring structures, the process is decomposed into several local models. To each model is associated a local diagnoser. Each diagnoser makes his decision based on local observation. In the case of global specifications, a protocol allows the communication between the different diagnosers in order to make a final decision. The distributed monitoring, based on chronicles, verifies that each observed event is consistent with the chronicles specified by the time constraints. If the chronicles represent an abnormal behavior model, chronicles recognition is used for monitoring purposes.


## 4.2. Distributed Monitoring Architecture

In the studied railway networks, every set of sensors providing useful information has its own monitoring sites. Then, these sites are connected in distributed architecture. According to Figure6, the chronicles are distributed into several subchronicles associated to the different monitoring sites Si. Each diagnoser monitors locally a subsystem and communicates with other diagnosers in order to get necessary information in order to take decisions relative to the distributed diagnosis; Figure 8. Therefore, the detection function of each site has to recognize an evolution of the monitored process by subchronicles which are defined as sets of patterns of observable events temporally constrained. The recognition is performed through the verification of the time constraints associated to the considered subchronicles. The distributed detection function monitors the system evolution through the recognition of chronicles; Figure 8.Figure 8
Distributed monitoring architecture based on chronicles.


## 4.3. Distributed Monitoring of the Railway Sousse-Monastir

In order to help the supervisor in charge of managing the studied railway networks (i.e., detecting traffic perturbations, alert traveler claims, and maintain stability and security of the networks), a monitoring task is needed. The overall aim of the proposed monitoring approach is to control that the railway traffic proceeded well to avoid undesirable situations. The distributed monitoring of the Sahel railway networks based on the chronicle method involves the following steps.


### 4.3.1. Step 1: Subsystems Identification

The first task is to subdivide the railway network system into sectors (subsystems).The railway network between Sousse and Monastir has been split into 10 sectors, as shown in Figure9. Each sector is indicated by the dashed lines. These sectors areFigure 9
Travel time from Sousse to Monastir stations.


### 4.3.2. Step 2: Sensor Recognition

This step consists of identifying the sensors needed to perform the system monitoring. In the proposed approach, each siteSi, Figure 9, monitors a specific zone of the railway from Sousse to Monastir, through the observation of a set of events generated by the process. In the railway network, the travelling times should be within two bounds. Any deviations (occurrence of a temporal disturbance) from the allowed lower (resp. upper) bounds will lead to a low service quality and can lead to disaster scenarios. Therefore, the proposed monitoring approach, based on chronicle formalism, uses the additional information provided by the knowledge of interval constraints and allows detecting disruptions when a constraint is violated. Figure 9 gives journey times between the railway stations and parking times of the metro in stations. These durations are represented by the static intervals “ISi”.


### 4.3.3. Step 3: Distributed Detection Function

To each event is associated a diagnoserSi, Figure 9, and the detection function is distributed on the different sites:eS1: Beginning journey from the Sousse Bab Jdid station (Sensor S1).eS2: Metro arrival at Sousse Med.5 (Sensor S2).eS3: Departure from the subway station Sousse Med.5 (Sensor S3).eS4: Metro arrival at Sousse Sud Station (Sensor S4).eS5: Departure from the subway station Sousse Sud (Sensor S5).eS6: Metro arrival at Sousse Z.Ind Station (Sensor S6).eS7: Departure from the subway station Sousse Z.Ind (Sensor S7).eS8: Metro arrival at Sahline Ville station (Sensor S8).eS9: Departure from the subway station Sahline Ville (Sensor S9).eS10: Metro arrival at Sahline Sabkha station (Sensor S10).eS11: Departure from the subway station Sahline Sabkha (Sensor S11).eS12: Metro arrival at the Les Hôtels station (Sensor S12).eS13: Departure from the subway station Les Hôtels (Sensor S13).eS14: Metro arrival at the L’Aeroport station (Sensor S14).eS15: Departure from the subway station L’Aeroport (Sensor S15).eS16: Metro arrival at the La faculté station (Sensor S16).eS17: Departure from the subway station La faculté (Sensor S17).eS18: Metro arrival at the Monastir station (Sensor S18).


### 4.3.4. Step 4: Scenario Generation

A scenario is defined for each possible deviation in each section bringing the system to an erroneous situation. Deviations can occur due to unavailability of transportation systems, temporal disruptions, or traffic management and disturbances.Let us detail the monitoring between the two events,eS18 (metro arrival at the Monastir station) and eS1 (departure from Sousse Bab Jdid); Figure 9. To monitor this duration, it is necessary to check the time constraint linking the occurrences of the two events eS18 and eS1. This timing constraint is a global one; therefore, the verification of this constraint can be performed through the measure of the travelling time between the station and parking time of a metro in stations. As long as these durations are included in the mentioned intervals, no disturbance is detected. Otherwise, a traffic disturbance is detected. As previously mentioned, the global constraint to compute is an interval constraint type defined as CS18,S1: dS18,S1 ≤ O(eS18) − O(eS1) ≤ fS18,S1, with(8)dS18,S1=∑i=4563Li,fS18,S1=∑i=4563Hi.According to the time intervals (Figure3), the minimum time (dS18,S1) of the travel between the Sousse and Monastir stations is 2315 s, whereas the maximum time (fS18,S1) is 3085 s.


### 4.3.5. Tested Scenario: Delay Departure of the Metro from the Sousse Bab Jdid Station

In the studied railway network, the detection function monitors the system evolution through the verification of time constraints. Let us suppose that(9)Φes1⟶es18=∑i=018Φi,where(10)Φ1=71s,Φ2=118s,Φ3=83s,Φ4=178s,Φ5=88s,Φ6=177s,Φ7=78s,Φ8=177s,Φ9=85s,Φ10=236s,Φ11=83s,Φ12=175s,Φ13=81s,Φ14=116s,Φ15=80s,Φ16=357s,Φ17=80s,Φ18=237s,Φ19=78s.From an experimental point of view, at the end of each journey, the current real-time values are collected and stored in an SNCFT database. The only data, from which the identification is computed is, therefore, a sequence of vectors indicating the planned schedule and the actual train departure times at each station. In order to illustrate our approach, a small part of the real data has been extracted from Table3. This table shows the planned and measured times associated to Sousse-Mahdia direction. These measurements have been recorded at the month of June, 2018.Table 3
Planned and measured times: Sousse to Mahdia direction.
StationPlanned timeReal timeSousse Bab Jdid05 : 40 : 0005 : 41 : 22Sousse Mohamed V05 : 42 : 0005 : 43 : 43Sousse Sud05 : 45 : 0005 : 48 : 28Sousse zone industrielle05 : 48 : 0005 : 52 : 34Sahline Ville05 : 51 : 0005 : 55 : 33Sahline Sabkha05 : 54 : 0005 : 57 : 33Les Hôtels05 : 58 : 0006 : 01 : 01Aeroport06 : 00 : 0006 : 02 : 01La faculté06 : 06 : 0006 : 08 : 00Monastir06 : 10 : 0006 : 21 : 25La faculté 206 : 24 : 0006 : 24 : 24Monastir zone industrielle06 : 26 : 0006 : 28 : 24Frina06 : 28 : 0006 : 30 : 05Khnis Bembla06 : 30 : 0006 : 33 : 24Ksibet Bennane06 : 34 : 0006 : 37 : 23Bouhadjar06 : 37 : 0006 : 40 : 23Lamta06 : 39 : 0006 : 43 : 22Sayada06 : 41 : 0006 : 44 : 43Ksar Hellal Z.I.06 : 43 : 0006 : 47 : 03Ksar Hellal06 : 45 : 0006 : 48 : 03Moknine Gribaa06 : 47 : 0006 : 50 : 02Moknine06 : 50 : 0006 : 52 : 03Téboulba Z.I.06 : 53 : 0006 : 56 : 20Téboulba06 : 56 : 0006 : 59 : 21Bekalta07 : 00 : 0007 : 06 : 02Baghdadi07 : 10 : 0007 : 13 : 04Mahdia Z.T.07 : 15 : 0007 : 18 : 03Sidi Massaoud07 : 18 : 0007 : 20 : 34Borj Arif07 : 21 : 0007 : 22 : 33Ezzahra07 : 24 : 0007 : 26 : 54Mahdia07 : 30 : 0007 : 28 : 54Let us suppose a late departure of the metro from the Sousse Bab Jdid station (departure at 05 : 41 : 22; see Table3), observed by the monitor S3: Φes1⟶es3 = 620 s. This delay may involve an illegal behavior and can lead to a degraded service. In fact, according to Figure 10, we deduce that the possibility of railway traffic delay is of 0.39 (possibility of violation of the global constraint CS18,S1). This delay can affect the stability of the studied railway network: according to Table 3, the metro arrives in the Monastir station with a delay of 685 s: violation of the constraint CS18,S1 (planned arrival time 06 : 10 : 00/measured arrival time at 06 : 21 : 25). Consequently, the distributed monitoring, based on chronicle, allows an early detection of traffic disturbance, to avoid catastrophic scenarios and preserve stability and security of the studied railway networks.Figure 10
Possibility function considering Φes1⟶es3.


## 4.3.1. Step 1: Subsystems Identification

The first task is to subdivide the railway network system into sectors (subsystems).The railway network between Sousse and Monastir has been split into 10 sectors, as shown in Figure9. Each sector is indicated by the dashed lines. These sectors areFigure 9
Travel time from Sousse to Monastir stations.


## 4.3.2. Step 2: Sensor Recognition

This step consists of identifying the sensors needed to perform the system monitoring. In the proposed approach, each siteSi, Figure 9, monitors a specific zone of the railway from Sousse to Monastir, through the observation of a set of events generated by the process. In the railway network, the travelling times should be within two bounds. Any deviations (occurrence of a temporal disturbance) from the allowed lower (resp. upper) bounds will lead to a low service quality and can lead to disaster scenarios. Therefore, the proposed monitoring approach, based on chronicle formalism, uses the additional information provided by the knowledge of interval constraints and allows detecting disruptions when a constraint is violated. Figure 9 gives journey times between the railway stations and parking times of the metro in stations. These durations are represented by the static intervals “ISi”.


## 4.3.3. Step 3: Distributed Detection Function

To each event is associated a diagnoserSi, Figure 9, and the detection function is distributed on the different sites:eS1: Beginning journey from the Sousse Bab Jdid station (Sensor S1).eS2: Metro arrival at Sousse Med.5 (Sensor S2).eS3: Departure from the subway station Sousse Med.5 (Sensor S3).eS4: Metro arrival at Sousse Sud Station (Sensor S4).eS5: Departure from the subway station Sousse Sud (Sensor S5).eS6: Metro arrival at Sousse Z.Ind Station (Sensor S6).eS7: Departure from the subway station Sousse Z.Ind (Sensor S7).eS8: Metro arrival at Sahline Ville station (Sensor S8).eS9: Departure from the subway station Sahline Ville (Sensor S9).eS10: Metro arrival at Sahline Sabkha station (Sensor S10).eS11: Departure from the subway station Sahline Sabkha (Sensor S11).eS12: Metro arrival at the Les Hôtels station (Sensor S12).eS13: Departure from the subway station Les Hôtels (Sensor S13).eS14: Metro arrival at the L’Aeroport station (Sensor S14).eS15: Departure from the subway station L’Aeroport (Sensor S15).eS16: Metro arrival at the La faculté station (Sensor S16).eS17: Departure from the subway station La faculté (Sensor S17).eS18: Metro arrival at the Monastir station (Sensor S18).


## 4.3.4. Step 4: Scenario Generation

A scenario is defined for each possible deviation in each section bringing the system to an erroneous situation. Deviations can occur due to unavailability of transportation systems, temporal disruptions, or traffic management and disturbances.Let us detail the monitoring between the two events,eS18 (metro arrival at the Monastir station) and eS1 (departure from Sousse Bab Jdid); Figure 9. To monitor this duration, it is necessary to check the time constraint linking the occurrences of the two events eS18 and eS1. This timing constraint is a global one; therefore, the verification of this constraint can be performed through the measure of the travelling time between the station and parking time of a metro in stations. As long as these durations are included in the mentioned intervals, no disturbance is detected. Otherwise, a traffic disturbance is detected. As previously mentioned, the global constraint to compute is an interval constraint type defined as CS18,S1: dS18,S1 ≤ O(eS18) − O(eS1) ≤ fS18,S1, with(8)dS18,S1=∑i=4563Li,fS18,S1=∑i=4563Hi.According to the time intervals (Figure3), the minimum time (dS18,S1) of the travel between the Sousse and Monastir stations is 2315 s, whereas the maximum time (fS18,S1) is 3085 s.


## 4.3.5. Tested Scenario: Delay Departure of the Metro from the Sousse Bab Jdid Station

In the studied railway network, the detection function monitors the system evolution through the verification of time constraints. Let us suppose that(9)Φes1⟶es18=∑i=018Φi,where(10)Φ1=71s,Φ2=118s,Φ3=83s,Φ4=178s,Φ5=88s,Φ6=177s,Φ7=78s,Φ8=177s,Φ9=85s,Φ10=236s,Φ11=83s,Φ12=175s,Φ13=81s,Φ14=116s,Φ15=80s,Φ16=357s,Φ17=80s,Φ18=237s,Φ19=78s.From an experimental point of view, at the end of each journey, the current real-time values are collected and stored in an SNCFT database. The only data, from which the identification is computed is, therefore, a sequence of vectors indicating the planned schedule and the actual train departure times at each station. In order to illustrate our approach, a small part of the real data has been extracted from Table3. This table shows the planned and measured times associated to Sousse-Mahdia direction. These measurements have been recorded at the month of June, 2018.Table 3
Planned and measured times: Sousse to Mahdia direction.
StationPlanned timeReal timeSousse Bab Jdid05 : 40 : 0005 : 41 : 22Sousse Mohamed V05 : 42 : 0005 : 43 : 43Sousse Sud05 : 45 : 0005 : 48 : 28Sousse zone industrielle05 : 48 : 0005 : 52 : 34Sahline Ville05 : 51 : 0005 : 55 : 33Sahline Sabkha05 : 54 : 0005 : 57 : 33Les Hôtels05 : 58 : 0006 : 01 : 01Aeroport06 : 00 : 0006 : 02 : 01La faculté06 : 06 : 0006 : 08 : 00Monastir06 : 10 : 0006 : 21 : 25La faculté 206 : 24 : 0006 : 24 : 24Monastir zone industrielle06 : 26 : 0006 : 28 : 24Frina06 : 28 : 0006 : 30 : 05Khnis Bembla06 : 30 : 0006 : 33 : 24Ksibet Bennane06 : 34 : 0006 : 37 : 23Bouhadjar06 : 37 : 0006 : 40 : 23Lamta06 : 39 : 0006 : 43 : 22Sayada06 : 41 : 0006 : 44 : 43Ksar Hellal Z.I.06 : 43 : 0006 : 47 : 03Ksar Hellal06 : 45 : 0006 : 48 : 03Moknine Gribaa06 : 47 : 0006 : 50 : 02Moknine06 : 50 : 0006 : 52 : 03Téboulba Z.I.06 : 53 : 0006 : 56 : 20Téboulba06 : 56 : 0006 : 59 : 21Bekalta07 : 00 : 0007 : 06 : 02Baghdadi07 : 10 : 0007 : 13 : 04Mahdia Z.T.07 : 15 : 0007 : 18 : 03Sidi Massaoud07 : 18 : 0007 : 20 : 34Borj Arif07 : 21 : 0007 : 22 : 33Ezzahra07 : 24 : 0007 : 26 : 54Mahdia07 : 30 : 0007 : 28 : 54Let us suppose a late departure of the metro from the Sousse Bab Jdid station (departure at 05 : 41 : 22; see Table3), observed by the monitor S3: Φes1⟶es3 = 620 s. This delay may involve an illegal behavior and can lead to a degraded service. In fact, according to Figure 10, we deduce that the possibility of railway traffic delay is of 0.39 (possibility of violation of the global constraint CS18,S1). This delay can affect the stability of the studied railway network: according to Table 3, the metro arrives in the Monastir station with a delay of 685 s: violation of the constraint CS18,S1 (planned arrival time 06 : 10 : 00/measured arrival time at 06 : 21 : 25). Consequently, the distributed monitoring, based on chronicle, allows an early detection of traffic disturbance, to avoid catastrophic scenarios and preserve stability and security of the studied railway networks.Figure 10
Possibility function considering Φes1⟶es3.


## 5. Conclusions

In this study, we investigate the monitoring of a railway network. The proposed monitoring architecture is a distributed architecture, based on monitoring sites. To each site is associated a subsystem receiving a set of events and providing a detection function. Abnormal behaviors and traffic disturbances are recognized with the cooperation of local detection functions. The verification of the global constraints supposes the existence of communication tools allowing event exchanges between monitoring sites. The chronicle recognition is exploited for the detection of traffic disturbances.The problem of the distributed recognition of subchronicles through the verification of local and global time constraints has been pointed out. This recognition is based on time constraints verification performed with a possibility evaluation. The results obtained for the illustrative example are promising. They show that the distributed monitoring improves the prevention of temporal disruption and traffic management by performing an early detection.In our future works, a comparative study based upon several cases should be developed. A comparison with the proposed monitoring architecture and results with the works based on Causal Time Signature (CTS) should also be considered. Likewise, it would be interesting to apply a robust control strategy facing disturbances in railway transport systems in order to compensate the disturbances once they have been detected [21].



---

*Source: 6265379-2020-08-25.xml*