Next Article in Journal
The Impact of Urban Recreation Environment on Residents’ Happiness—Based on a Case Study in China
Previous Article in Journal
Identifying Challenges and Solutions in Cultural Heritage Adaptive Reuse through the Historic Urban Landscape Approach in Amsterdam
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Holistic Multi-Objective Design Optimization Approach for Arctic Offshore Supply Vessels

by
Aleksander A. Kondratenko
1,*,
Martin Bergström
1,
Aleksander Reutskii
2 and
Pentti Kujala
1
1
Department of Mechanical Engineering, Aalto University, FI-00076 Aalto, Finland
2
Department of Fleet Operations, Central Marine Research and Design Institute, 191015 St. Petersburg, Russia
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(10), 5550; https://0-doi-org.brum.beds.ac.uk/10.3390/su13105550
Submission received: 22 March 2021 / Revised: 6 May 2021 / Accepted: 13 May 2021 / Published: 16 May 2021

Abstract

:
This article presents a new holistic multi-objective design approach for the optimization of Arctic Offshore Supply Vessels (OSVs) for cost- and eco-efficiency. The approach is intended to be used in the conceptual design phase of an Arctic OSV. It includes (a) a parametric design model of an Arctic OSV, (b) performance assessment models for independently operating and icebreaker-assisted Arctic OSVs, and (c) a novel adaptation of the Artificial Bee Colony (ABC) algorithm for multi-objective optimization of Arctic OSVs. To demonstrate the feasibility and viability of the proposed optimization approach, a series of case studies covering a wide range of operating scenarios are carried out. The results of the case studies indicate that the consideration of icebreaker assistance significantly extends the feasible design space of Arctic OSVs, enabling solutions with improved energy- and cost-efficiency. The results further indicate that the optimal amount of icebreaking assistance and optimal vessel speed differs for different vessels, highlighting the motivation for holistic design optimization. The applied adaptation of the ABC algorithm proved to be well suited and efficient for the multi-objective optimization problem considered.

1. Introduction

As per York and Bell [1], it took more than 50 years for oil and gas to replace oil and wood as the primary energy sources. Although there have been significant advances in science and technology since the 20th century, the ongoing transition from fossil fuel to renewable energy sources is expected to take decades to complete. To power this transition, a steady and as clean as possible supply of oil and gas is needed. In this context, the Arctic shelf, which is estimated to contain some 90 billion barrels of oil and 47 trillion cubic meters of natural gas, is expected to play a significant role [2].
A significant part of the Arctic oil and gas reserves is located under the seafloor in areas featuring thick sea ice for most of the year. The development of such locations requires a complex infrastructure, including an offshore installation, a coastal supply base, and an ice-going support fleet. Arctic Offshore Supply Vessels (OSVs) are a part of the support fleet providing year-round transportation of goods and supplies between the coastal supply base and the offshore installation.
An Arctic OSV is necessarily complex as it must serve both as a cargo ship and icebreaker. Because the acquisition and operation costs of such vessels tend to be high, cost-efficiency is important. Along the lines of Watson [3], a vessel’s cost efficiency can be measured by various metrics, including the Required Freight Rate (RFR), defined as the break-even revenue per transported cargo unit (USD/t). For sustainable operations, an Arctic OSV should not only be cost efficient but also eco-efficient. This is a challenge considering that Arctic OSVs, in comparison with equivalent non-Arctic ships, typically require significantly more propulsion power to enable safe and efficient operations both in open water and in ice. To promote the use of energy-efficient ships, the International Maritime Organization has introduced the Energy Efficiency Design Index (EEDI) regulations [4]. Specifically, the EEDI regulations regulate a ship’s allowable EEDI, indicating the maximum amount of CO2 emissions per unit of transport work (grams of CO2 per ton–nautical mile) [5]. Trivyza et al. [6] encourage the development of new models to calculate a vessel’s lifetime CO2 emissions for EEDI that consider its total lifecycle and non-linear power consumption.
As per Tarovik et al. [7], towards sustainable development of Arctic offshore resources, it is necessary to develop multi-objective optimization of OSVs, considering both internal factors (e.g., the technical characteristics of a ship and related constraints) and external factors (e.g., environmental, logistical, and market-driven factors such as the price of fuel). For the consideration of such factors, Papanikolaou et al. [8] present a holistic approach that contains (a) a parametric model of the vessel considered, (b) a model representing the vessel’s performance in the external operational context, and (c) an optimization algorithm.
Approaches for holistic vessel optimization have been applied successfully on open water ships. For instance, Gutsch et al. [9] and Marques et al. [10] present single-objective optimization approaches for offshore construction vessels and liquefied natural gas carriers. Multi-objective optimization approaches using genetic algorithms have been developed for open water containerships by Priftis et al. [11], for tankers and bulk carriers by Kanellopoulou et al. [12], for Ro-Ro vessels by Skoupas et al. [13], and for high-speed boats by Peri [14]. The approaches by Peri [14] and Martins et al. [15] represent an original adaptation of the genetic algorithm. Ray et al. [16] present applications of the Hooke and Jeeves method and the Rosenbrock method for multi-objective optimization of open water ships.
Sailing ice is associated with ice resistance and ice loads, which may significantly affect a ship’s speed and performance [17]. Therefore, spatial and temporal variations in the prevailing ice conditions are a challenge with regards to providing steady year-round transportation of cargo. The effects of sea ice can be mitigated by ship design measures (e.g., hull-ice strengthening, adaptations of the hull form, and added machinery power). However, increasing a ship’s icebreaking capabilities generally has drawbacks such as increased lightweight, reduced cargo capacity, added open water resistance, reduced seakeeping performance, and increased building and operating costs. An alternative way to mitigate the effects of sea ice is to use icebreaker assistance. However, the use of icebreaker assistance is associated with significant added costs in the form of icebreaker fees. Considering the above, when designing an Arctic OSV, it is necessary to make tradeoffs between multiple different design considerations. To the knowledge of the authors, presently there is no appropriate approach for this purpose.
Recent studies in Arctic maritime engineering provide a theoretical basis to develop a holistic optimization approach for Arctic ships. Bergström et al. [18] investigated the impact of stochastic factors and uncertainties on the design of ships and maritime transport systems for ice-infested waters. Topaj et al. [19] presented an ice routing optimization approach for a ship with icebreaker assistance. Dobrodeev and Sazonov [20] presented a method for estimating the impact of icebreaking assistance on the EEDI. Kondratenko and Tarovik [21] developed a parametric design model of an Arctic offshore support vessel.
It is noted that any new multi-objective optimization approach is an adaptation of some existing single-objective optimization algorithm [14,22]. Multi-objective optimization implies no single optimal solution but a set of Pareto-optimal solutions [8,22]. This means that in most cases, the improvement of one objective results in the deterioration of another. A Pareto front presents the optimization results in graphical form with the objectives plotted on the axis as an outcome of a Pareto optimization. When conducting a Pareto optimization, the main challenge is to obtain a well-informative Pareto front covering the whole design space.
An optimization technique that has been successfully applied to complex practical optimization problems is the so-called Artificial Bee Colony (ABC) metaheuristic, which is based on swarm intelligence, i.e., the collective behavior of a limited collection of interacting agents or individuals [23]. A study by Karaboda and Basturk [23] indicates that the ABC algorithm, when dealing with optimizing multi-variable functions, may outperform more commonly used techniques such as the genetic algorithm. Based on this finding, we think that there is motivation to assess whether the ABC algorithm would be well suited for holistic ship design optimization as an alternative to the genetic algorithm.
Based on this background, this study aims to develop a novel holistic multi-objective optimization approach for maximizing the cost- and eco-efficiency of Arctic OSVs. To this end, the approach includes (a) a parametric design model of an Arctic OSV, (b) performance assessment models for independently operating and icebreaker-assisted Arctic OSVs, and (c) a novel adaptation of the ABC algorithm for multi-objective optimization of Arctic OSVs. In addition, the study aims to analyze how such an optimization approach would affect the design of an Arctic OSV. In fulfilling these objectives, the study addresses the following research questions (RQ):
  • RQ 1: Is it feasible to apply the ABC algorithm for holistic multi-objective optimization of an Arctic vessel?
  • RQ 2: How does the consideration of icebreaker assistance affect the optimization of an Arctic OSV?
  • RQ 3: How does the assumed speed profile of an Arctic OSV affect the outcome of the vessel optimization process?
  • RQ 4: How do variations in the cargo capacity parameters (deadweight or cargo deck area), discount rate, or operation period affect the optimization results?

2. Materials and Methods

2.1. Optimization Process

The proposed optimization approach solves a mixed-integer non-linear programming (MINLP) problem in a multi-objective formulation, applying an adaptation of the single-objective ABC algorithm. The ABC algorithm is based on the foraging strategies of honeybees [23]. A general flow chart of the ABC algorithm is presented in Figure 1. A detailed description of the algorithm is presented in [23]. Each run of the ABC algorithm results in a Pareto front point representing a single optimal solution. A vessel evaluation model supplies the ABC algorithm with an estimation of the quality of a solution in terms of cost- and eco-efficiency. The ABC algorithm operates the model as a black box, which provides information about the solution’s feasibility and evaluates the value of the objective functions considered for each set of inputs.
The objective of the algorithm is to maximize the cost- and eco-efficiency of an Arctic OSV. The cost- and eco-efficiency of a vessel is measured in terms of a cost-efficiency key performance indicator (CKPI) and an eco-efficiency key performance indicator (EKPI), determined by adapting and simplifying the definitions of RFR [3] and EEDI [5], respectively (see Section 2.2). The optimization process results in a Pareto front consisting of a set of Pareto-optimal vessel designs, representing different compromises between cost- and eco-efficiency.
The quality of a Pareto front can be measured in terms of its homogeneity. A homogenous Pareto front is one in which the variations in the distances between consecutive points on the Pareto front are minor. This homogeneity measure is relevant because it represents the level of informativeness of a Pareto front for a decision-maker. In the case of complex scalarized optimization problems, the homogeneity often depends on the ability to find a Pareto front’s non-convex regions. This is because step-wise solution of an optimization problem using the weighted sum method can only specify the convex hull of the Pareto front [22]. In addition, the convex hull obtained will not be sufficiently homogeneous if the value ranges of the optimization objectives differ significantly from each other. This means that the concentration of the Pareto points would be unevenly distributed in favor of the objective with higher values. In Figure 2, we present a general flow chart of the algorithm to overcome these issues.
The proposed optimization process comprises two main phases. The first phase consists of a step-wise scalarized bi-objective optimization using a normalized objective function as per Equation (1). This ensures homogeneity for the convex parts of the Pareto front. The values of CKPImin and EKPImin are obtained in advance by two separate optimization runs. In the first run, CKPI is used as an optimization objective, whereas in the second run EKPI is used as an objective. This provides two boundary points of the Pareto front.
Objective = w · CKPI CKPI min + ( 1 w ) · EKPI EKPI min   ,
where w is a weighting factor defined in a range from 0 to 1. Starting from w = 0, the algorithm increases the value of w step-by-step by adding ∆w and finds an optimal solution, i.e., a new Pareto front point, while w < 1.
The second phase involves processing non-convex regions to provide missing points and homogeneity for the whole Pareto front. If the distance between two points exceeds a specific limit value (average distance + specific percentage × average distance), an additional constraint for either of the objectives is introduced to consider only the range between the points. Specifically, the algorithm divides the range by half and optimizes the unconstrained objective for each half of the range, resulting in two new Pareto front points. This process is repeated until the distance between all the points within the initial range is acceptable. Points that are too close to each other are dismissed.

2.2. Vessel Evaluation Model

The vessel evaluation model, presented in Figure 3, estimates a vessel’s EKPI and CKPI for a given set of inputs. Figure 4 provides a module decomposition of the model. First, the model calculates a vessel’s performance, considering ship design characteristics only. Second, the model calculates a vessel’s performance in a specific operational context. We define the operational context in terms of parameters presented in Figure 4.
The vessel module (see Figure 4) is based on a holistic parametric design model of an Arctic offshore support vessel, a detailed presentation of which is provided in Kondratenko et al. [24] and Kondratenko et al. [21]. The input parameters of the model are length between perpendiculars (Lpp), beam (B), draft (T), block coefficient (Cbwl), depth (H), cruising speed in open water (Vs), ice class (Arc), and icebreaking capability (hice). The ice classes considered are determined as per the Russian Maritime Register of Shipping [25]. The vessel module considers all of the essential ship qualities and related constraints, including the hull geometry, hydrostatics and stability, resistance and propulsion (in open water and ice), required power plant capacity, general arrangement criteria, estimations of lightweight and deadweight, cargo capacity, and freeboard criteria. This study also considers an additional constraint for the maximum tolerated vessel rolling to avoid dangerous situations in adverse weather and sea conditions [26]. The rolling constraint considers related environmental characteristics of the intended operating area.
The parametric model of the vessel is verified using real-life data on Arctic and non-Arctic OSVs, including vessels with installed machinery power of up to 22,500 kW [21]. It is noted that estimations of the required hull volume for vessels with an installed machinery power exceeding 22,500 kW are questionable as the layouts of very powerful vessels (e.g., specialized icebreakers) may differ from that of conventional OSVs.
Following the specification of a vessel’s design characteristics, we calculate the necessary parameters to estimate the values of the objective functions considering the vessel’s operational context. CKPI and EKPI are determined as per Formulas (2) and (3), which we adapted to fit Arctic OSVs. In the present model, Arctic OSVs are assumed to operate year-round with the option to call for icebreaker assistance if necessary. We assume that a vessel is bought at the beginning of the project. Equation (2) determines CKPI considering revenue corresponding to the assumed resale value of the considered vessel at the end of a project. The consideration of the resale value is important as an Arctic OSV is usually bought for a specific offshore project at a significant price, meaning that a vessel’s resale revenue at the end of a project may significantly affect its overall cost efficiency. We assume that unconsidered size-dependent cost categories such as insurance, maintenance, and port fees do not significantly affect the optimization results because of OSV’s size constraints (see Section 3.1.2).
CKPI = 1 n ( PW ( Operating   costs ) + PW ( Ship   acquisition   costs ) PW ( Ship   sale   income ) ) 1 n ( N v C cap )
where PW (present worth) is the value of future money flows discounted to the present, calculated as PW = 1 ( 1 + r ) n ; n is the operation period (years); r is the discount rate; Nv is the number of voyages per year; and Ccap is the cargo capacity parameter (tons). The operating costs are defined as the sum of the fuel and icebreaker assistance costs. The price of icebreaker assistance per nautical mile is calculated according to [27] considering an inflation adjustment. We assume that an icebreaker assists one ship at a time for the whole voyage, meaning that there is a maximum of one instance of icebreaker assistance per voyage. The formation of convoys is not relevant due to the peculiarities of offshore logistics; operating in a convoy suggests a high concentration of vessels at one location, which would be sub-optimal considering logistical efficiency. As a simplification, the model assumes that an icebreaker is always available. This assumption is reasonable as icebreaker operators along the Northern Sea Route tend to prioritize major Arctic offshore projects. We determine the acquisition price of a ship as a function of its lightweight as per [21]. The model calculates a ship’s market value at the end of a project, assuming an annual value loss of 2% due to aging.
The model considers two different cargo capacity parameters Ccap to account for the different possible combinations of different cargo types. If the cargo to be transported is dominated by deck cargo, then Ccap = Sd·mdc, where Sd (m2) is the cargo deck area and mdc is the average density of the deck cargo (tons/ m2). An Arctic offshore installation uses fuel oil for heating, which is typically delivered by a tanker. Alternatively, supply vessels deliver the fuel. In this case, the dominating type of cargo consists of cargo carried within the hull, meaning that the vessel’s cargo carrying capacity is limited by its deadweight (Dw) so that Ccap = Dw.
Because the EEDI regulations concern neither vessels with diesel-electric propulsion nor offshore supply vessels [5], we do not specify the maximum allowed EEDI index in the present study. Instead, the model uses the EKPI index as an optimization objective for eco-efficiency as per Equation (3). It is noted, this definition of EKPI is a simplification of the original EEDI formula by the IMO [5] (e.g., it does not consider ship size-specific coefficients). In addition, contrary to the original EEDI formula that calculates a ship’s fuel consumption as a linear function of the total installed power, the applied model calculates a ship’s fuel consumption based on an approach considering its total lifecycle. This approach, which is described below, makes it possible to define EKPI and CKPI across different vessel sizes.
EKPI = Fuel   consumption   per   hour   · C F C cap * Average   speed   ,
where CF is a conversion factor between fuel consumption and CO2 emission, and Ccap* is the cargo capacity parameter for EKPI (tons); Ccap* = Sd·mdc if the cargo is dominated by deck cargoes and Ccap* = fi·Dw if the cargo is dominated by cargos carried within the hull, where fi is the capacity correction factor for ice-classed ships [5].
We model the ice-going performance of an OSV using a dimensionless quadratic polynomial approximation of an h–v curve based on data provided in [28]. The h–v curve represents a vessel’s speed in level ice as a function of ice thickness at constant engine load. The prevailing ice conditions are modeled in terms of an equivalent ice thickness (heq), determined as per Equation (4) as a function of ice concentration (c), level ice thickness (hi), ice ridging (b), and snow cover thickness (hsn) [29]. The amount of ice ridging is quantified by integer values from 0 to 5, representing the share of ridged ice in a specific area, so that 0 indicates no ice ridging and 5 indicates 100% ridged ice.
h eq = c ( h i + 0.25 bh i + k sn h sn )   ,
where ksn is 0.5 if hsn ≥ 0.5 m and 0.33 otherwise.
A vessel’s h–v curve is typically specified for its maximum continuous rating (MCR) so that it indicates a vessels’ maximum attainable speed in different ice thicknesses. The assumption of operating at MCR in any ice condition is reasonable for low ice class vessels for which the power consumption at cruising speed in open water is close to the maximum continuous rating. High ice-class vessels such as Arctic OSVs, on the other hand, typically have a significant power reserve for icebreaking, resulting in a total installed machinery power that is significantly higher than what is needed to obtain a reasonable speed in open water. Therefore, for such vessels, it is not reasonable to assume continuous operations following an h–v curve calculated for the MCR.
The applied dimensionless approximation (see Equation (5)) makes it possible to draw h–v curves based on three parameters specified for a particular power output (j): attained speed in open water Vmax,j, icebreaking capability hice,j, and the speed at hice,j (V0,j, usually 2–3 knots). The methods from the vessel module [21,24] estimate the required power at various speeds. They also estimate the maximum icebreaking capability for a specific engine output. As a result, an h–v curve can be specified for any engine output.
V j = V 0 , j + ( 1 1.7306 ( h eq h ice , j ) + 0.7306 ( h eq h ice , j ) 2 )   ( V max , j   V 0 , j ) ,
where V j . is a vessel’s speed in ice with an equivalent ice thickness of h eq using power output (j).
As per Figure 5, the algorithm evaluates the maximum attainable speed during independent operation considering the prevailing equivalent ice thickness. If the attainable speed is lower than the set open water cruising speed, the vessel moves at the maximum attainable speed at MCR. Otherwise, the algorithm draws a new h–v curve for a lower engine load equal to the required power for operating at the set cruising speed in open water. Thereafter, the vessel sails at the attainable speed corresponding to the updated h–v curve considering the prevailing equivalent ice thickness. Controlled by the algorithm described above, the vessel operates at MCR only in heavy ice conditions. This prevents the vessel from operating at excessively high speed in ice and open water.
An OSV often operates at low power, for example, during standby or dynamic positioning near an offshore installation. For such operations, modern diesel-electric OSVs have auxiliary diesel generators and power distribution systems providing good fuel efficiency also at low engine load. Rao et al. [30] demonstrated that OSVs equipped with power distribution systems can operate at constant specific fuel consumption (g/kWh) for machinery loads above 2000 kW. Therefore, because all vessels considered in the present study have a minimum engine load above 2000 kW, we assume that all the vessels considered have a specific fuel consumption of 221 g/kW·h regardless of machinery loads.
The algorithm is assumed to provide realistic fuel consumption estimations, providing a fair comparison of vessels with different icebreaking capabilities. For a fair comparison, the vessels’ cruising speed in open water is also included among the optimization parameters. To extend the range of feasible solutions to include vessels with moderate icebreaking capabilities, it is assumed that Arctic OSVs may obtain icebreaker assistance when facing extreme ice conditions. Icebreakers are here modeled as resources that help OSVs to increase their speed without increasing their power output. Icebreaker’s h–v curves are defined for their MCR using the dimensionless approximation described above. The speed of an icebreaker-assisted vessel is estimated according to [31] as a function of the icebreaker’s attainable speed, the prevailing navigational season (summer–autumn/winter–spring), and the ice class of the assisted vessel. This methodology is based on the real-life experience of icebreaking operations along the Northern Sea Route.
The optimization of icebreaking assistance is essential because the optimal level of assistance differ for different vessels depending on their icebreaking capabilities. We optimize the use of icebreaker assistance for cost-efficiency separately considering real-life practice. Because the assisting icebreakers are assumed to be nuclear-powered, the eco-efficiency calculations do not consider their emissions. Separately for each voyage, the algorithm decides whether a ship calls for icebreaker assistance based on an assessment of the total operating costs, calculated as the sum of fuel and icebreaker assistance costs. It is noted that due to limited icebreaking resources, it is not reasonable to assume year-round operation with icebreaking assistance provided for all vessels.
The main outcome of the model is an estimate of a vessel’s total annual fuel consumption (see Equation (6)).
Fuel   consumption = m = 1 12 N m ( s = 1 s max i = 1 i max F m , s , i   p m , s , i + F CO ) ,
where m is the month number; Nm is the average number of voyages in month m; s is the route segment number; smax is the number of segments along the route; i is the ice condition type number corresponding to the equivalent ice thickness range; imax is the number of different types of ice condition occurring along the segment s; Fm,s,i is the fuel consumption in month m for segment s in ice condition i; pm,s,i is the probability of occurrence of ice conditions of type i along the route segment s in month m; and FCO—the fuel consumption for cargo operations. The values of Nm and Fm,s,i are evaluated considering the output of the algorithms for reduced fuel consumption and optimized icebreaking assistance. The total duration of offshore cargo operations is estimated as per [21]. For port and platform cargo handling operations, the fuel consumption is assumed to be two and ten tons per day, respectively. The fuel consumption per hour is calculated by dividing the annual fuel consumption by the total annual operating hours.

3. Results

3.1. Case Study Inputs

To test and demonstrate the merits of the approach, we carry out a series of case studies dealing with the optimization of an Arctic OSV for year-round operation in the northern part of the Kara Sea. Specifically, the considered offshore location to be supplied is the Pobeda field discovered by Rosneft and ExxonMobil in 2014, and the supply base considered is the port of Murmansk. The corresponding supply route is presented in Figure 6a.

3.1.1. Ice Condition

To facilitate the modeling of the ice conditions along the route, we divided the route into eight segments considering the spatial heterogeneity of the ice data. The distance between the waypoints is 100 nautical miles (n.m.) except for the open water segment between the supply base and waypoint 1 (170 n.m.), the segment between waypoints 4 and 5 (the end of the Barents Sea area, 70 n.m.), and the segment between waypoints 7 and 8 (the final segment, 70 n.m.). The ice conditions are modeled in terms of equivalent ice thickness. Ice thickness and concentration data are determined based on historical ice charts provided by AARI [32], an example of which is presented in Figure 6b. The snow cover thickness is modeled based on [33], whereas ice ridging characteristics are modeled based on [34]. The ice data cover the Barents Sea and the Kara Sea over the period 1998–2020 with a temporal resolution of one month. The ice data from the period are assumed to reflect the recent warming of the Arctic climate.
We determined the prevailing ice condition for each segment, assuming that vessels can make minor deviations from the planned route to avoid the most difficult local ice conditions. This assumption is implemented when extracting ice data from the ice charts. Based on the ice data considered, we determined ice thickness frequency distributions for each month and route segment. An example this distribution, representing segment 4–5 in March, is presented in Figure 7. As a simplification, we calculate a vessel’s performance based on the average ice thickness of a thickness range. The results were assembled into a database, which was used as input for the simulations.

3.1.2. Description of the Study Cases

Table 1 presents the parameters of seven different case studies carried out using the proposed multi-objective optimization approach. The outcomes of the case studies indicate the effects of variations in the cargo capacity parameter (Dw vs. Scargo), discount rate, operation period, and icebreaker assistance availability. A discount rate of 6% corresponds to a low-risk offshore project, whereas a discount rate of 12% corresponds to a high-risk project. We consider a fuel price of USD 400 and a 2% annual loss of the value of the considered vessel. The assisting icebreakers are assumed to be modern nuclear-powered icebreakers (hice = 2.8 m, maximum open water speed = 20.6 knots), similar to those currently operating along the Northern Sea Route. A ship’s icebreaking capability (hice) indicates the maximum equivalent ice thickness in which the ship can maintain a speed of 2 knots (V0,j = 2 knots).
Table 2 presents the range of vessel design parameter values used for all case studies. The values are determined considering a database [35] on existing open water and Arctic OSVs built within the period 1997–2017. The maximum vessel length (Lpp) considered is limited as per [7], according to which the performance and safety of an Arctic OSV suffer significantly if its Lpp exceeds the typical linear dimensions of an offshore installation (around 100 m). The maximum draft of 9 m is determined considering the minimum water depth in the operation region.
For some case studies, alternative estimates for EKPI and CKPI are calculated by deactivating the propulsion power control algorithm to demonstrate its effect on the absolute efficiency (EKPI and CKPI) and relative efficiency (compared with other solutions).

3.2. Case Study Outcomes

Figure 8 presents the optimization results for case study 1–2, which have the same discount rate (6%) but different operation periods (10 vs. 20 years). In this and all other cases, the Pareto front is smooth, generally non-convex, and close to equidistant (for the convex and non-convex regions of the Pareto front). For both charts, the EKPI values are in the same range. In the second case study, the CKPI values of the Pareto optimal vessels are considerably lower, as the increase in the amount of transported cargo outweighs the additional costs for the extended operation period. In addition, for a longer operating period (10 vs. 20 years) the effect of discounting is more significant.
Characteristics of the vessel designs on the Pareto fronts calculated for case study 1–2, optimized for CKPI, EKPI, and an intermediate solution, are presented in Table 3, where IButil is the ratio of the total distance covered with icebreaker assistance to the total covered distance. Waterline length (Lwl) is determined as a function of Lpp and a hull form. For each case study, a list of all Pareto optimized designs is presented in Appendix A. All vessels in Table 3 determined for case study 1 have similar values of Lwl, Cbwl, and ice class, while other characteristics are different. Increasing the weight of EKPI in the objective function results in solutions with lower fuel consumption per unit of transported cargo. The corresponding designs are significantly larger, measured in terms of B, T, H, Dw, and displacement (∆). An increase in vessel size also results in higher acquisition costs, something that the EKPI index does not account for. The designs also have lower cruising speeds, which, following the applied power control algorithm, also result in reduced emissions for most ice conditions. In addition to the data provided in Table 3, Figure 9 illustrates how the value of IButil, specified per month, is changing along the year for two specific optimized vessel designs (case study 1, Min CKPI and Min EKPI). Accordingly, for both vessels, the maximum utilization of icebreaker assistance is between February to May. The vessel optimized for CKPI has a higher ice-going capability (hice), which results in lower costs for icebreaker assistance. The vessel with a higher eco-efficiency has a more frequent need for icebreaker assistance to save fuel and limit its required power plant’s capacity (Npp), resulting in a lower lightweight and increased Dw.
The characteristics of the vessel designs optimized for minimum CKPI and EKPI presented in Table 3 are similar for case study 1 and 2. The EKPI optimized vessels must be similar because those parameters that were varied in case study 1–4 impact only the value of CKPI. The latter verifies the reliability of the optimization process as the algorithm demonstrates a stable solution despite its stochastic nature.
Even though the boundary solutions (Min CKPI, Min EKPI) are similar, the intermediary solutions (Table A1 in Appendix A) were affected by variations in the length of the operation period. Specifically, an increase in operation length (case study 2) results in a considerable shift to smaller, cheaper, faster, and less environmentally friendly vessels (the average rise of EKPI is about 6.4%). The average value of displacement and Dw decreased by 7.4% and 7.5% correspondingly, caused by reduced values of B, T, and H (a drop of 4.6%, 2.5%, 2.7% correspondingly). The second cost-reducing factor, besides the acquisition cost, relates to a reduced need for icebreaker assistance, as the average utilization of IB assistance (IButil) is 12.9% lower in case study 2 than in case study 1. The average value of the vessel open water speed is also reduced by 2.1%.
The main reason for this shift is that an Arctic OSV’s resale value at the end of an offshore project is decisive for its overall cost efficiency considering its high acquisition cost. A doubled operation period results in a 20% lower vessel resale value due to aging and discounting, which forces the algorithm to find more CKPI-oriented solutions.
Figure 10 illustrates the Pareto fronts for case study 3 and 4 (the operation periods are 10 and 20 years respectively), which differ from the first pair of cases in terms of discount rate (12%). We provided three vessels’ characteristics, optimized for CKPI, EKPI, and an intermediate solution in Table 4. Increasing the discount rate value from 6% to 12% for a 10-year operating period (moving from case study 1 to case study 3) reduces the operation cost, while the acquisition cost remains high. The vessel purchase occurs before the commissioning, while the discounting effect becomes more significant over time. In this case, the strategy of the algorithm is to minimize a vessel’s size without increasing its negative environmental impact due to higher utilization of icebreaker assistance. This results in a decrease in the values of B, T, H, displacement and hice (2.5%, 4%, 4.2%, 6.1%, 2.3% correspondingly) and an increase in the values of Vs (2.6%) and IButil (5.2%). An additional change in the length of the period of operation from 10 to 20 years for a constant discount rate (shifting from case study 3 to case study 4) results in an insignificant variation of the vessels’ characteristics (<1%).
A demonstration of the importance of the power control algorithm applied is provided for case study 4, representing a high-risk offshore project with duration of 20 years with icebreaker availability. To consider a situation where all vessels operate at the maximum attained speed regardless of the ice conditions (the power control algorithm was deactivated), we calculated the corresponding EKPI and CKPI values for the whole range of Pareto optimal vessel designs (black points in Figure 10b). The results are presented in Figure 10b) as yellow points. As per the figure, the yellow points are shifted upwards compared to the original Pareto front, which indicates a substantial efficiency reduction, and importantly the shift is not equidistant. This means that optimizing an Arctic OSV assuming that vessels are sailing at the maximum attainable speed in ice is not sensible as the optimal speed is significantly different for different vessels. Some vessels are less efficient than others while operating at maximum speed, meaning that they would be modeled unfairly without the power control algorithm.
Case 5–7 suggest a change in the cargo flow pattern, as a consequence of which the required cargo capacity is defined by Scargo instead of Dw. Figure 11 presents the Pareto fronts for Case 5–6. Characteristics of vessels optimized for CKPI, EKPI, and an intermediate solution are provided in Table 5. As the results indicate, the cargo capacity parameter has the most significant impact on the optimized characteristics of an Arctic OSV. The optimal vessels in Case 5–7 are significantly smaller in size than in Case 1–4, with a significantly smaller draft and hull main particulars. The cargo deck area Scargo is associated mainly with Lwl and B when the other particulars of a hull could be minimized if design constraints are satisfied. The optimal vessels in Case 5–7 have lower Dw and internal cargo capacity values than in Case 1–4. A visualization of the optimal hull developed by the model for the intermediate solutions of Case 4 (the cargo capacity parameter is Dw) and 6 (the cargo capacity parameter is Scargo) is presented in Figure 12.
Among the optimal solutions for Case 5–6, most vessels have ice class Arc 5. The exceptions are non-typical vessels optimized for Min EKPI, which were excluded as their Npp values were almost twice the maximum set considering potential general arrangement issues. These vessels have ice class Arc 8, a high icebreaking capability (hice) of almost 2.5 m, and operate at low power when Vs is near 6.5 knots. A constrained maximum speed saves fuel both in ice and open water. The latter is important considering the significant open-water resistance of ships optimized for icebreaking. A high ice class and icebreaking capability provide lower power consumption for independent operation and operation with icebreaker assistance. Although this strategy provides the best EKPI values, CKPI values for the Arc 8 vessels are very high. Thus, as per Figure 11, for case study 5–6, those designs form a small separate group standing far on the right from the Arc 5 vessels.
The yellow points in Figure 11b represent EKPI and CKPI (Case 6) values calculated for all the Pareto optimal solutions (black points in Figure 11b) without using the power control algorithm. The conclusions are the same as for Case 4: the efficiency of the obtained vessels is significantly reduced, but the reduction varies between different vessel designs. This is especially true for vessels of ice class Arc 8, which have a very high machinery power to provide the required icebreaking capability, resulting in high fuel consumption if the vessels operate at maximum speed.
Case 7 is similar to Case 6 except that no icebreaking assistance is provided. The results of the optimization for Case 7 are presented in Figure 13. The characteristics of several optimal vessels (Min CKPI, Min EKPI, Intermediate) are presented in Table 6. In the absence of icebreaking assistance, the optimal vessels have ice class Arc 7 or Arc 8, as vessels of lower ice classes do not guarantee safe year-round, independent operation in the assumed ice conditions. Compared with Case 6, the Pareto front determined for Case 7 is shifted into a less efficient zone, which highlights the importance of icebreaking assistance. Simultaneously, the results for Case 1–7 show that the optimal amount of icebreaking assistance is different for different vessels. The latter demonstrates the importance of icebreaking assistance optimization to ensure a fair treatment of all feasible vessel designs.

4. Discussion

This article presents a framework for holistic multi-objective optimization of Arctic OSVs for cost- and eco-efficiency. The framework is intended to be used in the conceptual design phase of a vessel and consists of three main components: (1) a parametric model of an Arctic OSV, (2) a model representing the performance of an Arctic OSV in a given operational context, and (3) an algorithm for design optimization consisting of an original adaptation of the ABC algorithm for a multi-objective, well-constrained optimization problem. The proposed approach is applied to an Arctic OSV, considering icebreaker assistance and propulsion power control, adapting a vessel’s propulsion power to the prevailing ice conditions.
The proposed adaptation of the ABC algorithm proved to be highly efficient for the multi-objective, mixed-integer non-linear programming (MINLP) problem with constraints. The algorithm treats the model as a black box when most constraints are encapsulated. The formulation represents a highly complex optimization problem and requires the algorithm’s maximum performance in calculation speed and efficiency. The Pareto fronts obtained are smooth, non-convex in general, and close to equidistant for a wide range of case studies.
As demonstrated through case studies, the consideration of icebreaker assistance significantly extends the feasible design space, indicating that disregarding icebreaking assistance in the optimization of an Arctic OSV may result in a suboptimal solution. This also applies to vessels with a high ice class and icebreaking capability, which might use icebreaker assistance for additional speed in difficult ice conditions, even if they would be able to operate independently. The case studies indicate that the overall best performing vessels have a moderate ice class and icebreaking capabilities, which requires that they call for icebreaker assistance in the most challenging ice conditions. The case studies further demonstrate that the optimal amount of icebreaking assistance is different for different vessels, which highlights the importance of considering icebreaker assistance as a part of the optimization to ensure that all relevant vessel operation strategies are given fair consideration.
The modeling of a vessel’s power consumption at reasonable speeds in different ice conditions is an essential component of the optimization problem considered. Propulsion power control may significantly reduce a vessel’s estimated fuel consumption and must therefore be considered in the optimization model. This is particularly important as the impact of propulsion power control is different for different vessels.
We analyzed the sensitivity of the proposed optimization model to variations in external factors, such as cargo capacity (deadweight or the effective cargo deck area), discount rate, and the length of the operation period. Variations in these factors proved to significantly influence the optimized vessel and operation characteristics (power consumption and icebreaker assistance).
In practice, following the actual optimization process, it is necessary to ensure that all the Pareto-optimal vessels suit the logistics of the offshore project considered. Due to the constrained size of OSVs, excess cargo capacity is typically not a major issue. A more important consideration is the cargo storage capacity (cargo deck area and cargo volumes) of the offshore installation, which should not be smaller than that of the OSVs to ensure enough space for the supplied cargo. Thus, OSVs with a cargo carrying capacity exceeding that of the cargo storage of the offshore installation should not be considered.

5. Conclusions

As demonstrated through case studies, the proposed approach is readily applicable for the conceptual design of an Arctic OSV. The approach is thought to be also well suited for other types of ice-going vessels, in particular for vessels that operate as a part of a complex maritime system, such as offshore wind farm vessels or Arctic shuttle tankers, whose efficiency depends significantly on their performance in offshore loading operations in ice-infested waters. The adaptation of the ABC algorithm presented could also be relevant for other types of multi-objective optimization problems utilizing black-box models, which are commonly used in the industry.
Limitations of the approach include: (1) the utility of the framework in its current form is limited to the conceptual design of Arctic OSVs; (2) the approach does not consider non-classic logistic solutions, such as the application of an intermediate floating storage between a supply base and an offshore installation; and (3) CKPI is a simplified version of RFR and does not include all revenue and cost categories.
Topics for future study include: (1) analysis of the impact of wind and waves parameters in open water and ice on the optimization outcome; (2) analysis of the impact of using RFR as a cost-efficiency indicator instead of CKPI; (3) consideration of potential financial losses due to vessel damages caused by ship-ice interactions and how these depend on a vessel’s ice class and the utilization of icebreaker assistance; (4) further development of the applied performance models to consider additional optimization objectives (e.g., safety, ergonomics, recyclability) and economic aspects in greater detail; and (5) adaptation of the approach to other types of vessels.

Author Contributions

A.A.K. (conceptualization, methodology, software, validation, formal analysis, investigation, data curation, writing—original draft preparation, writing—review and editing, visualization); M.B. (writing—review and editing); A.R. (visualization); P.K. (writing—review and editing, resources, supervision, project administration, funding acquisition). All authors have read and agreed to the published version of the manuscript.

Funding

This project has received funding from the Lloyd’s Register Foundation, a charitable foundation, helping to protect life and property by supporting engineering-related education, public engagement, and the application of research www.lrfoundation.org.uk (accessed on 15 May 2021).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We would like to thank the Arctic and Antarctic Research Institute for valuable information deposited in open sources, especially the ice charts and methodologies.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Table A1. Characteristics of the Pareto-optimal vessels for Cases 1–7.
Table A1. Characteristics of the Pareto-optimal vessels for Cases 1–7.
Lwl, mB, mT, mH, mCbwl∆, tVs, knArchice, mCKPIEKPIIButilDw, tSdeck, m2
Case 1:
110.818.46.88.20.635904113.550.8345.441.50.0344756851
110.418.47.59.00.635988513.050.8445.738.70.0365233847
110.818.57.99.50.63510,54112.550.8346.336.90.0375587856
110.718.48.510.20.63511,23212.050.8346.935.50.0385967850
110.818.58.910.90.63511,94012.050.8447.734.70.0396337855
110.821.98.410.20.63513,26612.550.8848.734.50.03970541104
109.920.88.710.60.62912,84312.050.8749.435.00.04067621039
110.122.19.010.90.63414,22812.050.8650.033.40.04175591107
110.318.37.69.20.631989112.050.7750.632.40.0585190842
110.818.58.310.10.63511,13712.050.7750.730.00.0605913854
110.818.58.910.80.63511,82212.050.7751.329.20.0626283852
110.918.59.010.90.63512,00211.550.7651.728.80.0636378854
110.321.08.810.70.63513,33011.550.7752.828.60.06571261050
110.521.69.010.90.63513,95912.050.7753.228.30.06574561082
110.822.58.910.80.63514,46412.050.7653.928.10.06677071135
110.622.59.010.90.63514,56111.550.7654.227.80.06777631130
110.824.88.510.40.63515,29510.950.7655.027.40.07081841349
110.824.79.010.90.63516,01011.550.7655.526.90.07085551342
110.825.09.010.90.63516,16211.050.7755.826.70.07186381358
110.825.09.010.90.63516,18911.050.7755.926.70.07186531360
Case 2:
110.218.56.98.30.635915213.550.8335.041.50.0344820851
110.818.47.08.40.635935212.550.8535.139.90.0354928852
110.518.56.57.70.635857613.550.8335.143.50.0344484851
110.418.47.69.20.63510,07112.550.8335.337.80.0365333845
110.818.58.09.60.63510,59712.550.8435.636.70.0375619852
110.718.48.210.00.63510,93512.550.8335.836.20.0375806850
110.918.58.510.30.63511,30212.550.8336.035.70.0386001854
110.718.48.610.40.63511,43912.550.8336.135.50.0386077851
110.618.48.810.70.63511,73912.050.8436.535.00.0396233849
109.221.08.810.60.63112,98312.550.8637.434.80.03969021039
110.321.78.910.80.63513,91512.550.8537.833.70.04074281088
110.422.89.010.90.63514,66512.050.8738.533.00.04278121142
106.019.28.310.00.63510,98312.050.7739.132.00.0605861848
110.918.58.19.80.63510,80011.550.7739.130.30.0615725854
110.818.58.810.70.63511,72512.050.7739.629.40.0616223853
110.818.68.910.80.63511,94611.050.7740.128.80.0646342858
110.821.38.910.80.63513,72212.050.7740.928.50.06573241072
110.822.48.910.80.63514,43111.550.7741.727.90.06776911129
110.624.68.810.70.63515,60111.550.7742.327.20.06983431335
110.925.09.010.90.63516,18411.050.7743.126.70.07186381361
Case 3:
110.718.56.57.80.635866013.550.8447.142.90.0344534853
110.618.56.98.30.635918913.050.8447.540.60.0354834851
110.818.57.38.80.635976313.050.8447.839.00.0365163852
110.818.57.99.50.63510,53613.050.8348.637.30.0365581853
110.618.57.28.70.635963613.050.7649.833.60.0565095853
110.618.47.48.80.635972413.050.7750.033.30.0565139847
110.518.37.99.60.63510,47212.950.7750.631.90.0585556843
110.818.47.99.50.63510,44112.050.7651.231.10.0595531850
110.818.48.19.80.63510,81412.050.7651.630.50.0605738852
110.618.58.410.20.63511,20912.050.7751.930.00.0605956850
110.818.59.010.90.63511,94112.550.7752.529.40.0616346853
110.818.59.010.90.63511,99712.050.7753.129.10.0626361854
110.719.38.910.80.63512,42611.450.7754.129.10.0646597894
110.721.69.010.90.63513,91712.550.7754.928.70.06474201084
110.821.99.010.90.63514,23112.050.7755.628.10.06675841103
110.824.68.410.20.63514,95812.050.7556.228.00.06779981339
110.824.98.410.20.63515,18611.550.7756.827.60.06881231357
110.624.68.910.80.63515,72411.550.7757.527.10.06984051337
110.824.68.910.90.63515,89011.050.7758.126.80.07184981338
110.925.09.010.90.63516,24011.050.7758.526.60.07186791361
Case 4:
110.518.47.08.40.635922613.450.8431.041.00.0344863849
110.718.57.28.70.635962213.550.8631.340.40.0355070853
110.418.27.69.10.632989613.050.8431.738.70.0365214840
110.818.47.79.30.63510,24212.950.8331.738.00.0365424848
110.918.56.57.80.635875113.050.7732.535.80.0544587856
110.618.47.28.60.635952813.050.7732.633.70.0565034848
110.518.57.38.70.635966012.550.7732.933.10.0575108852
110.318.38.09.60.63510,48512.950.7733.331.90.0585557842
110.318.48.09.60.63510,51011.950.7633.631.00.0595576844
110.318.48.410.20.6351112812.050.7734.030.10.0605912846
110.918.59.010.90.63511,98912.550.7734.529.40.0616370853
110.818.78.910.70.63511,92311.550.7634.929.00.0636337862
110.819.09.010.90.63512,32911.550.7635.329.00.0636548880
110.821.08.910.80.63513,53612.050.7635.928.60.06572211057
110.921.49.010.90.63513,85412.050.7636.228.40.06573831075
110.622.78.610.50.63514,13111.550.7736.728.40.06775271142
110.724.88.410.20.63514,99511.550.7737.127.70.06880191347
110.624.59.010.90.63315,83012.550.7737.527.60.06884391332
110.824.79.010.90.63515,97812.050.7637.827.20.06985261342
110.825.09.010.90.63516,21611.050.7738.426.60.07186651359
Case 5:
108.018.04.56.30.622564613.550.80122.1160.60.0302403814
108.818.34.66.20.624592212.950.84122.3157.20.0312647833
107.818.84.96.00.624636812.450.84123.8154.40.0323104846
109.220.75.26.60.621753113.950.86123.8145.70.03136371028
109.220.45.26.70.621744312.550.86125.0141.60.03235421017
109.421.25.56.60.625812212.550.87127.2139.70.03441151059
109.218.54.76.00.627610812.050.73132.3131.30.0502865843
110.518.44.96.10.632645312.050.76135.3128.30.0513131850
109.620.95.46.70.624794612.550.75138.1121.10.05339541045
108.523.55.97.60.616954312.050.76146.1113.90.05646281261
108.924.26.27.40.61910,28111.550.76148.8111.60.05952821301
108.824.16.57.80.61810,78412.050.77155.4111.10.05955791296
106.224.36.49.40.59810,12911.050.77159.5109.70.05641391380
105.623.76.79.70.59410,19710.450.76166.5109.60.05841611340
106.124.87.39.70.59611,69111.050.75172.9108.50.06053251412
103.825.07.911.10.57712,09911.550.76176.9106.60.05849611504
105.823.56.48.10.61099286.582.44227.699.30.10229551232
Case 6:
109.618.34.66.00.625593314.550.9279.3171.80.0272670841
109.018.64.85.90.628621814.050.8780.6161.20.0302961847
109.020.35.16.80.620723913.450.8682.9145.10.03133521010
109.018.74.96.00.626636212.050.8482.9151.80.0333088850
108.718.34.66.00.624584113.550.7683.1136.00.0472642832
110.118.34.86.00.629626412.550.7685.6129.20.0503011843
109.020.75.46.90.620767611.550.7692.2120.50.05436671028
107.622.85.98.00.610901711.950.7797.9114.70.05440811214
108.823.96.37.60.61810,39611.550.76101.6111.30.05953231287
106.224.46.59.50.59810,29511.050.77108.2109.90.05742071388
106.324.57.09.70.59711,17112.050.77112.2109.30.05748761394
106.024.87.39.70.59511,66410.550.77117.7107.90.06153021410
103.525.08.011.20.57512,15711.050.77121.5106.60.05949861499
107.524.36.98.40.61011,3236.482.50168.598.20.10637261295
105.325.07.810.00.59112,3466.282.47182.097.80.10836781413
Case 7:
108.820.05.07.00.617686715.081.8486.1240.701801991
107.721.95.58.20.608804215.571.8586.3230.7022931166
107.721.95.58.30.608807914.571.8387.0220.5023311165
108.022.45.68.00.611846414.581.8587.9208.7023541195
108.222.95.77.80.612892114.081.8389.9203.0027431225
108.022.85.98.00.611917514.081.8391.8199.7029181220
108.423.96.17.50.615990313.581.8393.4193.6035861282
106.425.06.69.40.59810,78813.581.8496.0184.7031551423
106.524.66.69.40.59910,64713.081.8397.4182.8031001404
106.424.76.79.50.59710,83313.081.8798.5182.7031871408
106.124.56.89.60.59610,77312.981.8499.2182.5031641392
106.424.76.99.60.59711,05613.081.86100.3181.8033201408
106.124.77.09.60.59611,26613.081.89101.7181.6034841407
106.025.07.39.60.59611,74213.081.83104.2181.0038861421
105.225.07.29.60.59511,56812.581.85105.1180.8037811411
102.925.08.311.10.57112,52113.081.83107.7178.2038471495
103.424.98.311.20.57412,55613.081.83108.1178.5037981497
102.925.08.211.10.57212,40712.081.83110.9177.9037741494

References

  1. York, R.; Bell, S.E. Energy transitions or additions?: Why a transition from fossil fuels requires more than the growth of renewable energy. Energy Res. Soc. Sci. 2019, 51, 40–43. [Google Scholar] [CrossRef]
  2. Panichkin, I. Arctic Oil and Gas Resource Development: Current Situation and Prospects; Russia International Affairs Council: Moscow, Russia, 2016; pp. 3–10. [Google Scholar]
  3. Watson, D.G. Practical Ship Design; Elsevier: Oxford, UK, 1998; pp. 495–497. [Google Scholar] [CrossRef]
  4. IMO. Energy Efficiency Measures. Available online: https://www.imo.org/en/OurWork/Environment/Pages/Technical-and-Operational-Measures.aspx. (accessed on 18 April 2021).
  5. IMO. Resolution MEPC.308(73) 2018 Guidelines on the Method of Calculation of the Attained Energy; International Maritime Organization: London, UK, 2018; pp. 1–36. [Google Scholar]
  6. Trivyza, N.L.; Rentizelas, A.; Theotokatos, G. A Comparative Analysis of EEDI Versus Lifetime CO2 Emissions. J. Mar. Sci. Eng. 2020, 8, 61. [Google Scholar] [CrossRef] [Green Version]
  7. Tarovik, O.V.; Topaj, A.G.; Krestyantsev, A.B.; Kondratenko, A.A.; Zaikin, D.A. Study on operation of Arctic offshore complex by means of multicomponent process-based simulation. JMSA 2018, 17, 471–497. [Google Scholar] [CrossRef]
  8. Papanikolaou, A. Holistic ship design optimization. Comput. Aided Des. 2010, 42, 1028–1044. [Google Scholar] [CrossRef]
  9. Gutsch, M.; Steen, S.; Sprenger, F. Operability robustness index as seakeeping performance criterion for offshore vessels. Ocean Eng. 2020, 217, 107931. [Google Scholar] [CrossRef]
  10. Marques, C.H.; Belchior, C.R.; Caprace, J. An early-stage approach to optimise a marine energy system for liquefied natural gas carriers: Part A—Developed approach. Ocean Eng. 2019, 181, 161–172. [Google Scholar] [CrossRef]
  11. Priftis, A.; Evangelos, B.; Osman, T.; Papanikolaou, A. Parametric design and multi-objective optimisation of containerships. Ocean Eng. 2018, 156, 347–357. [Google Scholar] [CrossRef] [Green Version]
  12. Kanellopoulou, A.; Kytariolou, A.; Papanikolaou, A.; Shigunov, V.; Zaraphonitis, G. Parametric ship design and optimisation of cargo vessels for efficiency and safe operation in adverse weather conditions. J. Mar. Sci. Technol. 2019, 24, 1223–1240. [Google Scholar] [CrossRef]
  13. Skoupas, S.; Zaraphonitis, G.; Papanikolaou, A. Parametric design and optimisation of high-speed Ro-Ro Passenger ships. Ocean Eng. 2019, 189, 10634. [Google Scholar] [CrossRef]
  14. Peri, D. Direct Tracking of the Pareto Front of a Multi-Objective Optimization Problem. J. Mar. Sci. Eng. 2020, 8, 699. [Google Scholar] [CrossRef]
  15. Martins, M.R.; Burgos, D.F. Multi-Objective Optimization Design of Tanker Ships via a Genetic Algorithm. J. Offshore Mech. Arct. Eng. 2011, 133, 041303. [Google Scholar] [CrossRef]
  16. Ray, T.; Gokarn, R.; Sha, O. A global optimization model for ship design. Comput Ind. 1995, 26, 175–192. [Google Scholar] [CrossRef]
  17. Kujala, P.; Goerlandt, F.; Way, B.; Smith, D.; Yang, M.; Khan, F.; Veitch, B. Review of risk-based design for ice-class ships. Mar. Struct. 2019, 63, 181–195. [Google Scholar] [CrossRef]
  18. Bergström, M.; Erikstad, S.; Ehlers, S. The Influence of model fidelity and uncertainties in the conceptual design of Arctic maritime transport systems. Ship Technol. Res. 2017, 64, 40–64. [Google Scholar] [CrossRef]
  19. Topaj, A.G.; Tarovik, O.V.; Bakharev, A.A.; Kondratenko, A.A. Optimal ice routing of a ship with icebreaker assistance. Appl. Ocean Res. 2019, 86, 177–187. [Google Scholar] [CrossRef]
  20. Dobrodeev, A.; Sazonov, K. The estimation of carbonic gas emission by ice-class large-size ships moving in ice using different escorting methods. In Proceedings of the ASME 2016 35th International Conference on Ocean, Offshore and Arctic Engineering OMAE 2016, Busan, Korea, 18–24 June 2016. [Google Scholar] [CrossRef]
  21. Kondratenko, A.A.; Tarovik, O.V. Analysis of the impact of arctic-related factors on offshore support vessels design and fleet composition performance. Ocean Eng. 2020, 203, 107201. [Google Scholar] [CrossRef]
  22. Keßler, T.; Logist, F.; Mangold, M. Use of predictor corrector methods for multi-objective optimization of dynamic systems. In Proceedings of the 26th European Symposium on Computer Aided Process Engineering—ESCAPE 26, Portorož, Slovenia, 12–15 June 2016. [Google Scholar] [CrossRef]
  23. Karaboga, D.; Basturk, B. A powerful and efficient algorithm for numerical function optimization: Artificial bee colony (ABC) algorithm. J. Glob. Optim. 2007, 39, 459–471. [Google Scholar] [CrossRef]
  24. Kondratenko, A.A.; Tarovik, O.V. Cargo-Flow-Oriented Design of supply vessel operating in ice conditions. In Proceedings of the ASME 2018 37th International Conference on Ocean, Offshore and Arctic Engineering OMAE2018, Madrid, Spain, 17–22 June 2018. [Google Scholar] [CrossRef]
  25. RS. International Association of Classification Societies. Symbols of the Classification Of Ships. Directory No 2-029901-002; The Russia Maritime Register of Shipping: St. Petersburg, Russia, 2015; pp. 28–35. [Google Scholar]
  26. IMO. Revised Guidance to the Master for Avoiding Dangerous Situations in Adverse Weather and Sea Conditions; International Maritime Organization: London, UK, 2007; pp. 1–8. [Google Scholar]
  27. The Federal Tariff Service. Decree No 45-T/1 on the Approval of Tariffs for Icebreaker Assistance of Ships, Provided by FSUE Atomflot on the Northern Sea Route; FST: Moscow, Russia, 2014; pp. 1–5. [Google Scholar]
  28. Brovin, A.I.; Klyachkin, S.V. Application of an empirical-statistical model of ship motion in ice to new types of icebreakers and ships. In Proceedings of the ASME 1997 16th International Conference on Ocean, Offshore and Arctic Engineering OMAE 1997, Yokohama, Japan, 13–17 April 1997. [Google Scholar]
  29. CNIIMF. Study of Special Features of Winter Navigations in GoF and Development Algorithms and Models for Calculations the Ship Movement in Ice with the Aim to Use Them in Risk Model When Evaluating the RCO Efficiency within the Frame of the WINOIL Project Works; Central Marine Research and Design Institute: St. Petersburg, Russia, 2014. [Google Scholar]
  30. Rao, S.K.; Chauhan, P.J.; Panda, S.K.; Wilson, G.; Liu, X.; Gupta, A.K. Optimal scheduling of diesel generators in offshore support vessels to minimize fuel consumption. In Proceedings of the IECON 2015—41st Annual Conference of the IEEE Industrial Electronics Society, Yokohama, Japan, 9–19 November 2015. [Google Scholar] [CrossRef]
  31. Buzuev, A.Y.; Dubovtsev, V.F.; Zakharov, V.F.; Smirnov, V.I. Conditions of Navigation in Sea Ice of the Northern Hemisphere; GUNIO USSR: Leningrad, USSR, 1988; pp. 154–157. (In Russian) [Google Scholar]
  32. Arctic and Antarctic Research Institute. Center of Ice & Hydrometeorological Information. Available online: http://www.aari.ru/main.php?lg=0&id=17 (accessed on 22 September 2020).
  33. Shalina, E.V.; Sandven, S. Snow depth on Arctic sea ice from historical in situ data. Cryosphere 2018, 12, 1867–1886. [Google Scholar] [CrossRef] [Green Version]
  34. Dumanskaya, I.O. Sea Ice Condition of European Russia; IG-SOTSIN: Obninsk, Russia, 2014; pp. 11–213. (In Russian) [Google Scholar]
  35. Kondratenko, A.; Tarovik, O. OSVs Characteristics Database. xlsx. Available online: https://www.researchgate.net/publication/340461340_OSVs_characteristics_databasexlsx_from_the_article_Analysis_of_the_impact_of_arctic-related_factors_on_offshore_support_vessels_design_and_fleet_composition_performance (accessed on 15 May 2021). [CrossRef]
Figure 1. A general flow chart of the applied ABC algorithm. A detailed description of the algorithm is presented in [23].
Figure 1. A general flow chart of the applied ABC algorithm. A detailed description of the algorithm is presented in [23].
Sustainability 13 05550 g001
Figure 2. A general flow chart of the optimization process.
Figure 2. A general flow chart of the optimization process.
Sustainability 13 05550 g002
Figure 3. A general flow chart of the vessel evaluation model.
Figure 3. A general flow chart of the vessel evaluation model.
Sustainability 13 05550 g003
Figure 4. Module decomposition of the vessel evaluation model. The modules that are the most important for this study are highlighted by a bold frame. The background color reflects the module category: gray—technical, red—safety, green—environment and emissions, blue—economics.
Figure 4. Module decomposition of the vessel evaluation model. The modules that are the most important for this study are highlighted by a bold frame. The background color reflects the module category: gray—technical, red—safety, green—environment and emissions, blue—economics.
Sustainability 13 05550 g004
Figure 5. The power control algorithm for the reduction of a vessel’s fuel consumption in a voyage.
Figure 5. The power control algorithm for the reduction of a vessel’s fuel consumption in a voyage.
Sustainability 13 05550 g005
Figure 6. (a) The route between the port Murmansk to the Pobeda field used in the case study; (b) an ice chart determined by AARI covering the Barents Sea [32].
Figure 6. (a) The route between the port Murmansk to the Pobeda field used in the case study; (b) an ice chart determined by AARI covering the Barents Sea [32].
Sustainability 13 05550 g006
Figure 7. Frequency distribution of different ranges of equivalent ice thickness along route segments 4 and 5 in March.
Figure 7. Frequency distribution of different ranges of equivalent ice thickness along route segments 4 and 5 in March.
Sustainability 13 05550 g007
Figure 8. Optimization results for case study 1 (a) and 2 (b). The black points highlight the Pareto front and blue points represent other feasible solutions.
Figure 8. Optimization results for case study 1 (a) and 2 (b). The black points highlight the Pareto front and blue points represent other feasible solutions.
Sustainability 13 05550 g008
Figure 9. Dynamics of the value of IButil during the year for vessel Min CKPI and Min EKPI (case study 1).
Figure 9. Dynamics of the value of IButil during the year for vessel Min CKPI and Min EKPI (case study 1).
Sustainability 13 05550 g009
Figure 10. Optimization results for case study 3 (a) and 4 (b). The Pareto front is shown as black points and blue points show other feasible solutions. The yellow points present the Pareto front vessels’ performance if vessels operate at maximum speed.
Figure 10. Optimization results for case study 3 (a) and 4 (b). The Pareto front is shown as black points and blue points show other feasible solutions. The yellow points present the Pareto front vessels’ performance if vessels operate at maximum speed.
Sustainability 13 05550 g010
Figure 11. Optimization results for case study 5 (a) and 6 (b). The Pareto front is shown as black points and blue points show other feasible solutions. The yellow points present the vessels’ performance from the Pareto front if vessels operate at maximum speed.
Figure 11. Optimization results for case study 5 (a) and 6 (b). The Pareto front is shown as black points and blue points show other feasible solutions. The yellow points present the vessels’ performance from the Pareto front if vessels operate at maximum speed.
Sustainability 13 05550 g011
Figure 12. Visualizations of the optimal hulls developed by the model for the intermediate Case 4 (a) and Case 6 (b).
Figure 12. Visualizations of the optimal hulls developed by the model for the intermediate Case 4 (a) and Case 6 (b).
Sustainability 13 05550 g012
Figure 13. Optimization results for Case 7 without icebreaking assistance. The black points represent the Pareto front and the blue points represent other feasible solutions.
Figure 13. Optimization results for Case 7 without icebreaking assistance. The black points represent the Pareto front and the blue points represent other feasible solutions.
Sustainability 13 05550 g013
Table 1. Case study parameters.
Table 1. Case study parameters.
Case StudyDiscount Rate (%)Cargo Capacity
Parameter
Operation Period (Years)Icebreaker Assistance (Yes/No)
16Dw10+
26Dw20+
312Dw10+
412Dw20+
56Scargo10+
612Scargo20+
712Scargo20-
Table 2. Applied design parameter ranges.
Table 2. Applied design parameter ranges.
ParameterMin.Max.
Length between perpendiculars, m50100
Beam, m1325
Draft, m3.59
Block coefficient0.570.78
Depth, m513.5
Icebreaking capability, m02.8
Table 3. Characteristics of selected vessels from the Pareto front, case study 1–2.
Table 3. Characteristics of selected vessels from the Pareto front, case study 1–2.
Case Study 1
CaptionLwl, mB, mT, m∆, tCbwlNpp, kWVs, knArchice, mIButilDw, tSdeck, m2
Min CKPI110.818.46.890410.635602513.550.830.0344756851
Intermediate110.122.19.0142280.634820912.050.860.04175591107
Min EKPI110.825.09.0161890.635955111.050.770.07186531360
Case Study 2
CaptionLwl, mB, mT, m∆, tCbwlNpp, kWVs, knArchice, mIButilDw, tSdeck, m2
Min CKPI110.218.56.991520.635607713.550.830.0344820851
Intermediate110.321.78.9139150.635806412.550.850.04074281088
Min EKPI110.925.09.0161840.635955111.050.770.07186381361
Table 4. Characteristics of selected vessels from the Pareto front, cases 3–4.
Table 4. Characteristics of selected vessels from the Pareto front, cases 3–4.
Case 3
CaptionLwl, mB, mT, m∆, tCbwlNpp, kWVs, knArchice, mIButilDw, tSdeck, m2
Min CKPI110.718.56.586600.635597013.550.840.0344534853
Intermediate110.618.58.411,2090.635651212.050.770.0605956850
Min EKPI110.925.09.016,2400.635956311.050.770.07186791361
Case 4
CaptionLwl, mB, mT, m∆, tCbwlNpp, kWVs, knArchice, mIButilDw, tSdeck, m2
Min CKPI110.518.47.092260.635607013.450.840.0344863849
Intermediate110.318.48.411,1280.635648312.050.770.0605912846
Min EKPI110.825.09.016,2160.635954711.050.770.07186651359
Table 5. Characteristics of selected vessels from the Pareto front, Cases 5–6.
Table 5. Characteristics of selected vessels from the Pareto front, Cases 5–6.
Case 5
CaptionLwl, mB, mT, m∆, tCbwlNpp, kWVs, knArchice, mIButilDw, tSdeck, m2
Min CKPI108.018.04.556460.622530613.550.800.0302403814
Intermediate109.620.95.479460.624649712.550.750.05339541045
Min EKPI105.823.56.499280.610397386.582.440.10229551232
Case 6
CaptionLwl, mB, mT, m∆, tCbwlNpp, kWVs, knArchice, mIButilDw, tSdeck, m2
Min CKPI109.618.34.659330.625652314.550.920.0272670841
Intermediate107.622.85.990170.610722711.950.770.05440811214
Min EKPI105.325.07.8123460.591419176.282.470.10836781413
Table 6. Characteristics of the selected vessels from the Pareto front for Case 7.
Table 6. Characteristics of the selected vessels from the Pareto front for Case 7.
Case 7
CaptionLwl, mB, mT, m∆, tCbwlNpp, kWVs, knArchice, mIButilDw, tSdeck, m2
Min CKPI108.820.05.068670.61721,08215.081.8401801991
Intermediate106.424.76.710,8330.59723,86913.081.87031871408
Min EKPI102.925.08.212,4070.57223,13912.081.83037741494
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kondratenko, A.A.; Bergström, M.; Reutskii, A.; Kujala, P. A Holistic Multi-Objective Design Optimization Approach for Arctic Offshore Supply Vessels. Sustainability 2021, 13, 5550. https://0-doi-org.brum.beds.ac.uk/10.3390/su13105550

AMA Style

Kondratenko AA, Bergström M, Reutskii A, Kujala P. A Holistic Multi-Objective Design Optimization Approach for Arctic Offshore Supply Vessels. Sustainability. 2021; 13(10):5550. https://0-doi-org.brum.beds.ac.uk/10.3390/su13105550

Chicago/Turabian Style

Kondratenko, Aleksander A., Martin Bergström, Aleksander Reutskii, and Pentti Kujala. 2021. "A Holistic Multi-Objective Design Optimization Approach for Arctic Offshore Supply Vessels" Sustainability 13, no. 10: 5550. https://0-doi-org.brum.beds.ac.uk/10.3390/su13105550

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop