Risk Modelling for Passages in Approach Channel

Methods of multivariate statistics, stochastic processes, and simulation methods are used to identify and assess the risk measures. This paper presents the use of generalized linear models and Markov models to study risks to ships along the approach channel. These models combined with simulation testing are used to determine the time required for continuous monitoring of endangered objects or period at which the level of risk should be verified.


Introduction
Logistic models are arguably one of the most widely used data analysis techniques.They are the most studied discriminative models [1][2][3].Logit models appear in a variety of forms in applications in biostatistics, epidemiology, economics, marketing research, and sociology.They are used to model the relationship between covariates and various types of discrete outcomes from the ubiquitous binary logit model for a two-level response to the conditional logit and multinomial (generalized) logit models concerning polytomous responses.Nested logit models allow for modeling the sequence of the decision process faced by the grouping alternatives at each stage into nests [4].
Logistic regression was applied to estimate the likelihood of mortality and severe injury in pedestrian casualties by considering the associations of such factors as demographic characteristics, injury characteristics, crash time, location, road environment, traffic control, and traffic conditions [5,6].Logit models were used for analyzing traffic crash severities aim at identifying and quantifying the effects of the factors which affect different crash injury severities [7,8].Ordered logit model were used in past studies about the great variability in conflict judgments by Air Traffic Controllers [9].Polynomial and logistic regression and maritime transportation simulation were basis to developed an oil outflow model for collision and grounding accidents of tankers [10].Navigation safety assessment can be carried out with the use of risk metrics analysis, multivariate statistical methods, stochastic models, and simulation methods, [11][12][13][14][15].
In this paper we present the possibility of using the category of ordered logit models to examine threats to vessels.Stochastic models were also proposed to allow, among other things, classification of ships with reference to the degree of risk of collision.Safe navigation requires knowledge of the whole picture of navigational, hydrometeorological situations and ship's maneuverability.An important factor affecting the safety of the vessel is the proper execution of maneuvers, [15].Weather conditions affect safe performance of ship's maneuvers during entering the port and berthing, especially the vessels with a big surface exposed to wind.Wind direction and force are very important factors.

Safety of Navigation in the Approach Channels
Safety of navigation is a state of a system related to performing certain maneuvers without collisions, accidents in restricted waters.Among all the hazards we can distinguish [12] (i) collision with another vessel or seamark in the fairway, (ii) running aground or hitting the bottom.

Mathematical Problems in Engineering
Probability of activating any of these risks can be represented by the following formula: (see [12]) where  is the index of hazard type,   are sea area parameters,   are vessel's parameters,   are hydrometeorological parameters,   are parameters of the performed maneuver,   are parameters of traffic density, and   are parameters of traffic control system.It is a dependent function; the variables   ,   ,   ,   ,   ,   , which reflect a number of factors, describe certain states of the system the vessel-water area-the environment.
Hydrometeorological conditions have a significant impact on the safety of navigation in restricted areas.You can, amongst other weather factors, include wind, waves, and currents.
Nature of the fluctuations in water level depends on what the cause is.They can be characterized as follows [16,17]: (i) short-term fluctuations in water level caused by (a) tides, (b) positive or negative surges, (ii) seasonal changes in water levels caused by long-term hydrometeorological changes in the sea area (region).
Water levels lowered or raised by wind are taken into consideration by adopting characteristic water level when determining fairways.Basing on these assumptions, procedures meeting the conditions of safe navigation are established to use a given sea area.Wind affects safe navigation in restricted waters and is always described by means of wind direction and force.
Sea currents influence the safety of the ship's maneuver as they act on the part of the hull which is below water.There are three important types of currents and streams in restricted waters [12]: (i) tidal streams, (ii) sea currents, (iii) wind current.
Wind sea is an important factor in determining the safety of maneuvering in restricted waters.Parameters of waves affecting the size of the safe maneuvering area for a vessel are the following: (i) height, length, and period of the waves, (ii) directions of waves flow in relation to the available navigable water, (iii) distributions of direction and wave height in available navigable waters.
In order to determine the effect of variables   ,   ,   ,   ,   ,   , and   on the safety of navigation, it is necessary to adopt quantitative indicators.They are different for different types of sea areas on which navigation takes place (or maneuver is performed).The limit values of these indicators are the basis for assessing the safety of navigation.
There are two types of evaluation criteria of safety of navigation in restricted waters: (i) basic criteria for assessing the safety of navigation determined by a single parameter, (ii) complex evaluation criteria of safety of navigation which take into account the effects of an accident and which are a function of a number of variables.
Navigational risk is a complex criterion of assessing safety of navigation, [13,14,18].In this work a deviation from the approach channel is assumed to be the hazard which could result in hitting the bottom or in collision.The level of risk depends on the distance from the centerline of the fairway.
Let us assume that the classification of vessel's states in terms of the degree of risk will be carried out basing on the value of a random variable  using vector of limits.Elements of limits (thresholds)   = [ 0 ,  1 , . . .,  −1 ] vectors are monotonic, nondecreasing sequence The different levels of risks and therefore states of the system are determined by (2)

Logit and Probit Models of Ordered States
Generalized linear models extend classical linear models.Logistic regression is an alternative to ordinary linear regression especially when we have discrete variables which describe categories in a given classification or states which may hold [19].Sequential logit and probit models are often sequences of binary outcome models [20].
Let us assume that for a vessel at the time  value of a variable   () deciding about assessing the degree of the th threat to be defined for the moment  +  depends linearly on the vector of features x  () [16], where  is a vector of unknown structural parameters of the model ( 3),  has a distribution independent of , and  of the expected zero value and of a constant variance and is defined with distribution function of the logistic distribution.Write the model as Let us divide the features vector into two "subvectors" where x 1 () is a vector of features defining the level of threat for a vessel  at a given moment  and x 2 () is a vector of quantifiable characteristics.
In this situation, it describes only the value of the unobservable variable  * , that if the random variable/component has a continuous distribution, (3), it can be written as where   is the cumulative distribution function of the random variable.An important issue is the way in which the various explanatory variables affect the dependent/explained variable and the extreme results in case of no unit variance of variable .
Let us consider a situation in which the ship enters the approach channel and we do not have any information about this channel.What is the probability that "the vessel makes an error" (navigational error)?Option 1.Let us assume that the vessel is typical and for which the value of an individual effect (resulting from a vessel typeoperating characteristics ( 5)) is zero. Then for a standard distribution , or generally Option 2. In normal and independent distribution  and  their sum V  =   +  is normally distributed with parameters (0;  2  +  2  ).This implies that we can write Of course, the parameter values and the variances are unknown, so you can take advantage of their assessment.
Let us consider a vessel about which we know how far it deviated from the centerline of the fairway in the approach channel in the past and what corrective maneuvers were performed.What is the probability that the vessel will make navigational error in the next period of time ()?
Knowing the vessel's history and being able to determine the current level of risk, we get information enabling to identify the individual effects.Let us write the probability in the form of conditional probability (3): We use the probability  calculated in accordance with Options 1 or 2, respectively, and then we take a threshold value (  ).If  >   we forecast an error and otherwise lack of an error in the next period of time, [16].
Let us consider probit model with delayed response variable as explanatory variable [4]: Since the probabilities in the subsequent periods for the same vessel are not independent, the model cannot be estimated with the presented method.However, conditional probabilities are independent with respect to the current value of the explanatory variable and to the individual effects.
To solve the initial conditions problem we can use the Heckman solution, [21], noting the model (11) in the form: In this model, for the first period an additional static equation defining the value of  1 was used.The explanatory variables in the equation for the first period need not be the same as in subsequent periods.
An assumption that the random variable has a logistic distribution leads to a logit model in which the probability of adopting the value 1 by the dependent variable with given values of explanatory variables is given as Elements of the vector of structural parameters of the model and vector of limit values can be estimated by maximum probability method using numerical procedures.Let us consider ordered model in which the observable dependent variable can take of  different values: The above model is not appropriate if the variable   has values which cannot be arranged according to a particular scale.
Probability of accepting by the variable values of  can be written as or by means of cumulative distribution of the random variable in the form of where   is the cumulative distribution function of the random variable.
In this way we obtain an ordered probit, and the estimation is based on maximizing the probability function.Estimation of the model requires the scale, taking a specific value of the variance of variable  and finding the maximum of the probability function.

Examples of Models Used for Examining Risks in the Approach Channel
The methodology and logic of ordered logit model can be presented using a following example.Let there be  ships in the channel, and assume that each ship's navigator "degree of risk acceptance" affects the probability of being in state .
Sometimes the outcomes in a response variable are perceived as a sequence with stages (2), states .The related probabilities can be written as [20]  0 =  (∑ . . .
In the logit model the  is logistic cumulative distribution function and subscripts , for  = 0, 1, . . .,  indicate the sets of x variables included in states  = 0,  = 1, . . .,  = , respectively.The parameters  can be estimated by dividing the sample into groups according to levels of risk.

Prioritizing Vessels due to the Risk That They Generate.
Let us examine the case of two vessels along the approach channel, and let us ask the question which of them generates the greater threat?Let us assume that the risk of the first ship is described by the random variable  and the other by random variable .Without any information, we assume that both these variables have the same probability distribution with cumulative distribution function () and independent (if ships are not interacting).In addition, without changing the generality, we can assume that both variables are continuous type, ( = ) = 0.
Let us calculate the probability ( > ) We obtained the result stating that without knowledge (for any distribution   ()) the probability of analyzed events for each vessel is the same and is 0.5.Let us consider the following algorithm, Figure 2: (i) let us assume a certain level of , (ii) assessing the ship for example, the first one, during the remaining period of time, we check if  < , (iii) if so, we predict that  < , (iv) otherwise, we have  > .
Let us calculate probability () of a correct assessment using the adopted assumptions about the type of random variables and probability properties we obtain: Since () is the cumulative distribution then for each level  value of () ∈ [0, 1], so let us examine the function () = − 2 +  + 0.5, Figure 1.When calculating we have   () = −2 + 1 = 0 for  = 0.5,   () = −2.This function has the following maximum value for  = 0.5.This means that an optimal level  opt is the median of the distribution (), ( opt ) = 0.5.Then () = 0.75.
The vessel's domain is defined by a rectangle of a length and width corresponding to the length and breadth of the vessel, where -heading, -maximum distance of the domain point from the centerline of the approach channel, Figure 3.The center of the domain, which is described around the ship, coincides with the geometric center of the vessel.A [] was assumed as the maximum distance which is the maximum distance measured from the domain vertices to the fairway centerline.
Simulation studies were carried out in the laboratories of Gdynia Maritime University, using navigationalmaneuvering simulator.This simulator in a very realistic way that represents sea areas and behaviour of vessels which are almost identical with real vessels.It is possible because these models use six degrees of freedom.The model of the vessel used was the model of a loaded LNG carrier of the following characteristics:  = 315 m,  = 50 m, and  = 12 m.The ship was on an even keel.Simulation of the ship's behaviour was carried out in the fairway and was influenced by interfering factors such as wind and swell-constant NE wind, period: 9 s., height of the swell: 1.5 m.
The value of 6 and 8 knots was taken as the nominal speed of the vessel.The wind direction adopted in simulations changed every 45 degrees starting from the .The simulation scenario of the behaviour of the model under way along the approach channel assumes that the initial position of the vessel is in the centerline of the existing fairway in its northern most part and that the line of symmetry of the model and the course over ground overlaps with the centerline of the fairway.The parameters of the generated waves were determined, and the simulation results were obtained for winds (no squalls) of speed up to 15,0 m.
For the assumed parameters of the fairway the distanced is given by where  LAT and  LONG Cartesian coordinates converted from geographic ones, [15].In the simulation, in particular, the course of the function of the distance from the model vessel to the centerline of the fairway with given parameters of influencing factors (hydrometeorological conditions) was examined.
Thanks to the research on the simulator some information was received according to which the following variables were determined: (i)   taking value of 1 ( = ) if the vessel at the time  was found in  class of the distance from the centerline of the fairway and 0 otherwise,  = 0, 1, . . ., 5 (ii)  1 vessel's speed, (iii)  2 wind speed, (iv)  3 wind direction.For  1 ,  2 , and  3 variables logistic models were obtained where eta is given in Table 1.
Test of compatibility  2 determines whether logistic function adequately fits the observed data.Since the value of  is greater than or equal to 0.05, there is no reason to reject the adequacy of the fitted model at least 95.0% of confidence level.

Markov Model of Leaving the Approach Channel.
Let us consider birth-death process (BD) with finite space of states {0, 1, 2, . . ., }.Graph of state of this process is shown in Figure 4.
Birth and death processes is a continuous-time Markov chain for which transitions may take place only between neighbouring states (i.e.,  to  + 1 or  − 1 only).

State transition probability:
(i)  01 = 1 (ii) birth before death: (iii) death before birth: Moreover let   () be the number of up state changes (from  =  to  =  + 1) during the interval [0, ] when the system is in state ; let   () be the number of down state changes (from  =  to  = −1) during the interval [0, ] when the system is in state ; and let   () be the total time during the interval [0, ] in which the system is in state .Then, [22] we have This estimation procedure need not produce an irreducible BD process, because there can be initial and final transient states.However, under the simplifying assumption of irreducibility, this estimated BD process has the unique stationary probability distribution.For example if there are some constants  1 ,  2 ,  1 , and  2 such that   () > 0 for  1 ≤  ≤  2 with   () = 0 otherwise and   () > 0 for  1 ≤  ≤  2 .From [22], this estimated BD process has the unique stationary probability distribution where  0 () = 1 and   () = ∏ The properly selected sets of explanatory variables in the model (3) allow you to define Markov chains with different transition probability matrix.If we assume that the set of variables was limited to two variables that as a result of the estimation of the logit model, we obtain the assessment of probabilities that can be interpreted as the transition probability of a one threat category to another category.The use of the full set of explanatory variables of the form (5) gives us the transition probabilities for any ship  and a time .

Conclusion
The study proves suitability of logit models used to estimate the hazards to navigation.They can be used in the development of guidelines for traffic safety management system and large ships manoeuvring in restricted waters, to ensure obtaining the assumed level of safety.The approach proposed in this work can be used in Markov models by using in estimations the transition matrix of other variables than just defining the transition between states.This will allow the use of semi-Markov models to estimate safety of navigation.Immediate work in simulation follows better evaluation measures and improvement of duration modeling, model system, and model confidence levels.The presented approach has to be further developed by more comprehensive experimental evaluations, examples of applications, and analytical models relating selected simulation responses with model parameters.
Discrete event models allow inclusion of individual variables without creating compound states, which could improve the precision of the model.Some interesting questions are still open.For example possible questions can relate to correlation between random variables.

Figure 2 :
Figure 2: The block diagram for the presented algorithm.

Figure 3 :
Figure3: Determination of the maximum distanced for the ship along the approach channel,[15].

Figure 4 :
Figure 4: Graph of the system states.

)
Consider a system in which arrivals and departures occur one at a time.Let () be the risk level of the system at time .We now consider fitting a BD model to data collected over an interval [0, ].Let   () and   () be direct estimates of the birth rates and death rates based on sample averages over the time interval [0, ].Similarly, let   () be estimates of the stationary distribution based on sample averages over the time interval [0, ].
System with Blocking.Let   be the probability that a ship will enter the higher class  + 1 if it is already at class . + ) ( [ −1 ] +  [  ]) 2 .