You are viewing a javascript disabled version of the site. Please enable Javascript for this site to function properly.
Go to headerGo to navigationGo to searchGo to contentsGo to footer
In content section. Select this link to jump to navigation

Matching method for emergency plans of highway traffic based on fuzzy sets and rough sets

Abstract

To deal with the problem in emergency plan matching of highway traffic that incident description is incomplete, incident properties are unclear, and plan matching is inaccurate, etc., a plan matching method is proposed based on fuzzy sets and rough sets. The property weight calculation method based on rough sets is used to reduce the dependency on prior knowledge; the structural similarity calculation is used to solve the problem of property missing and matching angle varying; the fuzzy set calculation method is adopted to solve the problem of fuzzy property similarity missing. The traffic emergency plan matching case for Changzhou section of Shanghai-Nanjing highway demonstrates that the proposed method can improve the accuracy and reliability of highway traffic emergency matching, reflecting the advantages of rough sets and fuzzy sets in emergency plan matching.

1Introduction

With the increase of vehicle traffic on highways, traffic emergencies continue to happen. It is an engineering problem urgently needed solving how to reduce the response time of traffic emergencies and increase the rescue efficiency. Practice has shown that the real-time acquisition of emergency plans is the key measure of reducing rescue decision time. Proposing automatic matching methods for emergency plans and developing corresponding application systems have important scientific and engineering value for improving emergency rescue decision efficiency.

Zografos [5] established a decision support system for emergency response to aid the dispatch decision of rescue resources. Logi, et al. designed a decision support system for traffic congestion incidents to implement the disposal decisions under congestion conditions. Alvear, et al. [3] designed a decision support system for highway tunnel emergency management to realize the command and decision of highway tunnel under incomplete information conditions. Konstantinos [6] developed a decision support system of the rescue resource scheduling for highway emergencies, which defines the rescue procedures for emergency incidents, and designs relevant databases and model bases to realize the optimal scheduling in rescue process. An emergency plan system was developed for Shanghai World Expo to improve the emergency response capabilities. Lu, et al. [9] adopted case-based reasoning (CBR) technology to develop an airport emergency plan management system. Yang, et al. [7] introduced CBR technology into highway emergency management and designed the retrieve, reuse, correct and save modules for cases, realizing plan matching based on nearest neighbor searching. There are still aspects that need improving in the above method and application system: 1) The property weights need manually determination in the plan matching process; 2) The information of historical cases is incomplete, the description of emergency incident information is incomplete, and the property missing affects the reliability of casematching; 3) The fuzziness of plan properties is not considered in the automatic retrieval of plans.

In this study, structural similarity is introduced to analyze the property structural characteristics of emergency incidents and historical cases. Different similarity calculation methods are adopted for different properties. The fuzziness of character properties is described by fuzzy sets and the similarity calculation methods of fuzzy properties are developed. The rough set theory is adopted to determine the property weights of emergency incidents, and design the algorithm process of emergency plan matching. We aim to propose an intelligent matching method for traffic emergency plan of freeway without prior knowledge, which can overcome the matching error caused by the uncertainty of emergency incidents to some degree, improving the accuracy and reliability of emergency plan matching.

2Overview of highway traffic emergency plan matching

2.1Incident character properties

According to the requirements of emergency rescue decision for highway traffic, emergency cases are divided into five parts, namely case character information, incident character properties, emergency command departments, decision plans and effect evaluation. The part playing a key role in rescue decision is incident character properties, including basic incident information and incident loss information. Basic incident information describes the traffic conditions and environmental conditions when the incident happens. Incident loss information describes the seriousness of emergency incidents, including fallouts, obstacles, the number of blocked lanes, the number and type of overturned vehicles, accident casualties, damaged condition of incident vehicles, etc. Incident character properties are divided into four classes [2]:

  • Numerical properties: expressed by determined numbers, such as death toll, number of injured, etc.;

  • Symbolic properties: expressed by clear terms, such as time of incident, site of incident, etc.;

  • Fuzzy conceptual properties: conceptual variables are the properties. A term set consists of all homogeneous fuzzy variable properties, with each fuzzy variable property corresponding to each term in the set, such as traffic congestion condition,etc.;

  • Fuzzy interval or fuzzy number properties: an interval without determined boundaries or an undetermined number. When the left and right boundaries of the fuzzy interval are equal, the fuzzy interval reduces to a fuzzy number. For example, accident duration, affecting scope, etc.

In this study, 24 highway incident character properties are extracted. Among them, there are four numerical properties shown in Table 1; fourteen symbolic properties shown in Table 2; four fuzzy conceptual properties shown in Table 3; and two fuzzy interval properties shown in Table 4.

2.2Matching method for traffic emergency plans

The traffic emergency plan matching is a calculation and comparing process of case similarity. By establishing a similarity measurement function, the current emergency incident is compared with the historical case in the case library. Then the case which is most similar to the current incident is retrieved and treated as the most suitable rescue plan for the current emergency incident. Therefore, choosing an appropriate similarity calculation method is the key factor of plan matching.

If we want to determine the similarity between an emergency incident and a historical case, we have to first determine the similarity between incident character properties. At present, commonly used similarity calculation method for character properties is nearest neighbor method [1], which is based on normalized distances between case properties. This method is inadequate for similarity calculation of qualitative case property description. For example, for qualitative properties “general” and “serious” levels of incidents, the similarity calculation result is 0 (which means dissimilar). However there is similarity between them, and only the extents of casualties and property losses are different. Other similarity calculation methods simply equalize fuzzy properties and determined properties, causing the loss of similarity calculation for qualitative description properties. Due to the uncertainty, complexity and fuzziness of highway traffic incidents, the character properties of incident sites are mainly fuzzy values. Therefore, using fuzzy sets to calculate similarity between properties can fully reveal the similarity between properties.

3Similarity calculation for emergency character properties

3.1Character weight calculation based on rough set theory

The importance of incident character properties in plan matching process is expressed by the weights in similarity calculation. Whether the property weights are reasonable directly affects the accuracy and reliability of plan matching. Current weigh assigning methods are mainly based on domain experts’ prior knowledge, which mainly depends on subjective decision and experience. In this study, rough set theory is used to determine each property weight of incidents. Only the data are compared and no other prior knowledge is needed. Therefore the method is objective [8]. Moreover, rough set theory has the ability to solve minimal invariant set (called kernel) and minimal rule sets (called reduction) in large-scale data. This ability can be used to simplify redundant properties and property values. The character properties of highway emergency incidents include quantitive properties and qualitative properties. The quantitive properties are firstly discretized. The discretization method is referred to literature [4].

According to the rough set theory, the knowledge expression system S = (U,  R,  F,  V) is designed, where, U is the incident historical case set; R = {x 1,  x 2, …, x n } is the incident feature property set; C = {c 1,  c 2, …,  c n-1} is the condition property set; D = {d} is the decision property set. Decision property is usually the “incident level” and R = C ∪ D. The dependency of decision property set D on condition property set C is given:

(1)
rc=Card(posc(D))Card(U)
where Card represents the cardinality of a set; pos c  (D) represents the C positive region of D.

The importance r c-i of the i-th property is given:

(2)
rc-i=Card(posc-i(D))Card(U)
where pos c-i  (D) represents the positive region of decision property set D in condition property set C with the condition property subset {c i } removed.

The objective weigh w i1 of the i-th property is given:

(3)
wi1=rc-rc-ii=1i=n-1rc-rc-i

The larger w i1 is, the more important the i-th property is to decisions. When property i is removed from condition property set, the positive region of class U/D is largely affected; when w i1 = 0, i.e., r c  - r c-i  = 0, which means removing the i-th property does not affect the positive region of class U/D. Therefore this property is redundant and should be removed.

The integrated weight w i of the i-th property is expressed by Equation (4):

(4)
wi=α×wi1+(1-α)×wi2
where, w i1 and w i2 are the objective weight and subjective weight, respectively, which are given by domain experts according to real situations; α is an empirical factor. The larger α is, the more importance is attached to objective weight, and the smaller α is, the more importance is attached to subjective weight.

3.2Calculating similarity between emergency incidents and historical case structures

The case properties corresponding to different incidents have different composition structures. There are numerous descriptive properties of highway emergency incidents. Therefore, it is inevitable that historical case information or case description is incomplete. Moreover, using structural similarity to choose historical cases not only increases the similarity between the chosen case and the emergency incident property, but also solves the property value missing problem, making the formed rescue plan more timely and useful [10].

Let the current emergency incident be P, and a certain historical case in case library be Q The structural similarity S (P,  Q) equals to the weight sum of the intersection set property of P and Q, divided by the weight sum of the union set property of P and Q, shown in Equation (5):

(5)
S(P,Q)=WPQWPQ=WGWH=a=1a=gwab=1b=hwb
where w a is the weight of the a–th property in set G = P ∩ Q; g is the number of elements in set G; w b is the weight of the b–th property in set H = P ∪ Q; h is the number of elements in set D.

3.3Property similarity calculation based on fuzzy sets

Fuzzy conceptual properties, fuzzy intervals and fuzzy numerical properties are collectively known as fuzzy properties. The similarities between fuzzy properties are special relation functions, such as trigonometric functions, trapezoidal function, Gaussian function, etc. The trapezoidal function is used in this study, the membership function S M  (x) of element x in fuzzy set is expressed [2]:

(6)
SM(x)={max(0,1-vl-xm)xvl1vlxvrmax(0,1-x-vrn)xvr
where M is the fuzzy set in discourse domain N; v l and v r are the lower and upper bounds of N, respectively; m and n are different with respect to properties. For fuzzy conceptual properties, m and n are usually determined by experts; for fuzzy intervals or fuzzy numbers, m and n are usually kv l and kv r , respectively, and k is usually 0.1.

The similarity between the fuzzy properties of incident P and historical case Q is calculated by summing areas. That is, to calculate the area overlapping ratio between the membership functions of two fuzzy properties, which is expressed by Equation (8):

(7)
sim(pi,qi)=A(piqi)A(pi)+A(qi)-A(piqi)
where p i , q i are the fuzzy sets corresponding to the properties of incident P and historical case Q, respectively; sim (p i ,  q i ) is the fuzzy similarity between two fuzzy sets p i and q i ; A represents the area of corresponding fuzzy set; p i  ∩ q i represents the intersection set of the two fuzzy sets p i and q i .

Point (pi*,qi*) is the intersection of the two fuzzy sets. pi* and qi* are shown in Equations (9) and (10), respectively [2]:

(8)
pi*=(npivlqi+mqivrpi)(npi+mqi)
(9)
qi*=1-(pi*vspi)/npi

If qi*<0, there is no intersection, sim (p i ,  q i ) =0.

For similarity between numbers, a linear function is adopted for normalization, which is expressed:

(10)
sim(pi,qi)=1-|pi-qimax(i)-min(i)|
where sim (p i ,  q i ) represents the similarity of the i–th numerical property for incident P and historical case Q. p i and q i represent the i–th property values of P and Q, respectively. max(i) and min(i) represent the maximum and minmum value of the i-th property, respectively.

Since there are only property values for symbolic properties and there is no real sense of similarity, the similarity measurement for symbolic properties only takes the value of 1 for identical properties and 0 for different properties, as shown:

(11)
sim(pi,qi)={1pi=qi0piqi
where sim (p j ,  q j ) is the similarity of the i–th symbolic property for incident P and historical case Q; p i and q i represent the i–th property value of incident P and historical case Q, respectively.

3.4Overall similarity calculation

Integrating structural similarity, property similarity and corresponding weights, we obtain the overall similarity shown in Equation (13):

(12)
sim(P,Q)=S(P,Q)×i=1nwisim(pi,qi)
where Sim (P,  Q) is the overall similarity; S (P,  Q) is the structural similarity which is obtained from Equation (5); sim (p i ,  q i ) is the similarity the similarity between i-th properties, calculated by Equations (4), (6) and (7), ∑w i  = 1, i = 1, …,  n, n is the total number of properties.

According to maximum similarity reuse strategy, select the case with maximum Sim (P,  Q) as the best matched emergency plan.

4Case analysis

According to the some cases in the emergency case library of Changzhou section of Shanghai-Nanjing highway, the proposed similarity calculation method is used to match exemplar plans, as shown in Table 5.

  • Determine the property weights

According to Equations (3) and (4), the property weights are:

W = {0.055, 0.034, 0.047, 0.047, 0.052, 0.050, 0.030, 0.028, 0.034, 0.036, 0.043, 0.048, 0.050, 0.053, 0.051, 0.033, 0.054, 0.056, 0.047, 0.038, 0.038, 0.037, 0.039}, and no redundant property exists.

  • - Calculate the structural similarity

According to Equation (5):

S(A,B1)=WAB1WAB1=a=1a=pwab=1b=pwb=0.9080.938=0.968

The calculation of S (A, B2) and S (A, B3) is similar. S (A, B2) = S (A, B3) =1.

It is clear that the property structure of emergency A is the same as historical cases B2 and B3, and is also similar to B1.

  • Calculate the similarity between properties

* The similarity calculation steps for fuzzy conceptual properties are as follows:

Use daily traffic congestion index to reflect the road network traffic congestion level. The road congestion is classified into 5 classes, with a value range of [0, 10], as shown in Table 6. The similarity of road congestion level property between incident A and historical case B1 is calculated. Since there are only two situations in Table 5, namely clear and minor congestion, only the similarity between (2,4] and (4,6] needs calculation.

Road congestion level of emergency A correspond to fuzzy set p 1, historical case B1 is q 1. Then v l p1  = 4,and v r p1  = 6. Since k = 0.1, m p1 = kv lp1 = 0.4, n p1 = kv rp1 = 0.6. Similarly, v lq 1  = 2, v rq 1  = 4, with m q1 = 0.2, n q1 = 0.4. According to Equations (8) and (9) p1*=4, q1*=1. The functions of p 1 and q 1 are shown in Fig. 1.

According to Equation (7):

sim(p1,q1)=A(p1q1)A(p1)+A(q1)-A(p1q1)=0.0909

Similarly, the similarity between the road congestion level of A and B1, B2, B3 is calculated. The results are shown in Table 7.

  • * The calculation of the similarity of fuzzy interval properties, is in accordance with the similarity calculation of fuzzy properties. The similarity results are shown in Table 8.

  • * The similarity calculation for numerical properties and symbolic properties is according to Equations (11) and (12), as shown in Tables 9 and 10.

  • * The overall similarity is calculated according to Equation (12) and shown in Table 11:

According to Table 11, the similarity between incident A and historical case B2 is larger than that of A and B1, A and B3, which indicates that the solution plan of historical case B2 should be the rescue plan of emergency incident A.

5Conclusions

The proposed method reduces the dependency on prior knowledge and eliminates redundant properties. The integrated structural similarity calculation for emergency incidents and historical cases improves the reliability of case structure matching and reduces the effect of property missing on similarity calculation. The traffic incidents are classified into different types according to their properties, and different calculation methods are adopted for different tupes. The fuzzy set is specially used to calculate the similarity between fuzzy properties, solving the previous problem of similarity loss. The proposed highway traffic emergency plan matching method is able to adapt to the uncertainty of incident description in plan matching, improving the accuracy of plan matching.

Considering it is difficult to collect the character of freeway emergency incident, we will focus on how to deal with emergency plan matching with a large number of event properties missing.

References

1 

Yan A, Qian L, Wang P (2014) A comparative study of attribute weights assignment for case-based reasoning Acta Automatica Sinica 40: 9 1896 1902

2 

Zhang B, Yu Y (2002) Hybrid similarity measure for retrieval in case-based reasoning system Systems Engineering-Theory & Practice 22: 3 131 136

3 

Alvear D, Abreu O, Cuesta A, Alonso V (2013) Decision support system for emergency management: Road tunnels Tunnelling and Underground Space Technology 34: 13 21

4 

Miao D (2001) A new method of discretization of continuous attributes in rough sets Acta Automatica Sinica 5: 296 302

5 

Zografos GK, Vasilakis MG, Giannouli MI (2000) Methodological framework for developing decision support systems for hazardous materials emergency response operations Journal of Hazardous Materrals 71: 503 521

6 

Konstantinos GZ, Konstantinos NA, George MV (2002) A real-time decision support system for roadway network incident response logistics Transportation Research Part C 10: 1 18

7 

Yang S, Xiong C, Ni F (2008) Freeway incident management system based on case-based reasoning Journal of Southeast University 38: 5 878 883

8 

Cao X, Liang J (2002) Determination method of character weight based on rough sets theory Chinese Journal of Management Science 10: 5 98 100

9 

Lu Z, Hu X, Yang H (2010) Airport emergency preparednesss management based on case-based reasoning Journal of Civil Aviation University of China 28: 1 5

10 

Zhang Z, Wang X, Dang J, Min Y (2012) Research on CBR decision method of railway emergency rescue based on integral similarity degree Journal of the China Railway Society 34: 11 49 53

Figures and Tables

Fig.1

Functions of fuzzy sets p 1 and q 1.

Functions of fuzzy sets p
1 and q
1.
Table 1

Numerical properties

Property nameProperty value
No. of damaged vehicles(determined number)
No. of injured(determined number)
Death toll(determined number)
No. of retained people(determined number)
Table 2

Symbolic properties

Property nameProperty value
Incident levelexceptional serious level I, serious level II,
   general level III, minor level IV
Incident data2008-12-31 (Specified date)
Incident sectionLianshui to Gaoyou (Specified section)
Incident timeworkday daytime (night), holiday daytime
   (night), activity day daytime (night)
Weather conditionsunny, foggy, rainy, snowy
Road conditiondry, wet, icy
Road occupancyone lane, two lanes, three lanes
Lane blockingbidirectional, unidirectional
   condition
Vehicle typemotorcycle, bus, sedan, van, truck,
   trailer, non-motor vehicle
Falloutyes, no
Road property lossyes, no
Hazardousyes, no
   chemical leak
Fireyes, no
Explosionyes, no
Table 3

Fuzzy conceptional properties

Property nameProperty value
Precipitationlight rain, moderate rain, heavy rain, storm
Snow amountlight snow, moderate snow,
   heavy snow, blizzard
Icing arealarge, small
Traffic congestionbusy, blocked, unblocked
Table 4

Fuzzy interval properties

Property nameProperty value
Visibility>500 m, 200–500 m, 100–200 m,
   50–100 m, <50m
Incident duration[1,3]h
Table 5

Property values of emergency incidents and historical cases

Property nameCase ACase B1Case B2Case B3
Incident levellevel IIIlevel IIIlevel IIIlevel IV
Incident date2015-01-202015-01-112014-12-022014-12-24
Incident sectionChangzhou-ZhenjiangChangshu-ChangzhouSuchou-ChangzhouDanyang-Changzhou
Incident timeworkday daytimeholiday daytimeholiday nightactivity day daytime
Weather conditionsunnysunnysunnyfoggy
Road conditiondrydrydrywet
Precipitation———-moderate———-———
Snow amount———-———-———-———
Icing area———-———-———-———-
Visibility300–400m200–300m300–400m100–200m
Incident duration[2,4]h[1,3]h[1,4]h[1,3]h
Road occupancyone laneone laneone laneone lane
Lane blocking conditionunidirectionalunidirectionalunidirectionalunidirectional
Traffic congestionminor congestionunblockedminor congestionunblocked
No. of damaged vehicles2230
Vehicle typesedansedansedanbus
No. of injured2421
Death toll0010
No. of retained people1120
Falloutnonoyesno
Road property lossnoyesnono
Firenononono
Explosionnononono
Table 6

Traffic congestion levels

Daily traffic[0,2](2,4](4,6](6,8](8,10]
congestion index
Traffic congestion levelVery clearclearminor congestionmedium congestionserious congestion
Table 7

Similarity of road congestion property

Property nameA and B1A and B2A and B3
Road congestion0.090910.0909
Table 8

Fuzzy property similarity

Property nameA and B1A and B2A and B3
Visibility0.130410
Incident duration0.38460.70770.3486
Table 9

Similarity of numeric properties

Property nameA and B1A and B2A and B3
No. of damaged vehicles10.6670.333
No. of injured0.33310.667
Death toll101
No. of retained people10.50.5
Table 10

Similarity of symbolic properties

Property nameA and B1A and B2A and B3
Incident level110
Incident date000
Incident section000
Incident time000
Weather condition110
Road conditions110
Road occupancy111
Lane blocking condition111
Vehicle type110
Fallout101
Road property loss111
Fire111
Explosion111
Table 11

Overall similarity

NameA and B1A and B2A and B3
Overall similarity0.60190.63290.4023