CC BY
20
MONTE CARLO METHOD OF SHIP'S UNDERKEEL CLEARANCE EVALUATION FOR SAFETY OF FERRY APPROACHING TO YSTAD PORT
DETERMINATION
L. Gucma, M. Schoeneich
Type the correct affiliation (Maritime Traffic Engineering Institute, Maritime University of Szczecin, Szczecin, Poland)
e-mail: lucek@am. szczecin.pl martas@am. szczecin.pl
ABSTRACT
The paper is concerned with the analysis of simulation research results of newly design Piast ferry entering to modernized Ystad Port. The ship simulation model is described. After execution of real time simulations the Monte Carlo method of underkeel clearance evaluation is applied to asses the probability of ferry touching the bottom. The results could be used in risk assessment of ships entering to the ports.
1 INTRODUCTION
Underkeel clearance (UKC) is most important factor which determines the possibility of ships hull touching the bottom [5] therefore it is one of the basic elements which decide of navigation safety in restricted waters. The basic navigator's responsibility is to keep safe underkeel clearance in any conditions.
The researches presented in this study are aimed towards determination of UKC of ferries on approach to Ystad Port. The UKC depends on squat which main factor is speed of ships. On approach to Ystad in unfavourable wind conditions speed of ferries cannot be decreased due to lack of maneuverability and this phenomenon is carefully investigated in the paper.
2 YSTAD PORT ENTERING SIMULATION RESEARCH
The main scope of the simulation research consisted of the following issues:
1. Determination of optimal parameters of:
- approach channels to reconstructed port of Ystad with respect to shape, width and depth,
- inner and outer port breakwaters with respect to its shape with respect to waving in port;
- turning places with respect to its shape and optimal depth;
- two new berthing places in inner port in respect to its shape, length, depth, maximal energy of ships contact, maximal speed of ships propeller and bowthruster streams in the bottom.
2. Determination of safety conditions of port operation in respect to:
- admissible meteorological conditions for given kind of ships and maneuvers;
- other navigational conditions and limitations like presence of other ships on berths, use of position fixing systems on approach, navigational markings, vessel traffic service.
3. Determination of maneuvering procedures during berthing and unberthing for different kind of ships and propulsion systems.
4. Determination of underkeel clearance by Monte Carlo method.
5. Determination of usage of main engine during entrance.
6. Determination of ferry distances to the most dangerous objects.
7. Carrying out most typical emergency runs and describe necessary emergency action for the captains.
3 SIMULATION MODEL
The model used in research is based on modular methodology where all influences like hull hydrodynamic forces, propeller drag and steering equipment forces and given external influences are modeled as separate forces and at the end summed as perpendicular, parallel and rotational ones (Gucma L., Gucma M., Przywarty, Tomczak 2006).
The model is operating in the loop where the input variables are calculated instantly (settings and disturbances) as the forces and moments acting on the hull and momentary accelerations are evaluated and speeds of movement surge, sway and yaw. The most important forces acting on the model are:
- thrust of propellers
- side force of propellers;
- sway and resistant force of propellers;
- bow and stern thrusters forces;
- current;
- wind;
- ice effects;
- moment and force of bank effect;
- shallow water forces;
- mooring and anchor forces;
- reaction of the fenders and friction between fender and ships hull;
- tugs forces;
- other depending of special characteristics of power and steering ships equipment.
The detailed diagram of model of the ship maneuvering simulation model is presented in Figure 1 .
Figure 1. The main diagram of simulation model
THE MONTE CARLO METHOD OF UNDERKEEL CLEARANCE ON FERRY APPROACH
The stochastic model of under keel clearance evaluation was presented in (Gucma 2005). It is based on Monte Carlo methodology where overall ships underkeel clearance is described by following mathematical model ( Figure 2):
UKC = (H 0 + - (T + Zst<) + (Aw + Z^) + Sn
where:
SHoi - the uncertainties concerned with depth and its determination, Sn - the uncertainties concerned with draught and its determination, Sswi - the uncertainties concerned with water level and its determination. SN - navigational and manoeuvring clearance.
(1)
Figure 2. Concept of probabilistic underkeel clearance of ships determination
The final model takes into account depth measurement uncertainty, uncertainty of draught determination in port, error of squat determination, bottom irregularity, tides and waves influence are deciding factors for underkeel clearance of ships. Program is capable to consider above mentioned uncertainties using distributions and their parameters. The following parameters are randomly selected from their distributions:
- depth - hi (averaged in following sections -100m, 0m, +100m, 200m, 300m from breakwater),
- sounding error- SBSi,
- mudding component clearance - SZt,
- draught determination error - St ,
- ship's heel error - SPi.
Random draught module
User-entered draught is corrected for draught determination error value and ship's heel error. Iterated draught (Ti) is calculated as follows:
T = T + S +SP (2)
4
where:
T - ships draught [m] assumed as 6.3m, St - draught determination error (assumed as +/-0.05m), SPt - ships heel error (assumed as +/-3 degrees).
Water level module
Water level PWi can be automatically load from online automatic gauges if such exists (Polish solution). In these researches the level was modelled as normal cut distribution with parameters (0, +/-0.1m).
Depth module
Depth hi was assumed as constant in given sections (it varies from 9 before and near breakwater to 8.5m inside the port).
Squat module
Squat (ship sinkage due to decrease of water pressure during movement) is calculated in three stages. First module calculates squat with analytical methods used to obtain moving vessel squat (Huusk, Milword 2, Turner, Hooft, Barrass 1, Barrass 2). Next standard errors of each method are applied. Squat model selection and their standard errors were verified by GPS-RTK experimental research (Gucma, Schoeneich 2006),( Research work 2006). As a result of the experiment uncertainty of each model was assessed and each squat method assigned weight factor. Method's weights and statistical resampling bootstrap method are used later on to calculate final ship's squat.
Underkeel clearance module
Underkeel clearance Zi is determined by using draught, depth, water level and squat results which were calculated before. Underkeel clearance is defined as:
Z = (h + S + SBSt) - (T + O + SN + SWPt + SF ) (3)
where:
hi - up-to-date depth in each iteration in sections (sounding from October 2007),
SZt - mudding component clearance (normal cut distribution with 0 and +/-0.1m),
SBSi - sounding error (normal cut distribution with 0 and +/- 0.1m),
Ti - ships draught with its uncertainty (0,+/-0.05m),
Oi - iterated squat (bootstrap model),
Sn - navigational clearance (0m),
SWPt - height of tide error (0m),
Sf - wave clearance (wave height assumed as h=0.4m before breakwater, h=0.2m in the breakwater and h=0m inside, direction from ferry traverse).
5 GENERAL SIMULATION RESEARCH ASSUMPTIONS
Simulation research is based on performance of maneuvering trial series (inbound/outbound) of significant number for detailed variants. These variants determines given problem. Comparing of results for each variant is done with use of navigational safety criteria.
Variants of researches were determined with following circumstances:
- exploitation conditions of given berth,
- previous researches results (m/f Polonia inYstad port),
- assumptions of analysis,
- investigated area,
- given vessels types,
- navigation conditions,
- manoeuvering tactics.
Simulations were conducted on PC based simulator described which interface is presented in Fig.3. It allows controlling and observing several ships parameters. Captains performed the simulations. In total 6 series were executed each consists of 15 trials. The simulation data were recorded and in the next step used as input data for Monte Carlo model.
Figure 3. Interface of simulation model ferry entering to Ystad
6 RESULT OF MONTE CARLO METHOD OF SHIP'S UNDERKEEL CLEARANCE EVALUATION
Figures 4a and 4b show speed which has the ferry on approach in particular simulation trials in two different wind conditions.
-»-Seriel Serie2 Serie3 Serie4 Serie5 Serie6 Serie7 Serie8 Serie9 Serie10 -•-Seriell Serie12 Serie13 Serie14 Serie15 Serie16 Serie17 Serie18 mean
Figure 4.a. Speed of ferry in knots on approach with E20m/s wind (x=0 outer breakwater)
-«-Serie1 Serie2 Serie3 -*- Serie4 -*- Serie5 Serie6 Serie7
— Serie8
— Serie9 Serie10 Serie11 Serie12 Serie13 Serie14 Serie15 Serie16 mean
Figure 4.b. Speed of ferry in knots on approach with W20m/s wind (x=0 outer breakwater)
The most important result is the probability that clearance is less than zero. This is the probability of accident due to insufficient water depth. The probabilistic model of underkeel clearance (Gucma, Schoeneich 2007) was used for these probability determinations. Histograms of simulated UKC by Monte Carlo method in different parts of Ystad Port are presented in Figure 5a and 5b.
PREDICTED UNPFRKFFI. CLEARANCE
UKC Lm]
Figure 5.a. Histogram UKC and squat value 130m before heads /sea wave =0.4 m/
PREDICTED UNDERKEEL CLEARANCE
UKC [rn]
Figure 5.b. Heads: Histogram UKC and squat value /sea wave =0.2 m/ Distributions show that probabilities of touching the bottom in all cases are near zero.
[m] 1wej1E20_2
breakwater
2wej1W20_2
________ ,2 ---^ --------*
•
—♦—UKC 95% squat
A ——__A
—A
-150 -100 -50 0 50 100 150 200 250 300 350
breakwater [m]
Figure 6. UKC on 95% level of confidence of Piast ferry approaching with E20m/s and W20m/s wind (x=0 means outer breakwater)
7 CONCLUSIONS
The presented methodology of UKC calculation can be used for safety determination and risk
analyses of ships in ports and approaches.
UKC calculated by Monte Carlo in Ystad Port by Monte Carlo model is above 1.6m in worst
assumed conditions.
REFERENCES
1. Gucma L., Gucma M., Przywarty M., Tomczak A. 2006, Eksperymentalne okreslanie rozkladow osiadania i trymu promu morskiego m/f Jan Sniadecki na podejsciu i w porcie Swinoujscie metodami RTK Proceedings of the XV-th International Scientific and Technical Conference "The Role of Navigation in Support of Human Activity on the Sea", Gdynia.
2. Gucma L. 2005, Risk Modelling of Ship Collisions Factors with Fixed Port and Offshore Structures, Maritime University of Szczecin, Szczecin.
3. Gucma L., Schoeneich M. 2006, Okreslanie niepewnosci metod osiadania na podstawie badan rzeczywistych przeprowadzonych na promie m/f Jan Sniadecki w porcie Swinoujscie, Proceedings of the XV-th International Scientific and Technical Conference "The Role of Navigation in Support of Human Activity on the Sea", Gdynia.
4. Gucma L., Schoeneich M. 2007, Probabilistic Model of Underkeel Clearance in Decision Making Process of Port Captain. Monograph- Advances in Marine Navigation and Safety of Sea Transportation TransNav'2007, Gdynia.
5. Gucma S., Jagniszczak I. 2006, Nawigacja morska dla kapitandw, Szczecin.
6. Research work 2006, Determination of squat of m/f Sniadecki by RTK method in Swinoujscie Port. Maritime University of Szczecin, Szczecin.
