Real-time, Adaptive Prediction of Incident Delay for Advanced Traffic Management Systems

Similar documents
OVERVIEW OF A SIMULATION STUDY OF THE HIGHWAY 401 FTMS

CHAPTER 5 DELAY ESTIMATION FOR OVERSATURATED SIGNALIZED APPROACHES

Impact of Day-to-Day Variability of Peak Hour Volumes on Signalized Intersection Performance

Signalized Intersection Delay Models

Assessing the uncertainty in micro-simulation model outputs

Incorporating the Effects of Traffic Signal Progression Into the Proposed Incremental Queue Accumulation (IQA) Method

FUZZY TRAFFIC SIGNAL CONTROL AND A NEW INFERENCE METHOD! MAXIMAL FUZZY SIMILARITY

Delay of Incidents Consequences of Stochastic Incident Duration

$QDO\]LQJ$UWHULDO6WUHHWVLQ1HDU&DSDFLW\ RU2YHUIORZ&RQGLWLRQV

A Speed-Delay Study of Highway 401 in Toronto, Ontario

Cumulative Count Curve and Queueing Analysis

Signalized Intersections

Capacity Drop. Relationship Between Speed in Congestion and the Queue Discharge Rate. Kai Yuan, Victor L. Knoop, and Serge P.

A STUDY OF TRAFFIC PERFORMANCE MODELS UNDER AN INCIDENT CONDITION

Performance Analysis of Delay Estimation Models for Signalized Intersection Networks

I. M. Schoeman North West University, South Africa. Abstract

Estimating the vehicle accumulation: Data-fusion of loop-detector flow and floating car speed data

Appendix BAL Baltimore, Maryland 2003 Annual Report on Freeway Mobility and Reliability

The Next Generation of Traffic Management Systems

Signalized Intersection Delay Models

Real-Time Travel Time Prediction Using Multi-level k-nearest Neighbor Algorithm and Data Fusion Method

Traffic Impact Study

Control Delay at Signalized Diamond Interchanges Considering Internal Queue Spillback Paper No

A Cellular Automaton Model for Heterogeneous and Incosistent Driver Behavior in Urban Traffic

Parking Regulations Dundas Street West, from Bathurst Street to Dovercourt Road

Accident Prediction Models for Freeways

Progression Factors in the HCM 2000 Queue and Delay Models for Traffic Signals

CONGESTION REPORT 1 st Quarter 2016

CE351 Transportation Systems: Planning and Design

NATHAN HALE HIGH SCHOOL PARKING AND TRAFFIC ANALYSIS. Table of Contents

Signalized Intersection Delay Models

CHAPTER 2. CAPACITY OF TWO-WAY STOP-CONTROLLED INTERSECTIONS

CHAPTER 3. CAPACITY OF SIGNALIZED INTERSECTIONS

Analytical Delay Models for Signalized Intersections

Individual speed variance in traffic flow: analysis of Bay Area radar measurements

Effect of Environmental Factors on Free-Flow Speed

Real-Time Crash Prediction Model for Application to Crash Prevention in Freeway Traffic

Real-Time Crash Prediction Model for the Application to Crash Prevention in Freeway Traffic

Empirical Analysis of Traffic Sensor Data Surrounding a Bottleneck on a German Autobahn

Travel Time Estimation with Correlation Analysis of Single Loop Detector Data

Traffic Flow Theory and Simulation

Appendixx C Travel Demand Model Development and Forecasting Lubbock Outer Route Study June 2014

Transportation and Road Weather

1.225 Transportation Flow Systems Quiz (December 17, 2001; Duration: 3 hours)

Uncertainty in the Yellow Change Interval

WRTM and Active Transportation and Demand Management

Anticipatory Pricing to Manage Flow Breakdown. Jonathan D. Hall University of Toronto and Ian Savage Northwestern University

TEEN DRIVER ELECTRONIC DEVICE OBSERVATION FORM

US 169/I-70 North Loop Planning & Environmental Linkages Study

Publication List PAPERS IN REFEREED JOURNALS. Submitted for publication

c) What are cumulative curves, and how are they constructed? (1 pt) A count of the number of vehicles over time at one location (1).

AN INVESTIGATION INTO THE IMPACT OF RAINFALL ON FREEWAY TRAFFIC FLOW

Market Street PDP. Nassau County, Florida. Transportation Impact Analysis. VHB/Vanasse Hangen Brustlin, Inc. Nassau County Growth Management

Using High-Resolution Detector and Signal Data to Support Crash Identification and Reconstruction. Indrajit Chatterjee Gary Davis May, 2011

Development and Evaluation of Online Estimation Methods for a Feedback-Based Freeway Ramp Metering Strategy

Traffic Flow Impact (TFI)

GIS-BASED VISUALIZATION FOR ESTIMATING LEVEL OF SERVICE Gozde BAKIOGLU 1 and Asli DOGRU 2

JEP John E. Jack Pflum, P.E. Consulting Engineering 7541 Hosbrook Road, Cincinnati, OH Telephone:

Estimation of Measures of Effectiveness Based on Connected Vehicle Data

TEEN DRIVER SEAT BELT OBSERVATION FORM

Variable Speed Approach for Congestion Alleviation on Boshporus Bridge Crossing

Traffic Flow Theory & Simulation

New Calculation Method for Existing and Extended HCM Delay Estimation Procedures

Automated Delay Estimation at Signalized Intersections: Phase I Concept and Algorithm Development

TREBALL FINAL DE GRAU

Yu (Marco) Nie. Appointment Northwestern University Assistant Professor, Department of Civil and Environmental Engineering, Fall present.

Existence, stability, and mitigation of gridlock in beltway networks

TRB Traffic Flow on Freeway Upgrades. By Werner Brilon and Andrea Bressler

Draft Final Report. SHRP2 Project L08 Incorporation of Travel Time Reliability into the Highway Capacity Manual. Prepared For: SHRP2

CONTINUING PLANNING PROGRAM LOCAL TRANSPORTATION AND TRAFFIC DATA PRODUCT REPORT [OH Corridors]

Investigating uncertainty in BPR formula parameters

Use of Performance Measurement System Data to Diagnose Freeway Bottleneck Locations Empirically in Orange County, California

Available online Journal of Scientific and Engineering Research, 2017, 4(4): Research Article

Texas A&M University

Path and travel time inference from GPS probe vehicle data

Traffic flow theory involves the development of mathematical relationships among

Validation of ATIS Journey Time and Traffic Speed Estimates by Floating Car Survey

Freeway Merge Junctions

Course Outline Introduction to Transportation Highway Users and their Performance Geometric Design Pavement Design

Prioritization of Freeway Segments Based on Travel Time Reliability and Congestion Severity

Traffic Surveillance from a Safety Perspective: An ITS Data Application

Estimation of Speeds from Single-Loop Freeway Flow and Occupancy Data Using Cusp Catastrophe Theory Model

Appendix C Final Methods and Assumptions for Forecasting Traffic Volumes

Airborne Traffic Flow Data and Traffic Management. Mark Hickman and Pitu Mirchandani

Conception of effective number of lanes as basis of traffic optimization

PATREC PERSPECTIVES Sensing Technology Innovations for Tracking Congestion

Reliability and Safety Prediction for Planning

Cipra D. Revised Submittal 1

Signalized Intersection Delay Models

2.1 Traffic Stream Characteristics. Time Space Diagram and Measurement Procedures Variables of Interest

Queuing and Statistical Analysis of Freeway Bottleneck Formation

EMPIRICAL STUDY OF TRAFFIC FEATURES AT A FREEWAY LANE DROP

Background and Hong Kong Statistics. Background. Estimation of Network Reliability under Traffic Incidents for ITS Applications

Measuring Wave Velocities on Highways during Congestion using Cross Spectral Analysis

Tahlequah Public Schools Inclement Weather Transportation Plan

April 10, Mr. Curt Van De Walle, City Manager City of Castle Hills 209 Lemonwood Drive Castle Hills, Texas 78213

Using PeMS Data to Empirically Diagnose Freeway Bottleneck Locations in Orange

TRAVEL TIME RELIABILITY ON EXPRESSWAY NETWORK UNDER UNCERTAIN ENVIRONMENT OF SNOWFALL AND TRAFFIC REGULATION

Chapter 5 Traffic Flow Characteristics

A Regression Model for Bus Running Times in Suburban Areas of Winnipeg

Transcription:

Real-time, Adaptive Prediction of Incident Delay for Advanced Traffic Management Systems Liping Fu and Bruce Hellinga Department of Civil Engineering, University of Waterloo, Waterloo, Canada Phone: 59 885-2 Ext. 3984; Fax: 59 888-697; email: lfu@uwaterloo.ca SUMMARY This paper presents a fuzzy queuing model that can be used to predict the possible delay or interval of delay that a vehicle will experience at an incident location based on real-time information on current queuing condition, future traffic arrival, lane closing and the vehicle s arrival time. Compared to most existing methods, the proposed model is unique in three aspects: first, it explicitly accounts for uncertainties involved in all influencing factors and thus allows easy incorporation of imprecise and vague information typically available in this type of prediction environment. Second, the model is adaptive in the way that it allows continuous update of estimates as new information is made available in real time. Third, delays obtained from the model are fuzzy numbers that can be conveniently mapped to linguistic terms for use in systems such as changeable message signs (CMS). A case study is presented to demonstrate the application of the proposed model in facilitating the composition of location-dependent delay messages for CMS. INTRODUCTION Provision of timely and reliable information on traffic incidents and subsequently induced congestion is a critical ability to the successful deployment of many envisioned Advanced Traffic Management and Information Systems (ATMIS). Ideally, anticipated and quantitative information such as time-dependent delay caused by an incident should be estimated and provided to drivers to maximize the effect of information provision (). Procurement of such information is however not a trivial task because of the complex interactions among various factors such as incident location and severity, incident response capacity, demand fluctuation and diverse driver responses to information. Moreover, most of these factors are subject to high uncertainty and information available to quantify them is often incomplete and subjective in nature. Consequently, provision of crisp values of expected delays to drivers through systems such as changeable message signs (CMS) would inevitably lower drivers trust in the accuracy of the provided information because the actual delays they would experience will be either larger or smaller than what were suggested. This underlying dilemma has become a major reason for many traffic management authorities to opt for less effective, but more credible alternatives such as providing qualitative information only (2). The goal of this paper is to demonstrate that, with an appropriate delay prediction model, it is possible to resolve this dilemma. Incident delay that a vehicle may experience at an incident location can be obtained using a deterministic model such as deterministic queuing model or shock wave theory if the future values of the associated parameters such as traffic arrival rate, capacity reduction, and incident duration can be identified exactly (3, 4, 5). In real-time application environments, however, these parameters are often subject to large variations and can not be predicted exactly. As a result, a deterministic model is not appropriate for use to model incident conditions (6).

Estimation and prediction of link travel time, which would include incident delay if an incident has occurred on the associated link, have been one of the major research interests in the area of ATMIS. Various methods have been proposed for demonstration ATMIS projects and simulation studies (7, 8, 9). The fundamental idea behind most of the proposed methods is a data fusion process that combines real-time travel time estimates from loop detectors and probe vehicles with historical travel times. Although impact of incident occurrences on travel times is considered through allocating higher weight to real-time estimates, these methods do not fully make use of available information such as lane closure, possible incident duration and future traffic volume. This paper proposed a framework based on fuzzy set theory to model the evolution of incident congestion or queue development. The rationale behind this approach is that information typically available under incident conditions is often in the form of linguistic descriptions characterized by imprecision and vagueness. The paper first describes various uncertainties involved during incident conditions and how they can be systematically modeled on the basis of fuzzy set theory. A fuzzy queuing model is subsequently presented for predicting the possible delay or interval of delay that a vehicle will experience at an incident location based on real-time information on current queuing condition, future traffic arrival, capacity reduction and the vehicle s arrival time. Lastly, a case study is presented to demonstrate the application of the proposed model in composing location-dependent delay messages for CMS. A FUZZY INCIDENT DELAY MODEL Delay that a vehicle will experience as a result of an incident depends on many factors including incident severity (capacity reduction), incident duration, traffic volume and the time when the vehicle arrives at the incident location. In a practical situation, each of these factors is subject to uncertainty. It would be a matter of a simple application of deterministic queuing theory or shock wave theory if we could predict the exact value of these factors. However, in a practical situation, each of these factors may be subject to a certain level of uncertainty and the information that may be available for estimating these factors is commonly imprecise or vague, as discussed below (0). Traffic Arrival Rate (V): Under normal traffic condition, the traffic arrival rate is usually stable and can be fairly accurately estimated based on historical traffic counts. However, during incident conditions, the prediction of traffic arrivals are no longer a trial task because some drivers may have be informed of the incident occurrence and decide to divert to other routes. How much traffic will divert and at what rate will depend on many factors such as traffic information coverage, drivers acceptance of the provided information and local network conditions. Currently, there is no dependable model available that can be used to model such complexity and to provide an accurate prediction of the dynamic traffic conditions. However, it can be reasonably expected that an approximate estimation of the traffic diversion may be available. For example, an experienced traffic manager may be able to give such estimation as About 20%~30% of the traffic will divert to avoid the incident congestion, Majority of the drivers will 2

still use the route even one lane is closed due to the incident. These linguistic terms can be properly modeled by a fuzzy set or a fuzzy number. Incident Duration (L): The time taken to remove an incident and recover the road capacity, or incident duration, is another key piece of information needed for predicting the incident delay. It has been observed that incident duration usually has a large variation depending on incident severity and location, traffic condition and the availability of incident management (). For example, Giuliano () showed that the mean duration is about 37 minutes with a standard deviation of 30 minutes. Therefore, it is nearly impossible to give a precise prediction on the incident duration even when there is a large amount of historical data available. However, it is not unusual for an experienced incident response team or highway police to give an estimation of the duration after they know the incident situation and location. For examples, they can present it as it would take about 30 to 40 minutes to remove the stalled car, it takes at least one hour or It shouldn t take longer than two hours. Such information presented in linguistic terms is commonly imprecise or vague and can be adequately represented using fuzzy numbers. Current Queue (Q): There is usually a time lag between the current time (the time to make a prediction) and the incident occurrence. Therefore, it is likely that a queue has formed at the incident location. The current queue can be estimated based on information such as the elapsed time from incident occurrence, traffic arrivals and reduced capacity. However, it is more likely that it can be directly obtained from various information sources such as observers, police or a special incident response team. This information is usually a linguistic description on the queuing status (e.g. the queue is about to backup to 2 th street ). It should be noted that in most cases, these description often gives the information on the current queue reach instead of the queue length. However, these two variables can be considered as the same before the incident is removed. For the same reason as for the previous parameters, it can be nicely represented using a fuzzy set. Incident Capacity (C): The capacity under normal traffic condition can be considered as constant for a given road section and estimated based on HCM (2). During incident condition, it has been observed that the departure rate from the queue (or capacity during incident condition) varies significantly because of the stop-and-go process and gawkers block. The actual value can be as low as 500 to as high as 2000 pcu/hour/lane, depending on the local condition and driver behavior (2). It has also been observed that the lane-blocking incidents have more than a proportional impact on capacity. For example, Urbanek and Rogers (3) indicated that the blockage of a signalized lane on a three-lane facility reduced freeway capacity by 40~50% (instead of 33% based on space reduction). Accordingly, it is desirable to use a fuzzy number to model the reduced capacity during incident. Vehicle Arrival Time (T): The incident delay that a vehicle may experience also depends on when the vehicle will arrive at the incident location. For example, if a Traffic Information Center (TIC) is to provide information to a vehicle currently at a known location, the prediction of incident delay also requires the estimation of travel time from the current location to the incident location. There is no doubt that this travel time involves uncertainty caused by many factors such as variation of traffic demands, traffic control and driving conditions. Although this travel time 3

can be modeled as a random variable, the underlying distribution may not follow a popular mathematical distribution such as normal, log-normal or Beta distribution. Therefore, it is also appropriate to use a fuzzy number to represent the travel time or arrival time at the incident location. A methodology that applies a traditional deterministic queuing model with fuzzy input parameters has recently been proposed (0). Given the input parameters described above, the functional relationship between the fuzzy incident delay and the input variables was established based on arithmetic operations of fuzzy numbers using the α-cut concept (4). Figure shows the α-cut representation of a fuzzy queuing model where the current queue (Q) is represented by its α-cut interval [q, q 2 ], the possible cumulative traffic arrival at the incident location is represented by two straight lines with rates of v and v 2 - the α-cut interval of the fuzzy traffic arrival rate V. The cumulative departure at the incident location is represented by two straight lines with rates of c and c 2, which intersect with lines representing the recovered full capacity (s) at the incident removal time, L, represented by l and l 2. The α-cut interval of the incident delay, D α, for a vehicle arriving at the incident location at a given time (T) - represented by t and t 2, can then be obtained by Equation : Cumulative Arrivals and Departures One possible incident delay, d Arrivals v 2 s v Current queue q 2 q c 2 c Departures t t 2 Arrival time l l 2 Incident removal time Time of Day Figure. An α-cut representation of a fuzzy incident delay queuing model ( α ) ( α ) Dα = [ d, d 2 ] () 4

where d and d 2 are known functions of v,v 2,q,q 2,c,c 2,c 2,s, l, l 2, t and t 2 (refer to (0) for the complete equations). The developed fuzzy incident delay models allows prediction of uncertainty associated with incident delay and can be readily used by ATMIS applications that need to consider uncertainty in information provision. Real-time information, such as queue length reported by traffic patrol and/or eyewitnesses, updated short-term prediction on traffic volume and estimate of incident duration provided by on-site incident response team can be translated into fuzzy numbers and fused as input to corresponding parameters. Prediction on delay can therefore be updated and improved adaptively in real time. CASE STUDY This section presents a hypothetical case to demonstrate the use of the proposed model. Consider the case that an accident was detected on a three-lane freeway section at 3:20pm, as shown in Figure 2. Two CMS were available for the traffic management center (TMC) to post incident delay information. The messages were intended for drivers who are just in the view of one of the CMS. In order to determine what message should be displayed, the TMC needed to predict the possible delays that would be experienced by vehicles if they were to continue to travel on the freeway section instead of diverting to other routes. It was assumed that some incomplete pieces of information were available which allowed representing the traffic arrival rate, incident duration, and capacity during incident and current queuing status as fuzzy trapezoidal numbers (TrFN) represented by {a, m, n, b}. The data used for analysis are summarized in Figure 2. Two prediction scenarios were considered. Scenario represented the estimation task at the time 3:20pm, that is, right at the time the incident was detected while Scenario 2 was set to model the prediction task at 3:40pm, at which the incident was expected to be removed soon. Note from the data that a larger value of current queue was used in Scenario 2 as compared to Scenario to reflect the likely congestion development. 0km - about 6 minutes of CMS 2 0km - about 6 minutes of CMS Incident detected at 3:20pm Figure 2. Freeway section with an incident Based on the given data, the fuzzy incident delay, as approximated using five levels of presumption with α ={0; 0.2; 0.4; 0.6; 0.8;.0}, can be calculated for further analysis. Figure 3 5

shows the predicted fuzzy incident delays under two given scenarios for vehicles at each CMS. The arrival times of the vehicles at a given CMS were generated based on the distance from the CMS to the incident location and a speed of 00km/h. The predicted delay values for both scenarios at each CMS are shown in Figure 3. The following findings are observed: There exists a significant amount of uncertainty in incident delay. This indicates the need to recognize it explicitly in information provision. For example, instead of displaying a single estimate of delay on a CMS, an interval of possible delay, such as incident delay between 5- CMS at different locations (distances from an incident spot) should display different delay information to account for differences in vehicles expected arrival time. Generally, during the time period that incident congestion starts to build up (Scenario ), CMS farther away from the incident spot (CMS 2) should display delay values higher than those on CMS nearer to the incident spot (CMS ). Conversely, during the time period that incident is soon to be removed (Scenario 2), CMS near the incident spot (CMS ) should display delay values higher than those on CMS farther away from the incident spot (CMS 2). CONCLUDING REMARKS This paper has presented a fuzzy queuing model that can be used to predict the delay that a vehicle would experience if traveling through an incident location. In contrast to the traditional deterministic models, the new model explicitly considers the uncertainties involved in future traffic arrivals, incident duration, departure rate during incident, current queue status and vehicle arrival time, allowing easy incorporation of imprecise and vague information on these variables. The proposed model does not require significant additional data or computational requirements over traditional methods and therefore may readily be adopted for ATMIS applications or simulation studies. The numerical example has demonstrated the necessity for providing drivers incident delay information that is location-varied and explicit in informing uncertainties. It should be pointed out that the methodology presented in this paper assumes that the input variables can be modeled as fuzzy numbers and the related membership functions are known a priori. As a result, the implementation of the proposed model requires an interface to generate membership functions of the input variables based on various sources of information in real-time. The next step of this research will focus on the development of this type of interface and the calibration and refinement of the proposed model for application in ATMIS. ACKNOWLEDGEMENT This research is funded by the Natural Science and Engineering Research Council of Canada. 6

Scenario Scenario 2 Prediction time: Current Queue (pcu): Incident duration (minute): Traffic volume (pcu/h): Capacity with incident (pcu/h): Full capacity (pcu/h): CMS 2: Arrival time (minutes) 3:20pm {4,8,22,26} {25,29,3,35} {3590,360,4390,440} {430,450,750,770} 5400 {4.9,7.,5.9,6.} 3:40pm {52,56,64,68} {4.5,5,5,5.5} {3590,360,4390,440} {430,450,750,770} 5400 {4.9,7.,5.9,6.} Predicted Delay Degree of membership 0.8 0.6 0.4 0.2 0 0 5 0 5 20 25 30 Incident delay (minutes) Degree of membership 0.8 0.6 0.4 0.2 0 0 2 4 6 8 0 Incident delay (minutes) Suggested Message CMS 2: Arrival time (minutes) Incident Delay 5-5 min. {0.3,.3,2.7,3.7} Incident Delay -4 min. {0.3,.3,2.7,3.7} Predicted Delay Degree of membership 0.8 0.6 0.4 0.2 0 Degree of membership 0.8 0.6 0.4 0.2 0 0 5 0 5 20 25 30 Incident delay (minutes) 0 2 4 6 8 0 Incident delay (minutes) Suggested Message Incident Delay 0-25 min. Incident Delay under 3 min. Figure 3. Estimation of fuzzy incident delay for CMS 7

REFERENCES [] Wardman, M, P. W. Bonsall and J. D. Shires (998) Driver Response to Variable Message Signs: A Stated Preference Investigation,, Vol. 5, No. 6, pp. 389-405. [2] Masters, P., G. Paesani and A. Coady (995) Congestion Management by Dynamic Guidance on Toronto s Dual Roadway Freeway, Proceedings of the Advanced Traffic Management Conference, St. Peterburg, Florida, October 3-8, 993, pp. 27-299. [3] Messer, J. C., C. L. Dudek and J. D. Friebele, (973) Method for Predicting Travel Time and Other Operational Measures in Real-time During Freeway Incident Conditions Highway Research Record 46, TRB, National Research Council, Washington, D. C., pp. ~6, 973. [4] Wirrasinghe, C. S. (978) Determination of Traffic Delays from Shock Wave Analysis, Transportation Research, Vol. 2, pp. 343~348, 978. [5] Morales, M. J. (986) Analytical Procedures for Estimating Freeway Traffic Congestion Public Road, Vol. 50, No. 2, September, pp. 55~6, 986. [6] Fu, L. and L. R. Rilett (997) Real-time Estimation of Delay due to Incidents in Dynamic Transportation Research Record 603, Transportation Research Board, National Research Council, Washington DC, 998. [7] Hoffman, C. and J. Janko (990) "Travel Time as a Basis of the LISB Guidance Strategy" Paper presented at the IEEE Road Traffic Control Conference, London, 990. [8] Koutsopoulos, H. N. and A. Yablonski (99) Design Parameters of Advanced Information System: the Case of Incident Congestion and Small Market Penetration, Proceedings of the Second Vehicle Navigation and Information System Conference, Dearborn, MI, 99. [9] Boyce, D., Rouphail, N., and A. Kirson (993) "Estimation and Measurement of Link Travel Times in the ADVANCE Project", IEEE-IEE Vehicle Navigation & Information Systems Conference, Ottawa, Ontario, 993. [0] Fu, L (998) A Fuzzy Queuing Model for Prediction of Incident Delay under Imperfect Information, submitted for publication in Transportation Research C. [] Giuliano, G. (989) Incident Characteristics, Frequency, and Duration on a High Volume Transportation Research -A, 23A, No. 5, pp. 387~396, 989. [2] TRB, Highway Capacity Manual, 994 Update, Special Report 209, Transportation Research Board, National Research Council, Washington, D.C., 994. [3] Urbanek, G. L. and R. W. Rogers (977) Alternative Surveillance Concepts and Methods for Freeway Incident Management, Volume : Executive Summary, Federal Highway Administration. Report # FHWA-RD-77-58, Washington, D.C. [4] Kaufmann, A. and M. M. Gupta, Introduction to Fuzzy Arithmetic: Theory and Applications, Van Nostrand Reinhold Company, New York, 985. 8