BaSSy Frequencies
BaSSy Frequencies
Page 2 of 59
TABLE OF CONTENTS
1. 2. 3. 3.1 4. 4.1 4.1.1 4.1.2 4.1.3 4.2 4.2.1 4.2.2 4.3 4.4 4.5 5. 5.1 5.2 5.2.1 5.2.2 5.2.3 5.3 6. 6.1 6.1.1 6.1.2 6.2 6.2.1 6.2.2 6.2.3 6.2.4 7. 8. 9. BACKGROUND .......................................................................................................................... 3 INTRODUCTION .......................................................................................................................... 3 PROBABILISTIC COLLISION AND GROUNDING ANALYSIS............................................................... 5 Risk models ............................................................................................................................ 5 PREDICTING COLLISION AND GROUNDING FREQUENCIES............................................................. 7 Frequency of collision............................................................................................................. 7 Head-on and overtaking collisions ......................................................................................... 8 Crossing collisions................................................................................................................ 10 Collision test cases............................................................................................................... 16 Probability of grounding........................................................................................................ 22 Powered grounding .............................................................................................................. 24 Drifting grounding ................................................................................................................. 26 Assessing the traffic spread across the route...................................................................... 28 Calculation procedure for estimating the collision frequency.............................................. 28 Combined causation factor PC i , j ........................................................................................ 29 CAUSATION PROBABILITY ........................................................................................................ 29 Causation probabilities from literature ................................................................................. 30 Risk model for obtaining the causation probability .............................................................. 32 Traditional approach............................................................................................................. 33 Using Bayesian Networks .................................................................................................... 33 Bayesian Network for ship-ship collisions............................................................................ 36 Default values used in GRISK.............................................................................................. 39 FACTORS THAT INFLUENCE THE CAUSATION PROBABILITY ........................................................ 39 Reported causes for grounding and collision ...................................................................... 39 Human and Organisational Errors ....................................................................................... 40 Human error evaluation........................................................................................................ 41 Aspects that the risk analysis should include ...................................................................... 42 Configuration of navigational area ....................................................................................... 43 Composition of ship traffic .................................................................................................... 45 Environmental conditions ..................................................................................................... 46 Configuration of considered vessel...................................................................................... 46 SHIP TYPES USED IN THE GRISK PROGRAM ............................................................................ 48 CONCLUSION .......................................................................................................................... 48 REFERENCES ..........................................................................................................................49
Date: 12/05/2009
Page 3 of 59
1. BACKGROUND
The objective of the present report is to describe the theoretical background for the collision and grounding frequency analysis that forms the basis of the developed program termed the BaSSy ToolBox or GRISK (GateHouse Risk). The BaSSy ToolBox (GRISK) was developed under the BaSSy-project from 2007 to 2009. The BaSSy project is a joint research project between Technical University of Denmark, GateHouse (Denmark), SSPA (Sweden), and VTT (Finland), which is funded in part by The Danish Maritime Foundation and Det Nordiske Ministerrd. The objective of the present report is to describe the theoretical foundation for the collision and grounding frequency analysis so that the interested user of BaSSy ToolBox may understand the fundament behind the program. It is assumed that the reader assumes a level of mathematical and probabilistic skills to gain full benefit of the report. At an early stage of the program development IALA (Internation Association of Marine Aids to Navigation and Lighthouse Authorities) found interest in the product to be developed. Therefore, a follow group of IALA members were established to assure that the final product could meet the needs of a potential large user group. The fruitful and successful collaboration with the IALA follow group resulted in a name change of the BaSSy ToolBox to Grisk, IWRAP MK II. IALA supports this version of the program and recommend its members to use this for assessing collision and grounding frequencies in navigational geographical areas.
Date: 12/05/2009
Page 4 of 59
The procedures developed during the development of GRACAT constitutes an essential part of the BaSSy project and thus also of the program with the GRISK, IWRAP MK II. Therefore, the theoretical foundation given in this document is to large extent routed in the basis established during the ISESO-project. The present document not only defines the theoretical background for the collision and grounding analysis, but it also summarises and discusses the background for the so-called causation probability. The document outlines a method for evaluating the collision and grounding frequency of vessels operating on a specified route. To identify the frequency of experiencing any collision or grounding in a given area involves first a specification of the routes and the associated traffic on the routes. Subsequently, the collision and ground frequency may be obtained by looping over all vessels operating on the route. The BaSSy program will contain tools for extracting the traffic distribution and traffic density functions from AIS data. Given that a collision or grounding has taken place the spatial distribution of the damages may further be calculated. Results of such analysis may in the BaSSy program be presented in terms of probability distributions, for indentation depth, length and height of the holes and for their location. Knowing the structural damage the resulting consequences in terms of bunker oil outflow and cargo outflow may subsequently be calculated. In future more consequence models may be implemented in the BaSSy program. One of the benefits of the formulated procedure is that it allows comparisons of various navigational routes by assessing the relative frequencies of collisions. Under the ISESO project the derived procedure was applied to different Ro-Ro passenger vessel routes (Great Belt, Dover-Calais, Turku-Stockholm). The results of the analyses were compared to registered data and good agreement was found in all cases. This constituted the validation of the software for frequency and damage distribution estimation by the GRACAT program. The GRACAT program has been completely re-programmed under the BaSSy project and several new facilities and tools have been implemented. Consult the User Manual for details on these. An independent validation of GRISK, IWRAP MK II program has been conducted for Baltic West, Sea of Aaland, and waters near Malaysia. The Baltic West validation constituted an independent assessment of the collision and grounding frequencies by COWI. The results of the analyses in Sea of Aaland, the Malaysian waters, as well as the Baltic West were compared to observed historically data. In all cases a good agreement was obtained. Further validation test cases may be found reported at http://www.iala-aism.org/. The applied model for calculating the frequency of grounding or collision accident involves the use of a so-called causation probability that is multiplied onto a theoretically obtained number of grounding or collision candidates. The causation factor models the probability of the officer on the watch not reacting in time given that he is on collision course with another vessel (or alternatively on grounding course). The numerical value of the causation probability is not a unified value but may often vary for different geographical locations. The applied value of the causation probability is therefore typically adjusted by a calibration to registered data. On the basis of a literature search this document summarises some of the causation probabilities that have been applied in different studies. The document also identifies some of the factors that are of importance when assessing the causation probability. Moreover, a Bayesian Network model for ship-ship collision is formulated for an analytical estimation of the causation probability. The obtained result agrees well with that obtained from statistical analyses of data.
Date: 12/05/2009
Page 5 of 59
Date: 12/05/2009
Page 6 of 59
An alternative risk analytical approach, the synthesis approach, see Gluver and Olsen [12], base the risk of grounding or collision on a set of scenarios where specific error situations or conditions are assumed to occur or exist in the vessel prior to or during the considered critical situation. Such an approach, however, requires that all significant accident scenarios are identified and analysed. Consequently, this also implies that the causation probability must be defined conditional on considered accident scenario. It therefore follows that the advantage of introducing the synthesis approach is that alternate risk-mitigating aspects more easily may be both identified and quantified. Examples of different accident scenarios could be rudder stuck, power failure, navigational error, etc. Each of these scenarios may be further sub-divided to describe the scenario in more details e.g., in what position the rudder is stuck and what other equipment is available to mitigate the problem.
Figure 1 Overall procedure for probabilistic prediction and spatial distribution of collision damages.
Date: 12/05/2009
Page 7 of 59
In the present work the scenario approach is applied and the procedure is schematically illustrated in Figure 1 for the collision analysis. A grounding analysis follows the same conceptual outline. Basically the procedure is as follows: First the relevant navigational area is described. This involves description of the traffic composition along all navigational routes and descriptions of all grounds in the vicinity of the route. Next the considered vessel (termed struck vessel in Figure 1) is defined to be operating on a specified route in the defined navigational area. All potential other vessels (striking vessel in Figure 1) or grounds is then identified and the probability of grounding and collision is calculated. Subsequently, the identified ground or striking vessel may further be used for calculating damage statistics. The ensuing consequence analysis (in terms of time to capsize, oil outflow, etc.) of the identified damages is not shown in the figure, but statistics for this may similarly be obtained. Although the procedure described above resembles the scenario approach the alternative synthesis approach may also be covered by careful application of the causation factors. Structured methods for this will be illustrated later by the application of Bayesian Networks for obtaining the causation factor.
The theoretical procedure laid out in this chapter represent the state-of-the-art framework that is applied for calculating the geometric number of collision and grounding candidates, N G . The values of the causation factor, PC , are typically in the range from . A prerequisite for the analysis is that the ship traffic has been grouped into a number of different ship classes according to vessel type, size, loaded or ballasted, with or without bulbous bow etc., and that the number of vessels per time unit have been registered for each waterway. It is noted that in the analysis presented below the time unit for the definition of number of vessels is in seconds-1 for dimension correctness.
Date: 12/05/2009
Page 8 of 59
The procedure for calculation of the number of collision candidates, N G , for the abovementioned two types are conceptually different since the geometric number of collision candidates first type becomes dependent on the lateral traffic spread on the route whereas the second is independent of the traffic spread. This can be seen by comparing Figure 2 and Figure 3. By inspecting Figure 2 it can be seen that the probability of the path of two meeting ships will overlap depends on the spreading of the lateral position where the vessels are sailing. The larger the -value the smaller becomes the probability of a collision. In Figure 3 it can be seen that although the risk area is affected by the spread of the traffic the probability of the ships meeting each other is not. In the following the head-on and overtaking collisions will first be treated, thereafter will the crossing collisions.
Figure 3 Crossing waterways with risk area of ship-ship collision. 4.1.1 Head-on and overtaking collisions Collisions along the route, see Figure 2, depends of
Date: 12/05/2009
Page 9 of 59
The length, LW , of the segment; The traffic composition, i.e. the number of passages per time unit for each ship type and size, Qi(1) and Q (j 2 ) , in each direction, (1) and (2), and their speed, Vi (1) and V j( 2 ) ; The geometrical probability distribution, f i (1) ( y ) and f j( 2 ) ( y ) , of the lateral traffic spread on the route. The traffic spread is typically defined by a Normal distribution but may in principle be of any type. The sign convention for the traffic distribution is measured from the centre of the channel and positive towards the right side in the sailing direction.
For head-on collisions the number of geometric collision candidates for ships sailing along the route segment in direction (1) and (2) can be expressed as,
head N G - on = LW PGhead - on i, j i, j
Vij Vi V
(1) ( 2) j
(Q
(1) i
Q (j 2 )
(4.2)
where Vij = Vi (1) + V j( 2 ) is the relative speed between the vessels and PG defines the probability that two ships will collide in a head on meeting situation. This probability is expressed as
B ( 2) B ( 2) B (1) B (1) PGhead - on = P yi(1) i < y (j2) + j I yi(1) + i > y (j2 ) j i, j 2 2 2 2 (1) (1) Bi(1) + B (j 2 ) Bi(1) + B (j 2) ( 2) ( 2) = P yi + y j < P yi + y j < 2 2
= =
Yi y j = yi B
yi + B
(4.3)
( yi ) fY j ( y j ) dyi dy j
Yi
( yi ) FY j ( yi + B) FY j ( yi B) dy j
It is noted that the random variable y (j2 ) is negative because of the positive sign convention in the sailing direction of the two vessels. In the last step it has been utilized that the two distributions are independent. It is possible to establish a closed form solution to Eq. (4.3) when the traffic distributions are normally distributed. In the general case Eq. (4.3) must in be solved by approximate procedures such as FORM/SORM or numerical integration. When f i (1) ( y ) and f j( 2 ) ( y ) both follow a normal distribution with
B ij B + ij ij PGhead - on = ij i, j ij ij
(4.4)
Date: 12/05/2009
Page 10 of 59
In which ( x ) is the standard normal distribution function, ij = i(1) + (j 2 ) is the mean sailing distance between the two vessels, ij = deviation of the joint distribution, and Bij =
( ) + ( )
(1) 2 i
( 2) 2 j
is the standard
Bi(1) + B (j 2) 2
head The frequency of head on collisions, Col - on , is obtained by multiplying the geometric
head number of collisions, N G - on , with the causation factor for head on collisions, PChead - on . In
the DROGDEN study a causation factor of 1.3 104 was applied for head-on and overtaking collisions. In the resent study Oil and Chemical spills in Danish waters [3] a factor of 3.0 104 was proposed. Based on collision statistics in Japanese waters, Fujii et al. [8] has estimated that for meeting ships in parallel waterways Pc = 4.910-5. For overtaking collisions the number of geometric collision candidates for ships sailing along the route segment in direction (1) is expressed by eq. (4.2) using the relative speed Vij = Vi (1) V j(1) , Vij > 0 . If Vij < 0 then vessel i will obviously not be able to overtake vessel j . In the practical implementation the absolute value of Vij is used and struck and striking vessel are registered. The geometric probability of meeting, eq. (4.3) becomes,
overtaking G i, j
(4.5)
For normally distributed variables the mean value in eq. (4.4) should be replaced by ij = i(1) (j1) to handle the overtaking situation.
4.1.2 Crossing collisions The frequency of crossing collisions depends on the angle between the two lanes. Figure 3 shows two crossing waterways for which the ship traffic also is given. The geometric number of crossing collisions candidates for crossing waterways can similarly to eq. (4.2) be expressed as,
crossing G
V
i, j
Qi(1)Q (j 2)
(1) i
( 2) j
D ij Vij
1 sin
for
(4.6)
where V ij =
(V ) + (V ) - 2 V
(1) 2 i ( 2) 2 j
(1) i
and Di , j defines the apparent collision diameter, see Figure 4. The sinus term stems from the variable transformation when integrating over the area of the joint probability distribution, see Figure 7. Note that contrary to head-on and overtaking collisions the distribution of the traffic spread is not relevant for crossing collisions, except for the sinus term of course. It is seen that when the crossing angle goes to zero the length of the crossing (or the time of the crossing) goes to infinity and hence does the number of collisions. For practical reasons it is therefore necessary to limit the crossing angle to an interval of, say, 10 to 170 degrees.
Date: 12/05/2009
Page 11 of 59
Figure 5 Calculation of the geometrical collision diameter Dij . As mentioned Dij is the geometrical collision diameter illustrated in Figure 4. If it is assumed that the ships can be approximated by rectangular shapes, then it can be shown, see Figure 5, that:
Date: 12/05/2009
Page 12 of 59
D ij =
(1) i
1/ 2
+ B i( 2) 1
where the relative velocity Vij is determined as
V - sin V ij
( 2) j
1/ 2
(4.7)
V ij =
(V ) + (V ) - 2 V
(1) 2 i ( 2) 2 j
(1) i
( 2) V j cos
(4.8)
Figure 6 Illustration of the apparent diameter for vessel 1 striking vessel 2. In the present work we are not only interested in the number of ship-ship collisions but also in the probability of the one or the other being the struck or striking vessel. To derive a simple expression for this event the apparent collision diameter is formulated, see Figure 6. The apparent collision diameter seen from vessel i , Di(1) , can be determined as
(1) L V sin 1 (2) (1) sin V i + Bj 1 - Di = 2 V ij V ij ( 2) j (1) i
1/ 2
(4.9)
( 2) 1 V + B i(1) 1 - sin j 2 V ij
1/ 2
Date: 12/05/2009
Page 13 of 59
Similarly, for the case where vessel j in waterway 2 is striking vessel i in waterway 1 the apparent collision diameter is:
( 2) Dj =
L V
(1) i
( 2) j
sin
V ij
( 2) 1 V + B i(1) 1 - sin j 2 V ij
1/ 2
(3.10)
1/ 2
It is seen that the total collision diameter Dij is the sum of the two apparent collision diameters, i.e.:
(1) (2) Dij = Di + D j
The probability of vessel i in waterway 1 striking vessel j in waterway 2 given a collision may then be determined as
(4.12)
Dj D ij
( 2)
(4.13)
The frequency, ship ship , of ship-ship collision per time unit is then determined as
ship.ship = P C i , j N G i , j
(4.14)
Due to the fact that both of the involved two ships have the possibility of making aversive manoeuvres, the causation probability, PC , for ship-ship collision is smaller than the one given for grounding and collision against fixed objects, see Section 5.1. Based on collision statistics in Japanese waters, Fujii et al. [8] has estimated that for crossing ships Pc = 1.210-4 and for meeting ships in parallel waterways Pc = 4.910-5. Given the frequency of the (annual) number of collisions ship ship the probability of having a collision during time interval t can be estimated on the assumption of arrivals of the collisions as points in a Poisson process:
ship - ship 0
Date: 12/05/2009
Page 14 of 59
Figure 7 Basic layout of the simulated crossings with results presented in Figure 8 and Figure 9.
Figure 8 Comparison between simulated results assuming a Poisson distribution of ships in the two waterways and analytical results using Eq. (3.6)
Date: 12/05/2009
Page 15 of 59
Figure 9 Simulated and analytical results for the probability of vessel 2 colliding with vessel 1 given a collision. In order to verify the established analytical model and at the same time gain insight in the sensitivity of the number of collision candidates with respect to parameters such as crossing angles , ship dimensions, and ship speeds a program has been written which is based on time simulation. What has been considered is vessels of the same type in waterway 1 crossing the axis , see Figure 7, as points in a Poisson process with intensity Q1 = 20000 vessels per year. The vessels in waterway 1 are assumed normally distributed over the width of the channel with = 100 m and = 45 m. Similar assumptions are made for vessels in waterway 2 crossing the axis in Figure 7. Here Q2 = 50000 vessels per year and = 100 m and = 45 m. The vessels in the two waterways are moving with speed V1 and V2, respectively. A time history of +/- one hour of simulated vessels in waterway 1 is kept for matching against simulated vessels in waterway 2, i.e. t1 - 3600 s < t2 < 3600 s + t1. During the simulation it is calculated whether or not the two vessels are colliding. If they are observed to collide, it is further identified which of the vessels is the struck vessel and which is the striking vessel. Figure 8 shows the number of collision candidates during a period of forty years determined by time simulation and determined by the analytical expression Eq. (3.6) as functions of the angle between the two waterways.
Date: 12/05/2009
Page 16 of 59
4.1.3 Collisions with area traffic By area traffic is understood that traffic that do not follow ordinary routes. The area traffic is composed of leisure crafts and fishing vessels. These vessels will cross the routes at which the line traffic operates at random angles. The number of collisions between the area traffic and the line traffic is calculated by assuming that the area traffic crosses the route the line traffic operates on at eight different directions. The number of collisions from each direction is calculated by use of the collision algorithm describe in the previous section. 4.1.4 Collision test cases This subsection describes the result of a series of selected test cases that were analysed by use of the GRISK program and by hand calculation. Only the number of collision candidates is calculated. 4.1.4.1 Test 1, Head on collision Test case 1 calculates the number of head on collisions per year. The scenario is: Length of leg Ships in each direction Length of ships Breadth of ships Speed of ships Traffic distribution Mean position from leg Standard deviation Causation factors 35,046 m 10,000 214 m 33.4 m 14.7 knots Normal dist 300 m 150 m 1.0
To calculate this scenario in the BaSSy toolbox, GRISK, do the following steps:
1. Define leg
Date: 12/05/2009
Page 17 of 59
5. Run the job and read results Note that we here use 10000 crude oil tankers in the length interval 200-225 m. From this length interval and ship type GRISK looks up the breadth and speed using the predefined dimension tables described in the appendix Hand calculation of the head-on collision scenario
head N G - on = LW PG i , j (head - on ) i, j
Vij Vi V
(1) ( 2) j
(Q
(1) i
Q (j 2)
33.4m 600m 33.4 + 600m PGhead - on = = (2.67) (2.99) = 0.00237 i, j 212.1m 212.1m head N G - on = 35,046m 0.00237 2.646 108 (360 24 3600) = 70.5
i, j
Date: 12/05/2009
Page 18 of 59
4.1.4.2 Test 2: Overtaking collision Calculates the number of collisions per year on a leg where ships sail in the same direction but at different speeds. The scenario is: Length of leg Number of ship 1 Length of ship 1 Breadth of ship 1 Speed of ship 1 Number of ship 2 Length of ship 2 Breadth of ship 2 Speed of ship 2 Mean position from leg Standard deviation Causation factors 35,046 m 10,000 214 m 33.2 m 14.7 knots 10,000 162 m 25.0 m 18.9 knots 300 m 150 m 1.0
1. Define leg
Date: 12/05/2009
Page 19 of 59
Hand calculation of the head-on collision scenario This only difference to the head on collision calculation is the sign of the two speeds.
overtaking = LW PG i , j (head - on ) NG i, j
Vij Vi V
(1) ( 2) j
(Q
(1) i
Q (j 2 )
B + ij B ij overtaking ij PG i , j = ij ij ij
Date: 12/05/2009
Page 20 of 59
4.1.4.3 Test 3: Crossing collision Test case 3 calculates the number of crossing collisions per year. The scenario is: Ships North going Ships East going Length of ships Breadth of ships Speed of ships Angle between legs Causation factors 10 000 10 000 200 m 33.4 m 14.7 knots 88.8 deg 1.0
4. Define how much traffic sails from one Run the program and read the results leg to the others
Date: 12/05/2009
Page 21 of 59
The geometric number of crossing collisions candidates for crossing waterways can similarly to Eq. (3.2) be expressed as,
crossing NG =
V
i, j
Qi(1)Q (j 2)
(1) i
( 2) j
D ij Vij
1 sin
for
(3.6)
where V ij =
(V ) + (V ) - 2 V
(1) 2 i ( 2) 2 j
(1) i
V ij =
2 200 7.56 + 200 7.56 7.56 = + 33.2 1 - sin(88.8) sin(88.8) Dij 10.58 10.58 2 7.56 + 33.2 1 - sin(88.8) 10.58
1/ 2
= 332.5
crossing NG =
10 000 10 000
= 197.8
Date: 12/05/2009
Page 22 of 59
4.1.4.4 Test 4: Merging collision Collision due to merging traffic is calculated as crossing collisions
Merging traffic
Date: 12/05/2009
Page 23 of 59
Figure 10 Observed grounding events over a 15-year-period in a Danish Strait, from [25]. In the following, expressions are presented for predicting the expected annual number of grounding events of category I and II accidents. The probability of category III and IV grounding events are today normally found by modification of the traffic distribution along the route. In the present work drifting ships (category IV) are modelled by assuming a drifting direction according to a user specified wind rose. Evasive manoeuvres (category III are not explicitly dealt with in the present version. We revert to these issues later.
Date: 12/05/2009
Page 24 of 59
Ships in category I and II, following an ordinary route, are distributed over a transverse section of the waterway with some probability density function, f i (z ) , where index i refers to a ship class and z is the transverse coordinate, see Figure 11. The shape of fi is a strong function of the considered waterway so a major challenge of the present approach is to define rationally f i (z ) along a given route. Given fi the number of candidates of grounding events can be calculated as an integral of fi over the width, zmin to zmax , of the obstacle. The hatched area in Figure 11 illustrates this. Most of these candidates will be aware of the danger and take the necessary aversive actions before they hit the obstacle. However, a fraction, Pc , of the candidates will fail to avoid the obstacle, due to for example human and technical errors. The fraction Pc is normally referred to as the causation probability, and it will be shown later how it can be estimated. Groundings that are caused by a meeting situations where ships may feel forced to give way, which then subsequent may result in grounding, has not been considered. Such a model requires much more advanced modelling than what is implemented at the present stage. Neither have groundings that are caused by a rudder stuck failure. In this case the rudder may either go the extreme starboard/port side or it will get stuck in a central position and cause the vessel either to turn in circles or to follow its path. The model requires more data information before it will be implemented.
Figure 11 Illustration of model for predicting the expected number of grounding events or collisions with fixed objects on a given ship route, from [25]. 4.2.1 Powered grounding According to the model described above, the expected number of grounding events in Category I and II can be calculated as
NI =
Ship class, i
z max
Pc,iQi
fi ( z )dz
z max
(4.15)
z min
N II =
Ship class, i
Pc ,iQi exp(d / ai )
z min
f ( z )dz
i
(4.16)
Date: 12/05/2009
Page 25 of 59
ai
d
Average distance between position checks by the navigator. Distance from the obstacle to the bend in the navigation route varying with the lateral position, s, of the ship. Index for ship class, categorised after vessel type and dead weight or length. Probability density function for the ship traffic Expected number of category I grounding events per year. Expected number of category II grounding events per year. Causation probability, i.e. ratio between ships grounding and ships on a grounding course. Number of ships in class i passing a cross section of the route per year. Coordinate in the direction perpendicular to the route. Transverse coordinates for an obstacle.
i
fi (z )
NI N II
Pc,i Qi
z
zmin , zmax
In the above it is assumed that the event of checking the position of the ship can be described as a Poisson process. Thus, the factor exp( d / ai ) represents the probability of the navigator not checking the position from the bend to the obstacle. The average distance between position checks is conveniently expressed in terms of the expected value of the time between position checks, , (approximately equal to 3 minutes) whereby the factor exp( d / ai ) becomes a function of the ship speed, ai = V . Eq. (4.16) is only correct for the case when the ground is orthogonal to the sailing route, which rarely is the case. In the event of the ground not being perpendicular to the sailing direction, but inclined so that the distance d from the bend to the ground may be expressed as d = az + b . In the event that the two traffic spread distributions follow a normal distribution then Eq. (4.16) can be simplified to
max 1 a 2 2 2bV zV + a 2 2 1 = Pc,iQi 2 exp 2 2 V V 2 Ship class, i z min z
N II
(4.17)
For other type of distributions, especially for mixed distributions, it is not a straight forward task to formulate a closed form solution to the number of grounding candidates. For such distributions the integral can be effectively solved by solving part of the integral analytical and then performing a numerical integration over the remaining variables. With the formulation above the expected number of annual grounding events becomes a function of traffic distributions, bottom topology, route layout etc. It is seen from Eq. (4.15), (4.16) and (4.17) that another important parameter is the causation factor, Pc , determining how large a fraction of the accident candidates actually run aground or hit the obstacle. Chapter 5 gives a thorough presentation of the causation factor. Often the causation probability is selected to be in the vicinity to Pc = 2 10 4 .
Date: 12/05/2009
Page 26 of 59
The calculated expected number of yearly (powered) grounding events, N g =N I + N II , can be considered as the intensity in a Poisson process. The probability of no grounding events in one year is then
P[Grounding ] = 1 exp( N g )
(4.18)
4.2.2 Drifting grounding The probability of category III (evasive manoeuvres) and category IV (drifting ships) grounding events are today typically found by modification of the traffic distribution along the route. Combining a 98% Gaussian distribution and a 2% uniform distribution performs the modification of the traffic distribution; see e.g. Gluver and Olsen [12] or Karlson et al. [19]. The value of 2% is based on engineering judgement and the results are dependent on the value especially in narrow restricted waters. Although this approach is very fast and easy to implement, it is considered to be too coarse a model that does not properly account for the physical effects that governs the drifting problem. In the section the implemented drifting model is defined. The two main causes for a ship to be not under command are rudder stuck and blackout of the main engine. Rudder stuck will not be dealt with in this study. Most ships experiences of the order of one black out of the main engine per year. The number of any blackout for a given ship will typically lie in the interval from 0.1 to 2 blackouts per year. The actual frequency of blackouts depends on the degree of redundancy and the general maintenance level of the ships. Ferries and ro/ro vessels generally have a high degree of built-in redundancy into the engine room (2 to 4 propulsion units) and hence they have a low frequency of blackouts. For many other single propulsion unit ships the frequency of blackouts are higher. In the present study the following blackout frequencies are selected as base values: Vessel type Passenger / Ro-Ro Other vessels Annual frequency 0.1 y-1 0.75 y-1 Hourly frequency 1.1510-5 h-1 8.5610-5 h-1
A blackout may be caused by contaminated fuel, internal fault in the main engine, or failure of the electrical system. The seriousness of the incident depends on the location at which the blackout occurs, the wind direction, wind speed, and of course the duration of the blackout (that is the drifting time). If a high degree of redundancy has been built into the engine room then the command over vessel may be regained in relative short time. In other situations, the drifting time may be of order of hours. The drifting ship will drift side ways and it will drift (approximately) in the direction of the wind. The drifting scenario may be remediated either by repairing the problem, by anchoring the vessel or by calling a tug boat. Failure of propulsion machinery may occur at any location along the waterway. Assuming that blackouts occur as points in a Poisson process then the probability of having a blackout along a leg segment of length Lsegment is:
Date: 12/05/2009
Page 27 of 59
In which blackout is the frequency of blackout and vvessel is the operational speed of the
drift vessel. The number of drifting groundings, N grounding , out of the N ship candidates of a
P
=0
360
wind
( )
Lsegment
( Lsegment )
x = 0 All v drift
In which Pwind ( ) defines the probability of having wind coming from direction . The probability of no repair is defined by the complementary distribution function of the repair time distribution. The default repair time distribution is modeled as a Weibull distribution,
N g =N I + N II + N III + N IV
Date: 12/05/2009
Page 28 of 59
Date: 12/05/2009
Page 29 of 59
Finally, when more than two legs meet at a waypoint the model calculates the probability for crossing collisions, as indicated in Figure 3. Dependent on how the vessel traffic on the legs meets each other, the scenario will be characterized either as either a crossing collision, a merging collision, or a bend collision. For each leg the identified number of collision candidates related to head-on, bends, and crossings is calculated for each vessel type and is subsequent multiplied by a causation factor. The following causation factors inspired by Fujii et al. [8]: P[head on] P[bend] P[crossing] = 4.910-5 = 1.310-4 = 1.310-4
These values for the causation factor are typical values for well regulated ship traffic in Japan. The causation factor will be a function of visibility, darkness, current and wind in the actual geographical area. All these factors suggest that larger values should be used around the Nordic countries. However, Fujii has also observed that passenger ferries have smaller collision probabilities than ordinary merchant vessels. This is due to the navigator awareness of the area and the fact that there are two navigator onboard passenger ferries. The following chapter discusses the assessment of the causation probability in more detail. The causation factors suggested used in the present study are presented in section 5.3.
PC i , j = PC i PC j
(3.15)
This is a justifiable pragmatic approach that assures a balance between failures of the navigational watch keeping on the two vessels.
5. CAUSATION PROBABILITY
This chapter presents a comprehensive collection of causation probabilities that have been proposed in literature. Further, a risk model that may be used for evaluating the causation probability is presented. Inadequacies of a frequently cited risk model are discussed and instead we propose to apply Bayesian Networks. A Bayesian Network for obtaining the causation factor for ship-ship collision is established, and the results are compared to available statistics, where good agreement was found.
Date: 12/05/2009
Page 30 of 59
Vessel grounding
Location Pc [10 ]
Japanese Straits Japanese Straits Japanese Straits Dover Strait Dover Strait Strait of Gibraltar resund, Denmark [1.0; 6.3] 1.58 [0.8; 4.3] 1.55 1.41 2.2 2.0 No traffic separation With traffic separation Collisions and grounding
-4
Comment
Ship-ship collisions
Location Pc [10 ]
Dover Strait Dover Strait resund, Denmark Japanese Straits Japanese Straits Dover Strait Dover Strait Strait of Gibraltar Japanese Straits Great Belt, Denmark Danish waters 5.18 3.15 0.27 0.49 1.23 1.11 0.95 1.2 1.10 1.30 3.0 Overtaking At bends in lanes Head-on and overtaking Crossings also? Head-on, no traffic separation Head-on, with traffic separation Head on Head on Crossings Crossings, no traffic separation Crossings, with traffic separation
-4
Comment
Date: 12/05/2009
Page 31 of 59
Ship-bridge collisions
Location Pc [10-4]
Great Belt Great Belt East and West Bridge Great Belt East and West Bridge Tasman Bridge 0.4 1.1 3.2 [0.7; 1.0] Traffic regulations, marking of route, detectability Having pilot on board Without pilot on board Visibility, env. conditions, human error, mechanical failure, traffic intensity Traffic density, use of pilots, traffic restrictions
Comment
Sunshine Skyway Bridge, Florida Annacis Island Bridge, Fraser River, British Columbia Sunshine Skyway Bridge, Florida Sunshine Skyway Bridge, Florida Francis Scott Key Bridge Francis Scott Key Bridge Wm Preston Lane, Jr. Men. Bridge, Maryland Wm Preston Lane, Jr. Men. Bridge, Maryland Chesapeake Bay Bridges and Tunnels, Virginia Dames Point Bridge, Florida Dames Point Bridge, Florida Vicksburg Bridge, Mississippi River Huey P. Long Bridge, Mississippi River Greater New Orleans Bridge, Mississippi River Strait of Gibraltar Japanese Straits
0.5 3.6
1.3 2.0 1.0 2.0 1.0 2.0 0.7 1.3 4.1 5.4 2.5 1.3 0.6 1.86
Ships only Barges only Ships only Barges only Ships only Barges only
Knott et al. Knott et al. Knott et al. Knott et al. Knott et al. Knott et al. Knott et al.
Knott et al. Knott et al. Modjeski & Masters Modjeski & Masters Modjeski & Masters
The values of the causation probabilities by Fujii and Mizuki [9] given in the tables above are mean values. Fujii and Mizuki [9] have given the following ranges:
Date: 12/05/2009
Page 32 of 59
log Pc = -4.31 0.35 log Pc = -3.96 0.36 log Pc = -3.89 0.34 log Pc = -3.80 0.26 log Pc = -3.73 0.36
for head-on collisions for collisions in overtaking for collisions in crossing for grounding for collisions with objects
Further, Fujii and Mizuki [9] states that the above given causation probabilities are obtained for a frequency of visibility less than 1 km that is equal to 263 hours pr. year (i.e. 3%). They further state that the influence of low visibility on the causation probability is approximately proportional to the inverse of the visibility. Finally, they suggest to multiply the above given causation probabilities with a factor of 2 if the frequency of visibility less than 1 km. is in the range of 3% to 10%, and a factor of 8 if it is in the range of 10% to 30%.
Figure 13 Fault tree for causation probability Pc for collision against fixed object.
Date: 12/05/2009
Page 33 of 59
5.2.1 Traditional approach The traditional approach for calculation of Pc (i.e. analysing the cause leading to human inaction or external failures) is to formulate a fault tree or an event tree analysis, see Haugen [13], as shown in Figure 13. From this fault tree it is found that the causation probability Pc can be expressed as
PC = X A + (1 X A ) X C1 X C 2
where XA XC1 XC2 is the probability of human failure is the probability of radar failure, which will depend on vessel size, age, nationality, etc. is the fraction of the year with low visibility.
By application of such fault tree analyses for estimation of the causation probability, it is possible to examine the beneficial effect of new bridge procedures, of having a pilot on board, or of introducing a VTS system in certain geographical areas. Olsen et al. [22] studied the effect of a VTS system by an event tree analysis, see Section 6.2.1. When inspecting the above fault tree it is questionable whether the modelling actually captures any of the important failure mechanism relevant for the considered critical situation. Factors that relate to navigational complications are not included in the analysis, although these are of importance for the relevant set of human errors. Moreover, it is seen that human failure contributes with 75% (2.610-4) to the causation probability. The dominance of the human failure is in agreement with observations. However, the Asleep node is the dominant contributor (2.010-4) and it accounts for 60% of the causation probability. Although the dominating cause may be attributed to human errors this does not seem to be correct as high vigilance is expected in confined navigational areas. An important concern of the fault tree modelling is that the human factor model does not capture the relevant tasks that must be considered in the considered critical situation. 5.2.2 Using Bayesian Networks Most practical risk analysis problems are characterised by a large set of interrelated uncertain quantities and alternatives. Within the conventional risk analysis different methods such as fault tree analysis and event tree analysis have been developed to address these problems. A fault tree analysis seeks the causes of a given event, and an event tree analysis seeks the consequences of a given event. The two analysis techniques are supplementary methods, and when applied correctly the formulated model may reveal the entire probability structure of the model. Both fault tree analysis and event tree analysis applied separately and combined have in the past with success been used in the evaluation of the risk of various hazardous activities. Unfortunately, both fault tree and event tree analyses do have their drawbacks. Firstly, it is difficult to include conditional dependencies and mutually exclusive events in a fault tree analysis (a conditional dependency is, for example, the dependence of the visibility on the weather; mutually exclusive events are, for example, good weather and storm). If conditional dependencies and mutually exclusive events are included in a fault tree analysis the implementation and the pursuing analysis must be performed with utmost care. Secondly, the size of an event tree increases exponentially in the number of variables. Thirdly, if the analysis should capture the primary failure mechanism, the global model, which is combined fault trees and event trees, generally becomes so big that it is virtually impossible for third parties (and sometimes even for first parties) to validate the model.
Date: 12/05/2009
Page 34 of 59
Here we advocate for using Bayesian Networks as the risk modelling and analysis tool. A Bayesian Network is a graphical representation of uncertain quantities (and decisions) that explicitly reveals the probabilistic dependence between the set of variables and the flow of information in the model. A Bayesian Network is designed as a knowledge representation of the considered problem and may therefore be considered as the proper vehicle to bridge the gap between analysis and formulation.
Date: 12/05/2009
Page 35 of 59
Weather
Radar status
Day light
Visual dist.
Speed reducti
Maneuv. time
Other alarms
OOW training
Stress level
Bridge
OOW Task OOW radar OOW visual Looking freq Radar freq
Alarm transfer
OOW acts
Figure 14 Example for Bayesian Network for a navigating officer reacts in the event of being on collision course with an object, from Friis Hansen and Pedersen [10].
Date: 12/05/2009
Page 36 of 59
A Bayesian Network is a network with directed arcs and no cycles. The nodes (to which the arcs point) represent random variables and decisions. Arcs into random variables indicate probabilistic dependence, while arcs into decisions specify the information available at the time of the decision. As an example, one node in the network may represent the weather, whereas another may represent the visibility. An arc from weather to visibility indicates that visibility is conditionally dependent on weather; see Figure 14. The diagram is compact and intuitive, emphasising the relationship among the variables, and yet it represents a complete probabilistic description of the problem. For example, it is easy to convert any event tree or fault tree into a Bayesian Network. Conversely, it may not always be an easy task to convert a Bayesian Network into a combined fault tree and event tree, although theoretically possible. A drawback of Bayesian Network is that they require the state space of the random variables (the nodes) to be defined as discrete states. In our above-mentioned example of weather and visibility, the state space of weather may easily be discretised into states as good weather, storm, etc., whereas the state for visibility more naturally would have been defined as a continuous state space. The Bayesian Network modelling does, unfortunately, require the state space of visibility to be discretised in ranges as for example, 0 to 1 km, 1 to 2 km, etc. Although this is mentioned as a drawback, neither fault trees nor event trees offer any better alternatives. A consequence of the discretisation is partly that the result of the Bayesian Network may be sensitive to the selected discretisation, and partly that the calculations involved in the evaluation of the Bayesian Network grow almost exponentially in the number of states of the nodes. The latter is a consequence of Bayesian Networks encodes the entire probabilistic structure of the problem. A focus on the causal relationship among the variables most effectively does the building of a Bayesian Network. This implies that a Bayesian Network becomes a reasonably realistic model of the problem domain, which is useful when we try to get an understanding about a problem domain. In addition, knowledge of causal relationships allows us to make predictions in the presence of interventions. Last, but not least, the model building through causal relationship makes it much easier to validate and convey the model to third parties. We will not give any details here on how Bayesian Networks are analysed. Instead reference is left to Jensen [18] and Pearl [23]. The Bayesian Network described above is taken from Friis Hansen and Pedersen [10] where a comparative risk evaluation of traditional watch keeping and one-watch keeping has taken pace. The results of the modelling were compared to observations, and good agreement was obtained. Here we extend the modelling to also cover ship-ship collisions. 5.2.3 Bayesian Network for ship-ship collisions The network for predicting the causation factor for ship-ship collisions is rooted in the network shown in Figure 14. The Bayesian Network was extended to model two ships, i.e. ship-ship collision situations. The network used for this analysis is presented as Figure 15. It is seen that this Bayesian Network take into account the correlation between the two vessels, that is, they have to detect each other under the same conditions. Although the network appears complicated, the elements from the basic network in Figure 14 are recognised. It is noted that the to more isolated groups in the lower part of the network models the behaviour on the two bridges, whereas the central group in the upper part of the figure models the two vessel that the two vessels has to be detected by each other.
Date: 12/05/2009
Page 37 of 59
Weather
Visibility
Radar status
Radar status B
Speed red A
Speed red B
Radar dist. A
Vessel B
Vessel A
Radar dist. B
Visual dist. A
Visual dist. B
Rel. Speed
Visual time A
Visual time B
Maneuv. time
Traffic intensit
Maneuv. time
OOW train A
Alarms A
Alarms B
OOW train B
Stress A
Stress B
Bridge A
Bridge B
OOW Task A
OOW Task B
OOW radar A
Radar freq A
Look freq A
OOW vis A
OOW vis B
Look freq B
Radar freq B
OOW radar B
Alarm trf A
Alarm trf B
OOW A acts
OOW B acts
Collision
Figure 15 Bayesian Network model for ship-ship collisions accounting for the correlation between the two vessels.
Date: 12/05/2009
Page 38 of 59
Table 1 shows the calculated causation factors for all the combinations of meetings between large vessels with conventional bridge layout and vessels equipped for sole look-out. It is seen that the calculated causation factor for meetings between conventional vessels is found to be, [10], Pc = 9.0010-5. This value can be compared with observed causation probabilities determined from large data sets published by Fujii et al. [9]. These observed values are given in Table 2. In Table 3 the different headings have been weighed to obtain one global causation factor. The result is that the observations indicate that the causation factor is close to Pc = 8.4110-5 That is a factor which is very close to the causation factor Pc = 9.00 10-5 calculated by the Bayesian Network procedure for conventional vessels operating in geographical areas where the frequency of visibility less than 1 km is 3%. The modelling illustrates that it indeed is possible to establish a realistic modelling of the causation probability. Table 1 Causation factors determined by Bayesian Network
Table 2 Causation Probabilities from Fujii and Mizukis observations, Ref. [9].
Log P Head-on Overtaking Crossings Grounding Object -4.31 -3.96 -3.89 -3.80 -3.73
Date: 12/05/2009
Page 39 of 59
This value setting is mainly rooted in the observations Fujii and Mizukis, Ref. [9]
Date: 12/05/2009
Page 40 of 59
Group 1 and 2 in the list above represent the contribution from human errors. Unquestionable, human error is an important cause to navigational accidents perhaps dominant, as it is quoted that human errors account for at least 1 80% of all accidents. More precisely it could be stated that approximately 80% of navigational accidents involves at least some human errors or questionable judgements rounded in organisational factors. What complicates the assessment is that the blame (or cause) for an accident can be allocated in different ways according to the perspective of the investigator. Typically a serious accidents start from basic human errors but the seriousness of the accident is rather a compound of a set of technical failure, operators error, fundamental design errors, and management errors. Therefore, any realistic modelling must provide a detailed representation of human error in order to be successful. Unfortunately, the human error mechanisms differ from technical or environmental cause (viz. the remaining 20%), and are in fact not yet well understood. A major problem in this respect is that there exists no such thing as a recipe for doing a specific task in the right way (e.g. performing a turn). In an examination of a series of manoeuvring simulation that have led to a grounding accident, Thau [29] found that the primary human error leading to the accident often occurred more than 10 minutes prior to the accident. Contrary, technical or environmental causes are generally simpler to model and understand. 6.1.1 Human and Organisational Errors Human errors can be described as actions taken by individuals that can lead an activity (design, construction, and operation) to realise a quality lower than intended. Human errors also include actions not taken, as these also may lead an activity to realise a quality lower than intended. Many people typically think of human error as operator error or cockpit error, in which the operator makes a slip or mistake due to misperceptions, faulty reasoning, inattention, or debilitating attributes such as sickness, drugs, or fatigue. However, there are many other important sources of human error. These includes factors such as management policies which pressure shipmasters to stay on schedule at all costs, poor equipment design which impedes the operators ability to perform a task, improper or lack of maintenance, improper or lack of training, and inadequate number of crew to perform a task. The human error factors range from those of judgement to ignorance, folly, and mischief. Inadequate training is the primary contributor to many of the past failures in marine structures. Also boredom has played a major role in many accidents. Based on a study by Bea [1] of human error factors in marine engineering the following primary factors were identified: Inadequate training Physical limitations Inadequate communication Bad judgement Fatigue Boredom Carelessness Wishful thinking Ignorance Negligence Folly Panic Ego Laziness Greed Alcohol Mischief Violations
Some researchers even argue that 100% of all accidents are due to human error, since poor man-machine interface, failure of instrumentation (should have been checked more properly), under design, etc. all may be attributed as the result of some sort of human error. Any design is the consequence of human decisions.
Date: 12/05/2009
Page 41 of 59
Organisation errors are a departure from acceptable or desirable practice on the part of a group of individuals that results in unacceptable or undesirable results. Primary organisational error factors includes, [1]: Ineffective regulatory requirements Poor planning / training Poor communications Low quality culture Low worker morale Production orientation Cost-profit incentives Time pressures Rejection of information Complex structure Inequitable promotion / recognition Ineffective monitoring Ego Negative incentives Violations
For example, the goals set by the organisation may lead rational individuals to conduct certain operations in manner that the corporate management would not approve if they were aware of their reliability implications. Similarly, corporate management, under pressures to reduce costs and maintain schedules, may not provide the necessary resources required allowing adequately safe operations. Other types of organisation and management procedure that affect the system reliability include, for example, parallel processing such as developing design criteria at the same time as the structure is being designed a procedure that may not be appropriate in economic terms according to the costs and uncertainties. 6.1.2 Human error evaluation To date, four methodologies have been developed or adapted for maritime use. These are: 1. 2. 3. 4. The operator function model (OFM) type of task analysis Cognitive task analysis Skill assessment Error analysis
The OFM task analysis, developed in1986 by Mitchell and Miller, see Rasmussen [26], provides a breakdown of a function (such as avoiding collisions with neighbouring vessels) into the tasks that must be performed. This also includes the information needed to perform each task, and the decisions that direct the sequence of tasks. This type of task description is independent of the automation; that is, the same tasks, information, and decisions are required, regardless of whether they are performed by a human or by a machine. For example, in collision avoidance, other vessels must be detected, their relative motions analysed to determine whether there is a threat of collision, and a decision made regarding how to change own ships course or speed in order to avoid a potential collision. These tasks must be performed regardless of who (human or machine) executes them.
Date: 12/05/2009
Page 42 of 59
The cognitive task analysis method extends the OFM by considering the mental demands that would be placed on a human operator while performing tasks. For example, in order for a human to detect a new ship as soon as it appears, vigilance (sustained attention) and discrimination (the ability to spot a target against the background) are required. The mental demands of analysing the relative motion of the target vessel include plotting a series of target ranges (distance) and bearings (its angular position relative to own ship) and evaluating the ratio of change over time. Hollnagel [14] introduced a task transaction vocabulary that categorises mental demands, such as search, detect, code, interpret, and decide/select. Assigning the appropriate OFM tasks to humans or machines can thereby represent different levels of automation. Then the cognitive impact of automation can be identified by comparing the number and types of cognitive demands placed on the human operator under the different levels of automation. For example, Froese et al. [4] found that when collision avoidance by manual methods was compared to the use of ARPA radar, then virtually all of the computational demands of the manual method had been eliminated through automation. In order to evaluate the impact of automation on training requirements, a skill assessment technique was developed at US Coast Guard [30] by combining the OFM and cognitive task analyses with the Knowledge, Skills, and Abilities (KSA) analysis. The skill assessment is performed by taking each cognitive task (from the OFM/cognitive task analysis) and determining what types of knowledge or skill that is required for the proper performance of a task. The hybrid analysis thereby focuses the knowledge and skill assessment on the task level. For example, when comparing the manual task in collision avoidance of plotting target range and bearing to the automated scenario that displays target information on the ARPA, then the basic knowledge requirements of collision avoidance do not change with automation. However, the procedural requirements change radically. That is, the mariner has to understand the theory behind collision avoidance regardless of the level of automation, but the specific set of procedural knowledge and skills the mariner needs is dependent on the level and type of automation. Application of the described skill assessment technique has allowed both US Coast Guard [30] and Schraagen et al. [28] to distinguish changes in skill level as a result of automation. The studies by Froese et al. [4] and by Scraagen et al. [28] concludes that the way an automated system is designed can also affect the mariners performance. Some automation hides information from the mariner, presenting only what the designer thought was needed. Unfortunately, many system designers do not fully understand the users task, and consequently we end up with less-than-perfect, error inducing designs. By studying the types of errors commonly made by operators, and by understanding the ramifications of these errors (i.e., are they just nuisance errors or can they cause an accident?), important information is gained that further can be used in training and system redesign. Both error analyses adopted in [4] and [28] consisted of interviewing mariners and instructors, and observing the use of automation during routine shipboard operations.
Date: 12/05/2009
Page 43 of 59
Configuration of the considered navigational area Composition of the ship traffic in the area Environmental conditions, such as weather, visibility, current, etc. Configuration of considered vessel, such as main particulars, manoeuvrability, bridge layout and procedures.
Today, many ships have periodically unmanned engine rooms connected by computerised alarm systems to the bridge. Further, microcomputers for accounting, general record-keeping, and e-mail to land-based operations, automated satellite positioning systems (e.g., the global positioning system or GPS), navigation and collision-avoidance systems like electronic charts (ECDIS) and automated radar plotting aids (ARPA). With this boom in technology comes the concern that not all mariners understand how to use the automation effectively and safely. Indeed, there have been several automation-assisted accidents in recent years in which otherwise experienced mariners either did not know how to use the automated system or had trouble using it because of poor system design, Rothblum and Carvalhais [27]. The related human error modelling is best analysed using the cognitive task analysis. In a subsequent subsection the technical aspects of the different electronic systems is described. 6.2.1 Configuration of navigational area System knowledge of the configuration of the navigational area concerns the arrangement of the route in the vicinity of the considered area and identification of all difficulties in following the route before the considered location. Routes in the considered region that crosses the route prior to the considered location may have influence on the navigational safety and may thus indirectly have influence on faults at the considered location. The navigational markings, such as type of buoys that constitutes the routing system, must be identified. Further, presence and configuration of VTS system in the area as well as requirements for having pilot on board is part of the routing system. The purpose of the routing system is to improve the safety of navigation in converging areas and in areas where the density of traffic is great or where freedom of movement of shipping is inhibited by restricted sea-room, the existence of obstructions to navigation, limited depths and unfavourable meteorological conditions. This subsection describes some relevant aspects of the routing system. Navigational route, markings, aids, and restrictions Traffic lane - An area within defined limits in which one-way traffic is established. Natural obstacles, including those forming separation zones, may constitute a boundary. Traffic Separation Scheme - A routing measure aimed at the separation of opposing streams of traffic by appropriate means and by establishment of traffic lanes. Separation zone and lines - A zone or line separating the traffic lanes in which ships are proceeding in opposite or nearly opposite directions; or separating a traffic lane from the adjacent sea area; or separating traffic lanes designated for particular classes of ship proceeding in the same direction. Inshore traffic zone - A routing measure comprising a designated area between the landward boundary of a traffic separation scheme and the adjacent coast, to be used in accordance with the provision of amendment to International Regulations for Preventing Collision at Sea, 1972 (Collision Regulations).
Date: 12/05/2009
Page 44 of 59
Deep-water route - A route within defined limits that have been accurately surveyed for clearance of sea bottom and submerged obstacles as indicated on the charts. Precautionary areas - A routing measure comprising an area within defined limits where ships must navigate with particular caution and within which the direction of traffic flow may be recommended. Navigational complications Complications that may impinge on the operational safety, e.g. bridges, multiple routes, crossing traffic, etc. Local and regional bathymetry Have an influence on the vessel sizes that are able to operate in the area or collide with specific obstacles. The distance from the route to the ground affects Pc for both collision and grounding. VTS system A VTS system is typically present in areas of high navigational complexity where an accurate monitoring or guidance of the vessels in the area is of importance. Typical such areas may be near location of large bridges, areas with high rate of icebergs, or highly trafficked areas. The main effect of a VTS system, for a ship in contact with the VTS system, will be on the selection of route and distribution of ships across the routes. Reportedly, Olsen et al. [22] found that the effect of the presence of a VTS system might reduce the causation probability for ship-bridge collisions by a factor of 2 to 3. VTS systems may consist of the following equipment in different configurations, Olsen et al. [22]:
Radar installations VHF radio and VHF direction finder Closed Circuit Television Infrared Television Presence of a guard ship
A VTS system consisting of only radar, VHF radio and VHF direction finder constitutes the basic system. Closed Circuit Television and Infrared Television are additional equipment. In some areas a guard ship may be attached to the VTS system. Ships participating in the VTS system must if mandatory when entering the VTS area report to the VTS centre via the VHF radio. Local authorities define the requirement to the ship sizes that should participate in the VTS system. According to the IMO regulation it is mandatory for vessels above 300 GRT to have VHF radio onboard. Some of the benefits of a VTS system is that the radar can detect navigational errors and thereby be corrected via the VHF communication. For ships violating the navigational regulations for the area, attempts can be made to establish contact with these over the VHF radio. Presence of Closed Circuit Television or Infrared Television allows for an improved surveillance of the navigation in the approach channels, for instance detecting a ship that omits to turn at a sea buoy or navigates of the channel. Presence of a guard ship may be able to help wandering vessels or like. This, of course, is highly dependent on the location of the guard ship and on the weather conditions.
Date: 12/05/2009
Page 45 of 59
The degree of vessel participating in the VTS system varies considerably for different locations and is highly dependent on the presence of identifiable hazards in the waterway (e.g. fishing boats, icebergs, bridges, etc.). Presence of identifiable hazards increases the degree of participation. It should be noted that ship owners might obtain lower insurance premiums if their vessels participate in local VTS systems. This aspect therefore presents incitement to participate in the VTS system. The following probabilities of a vessel not participating in the VTS system have been extracted from [22].
Participation conditions Mandatory Mandatory in domestic waters Voluntary in domestic waters Voluntary in the Dover Strait Probability of not participating 110-4 110-3 0.01 to 0.4 0.2 Reported by Japanese studies Canadian Coast Guard U.S. Coast Guard U.K. Department of Transport
Moreover, in the event of a vessel not reporting to the VTS system, then almost all (>99%, [22]) vessels respond to a direct call if the VTS system broadcast position, speed and course of the vessel. Some vessels, however, have proved impossible to contact by VHF or from a guard ship. When receiving an advice by VHF from the VTS centre, Olsen et al. [22] also reports that an average of 90% to 95% comply with the VTS advice. It is noted that the compliance is dependent on the nature of the advice and on the credibility of the system with the mariner. Local conditions near the vessel and unknown to the VTS may prevent the ship operator from following the advice. Requirements for pilot on board In some navigational areas it is required that vessels above a specific size must take a pilot on board. Aspect that must be addressed relates to how well does the pilot inform the master of the vessel of navigational plans? What are standard procedures? Are there requirements to the pilot of specific knowledge of the manoeuvrability of the vessel? Etc. 6.2.2 Composition of ship traffic The vessels that operate at general international routes range from traditional sailing ships, leisure crafts and fishing vessels (whose courses are unpredictable) to large tankers that are confined to deep-water routes only. The large diversity in the vessel traffic composition must be taken properly into account. This concerns bulk carriers and tankers in ballast having poor manoeuvrability; container ships with high cruising speed, hard pressed to arrive at their designated terminals just in time. Smaller petroleum, chemical and gas tankers feeding depots around the region, tow-boats and barges requiring plenty of sea-room to manoeuvre, and passenger ferries crossing the considered operational route. Among the shipmasters of these vessels there is a wide variance in the interpretation of safety and the choice of accepting a particular standard, which varies from criteria used, the circumstances and in most cases opinion.
Date: 12/05/2009
Page 46 of 59
In gathering information on the ship traffic, focus will normally be on the commercial traffic, since these always will represent the primary threat to the navigational area. Leisure traffic and local fishing activity, however, can disturb the commercial traffic and thus be a source of errors. The extent and pattern of this type of traffic should be quantified. Type, size, and frequency of vessels operating in the area should be registered. When combined with information of the configuration of the navigational area this information provides guidance of the possibility for performing evasive manoeuvres. In essence, more ships mean more risks! For long-term design purposes forecasting of traffic intensity and composition is important. In this respect local bathymetry provides guidance for limiting vessel sizes, at least with respect to draft. 6.2.3 Environmental conditions The annual conditions for Weather condition, wind variations, cross wind and in sailing direction waves, visibility (fog, precipitation) current variations, cross current and in sailing direction ice conditions Major parts of these aspects were addressed in Friis Hansen and Pedersen [10]. 6.2.4 Configuration of considered vessel Aspects that should be described Vessel type and particulars: speed, profile. Manoeuvrability of considered vessel Layout of Man-Machine interface Number of officers on the bridge Instrumentation: ARPA, ECDIS, GPS, collision avoidance, and track keeping, etc. In the last few years, the Electronic Chart Display and Information System (ECDIS) has emerged as a powerful addition to the modern bridge. ECDIS offers the possibility for major changes in the navigation process and improves the safety and efficiency of maritime operations. By superimposing three items: a chart, the ship's real-time position, and radar on one display, ECDIS has the potential to improve the accuracy of navigation, increase awareness of dangerous conditions, and reduce the mariner's workload. At US Coast Guard, [30], the potential effects of these systems on bridge operations were examined, using the controlled conditions possible on a full-mission simulator. Four issues were examined: the potential of ECDIS to contribute to navigational precision, its potential to reduce navigation workload, the chart features and navigation functions required by the mariner, and the potential contribution of the integration of radar features on ECDIS. The results provided support to the U.S. position on the International Maritime Organization's (IMO) Standards for ECDIS and recommendations for future system design and for the incorporation of the system in bridge operations, [30].
Date: 12/05/2009
Page 47 of 59
Figure 16 Illustration of an ECDIS display, from [30]. Other relevant concerns relates to: Navigational procedures and practice: voyage planning, pre-planning of actions and procedures in the event of evasive manoeuvres. Communication on the bridge Human failures: no action: absence, present but not attentive, attentive but problem not realised Unintended wrong action: situation misunderstood, wrong action chosen, communication problems Intended wrong action: navigational basis (charts) not updated, confusions of buoys and/or landmarks, manoeuvring capabilities overestimated, clearance requirements underestimated (relevant for ship-bridge collisions) Technical failures: loss of propulsion steering system failures radar failure GPS failure
Influence of the effects of automation (ARPA and ECDIS) on navigational functions: voyage planning, collision avoidance, and track keeping. What is the management attitude towards level of detail in voyage planning? Concerns should also be given on how training may affect the situation? Changes in training: less on computation and more on interpretation is needed given the wide usage of ARPA. At USCG [30] the skills assessment and error analysis techniques identified several important types of skill and knowledge that were not fully covered in current internationally recommended training course objectives for ARPA. These same techniques also allowed the development of training course objectives for ECDIS, a relatively new piece of equipment for which no formal training courses exists. In Froese et al. [4] and in Schraagen et al. [28] the cognitive task analysis and error analysis also proved valuable in identifying aspects of the user interface and equipment functionality which were inconsistent with the needs of the crew in the performance of the automated tasks. Taken together, these tools provide a powerful and comprehensive method of identifying the impact of automation on task and training requirements.
Date: 12/05/2009
Page 48 of 59
Chemical tanker Gas tanker Container ship General cargo ship Bulk carrier Ro-Ro cargo ship Passenger ship Fast ferry Support ship
Fishing ship
8. CONCLUSION
A complete procedure has been presented for the analysis of grounding and ship-ship collision rates and the associated damage caused by collisions. The procedure has been applied to analysis of grounding and collision risks for selected geographical areas (Sea of Aaland and Bornholms Gat / Baltic West). These studies have not yet been reported, but reference will be added when these reports are completed. The document also presents a risk-based framework for calculating the causation probability for grounding and collision. The causation factor for ship-ship collision has been calculated using a Bayesian Network model. The result of the analysis was compared to reported causation probabilities and surprisingly good agreement was obtained. In order to be fully complete, however, the modelling needs to be extended. The concern that needs to be addressed was also described in the present document.
Date: 12/05/2009
Page 49 of 59
9. REFERENCES
[1] [2] [3] [4] [5] Bea, R.: The Role of Human Error in Design, Construction, and Reliability of Marine Structures. Technical report, Ship Structure Committee, 1994. SSC-378. Comstock, J.P. and Robertson, J.B.: "Survival of Collision Damage Versus the 1960 Convention of Safety of Life at Sea", SNAME, pp. 461-522. 1961. Cowi: "Risk analysis of oil and chemical spills in Danish waters", Minister of Defence, 147 pages. 2007. (In Danish) Froese, J, Hartung, R. and Schack, C.: Risk Assessment for Solo Watchkeeping at Night under Defined Conditions. ISSUS, Germany. March, 1996. Fujii, Y. Yamanouchi, H and Mizuki, N.: Some Factors Affecting the Frequency of Accidents in Marine Traffic. II: The probability of Stranding, III: The Effect of Darkness on the Probability of Stranding. Journal of Navigation, Vol. 27, 1974. Fujii, Y. and Yamanouchi, H.: Visual range and the Degree of Risk, Journal Of Navigation Vol. 27, No. 2, pp 248- 252, 1974. Fujii, Y.: "Integrated Study on Marine Traffic Accidents", IABSE Colloquium on Ship Collision with Bridges and Offshore Structures, Copenhagen, Vol. 42, pp. 91-98. 1983. Fujii, Y., Yamanouchi, H. & Matui, T.: "Survey on Vessel Traffic Management Systems and Brief Introduction to Marine Traffic Studies", Electronic Navigation Research Institute Papers No. 45. 1984. Fujii, Y. and Mizuki, N.: Design of VTS systems for water with bridges. Proc. of the International Symposium on Advances in Ship Collision Analysis. Gluver & Olsen eds. Copenhagen, Denmark, 10-13 May, 1998. pp. 177-190. Friis Hansen, P and Pedersen, P.T.: "Risk Analysis of Conventional and Solo Watch Keeping Submitted to Int. Maritime Organisation (IMO) Maritime Safety Committee by Denmark at the 69th Session. 1998. Friis-Hansen, P and Cerup Simonsen, B.: "Dogden Feasibility Studie 2001, Aktivitet 3.9, Suheld, December 2001, HLD Joint Venture Gluver, H. and Olsen, D.: Current practice in risk analysis of ship collisions to bridges. Proc. of the International Symposium on Advances in Ship Collision Analysis. Gluver & Olsen eds. Copenhagen, Denmark, 10-13 May, 1998. pp. 85-96. Haugen, S.: Probabilistic Evaluation of Frequency of Collision Between Ships and Offshore Platforms. Ph.D. thesis, Marine Structures, University of Trondheim. Aug. 1991. Hollnagel, E.: A Cognitive Task Analysis of the STGR Scenario. Nordic nuclear safety research (NKS). NKS/RAK-1(96)R3. April 1996. IMO: "Resolution and Other Decisions", Resolution 680-732, London. 1992. Inoue, K.: "On the Separation of Traffic at Straight Waterway by Distribution Model of Ship Paters", J. Nautical Society of Japan, No. 5. 1972. ISESO: "Information Technology for Enhanced Safety and Efficiency in Ship Design and Operation", Danish Maritime Authority. http://www.sofartsstyrelsen.dk/sw1161.asp Jensen, F.V.: An Introduction to Bayesian Networks. UCL Press. 1996. Karlson, M. Rasmussen, F. and Frisk, L: Verification of ship collision frequency model. Proc. of the International Symposium on Advances in Ship Collision Analysis. Gluver & Olsen eds. Copenhagen, Denmark, 10-13 May, 1998. pp. 117121. Larsen, O. Damgaard: Ship Collisions with Bridges The interaction between Vessel Traffic and Bridge Structures. Structural Engineering Documents 4. International Association for Bridge and Structural Engineering. 1993.
[6] [7]
[8]
[9]
[10]
[11] [12]
[13]
[18] [19]
[20]
Date: 12/05/2009
Page 50 of 59
[21] MacDuff, T.: The Probability of Vessel Collisions. Ocean Industry, September 1974. pp. 144-148. [22] Olsen, D. Gotfredsen, H.H. and Fujii, Y.: Risk Reducing Effects of the Great Belt VTS System. 7th International VTS Symposium, Vancouver, Canada, June 1992. [23] Pearl, J.: Probabilistic Reasoning in Intelligent Systems: Networks of Plausible Inference. Morgan Kaufmann Publishers, Inc. 1988. [24] Pedersen, P. Terndrup, Hansen, P. Friis, and Nielsen, L.: Probabilistic Analysis of Collision Damages With Application to Passenger Ro-Ro Vessels. Safety of Passenger Ro-Ro Vessels. Dept. of Naval Architecture and Ocean Eng. Doc. pac001. 1995. [25] Pedersen, P. Terndrup: "Collision and Grounding Mechanics". Proc. WEMT 1995, Copenhagen, Volume 1, pp.125-157. 1995. [26] Rasmussen, B. and Whetton, C.: Hazard Identification Based on Plant Functional Modelling. Ris National Laboratory, Roskilde, Denmark. October 1993. [27] Rothblum, A. M. and Carvalhais, A.B.: Maritime Applications of Human Factors Test and Evaluation. Chapter 15 in the book, Handbook of Human Factors Testing and Evaluation, edited by T. G. O'Brien and S. G. Charlton, Mahwah, NJ: Lawrence Erlbaum Assoc., 1996. [28] Schraagen, J.M.C, van Breda, L., and Rasker, P.C.: Sole look-out during periods of darkness. TNO Human Factors Rresearch Institute. August 22, 1997. [29] Thau, J. Personal communication. Danish Maritime Institute, Denmark. 1999. [30] U.S. Coast Guard: Homepage of Research and Development Center. http://www.rdc.uscg.mil
Date: 12/05/2009
Page 51 of 59
Appendix: Used ship types compared to the types defined in the Lloyds data base The GRISK-program uses internally the following 14 ship types
GRISK Shiptype Crude oil tanker Oil products tanker Ship Code This seems as a large and homogenous group 1 This is chosen because the oil products carried by this type, has 2 different properties than crude oil. Chemical tankers have generally more separate tanks than other tankers. They carry chemicals which in many cases are dissolved 3 in the ocean when spilled. This is why marginal types such as wine and juice tanks are included here. Present a different risk (explosion) than oil and chemical tankers 4 Often fast ships. It could be categorized as just cargo 5 Often older and slower ships. It could be categorized as just cargo 6 This seems as a large and homogenous group 7 This could also be classified as just cargo. It has been chosen 8 because of its special stability problems All ships carrying more than 12 passengers sailing less than 30 9 knots All passenger ships sailing faster than 30 knots. Lloyds do not 10 have a category for this A large group consisting mainly of small and slow work related 11 crafts. However it also includes supply ships, tugs and pilots. They typical sail more randomly than larger ships Most fishing ships do not carry an AIS transponder but from a collision analysis point of view there present could be important. 12 There is of course also the question of whether the ship is fishing or just sailing. This I believe is included in the activity part of the AIS data 13 All other. Includes naval ships Is not relevant for AIS, but from a collision analysis point of view 14 they would be nice to include.
Chemical tanker Gas tanker Container ship General cargo ship Bulk carrier Ro-Ro cargo ship Passenger ship Fast ferry Support ship
Fishing ship
Date: 12/05/2009
Page 52 of 59
N=Number of ships E() = average L=Lpp=perpendicular B=Breadth moulded D=Depth T=Draught Cb=Block coefficient V=Speed Ships are from 1980-
11.3 12.4 12.1 13.9 14.4 14.6 14.5 14.7 14.9 14.9 15.4 14.7 16.1
0.89 0.88 1.00 0.30 0.99 0.96 0.97 1.00 1.00 1.00 0.96 1.00 1.00
Lpp [a;b[ -1-425 0-25 25-50 50-75 75-100 100-125 125-150 150-175 175-200 200-225 225-250 250-275 275-300 300-325 325-350 350-375 375-400 400-425
Ntotal 5207 122 477 900 1079 738 391 968 283 178 68 2 1
E(L) 115 24 39 65 89 111 137 169 178 218 235 264 279
E(L/B) 5.87 3.43 4.81 5.72 5.97 6.24 6.42 5.67 5.97 6.70 5.65 5.28 6.27
E(B/D) 2.03 2.31 2.45 2.32 2.03 2.00 1.96 1.79 1.77 1.62 1.97 2.16 1.87
E(D/T) 1.32 1.25 1.15 1.17 1.25 1.33 1.38 1.49 1.47 1.46 1.48 1.40 1.82
E(Cb) 0.42 0.17 0.19 0.32 0.46 0.56 0.63 0.61 0.70 0.63 0.85 0.94
E(V) 13.2 11.2 9.6 11.2 12.7 13.1 13.8 14.7 14.9 14.9 14.9 15.6 15.0
Bulb pct 0.84 0.44 0.30 0.58 0.89 0.85 0.82 0.98 0.77 1.00 0.98 1.00
Lpp [a;b[ -1-425 0-25 25-50 50-75 75-100 100-125 125-150 150-175 175-200 200-225 225-250 250-275 275-300 300-325 325-350 350-375 375-400 400-425
E(L/B) 6.10 3.57 5.15 5.85 6.40 7.28 6.48 5.84 5.67 5.52
E(B/D) 2.14 2.41 2.31 2.22 2.12 2.02 1.93 1.89 1.97 1.98
E(D/T) 1.25 1.32 1.13 1.12 1.30 1.38 1.36 1.42 1.47 1.45
E(Cb) 0.33 0.21 0.29 0.33 0.34 0.47 0.51 0.64 0.82
E(V) 12.6 11.3 9.7 11.4 12.5 13.4 14.6 15.1 16.8 15.4
Bulb pct 0.68 0.29 0.30 0.57 0.77 0.67 0.93 1.00 1.00 1.00
Date: 12/05/2009
Page 53 of 59
Lpp [a;b[ -1-425 0-25 25-50 50-75 75-100 100-125 125-150 150-175 175-200 200-225 225-250 250-275 275-300 300-325 325-350 350-375 375-400 400-425
E(L) 139 44 63 90 110 142 163 191 215 230 268 278
E(L/B) 5.82 4.63 5.36 5.74 6.14 6.14 6.02 6.17 6.01 5.99 6.07 6.04
E(B/D) 1.95 2.13 2.83 2.21 2.12 1.83 1.66 1.57 1.60 1.73 1.51 1.72 1.78
E(D/T) 1.56 1.74 1.31 1.23 1.37 1.40 1.55 1.66 1.73 1.81 2.22 2.17 2.22
E(Cb) 0.41 0.18 0.23 0.40 0.44 0.56 0.58 0.66 0.57 0.47 0.44 0.32
E(V) 15.1 12.2 10.2 12.2 13.3 14.9 16.0 16.2 16.4 16.6 17.7 19.7 19.6
Bulb pct 0.92 1.00 0.80 0.82 0.94 0.90 0.82 1.00 1.00 0.93 1.00 1.00 1.00
Lpp [a;b[ -1-425 0-25 25-50 50-75 75-100 100-125 125-150 150-175 175-200 200-225 225-250 250-275 275-300 300-325 325-350 350-375 375-400 400-425
Ntotal 4678 101 25 94 370 615 886 669 517 249 264 409 327 105 36 10 1
E(L) 176 23 40 65 89 114 137 162 190 212 237 263 283 315 333 361 376
E(L/B) 6.48 3.47 4.83 5.43 5.67 6.01 6.17 6.09 6.34 6.76 7.37 7.33 8.09 7.29 7.78 7.58 6.67
E(B/D) 1.84 1.76 2.06 2.25 2.08 1.93 1.89 1.89 1.80 1.75 1.68 1.62 1.58 1.76 1.71 1.81 1.87
E(D/T) 1.48 1.54 1.35 1.36 1.37 1.37 1.41 1.44 1.50 1.52 1.59 1.69 1.66 1.72 1.71 1.72 1.89
E(Cb) 0.27 0.26 0.11 0.10 0.24 0.26 0.26 0.31 0.21 0.27 0.32 0.28 0.32 0.18 0.50 0.66 0.63
E(V) 19.8 22.2 11.5 12.5 14.2 16.3 18.6 19.0 20.8 21.4 22.6 24.3 24.4 25.0 25.0 25.0 25.0
Bulb pct 0.99 1.00 0.85 0.96 0.99 0.99 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00
Lpp [a;b[ -1-425 0-25 25-50 50-75 75-100 100-125 125-150 150-175 175-200 200-225 225-250 250-275 275-300 300-325 325-350 350-375 375-400 400-425
Ntotal 9457 413 1101 2690 2894 1180 668 300 208 2 1
E(L/B) 5.89 3.11 4.88 5.62 5.95 6.55 6.73 6.35 6.20 6.53 7.35
E(B/D) 2.03 2.29 2.19 2.02 2.00 2.07 1.94 1.80 1.79 1.75 1.76
E(D/T) 1.40 1.38 1.38 1.46 1.35 1.39 1.40 1.41 1.42 1.51 1.55
E(Cb) 0.22 0.13 0.12 0.10 0.25 0.34 0.36 0.33 0.31 0.76 0.67
E(V) 12.5 11.3 10.0 11.0 12.3 13.3 14.8 15.6 15.7 15.8 20.5
Bulb pct 0.80 0.56 0.64 0.82 0.79 0.69 0.83 0.99 0.98 1.00
Date: 12/05/2009
Page 54 of 59
Lpp [a;b[ -1-425 0-25 25-50 50-75 75-100 100-125 125-150 150-175 175-200 200-225 225-250 250-275 275-300 300-325 325-350 350-375 375-400 400-425
Ntotal 6090 63 160 263 121 212 396 1032 1556 1376 174 255 431 47 4
E(L) 185 44 62 88 111 141 165 182 216 235 262 281 307 327
E(L/B) 6.14 4.07 4.92 5.54 6.22 5.91 6.22 6.03 6.70 6.57 6.02 6.12 6.16 5.64
E(B/D) 1.87 2.18 2.46 2.13 2.19 1.98 1.91 1.88 1.83 1.72 1.90 1.85 1.89 2.02 1.87
E(D/T) 1.41 1.35 1.41 1.42 1.35 1.37 1.39 1.41 1.45 1.40 1.43 1.38 1.36 1.40 1.39
E(Cb) 0.30 0.19 0.18 0.27 0.32 0.20 0.27 0.27 0.33 0.45 0.53 0.32 0.39 0.63
E(V) 14.1 13.5 10.1 10.9 12.1 13.1 13.9 14.3 14.4 14.5 14.3 14.2 14.6 13.9 13.8
Bulb pct 0.96 0.80 0.64 0.78 0.75 0.77 0.96 0.98 0.99 0.98 0.89 0.98 0.99 1.00 1.00
Lpp [a;b[ -1-425 0-25 25-50 50-75 75-100 100-125 125-150 150-175 175-200 200-225 225-250 250-275 275-300 300-325 325-350 350-375 375-400 400-425
E(L/B) 5.24 3.23 4.11 4.52 5.13 5.71 6.00 5.55 5.93 6.70 7.31
E(B/D) 2.42 3.38 3.57 3.63 2.28 2.23 2.05 1.82 1.58 1.38 1.79
E(D/T) 1.91 1.63 1.35 1.44 1.67 1.66 1.90 2.21 2.43 2.60 1.56
E(Cb) 0.16 0.18 0.05 0.10 0.23 0.23 0.18 0.21 0.20 0.25 0.48
E(V) 16.7 10.4 9.5 9.9 14.5 16.5 18.1 18.8 19.3 20.2 20.0
Bulb pct 0.87 0.09 0.23 0.93 1.00 0.96 0.99 1.00 1.00 0.88
Lpp [a;b[ -1-425 0-25 25-50 50-75 75-100 100-125 125-150 150-175 175-200 200-225 225-250 250-275 275-300 300-325 325-350 350-375 375-400 400-425
Ntotal 2999 489 1134 406 225 147 190 175 122 38 33 36 1 3
E(L) 80 22 35 61 86 112 136 162 184 215 238 263 275 303
E(L/B) 4.78 2.96 3.93 4.59 5.21 5.76 5.89 6.26 6.73 6.89 6.96 8.01 7.13 7.69
E(B/D) 2.72 2.90 2.95 3.00 2.60 2.27 2.12 2.37 2.09 2.09 2.54 2.33 3.30 2.64
E(D/T) 1.81 1.83 1.76 1.64 1.83 1.82 2.02 1.89 2.05 2.28 1.82 2.10 1.33 1.90
E(Cb) 0.16 0.09 0.09 0.15 0.19 0.27 0.27 0.25 0.22 0.23 0.31 0.26
E(V) 17.9 18.5 16.8 14.1 15.5 18.0 19.8 21.6 23.4 20.8 21.6 22.9 22.1 24.5
Bulb pct 0.61 0.11 0.19 0.53 0.56 0.81 0.95 0.95 0.99 1.00 1.00 1.00 1.00 1.00
Date: 12/05/2009
Page 55 of 59
Lpp [a;b[ -1-425 0-25 25-50 50-75 75-100 100-125 125-150 150-175 175-200 200-225 225-250 250-275 275-300 300-325 325-350 350-375 375-400 400-425
E(L/B) 3.67 2.82 3.57 3.35 4.45 4.94 5.57 7.24 8.00
E(B/D) 3.00 3.00 2.88 3.71 3.14 2.76 2.22 1.68 1.40
E(D/T) 2.38 1.99 2.44 2.26 2.57 2.77 2.73 2.15 2.51
E(V) 36.2 36.5 35.5 37.3 38.2 38.9 40.5 32.0 30.5
1.00 1.00
Lpp [a;b[ -1-425 0-25 25-50 50-75 75-100 100-125 125-150 150-175 175-200 200-225 225-250 250-275 275-300 300-325 325-350 350-375 375-400 400-425
E(L) 42 22 32 60 83 112 134 160 189 213 235 263 284 308 329
E(L/B) 3.63 2.94 3.39 4.34 4.68 5.25 5.25 5.44 5.79 5.74 5.59 5.64 5.23 5.46 5.29
E(B/D) 2.32 2.27 2.28 2.45 2.41 2.30 2.29 2.46 1.86 1.93 1.89 1.86 1.97 1.95 2.05
E(D/T) 1.33 1.35 1.33 1.28 1.34 1.53 1.57 1.49 1.61 1.48 1.46 1.45 1.39 1.38 1.33
E(Cb) 0.19 0.37 0.12 0.15 0.22 0.22 0.32 0.20 0.42 0.40 0.49 0.37 0.34 0.41 0.82
E(V) 13.0 12.4 13.1 13.0 13.9 13.7 14.7 14.8 12.3 12.7 13.3 13.4 14.5 14.5 15.1
Bulb pct 0.13 0.02 0.05 0.28 0.46 0.60 0.52 0.47 0.80 0.82 0.75 1.00 0.33 1.00 1.00
Lpp [a;b[ -1-425 0-25 25-50 50-75 75-100 100-125 125-150 150-175 175-200 200-225 225-250 250-275 275-300 300-325 325-350 350-375 375-400 400-425
Date: 12/05/2009
Page 56 of 59
Lpp [a;b[ -1-425 0-25 25-50 50-75 75-100 100-125 125-150 150-175 175-200 200-225 225-250 250-275 275-300 300-325 325-350 350-375 375-400 400-425
E(L/B) 4.05 3.09 3.59 4.19 4.10 4.35 5.52 5.22 6.23 5.75 4.60 4.28 4.74
E(B/D) 3.06 3.40 3.12 2.91 3.05 3.14 2.57 2.47 1.94 2.64 2.79 2.00 1.93
E(D/T) 1.61 1.81 1.52 1.48 1.60 1.81 1.39 1.99 1.62 1.62 2.67 1.43 1.46
E(Cb) 0.16 0.15 0.21 0.13 0.13 0.13 0.08 0.32 0.25 0.46 0.31
E(V) 15.4 22.0 14.2 13.3 15.0 16.2 13.8 15.4 18.4 14.3
Bulb pct 0.44 0.11 0.39 0.58 0.54 0.44 0.50 1.00 1.00 0.67
Lpp [a;b[ -1-425 0-25 25-50 50-75 75-100 100-125 125-150 150-175 175-200 200-225 225-250 250-275 275-300 300-325 325-350 350-375 375-400 400-425
Date: 12/05/2009
Page 57 of 59
Date: 12/05/2009
Page 58 of 59
Date: 12/05/2009
Page 59 of 59
Date: 12/05/2009