The Model of Optimal Allocation of Maritime Oil Spill Combat Ships
Abstract
:1. Introduction
2. Problems of Combating Oil Spills
2.1. Behavior of Oil in the Sea Water
2.2. Problems of Optimizing the Allocation of Response Resources. Modeling of Spills
3. Allocation Optimiation Methods; Solving Allocation Problems
Evolutionary Algorithms
- Individual—a basic unit of evolution, arriving in a certain environment, to which it should be more or less adapted (example solution of the task, one of the possible configurations of ship deployment in ports);
- population—a pool of individuals living in the common environment and competing for its resources (collection of possible configurations for ship deployment in ports);
- phenotype—external characteristics of the individual. In evolutionary algorithms, these are the parameters of the solution to be assessed (cost of pollution combat for a given ship allocation);
- genotype—unambiguous description of an individual contained in its genes (assignment of ships to specific ports);
- chromosome—a place where the genotype of an individual is stored;
- encoding solutions—a way to save any acceptable solution of the problem in the form of an individual genotype; each solution must be able to be recorded in the form of a genotype.
- Initiation P0—the base population is filled with randomly generated individuals. A fitness value is calculated for each of them. After the basic population has been prepared, the main program loop takes place, in which an artificial evolutionary process is defined.
- Reproduction—copying to temporary population Tt of randomly selected individuals from the basic population. Individuals with a higher adaptation function are more likely to reproduce.
- Genetic operations—Tt temporary population individuals are subjected to genetic operations (cross and mutation).
- Randomly selected individuals are selected in pairs, and a decision is made (randomly and with a certain probability) about crossbreeding (Figure 1). The resulting offspring replace their parents. The probability of crossing pc is an important parameter of the algorithm.
- Then, with a certain probability of mutation, which is an algorithm parameter, mutations can be performed on individuals (Figure 2). A randomly selected gene is selected from the population, in which the perturbation of the genotype takes place, and a randomly selected gene is replaced by another.
- In the next step, the resulting solution is assessed using a target function.
- Selection procedure. On the basis of the calculated value of the target function, the adjustment value is determined (very often the same function). Then, N individuals of the temporary population are randomly selected from among the N individuals of the population of the progeny (with repetitions) by means of the algorithm “roulette circle”. Individuals with better fitness are given a higher weightage in the roulette wheel (Figure 3).
- In the next loop cycle, the resulting population becomes the base population (return to reproduction).
4. Model of Optimal Allocation of Oil Spill Response Resources
4.1. General Method for Optimizing the Allocation of Oil Pollution Response Resources
4.1.1. Assumptions and Problem Definition
4.1.2. Input Data
4.1.3. Formulation of Objective Functions
- j—the index of ports (j = 1, 2, …, J);
- i—the index of ships (i = 1, 2, …, I);
- k—the index of spills (k = 1, 2, …, K).
- J—total number of ports;
- I—total number of ships;
- K—total number of spills;
- pj—capacity of the i-th port [ships];
- cpj—the cost of berthing of the ship in port [PLN/day];
- vj—speed of the ship [kn];
- vrk—the size of the k-th spill [t];
- tijk—the time of arrival of the ship to the spill [h];
- tik—operating time of the ship at the action scene [h];
- cik—the labor cost of the ship at the action site [PLN/h];
- djk—distance between ports and spill sites.
- total (total) cost of arrival at spill from berth;
- the total cost of oil slick clean-up by the vessels involved;
- the cost of environmental damage;
- the cost of maintaining ships in port:
- tijk—the time of arrival of the ship to the spill [h];
- cik—the labor cost of the ship at the action site [PLN/h];
- xij—decision variable (if the i-th ship is assigned to a j-th port xij = 1, otherwise xij = 0);
- yik—decision variable (if the i-th ship is assigned to a k-th spill yik = 1, otherwise yik = 0).
4.2. Method for Optimal Allocation of Oil Spill Response Resources
- port positions (geographical coordinates);
- the current distribution of ships in ports;
- technical characteristics of the ships;
- cost of using ships to combat oil spills [PLN/h];
- geographical coordinates of spill locations;
- type of spilled substance;
- hydrometeorological conditions.
4.3. Optimizing Allocation Using Evolutionary Algorithms
- An allele was selected by uniform distribution, in our case, a ship in the form of number from one to chromosome length (number of ships),
- A value from one to the number of ports was randomized (also evenly distributed) and assigned to a given point of the chromosome.
4.4. Application Interface
5. Experiment Planning
5.1. Hydrometeorological Conditions
5.2. Response Resources for Combat Oil Spill
5.3. Estimated Probable Spill Locations for Simulations
6. Experiment and Analysis of the Results Obtained
6.1. Performing the Experiment
- Input data:
- number (p = 3) and port’s location (geographical coordinates);
- number (r = 3) and spill’s location (geographical coordinates);
- list of available ships for combating oil spill and their parameters;
- the current allocation of ships in ports.
- Analysis of input data:
- determination of the total distance and time at which ships arrive at the spill site for subsequent allocations;
- selection of allocation with the shortest total time of arrival of ships to the place of action;
- Selection of the solution (allocation) with the lowest value of the aim function.
- Copy the best (or selected) Stage I solution (allocation);
- Performance of oil spill combat action simulation for selected allocation in the PISCES II simulator for spill locations 1, 2, and 3 (r = 1, 2, 3) and:
- good hydrometeorological conditions;
- spill size 3500 t;
- spill size 7000 t;
- spill size 14,000t;
- bad hydrometeorological conditions;
- spill size 3500 t;
- spill size 7000 t;
- spill size 14,000t;
- Copy the next solution from the first step (return to stage II, first point); (the second stage can be carried out any number of times, as long as the solutions obtained in the first stage are exhausted, by selecting the subsequent allocations obtained in the first stage of the model).
- Summary of results (costs) of oil pollution combat actions carried out for all analyzed allocations.
6.2. Analysis of Costs of Combating Spills Generated by Simulation Scenarios
6.3. Statistical Analysis of Research Results
- Allocation of resources taking values according to location of ships as: A = {0; 1; 2; 3}.
- Spillage site with values M = {1; 2; 3}.
- Spillage size taking the values W = {3; 5; 7; 14}.
- Hydrometeorological conditions with two values of H = {bad, good}.
6.3.1. Application of a Multifactor Variance Analysis Model
- samples shall be taken from independent observations,
- normality of the distributions from which the samples come (the Shapiro–Wilk test was used),
- homogeneity of variance between groups (the Levene test was used).
- If a > p, then at the materiality level a, the zero hypothesis (H0) should be rejected.
- Otherwise—i.e., when a < p—at importance level, there are no reasons to reject the zero hypothesis (H0).
6.3.2. Variation Analysis
- the conditions influence the cost of spills removal,
- spillage size clearly influences the cost of disposal,
- the allocation does not show a statistically significant impact on the cost of spill removal.
6.3.3. Summary of Statistical Data Analysis
- The hypotheses about the equal average cost of spillage removal for all of the examined factors cannot be accepted.
- The allocation does not have a statistically significant impact on the cost of the action.
- The size, conditions, and location of the bottling operations have a significant impact on the cost of the combat action.
- Costs rise significantly with any increase in the volume and deterioration of conditions.
- There is no statistically significant difference between average costs for the five allocations tested, including the current and suboptimal ones based on the EA.
- In spite of this, the lowest average cost of spillage combat action was observed for the 0 and 1 allocations, which were significantly different from the other allocations.
7. Conclusions
- application of any base of resources;
- optimizing the distribution of any number of ships in any number of ports at any number of spill locations;
- an assessment of the costs arising from any allocation of ships;
- support for the establishment of contingency plans.
Author Contributions
Funding
Conflicts of Interest
References
- Cai, L.; Yan, L.; Ni, J.; Wang, I.C. Assessment of Ecological Vulnerability under Oil Spill Stress. Sustainability 2015, 10, 13073–13084. [Google Scholar] [CrossRef]
- Kang, J.; Zhang, J.; Bai, Y. Modeling and evaluation of the oil-spill emergency response capability based on linguistic variables. Mar. Pollut. Bull. 2016, 113, 293–301. [Google Scholar] [CrossRef] [PubMed]
- Cocârţă, D.; Stoian, M.; Karademir, A. Crude Oil Contaminated Sites: Evaluation by Using Risk Assessment Approach. Sustainability 2017, 9, 1365. [Google Scholar]
- Fingas, M. The Basics of Oil Spill Cleanup, 2rd ed.; Lewis Publishers: Boca Raton, FL, USA, 2001. [Google Scholar]
- Fingas, M.F. Oil Spill Science and Technology Prevention, Response, and Clean up; Gulf Professional Publishing: Burlington, MA, USA, 2011. [Google Scholar]
- Ornitz, B.; Champ, M. Oil Spills First Principles Prevention and Best Response; Elsevier: Amsterdam, The Netherlands; New York, NY, USA, 2002. [Google Scholar]
- National Research Council (U.S.), Red. Oil in the Sea III: Inputs, Fates, and Effects; National Academy Press: Washington, DC, USA, 2003. [Google Scholar]
- Psaraftis, H.N.; Tharakan, G.G.; Ceder, A. Optimal response to oil spills. Oper. Res. 1986, 34, 203–217. [Google Scholar] [CrossRef]
- Iakovou, E.; Ip, C.M.; Douligeris, C.; Korde, A. Optimal location and capacity of emergency cleanup equipment for oil spill response. Eur. J. Oper. Res. 1996, 96, 72–80. [Google Scholar] [CrossRef]
- Lehikoinen, A.; Luoma, E.; Mäntyniemi, S.; Kuikka, S. Optimizing the Recovery Efficiency of Finnish Oil Combating Vessels in the Gulf of Finland Using Bayesian Networks. Environ. Sci. Technol. 2013, 47, 1792–1799. [Google Scholar] [CrossRef] [PubMed]
- Verma, M.; Gendreau, M.; Laporte, G. Optimal location and capability of oil-spill response facilities for the south coast of Newfoundland. Omega 2013, 41, 856–867. [Google Scholar] [CrossRef]
- Spaulding, M.L. State of the art review and future directions in oil spill modeling. Mar. Pollut. Bull. 2017, 115, 7–19. [Google Scholar] [CrossRef] [PubMed]
- Transas. PISCES II Potential Incidents, Control and Evaluation System. 2008. Available online: http://www.transas.com/about/press/news/16531 (accessed on 29 June 2018).
- ASA. OILMAP: Oil Spill Model and Response System Overview. Available online: http://asascience.com/software/PDF/OILMAP_Suite.pdf (accessed on 29 June 2018).
- Marta-Almeida, M. Efficient tools for marine operational forecast and oil spill tracking. Mar. Pollut. Bull. 2013, 71, 139–151. [Google Scholar] [CrossRef] [PubMed]
- SMHI’s Contingency Centre. Seatrack Web Manual. 2012. Available online: https://www.smhi.se/polopoly_fs/1.15599!/Seatrack%20Web%20manual.pdf (accessed on 29 June 2018).
- Łazuga, K.; Gucma, L.; Perkovic, M. M/t “Baltic Carrier” accident. The reconstruction of oil spill with PISCES II simulator application. Zesz. Nauk. Akad. Morska W Szczecinie 2013, 36, 110–115. [Google Scholar]
- Amborski, K. Podstawy Metod Optymalizacji; Oficyna Wydawnicza Politechniki Warszawskiej: Warszawa, Poland, 2009. [Google Scholar]
- Dąbrowski, A. Demokracja I NP-Trudne Problemy; Instytut Matematyczny, Uniwersytet Wrocławski: Wroclaw, Poland, 2015. [Google Scholar]
- Rabbani, M.; Yousefnejad, H. A novel approach for solving a capacitated location allocation problem. Int. J. Ind. Eng. Comput. 2013, 4, 203–214. [Google Scholar] [CrossRef]
- Stachurski, A.; Wierzbicki, A. Podstawy Optymalizacji; Oficyna Wydawnicza Politechniki Warszawskiej: Warszawa, Poland, 2001. [Google Scholar]
- Goldberg, D.E. Algorytmy Genetyczne I ich Zastosowania; Wydawnictwa Naukowo-Techniczne: Warszawa, Poland, 1995. [Google Scholar]
- Wang, S.; Schonfeld, P. Genetic Alghorithms for Selecting and Scheduling Waterway Proje; University of Meryland: College Park, MD, USA, 2006. [Google Scholar]
- Gen, M.; Cheng, R.; Lin, L. Network Models and Optimization; Springer: London, UK, 2008. [Google Scholar]
- Holland, J.H. Adaptation in Natural and Artificial Systems, 1992 ed.; MIT Press: Cambridge, MA, USA, 1992. [Google Scholar]
- Łazuga, K.; Gucma, L.; Juszkiewicz, W. Optimal Planning of Pollution Emergency Response with Application of Navigational Risk Management. Annu. Navig. 2012, 19, 67–77. [Google Scholar] [Green Version]
- EMSA. Inventory of EU Member States Oil Pollution Response Vessels. 2012. Available online: http://www.emsa.europa.eu/emsa-documents/latest/download/2292/487/23.html (accessed on 29 June 2018).
- LAMOR. Lamor Product Reel. Solutions for Selected Scenarios; Lamor Co.: Porvoo, Finland, 2012. [Google Scholar]
- EMSA. Network of Standby Oil Spill Response Vessels and Equipment. Available online: http://www.emsa.europa.eu/news-a-press-centre/external-news/item/1439-network-of-stand-by-oil-spill-response-vessels-and-equipment-handbook-2014.html (accessed on 29 June 2018).
- BRISK. Additional Study. Polish Marine Areas; HELCOM: Helsinki, Finland, 2011. [Google Scholar]
- Gucma, L.; Przywarty, M.; Gralak, R. Evaluation of failure probability of ships on the baltic sea by means of simulation model and statistical data. J. KONBiN 2008, 5, 5–23. [Google Scholar]
- Netka, K. Rozlewy olejowe- jak je ujarzmiać? Środowisko 2010, 13–14, 34–36. [Google Scholar]
- Gucma, L.; Przywarty, M. The Model of Oil Spills Due to Ships Collisions in Southern Baltic area. Available online: http://www.transnav.eu/Article_The_Model_of_Oil_Spills_Due_Gucma,8,124.html (accessed on 29 June 2018).
- Admiralty Sailing Directions: Baltic Pilot Vol. II—NP19. Available online: https://www.amnautical.com/products/ba-sailing-directions-np19-baltic-pilot-vol-2#.WzyEtKsRWUl (accessed on 29 June 2018).
- Hasselström, L.; Consultancy, E.E.E. Costs associated with a major oil spill in the Baltic Sea; Baltic Master II 2011. Available online: https://docplayer.net/6337753-Costs-associated-with-a-major-oil-spill-in-the-baltic-sea.html (accessed on 29 June 2018).
- Śmierzchalski, R.; Michalewicz, Z. Modeling of ship trajectory in collision situations by an evolutionary algorithm. IEEE Trans. Evolut. Comput. 2000, 4, 227–241. [Google Scholar] [CrossRef]
Ship | 1 | 2 | 3 | 4 | 5 | … | I |
---|---|---|---|---|---|---|---|
Port | 1 | 1 | 3 | 3 | 3 | … | … |
Spill Site | Weather | Current | Wind | Water Temp. (°C) | Air Temp. (°C) | Significant Wave Heights (m) | Density (kg/m3) | ||
---|---|---|---|---|---|---|---|---|---|
V (kn) | K (°) | V (m/s) | K (°) | ||||||
Świnoujście | Good | 0.25 | SE | 2.0 | N | 17 | 12–20 | 0.3 | 1006 |
Świnoujście | Bad | 0.25 | SE | 5.7 | N | 17 | 12–20 | 1.0 | 1006 |
Kołobrzeg | Good | 0.25 | SE | 2.0 | N | 17 | 12–20 | 0.3 | 1006 |
Kołobrzeg | Bad | 0.25 | SE | 5.7 | N | 17 | 12–20 | 1.0 | 1006 |
Łeba | Good | 0.25 | SW | 2.0 | N | 17 | 12–20 | 0.3 | 1006 |
Łeba | Bad | 0.25 | SW | 5.7 | N | 17 | 12–20 | 1.0 | 1006 |
No. | Ship’s name | Allocation 0 (Actual) | Allocation 1 (the Best) | Allocation 2 (Middle) | Allocation 3 (Middle) | Allocation 4 (Worst) |
---|---|---|---|---|---|---|
1 | Kapitan Poinc | Gdynia | Kołobrzeg | Świnoujście | Świnoujście | Świnoujście |
2 | Czeslaw II | Świnoujście | Kołobrzeg | Kołobrzeg | Gdynia | Świnoujście |
3 | Zodiak | Gdynia | Kołobrzeg | Kołobrzeg | Gdynia | Świnoujście |
4 | Orkan | Gdynia | Kołobrzeg | Kołobrzeg | Kołobrzeg | Świnoujście |
5 | Passat | Świnoujście | Kołobrzeg | Kołobrzeg | Świnoujście | Świnoujście |
6 | Santa Barbara | Gdynia | Kołobrzeg | Świnoujście | Świnoujście | Świnoujście |
7 | Kambr | Gdynia | Kołobrzeg | Świnoujście | Świnoujście | Świnoujście |
8 | Bazalt | Gdynia | Kołobrzeg | Gdynia | Gdynia | Świnoujście |
9 | Aphrodite | Gdynia | Kołobrzeg | Gdynia | Kołobrzeg | Świnoujście |
1 | spill site number (location) |
3500_IFO180 | spill size_substance type |
dobra | weather conditions |
© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Share and Cite
Łazuga, K.; Gucma, L.; Perkovic, M. The Model of Optimal Allocation of Maritime Oil Spill Combat Ships. Sustainability 2018, 10, 2321. https://doi.org/10.3390/su10072321
Łazuga K, Gucma L, Perkovic M. The Model of Optimal Allocation of Maritime Oil Spill Combat Ships. Sustainability. 2018; 10(7):2321. https://doi.org/10.3390/su10072321
Chicago/Turabian StyleŁazuga, Kinga, Lucjan Gucma, and Marko Perkovic. 2018. "The Model of Optimal Allocation of Maritime Oil Spill Combat Ships" Sustainability 10, no. 7: 2321. https://doi.org/10.3390/su10072321
APA StyleŁazuga, K., Gucma, L., & Perkovic, M. (2018). The Model of Optimal Allocation of Maritime Oil Spill Combat Ships. Sustainability, 10(7), 2321. https://doi.org/10.3390/su10072321