Next Article in Journal
Development of a MELCOR Model for LVR-15 Severe Accidents Assessment
Next Article in Special Issue
Analysis of Non-Road Mobile Machinery Homologation Standards in Relation to Actual Exhaust Emissions
Previous Article in Journal
Investigation of Split Diesel Injections in Methanol/Diesel Dual-Fuel Combustion in an Optical Engine
Previous Article in Special Issue
Visualisation Testing of the Vertex Angle of the Spray Formed by Injected Diesel–Ethanol Fuel Blends
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reducing CO2 Emissions through the Strategic Optimization of a Bulk Carrier Fleet for Loading and Transporting Polymetallic Nodules from the Clarion-Clipperton Zone

Faculty of Navigation, Maritime University of Szczecin, 1-2 Wały Chrobrego St., 70-500 Szczecin, Poland
*
Author to whom correspondence should be addressed.
Energies 2024, 17(14), 3383; https://doi.org/10.3390/en17143383
Submission received: 13 June 2024 / Revised: 5 July 2024 / Accepted: 8 July 2024 / Published: 10 July 2024
(This article belongs to the Special Issue CO2 Emissions from Vehicles (Volume II))

Abstract

:
As global maritime cargo transportation intensifies, managing CO2 emissions from ships becomes increasingly crucial. This article explores optimizing bulk carrier fleets for transporting polymetallic nodules (PMNs) from the Clarion-Clipperton Zone (CCZ) to reduce CO2 emissions. Our analysis shows that larger bulk carriers, despite greater drifting forces from environmental conditions, emit less CO2 over the entire transport mission, including loading and transit. Deploying large ships in global maritime trade could significantly reduce CO2 emissions. This study also introduces a novel artificial neural network (ANN) model to estimate drifting forces during loading operations and proposes a new method for estimating CO2 emissions, considering environmental conditions and ship seakeeping properties. These findings highlight the importance of fleet size optimization and effective operational planning in achieving environmental sustainability in maritime transport.

Graphical Abstract

1. Introduction

The 2021 United Nations Conference on Trade and Development report [1] highlights maritime transport as the cornerstone of international trade and the global economic framework. The report indicates that over 80% of global trade volume, particularly in goods, occurs via sea routes. This significant reliance on maritime transportation, in the context of expanding international trade and coupled with numerous pollution control measures targeting other sources, has led to an increasingly substantial contribution of ship emissions to total anthropogenic emissions [2]. Nunes et al. [3] observed that the maritime sector is likely to continue its significant growth alongside global trade. However, the full extent of its impact on the environment, societal aspects, and human health remains unclear. Moreno-Gutiérrez et al. [4] identified the combustion of fuel in ship engines as the primary source of ship-emitted pollutants. Hoang et al. [5] argue that the pursuit of intelligent strategies through the utilization of renewable energy sources, clean fuels, smart grids, and measures for efficient energy use is beneficial for achieving the main IMO objectives, particularly reductions in future CO2 emissions. Fan et al. [6], in their analysis of ship energy management, justified that managing energy savings and emission reductions is a systemic issue that should be comprehensively considered from multiple perspectives, such as ship design, operational management, and performance assessment. The technical modernization of the ship can also contribute to reducing CO2 emissions. Barone et al. [7] conducted simulations on a ship equipped with five electric chillers, two auxiliary boilers, two reverse osmosis devices, and two multi-stage flash desalination devices. The simulation results show remarkable energy savings obtained through the proposed optimization approach, corresponding to a fuel consumption reduction of 2.5 kt/y (−1.6 M$/y) and to avoided CO2 emissions of 8.0 kt/y. The application of solar and wind energy devices affects the reduction in CO2 emissions. Nyanya et al. [8], in their analysis using a bulk carrier as an example, demonstrated that considering the optimal sail angle and optimizing the available deck area through the installation of solar and wind systems allowed for the maximum utilization of renewable energy, which contributed to a 36% reduction in carbon dioxide emissions when compared to the same ship without innovative technologies. The conclusions showed that if the ship’s speed were reduced to 56% of its original speed, the ship could sail exclusively on renewable energy collected on board. Ytreberg et al. [9] reported that maritime shipping contributes to several environmental challenges, including deteriorations in air quality and human well-being, leading to numerous quantifiable and unquantifiable harms. Their findings demonstrate that transporting goods solely in the Baltic Sea incurs costs due to air quality deterioration, which exceed those associated with climate change impacts.
It is worth noting that oceans and seas serve not only as crucial transport routes for global trade, but also as valuable sources of natural raw materials upon which many industries depend. Recent studies suggest that the deep-sea floor at depths of 4000–6000 m represents a significant source of polymetallic nodules (PMNs), comprising minerals such as manganese (Mn), copper (Cu), cobalt (Co), nickel (Ni), and some rare earth elements [10]. These nodules are considered a promising source of metallic raw materials, with diverse mining systems currently being developed for their extraction [11]. Cunningham [12] concluded that the extraction of nodules from the ocean floor could provide a more stable and decentralized alternative to the current cobalt supplies. The increased cobalt supply from nodule extraction may offset the harmful impacts on the marine environment. Studies by Gollner et al. and Bonifácio et al. [13,14] noted that the extraction process could have adverse environmental effects, impacting biodiversity in the area. Nodule extraction disrupts the composition of seabed sediments, which may lead to long-term negative consequences considering the slow pace of natural habitat recovery. The CC zone features 300 morphotypes, including numerous groups of invertebrates such as corals, crustaceans, and fish [15]. Stratmann [16] argues that conducting extraction operations could lead to an 18% loss of organisms (taxa). B. Li et al. [17] suggests the application of ore separation technology on extraction ships and transporting the separated sediments to land for comprehensive utilization on a separate ship. However, Katona et al. [18] argues that the long-term environmental impact of nodule extraction on marine ecosystems might be lesser when compared to terrestrial mining.
According to recent findings [19], it is anticipated that nodules extracted from the ocean floor will be directly loaded onto high-capacity extraction vessels and transferred to transport ships that will deliver the cargo to destination ports.
During the transfer of nodules from the extraction unit to the transport ship on the ocean surface, the marine environment, i.e., wind and waves, interacts with the transport ship. This can lead to excessive ship swaying, bending, or a loss of stability, as indicated by Kacprzak [20]. In such events, it becomes necessary to suspend the transfer, potentially disrupting the extraction of the nodules from the ocean floor and the operation of the extraction system. Therefore, the transport ship should have the appropriate maritime characteristics to ensure high-efficiency loading and the ability to operate in challenging weather conditions for as long as possible over the mining field, thus minimizing the risk of downtime. An inappropriate selection of transport ships might disrupt the entire extraction system.
According to Kacprzak [21], the transport of PMNs to the destination port can be achieved by both small and large bulk carriers. Research by this author shows that loading onto smaller bulk carriers would lead to shorter loading times and good longitudinal strength, which are significant advantages. However, the disadvantage of these units is their smaller transport capacity, increased ship movement, and poorer stability. On the other hand, large bulk carriers have the advantages of high transport capacity, minimal ship movement, and good stability, with the downside being longer loading times and reduced longitudinal strength.
In addition to the technical limitations with regard to safe loading, a key aspect is the environmental pollution caused by the transporting ship during loading and transport operations. The total CO2 emissions into the marine environment by the transporting ship depends on the amount of fuel being burned [22]. This amount mainly depends on the power of the dynamic positioning (DP) system which stabilizes the bulk carrier at the loading site, as well as environmental conditions. The required power of the DP system depends on many factors, including the technical parameters of the ship and environmental conditions. Among the technical properties, the size of the ship is most important; according to Jurdziński [23], larger bulk carriers are characterized by a larger windage area, exposed to wind action, and a larger wetted surface, exposed to wave action, which may result in greater drift force. An increase in drift force leads to an increase in the energy needed to maintain the ship position, and thus an increase in pollutant emissions.
In turn, the amount of fuel burned during transport depends on the technical–operational parameters and environmental conditions. The most important technical–operational parameters include ship speed, capacity, the hydro-mechanical properties of the hull, engine efficiency, propulsion system and propeller, type and quality of fuel, draught, trim and hull, and the propeller condition from the fouling process, as mentioned in [24,25,26]. Environmental conditions affecting fuel consumption include waves, wind, ocean currents, air, and water temperature, as noted in [27]. The reloading and transport of PMNs can be carried out by a fleet of small, medium, and large bulk carriers. Research by Kacprzak [21] has shown that transporting a specified amount of PMNs from the Clarion-Clipperton Zone (CCZ) to the destination requires the consideration of the environmental impact on the technical properties of ships, including stability, motions, and longitudinal strength. Therefore, the bulk carrier fleet should consist of appropriately selected units of various sizes. The key concern is what impact the use of the entire fleet to transport the specified amount of PMNs to different destinations, on routes of varying distances, will have on emissions to the marine environment. Will the use of a fleet of small, medium, or large bulk carriers ensure lower exhaust emissions? How does the distance to the destination port from the CCZ affect this choice? Therefore, the main goal of this study is to answer these questions, i.e., to determine which fleet of bulk carriers (small, medium, large) will be most efficient in terms of environmental CO2 emissions during the loading and transport of PMNs from the CCZ to various ports around the world.
Studies [28,29] have shown that a large ship can transport a larger amount of cargo at once, burning relatively less fuel during transport to the destination, which translates into less environmental pollution. Therefore, to assess the energy efficiency of the ship, many indicators have been developed, such as
  • the energy efficiency operational indicator (EEOI), which is calculated as the ratio of CO2 emissions to the amount of cargo carried and the distance traveled [30],
  • the energy efficiency design index (EEDI), which is expressed as the amount of carbon dioxide CO2 emitted into the environment based on the estimated fuel consumption in relation to ship capacity and speed [31],
  • the annual efficiency ratio (AER), which was created to measure the energy efficiency of a ship based on its annual activity, taking into account the total fuel consumption and distance traveled [32],
  • the deadweight tonnage (DWT) efficiency, which determines the efficiency of a ship by the ratio of its cargo capacity (carrying capacity) to fuel consumption,
  • the carbon intensity indicator (CII), which measures CO2 emissions per unit of transport and distance [33].
Chen et al. [34] argue that the ship’s EEOI calculated at full load is naturally lower than the equivalent at partial load. The energy efficiency of the fleet shows a slight increase (at least by 1%) due to the implementation of the Ship Energy Efficiency Management Plan. Lee [35] evaluated the possibility of achieving the 40% CO2 emission reduction target by 2030 through the implementation of the EEDI and the EEXI for three types of ships as follows: bulk carriers, container ships, and tankers. The results show that achieving the 40% CO2 emission reduction target by 2030 is possible for ships regulated by the EEDI, while the calculated CO2 emissions for ships regulated by the EEXI range from 17.4% to 24.6%. This difference arises from the fact that these calculations are based on the maximum speed limit of the EPL’s engine output, which exceeds the actual operational speed.
These indicators can be used to assess CO2 emissions during the transport of nodules. However, they do not consider the energy supplied to the system stabilizing ship movement during reloading operations. During this operation, the transport ship is stationary or moves at a low speed along with the extraction ship, and the DP system stabilizes any ship movement caused by the drift force resulting from wave and wind action. There is no appropriate indicator which assesses the energy efficiency of this process. Therefore, the second goal of the research is to develop a method for estimating CO2 emissions that considers the impact of environmental conditions on ship loading operations at sea with a cargo of nodules.

2. Materials and Methods

2.1. Data

The subject of the study was a fleet of bulk carriers comprising small, medium, and large ships. Representatives of the three groups of ships were selected for the study. The parameters of these ships are presented in Table 1.
To compare the efficiency of the different ship types, assumptions regarding the mining system and the locations of the destination ports to which the mined PMNs could be transported were made. It was assumed that the loading of PMNs would take place in the CCZ and would be carried out by the mining ship Hidden Gem. The technical characteristics of this ship are presented in Table 2.
Assumptions regarding the mining system were adopted in accordance with [36] and are presented below:
  • the number of planned days for transshipment operations—250 days,
  • the assumed annual transport of PMNs—2,000,000 tons,
  • the mining capacity of the mining ship—333 tons/h,
  • the cargo carrying capacity of the mining ship Hidden Gem—20,000–25,000 t.
It was assumed that the transportation of PMNs from the CCZ would be carried out to the destination ports along the routes described in Figure 1.
The distribution of wave parameters in the CCZ was developed based on data from [37] and is presented in Table 3. These distributions relate to the occurrence frequency of these parameters throughout the year.
The distribution of wind parameters in the CCZ was adopted in accordance with data from [38] and is presented in Table 4.
To estimate the drift force (Fwy), results from ship model tests, the dimensions of which are presented in Table 5, were utilized. These results are based on the studies [39,40,41,42,43,44] and include the nondimensional drift force coefficient (Cwy), calculated using the following equation:
C w y = F w y ρ · g · ς a 2 · L
where
Cwy—nondimensional drift force coefficient,
Fwy—mean drift force in a regular wave,
ρ—seawater density,
g —gravity acceleration,
ς a —amplitude of a regular wave,
L—length of the ship.
All data from these model tests are publicly available, and the references for each dataset are provided in Table 5. Figure 2 shows the values of the drift force coefficient (Cwy) in regular waves for these ship models at zero speed as a function of the ratio of the wavelength (λ) to the length of the ship (L). Based on the nondimensional coefficient Cwy, calculated using Equation (1) for the ship model scale, it is possible to determine the drift force Fwy for a full-scale ship of length L.

2.2. General Research Framework

To achieve the research objective, the estimation of CO2 emissions during the loading and transportation of PMNs to various long-range destinations was identified as a key step. This estimation was conducted for each type of bulk carrier fleet, designated as A, B, and C, as described in Table 1, depending on the location of the destination port to which the PMNs were to be delivered.
Accordingly, the total CO2 emission is made up of emissions during loading and transportation as follows:
CO2 = CO2load + CO2transp,
where
CO2—total CO2 emissions by bulk carrier fleets A/B/C,
CO2load—CO2 emissions by bulk carrier fleets A/B/C during loading,
CO2transp—CO2 emissions during transportation to the destination port using the bulk carrier fleets A/B/C.
In this study, it was assumed that loading from the mining ship could be executed concurrently with the extraction of PMNs and continuously without breaks for nodule reloading.

2.2.1. Estimation of CO2 Emissions during Polymetallic Nodules Loading in the Clarion-Clipperton Zone

Figure 3 illustrates a general block diagram for the estimation of CO2 emissions while nodule loading in the CCZ. The calculation of CO2 emissions during the loading process is based on the fuel consumption of the transporting bulk carrier while it loads PMNs from the extraction ship as follows:
CO2 = EF ∙ FC [t]
where
FC—fuel consumption,
EF—emission factor, whose value for a marine engine can be assumed as approximately 3.206 kg CO2/kg for marine diesel oil (MDO) or approximately 3.114 kg CO2/kg for heavy fuel oil (HFO).
The fuel consumption during this operation stems from the necessity to generate thrust via the propeller of the DP system to balance out environmental forces, thereby maintaining the ship position during loading under the specified environmental conditions. This fuel consumption (FC) was calculated based on the general relation as follows:
FC = P · t · FOC
where
P—power supplied to the DP system [kW],
t—time under specified environmental conditions [s],
FOC—fuel consumption rate.
Equation (4) finds extensive application in the analysis and estimation of fuel consumption, both under operational conditions and in the design of ship propulsion systems. A crucial component of this Equation is the fuel consumption rate (FOC), whose values can be derived from the nominal fuel consumption of the marine engine or from real-time fuel flow measurements under operational conditions. The typical FOC values range from 0.170 to 0.185 kg/kWh for marine diesel oil (MDO) and marine gas oil (MGO), and from 0.195 to 0.230 kg/kWh for heavy fuel oil (HFO), depending on the engine type and the operational conditions.
The required power supplied to the propellers of the DP system can be calculated based on the following formula [45]:
P = Tp/Cp
where
P—propeller power [kW],
Tp—thrust force of the DP system propeller required to overcome environmental forces [kN],
Cp—power-to-thrust efficiency ratio for the given propeller; for tunnel thrusters,
Cp = 0.12–0.16.
The calculation of propeller thrust force generally considers the wind drift force, wave drift force, and current drift force. Given that the CCZ, located in the Pacific between Hawaii and Mexico, is known for its slow-moving ocean currents with a maximum speed of approximately 0.3 m/s [37] and given the negligible speed of the ship, the influence of ocean currents on the ship was omitted in this study. Therefore, the thrust force (Tp) was calculated as follows:
Tp = Fw + Fwy
where
Fw—force from wind pressure,
Fwy—force from wave drift.
The magnitude of forces from wind and waves depends on their direction relative to the ship. Given the lack of detailed information about wave and wind directions in the CCZ, this study assumes that waves and wind impact the ship unfavorably, i.e., perpendicular to the beam. The force from wind pressure is calculated as follows:
Fw = ½ · vW2 · ρ · Aw
where
vW—wind speed according to Table 4,
ρ—air density,
Aw—the lateral plane of the vessel above the plane of draught according to Table 1.
The wave drift force was calculated based on the following relation:
F w y = L · ρ · g · S ς ς ω · C w y ω d ω
where
Fwy—mean wave drift force from irregular waves [kN],
L—ship length [m],
ρ—seawater density,
g—acceleration due to gravity,
Sζζ—wave spectral function,
Cwy—mean nondimensional sway force in regular wave.
The wave spectrum function was implemented according to the ITTC guidelines [46] as follows:
S ζ ζ ( ω ) = A ω 5 exp ( B ω 4 )
where
A, B—coefficients calculated as follows:
A = 173   H s 2 T z 4
B = 691 T z 4
HS—significant wave height,
Tz—zero-up-crossing wave period,
ω—circular frequency.
To calculate the drift force, numerical methods, such as the far-field and near-field pressure methods, as well as experimental studies, are used. Far-field methods, initiated by Maruo [47], further developed by Newman, Salvesen, Gerritsma, and Beukelman [48,49,50] and most recently by Kashiwagi [51], are based on considerations of the energy of diffracted (reflected and transmitted) and radiated waves, momentum flux at infinity, and work conducted by the body in the near field, leading to a constant force acting through the total rate of momentum change. Alternatively, near-field methods, initiated by Boese [52], calculate second-order constant forces/moments through the direct integration of hydrodynamic, steady second-order pressure acting on the wetted surface of the body. More accurate numerical implementations were later introduced by Faltinsen et al. [53], and then by Papanikolaou and Nowacki [54], which led to an increase in precision. These methods have evolved and been validated based on experimental studies, mainly in terms of the longitudinal component, referred to as additional wave resistance. However, the validation of these methods for a transverse component is limited, mainly due to the scarcity of experimental research. Recently, Liu S and Papanikolaou [55] approximated the mean sway force with an empirical formula, only utilizing main ship particulars and wave parameters.

2.2.2. Application of Artificial Neural Networks to Estimate Mean Nondimensional Sway Force in Regular Waves

The application of artificial neural networks (ANNs) to estimate the nondimensional drift force coefficient in regular waves (Cwy) for bulk carriers A, B, and C was assessed using the results of model studies described in Table 5, which are presented in Figure 2. An ANN was used for the estimation process. An additional objective of this study was to examine the effectiveness of ANN use to predict drift force based on the basic dimensions of the ship.
In the field of maritime engineering, the application of ANNs has been increasingly acknowledged for its critical role in enhancing the predictability of the experimental measurements pertinent to ship theory. The development and implementation of ANNs have enabled scientists and engineers to process and interpret extensive datasets efficiently and accurately, significantly improving the accuracy of predictive models in ship design and operational strategies. The integration of ANNs in naval architecture has notably enhanced the accuracy of prognostic models, facilitating substantial advancements in various aspects of ship construction and functionality.
For example, research by Cepowski [56,57] utilized Multilayer Perceptron (MLP) structured neural networks to estimate the additional wave resistance of ships based on fundamental design parameters. This approach resulted in remarkably accurate estimations when compared to outcomes from model experiments. Similarly, Yangjun and Yonghwan [58] applied an ANN to predict sloshing effects in model tests of floating structures larger than conventional LNG carriers. The development of this ANN was based on a thorough analysis of over 540 terabytes of experimental data, covering a wide range of variables, such as cargo hold dimensions, vessel types, and environmental and operational conditions. Furthermore, Dyer et al. [59] employed an ANN to estimate the remaining service life of offshore oil and gas platforms, incorporating a comprehensive dataset representing both natural and engineered aspects of offshore systems. Their findings highlighted the capacity of ANNs to provide highly accurate predictions, indicating life extension, maintenance, and risk minimization strategies for offshore platforms.
In this study, the development of an ANN to estimate the Cwy coefficient utilized ship data from Table 5 and model study results presented in Figure 2. Following [56,57], a Multilayer Perceptron (MLP)-structured ANN was adopted.
Figure 4 illustrates the procedure for developing the ANN. The first step involved normalizing the dataset to the [0, 1] range to accelerate learning, avoid numerical issues such as overflow or underflow, prevent feature domination, facilitate weight initialization, and enhance the performance of algorithms based on gradients or distances. Following this procedure, a set of alternative ANNs was developed, differing in the number of neurons in the hidden layer and activation functions. In each ANN, to detect overfitting, the hold-out method with an additional validation set was applied, consisting of the following steps:
  • randomly splitting the data into three parts: training set (50%), testing set (25%), and validation set (25%),
  • training the model on the training set,
  • applying early stopping to prevent overfitting by monitoring the error on the testing set,
  • testing the model on the testing set after each epoch and calculating performance metrics,
  • validating the model on the validation set after training is complete,
  • calculating final performance metrics and comparing them with the results from the testing set to assess the model’s generalization capability,
  • analyzing the results from the validation set to evaluate the model’s ability to generalize to new data.
In the hold-out procedure, the training and testing process involved early stopping to prevent overfitting by halting the training process if errors increased. The early stopping technique consisted of the following steps:
  • input data: feeding the network with training dataset inputs,
  • prediction: calculating the network’s predictions (outputs),
  • error calculation: computing the difference between the predictions and the actual output values using an error function,
  • repetition: repeating steps 1 and 2 for all input–output pairs in the training dataset,
  • weight adjustment: using a learning algorithm to adjust the weights of the neurons to minimize the prediction error,
  • testing: feeding the network with all cases from the testing dataset, obtaining predictions, and comparing them with the actual output values from the dataset to calculate the network error,
  • error comparison: comparing the current network error with the error from the previous cycle. If the error has decreased, the training continues; otherwise, the process is stopped.
The mean absolute error (MAE) was used as the error function in this algorithm, which is defined as follows:
MAE = i = 1 n | o u t ( i ) o u t p ( i ) | n
where out is the output value, outp is the predicted output value, and n is the number of instances.
The efficacy of the proposed ANN was assessed using the following established performance metrics: Pearson correlation coefficient (PCC), determination coefficient (R2), RMSE, and relative RMSE (RRMSE), as delineated in [60,61]. Values of RMSE and RRMSE approaching zero, along with R2 and PCC values nearing one, indicated a robust alignment between forecast outcomes and empirical data. According to a previous study [62,63], THE estimation accuracy was divided into THE FOLLOWING four categories based on the RRMSE criterion: excellent (below 10%), good (10% to 20%), acceptable (20% to 30%), and unacceptable (above 30%). The performance metrics were calculated as follows:
PCC = ( y e   y ¯ ) ( y e y e ¯ ) ( y y ) 2 ( y e y e ¯ ) 2
-
Determination coefficient (R2):
R 2 = ( y e   y ¯ ) 2 ( y   y ¯ ) 2
-
Root mean square error (RMSE):
RMSE = ( y y e ) 2 n
-
Relative root mean square error (RRMSE):
RRMSE = RMSE y e ¯ × 100
In the final stage, the best ANN was selected based on two criteria as follows:
  • highest accuracy of prediction relative to measured data,
  • smallest difference between the performance metrics calculated for the training, testing, and validation sets.

2.2.3. Application of Artificial Neural Networks to Estimate Mean Nondimensional Sway Force in Regular Waves

The general Equation (3) presented above for calculating CO2 emissions is based on Equation (4) for estimating fuel consumption. A key element in calculating fuel consumption is the power supplied to the DP system over time under specified environmental conditions, such as wind and waves. To comprehensively include the impact of wave and wind interaction on fuel consumption, this article expands Equation (4) and proposes a new approach to calculating fuel consumption by considering the impact of wind and waves as follows:
F C = F O C · t · ( i = 1 8 j = 1 11 p w a v e ( H S ( i ) , T z ( j ) ) · P w a v e ( H S ( i ) , T z ( j ) ) + i = 1 8 p w i n d ( i ) · P w i n d ( i ) )
where
Pwave—power needed to counteract the drift force caused by the wave height Hs and period Tz [kW],
pwave—annual frequency of the wave height (Hs) and period (Tz) in the CCZ according to Table 3 [-],
Pwind—power needed to counteract the drift force caused by wind [kW],
pwind—annual frequency of wind occurrence at a given speed according to Table 4 [-].
The above formula takes fuel consumption calculations into account in the case of the wave interaction time at a given height Hs and period Tz and the time of wind interaction at a given speed. These times are calculated based on the wave occurrence frequency and wind parameters presented in Table 3 and Table 4.
Equation (17) can be used to calculate fuel consumption and considers the full range of wave and wind parameters, as well as taking a limited range of wave and wind conditions into account, e.g., for waves up to 3 m and corresponding wind speeds.

2.2.4. CO2 Emissions during Transport

The estimation of CO2 emissions during maritime transport is based on the energy efficiency design index (EEDI), which, according to the IMO guidelines [31], is defined as the ratio of the CO2 emitted by a ship per nautical mile to its transport capacity and speed. The calculations considered the threshold values of the EEDI that must be met by each bulk carrier.
Figure 5 illustrates the block diagram for estimating CO2 emissions during the transport of PMNs to the destination port. The CO2 emission was calculated using the following formula:
CO2transport = EEDIc · DWT · dT
where
DWT—deadweight tonnage of the bulk carrier from fleet A/B/C [t],
dT—total distance [NM],
EEDIc—the required energy efficiency design index for bulk carriers A/B/C according to the International Maritime Organization [31] guidelines, calculated as follows:
EEDI c = 1 X 100 · 961.79 · DWT 0.477
where
EEDIc—the required energy efficiency design index, expressed as the number of grams of CO2 emissions per nautical mile per ton of cargo,
X—is the reduction factor, whose value for phase 3, according to [31], which starts from 2025, is 30.
The EEDIc index, calculated using Equation (19), serves as the reference line for determining the required emission reductions within the framework of the energy efficiency design index (EEDI). In calculating the EEDIc, various technical and operational factors were considered, including deadweight tonnage (DWT), ship type, technical specifications, and operational speed. Weather conditions were also indirectly accounted for through the inclusion of the sea margin and historical data on the actual fuel consumption and CO2 emissions of ships, which formed the basis for developing the EEDIc. Figure 6 shows the values of EEDIc for different types of ships depending on the DWT, assuming a reduction factor X, applicable from 2025.
The total distance dT was calculated based on the formula:
dT = nload · d
where
d—distance from the CCZ to the destination port,
nload—number of transfers from the extraction ship to the transporting bulk carrier.
The number of PMN transfers from the extraction ship to the transporting bulk carrier was calculated based on the total assumed mass of nodules to be transferred and the deadweight tonnage of carriers from fleets A, B, and C as follows:
nload = M/DWT
where
M—total mass of PMNs to be transported, M = 2,000,000 t.

3. Results

3.1. Prediction of a Nondimensional Sway Force Coefficient Using an ANN

A nondimensional sway force coefficient, Cwy, was predicted using a Multilayer Perceptron (MLP) neural network structure. The model utilized the following ratios and coefficients that impact the sway force coefficient, Cwy:
  • L/B (length-to-breadth ratio),
  • B/T (breadth-to-draught ratio),
  • L/T (length-to-draught ratio, as a combination of L/B and B/T),
  • CB (block coefficient),
  • λ/L (wavelength-to-ship length ratio).
These features were chosen based on their physical significance and impact on ship hydrodynamics, as supported by the literature [39,40,41,42,43,44]. This selection ensured that the model utilized the most informative and relevant variables to accurately predict the Cwy coefficient.
The training dataset was sourced from publicly available ship model tests, as reported in the literature [39,40,41,42,43,44]. The design characteristics of the ship models used in these tests are presented in Table 5. The dataset consisted of 32 cases. To evaluate the model’s performance and prevent overfitting, early stopping was employed during the network’s training. Additionally, the dataset was randomly divided into the following three parts using the hold-out method: training (sixteen cases), testing (eight cases), and validation (eight cases). This approach allowed for the detection and prevention of overfitting by providing separate datasets for model training, validation, and final performance evaluation.
To optimize the learning process, the data were normalized to a range from 0 to 1. This normalization aimed to accelerate learning, avoid numerical issues such as overflow or underflow, prevent feature domination, facilitate weight initialization, and improve the performance of gradient-based algorithms.
During the training of the ANN, additional hyperparameters were used as follows: a learning rate of 0.1, a batch size of 32, and a maximum of 100 epochs. The BFGS optimization algorithm [64,65,66,67] was employed, which has proven to be highly effective for training ANNs to predict additional wave resistance [57]. The most accurate neural network identified through this research contained five neurons in the input layer, six neurons in the hidden layer, and one neuron in the output layer. All neurons were activated using a logistic function as follows:
σ ( x ) = 1 1 + e x
where σ(x) is the value of the activation function for the input x.
Graphically, the general structure of this ANN is illustrated in Figure 7. Normalization coefficients and weights are provided in Appendix A.
Additionally, during training, early stopping was utilized to prevent overfitting by halting the training process if errors increased for the training and testing data. This process was stopped after 54 epochs, with the stopping condition being a change in error of 1 × 10−5.
To further increase the reliability of the results, four different performance evaluation metrics were applied as follows: Pearson correlation coefficient (PCC), determination coefficient (R2), root mean square error (RMSE), and relative root mean square error (RRMSE). Table 6 presents these model evaluation metrics for the ANN. This table shows that the PCC values are close to 1, indicating a strong linear correlation between the predicted and actual values for all datasets. High PCC values suggest that the model predictions are well correlated with the actual data. For the training set, the PCC value is 0.98, for the validation set, it is 0.99, and for the test set, it is 0.97. The R2 values are also high, indicating that the model explains the variability of the data well. High R2 values mean that the model fits the data well for the training, validation, and test sets, where the R² values are 0.95, 0.99, and 0.95, respectively. Low RMSE values indicate small root mean squared errors, meaning that the model predictions are close to the actual values. The RMSE for all datasets is 0.05. RRMSE expressed as a percentage also shows relatively small errors relative to the range of data, with values of 19.58% for the training set, 16.36% for the validation set, and 14.44% for the test set. Lower RRMSE values for the validation and test sets suggest good overall model performance and no signs of overfitting. Despite the small number of measurements, with only 33 data points, the network exhibits good metrics. Based on these parameters, it can be concluded that the neural network model performs very well. High PCC and R² values, along with low RMSE and RRMSE values, indicate accurate and stable predictions for the training, validation, and test sets. The relatively low RMSE indicates that the average deviation of the model predictions from the actual values is small. According to a previous study [48,49], the estimation accuracy can be divided into four categories based on the RRMSE criterion as follows: excellent (below 10%), good (10% to 20%), acceptable (20% to 30%), and unacceptable (above 30%). Thus, the RRMSE values indicate good precision within the model. The model appears to be well-fitted without signs of overfitting. These model evaluation metrics for the test and validation sets demonstrate the sufficient generalization capability of the ANN for new, unseen data, which is crucial for the practical application of the network.
Figure 8 presents a residuals plot against the predicted values, indicating that the largest portion of data yields estimation errors ranging from −0.15 to 0.1.
Figure 9 shows scatter plots which indicate that the trend line between the measured and predicted Cwy values follows the linear y = x line (45°) relatively closely, with slopes of 0.959, 1.0345, and 1.096, and low constants of 0.0097, 0.0003, and 0.0068 for the training, test, and validation data, respectively. This means that predictions using the ANN are characterized by a fairly accurate prediction ability.
Figure 10 compares predictions made by the ANN with results from a model test. These plots show that the predictions obtained with the ANN are consistent with the test results and exhibit corresponding trends.
Figure 11 presents predictions for fleets A, B, and C in relation to the data used for training the neural network. This Figure indicates that these predictions exhibit trends similar to those observed in the model test.

3.2. CO2 Emissions during Nodule Loading in the CCZ

Using the developed ANN and Equations (8)–(11), wave-induced drift force values for bulk carriers A, B, and C were estimated by considering the significant wave height (Hs) and the average zero-up-crossing wave period (Tz) in the CCZ, as shown in Table 3. The results are presented in Table A4, Table A5 and Table A6 in Appendix B. By applying Equation (7) and considering the wind speed distribution shown in Table 4, along with the characteristics of bulk carriers A, B, and C, the wind-induced drift forces were estimated. These are depicted in Figure 12.
When considering the most frequently occurring wave heights in the CCZ, which are 4 m according to Table 3, it is suggested that the largest drift forces, as shown in the figures and tables, are as follows:
  • 775 kN, 1287 kN, and 1688 kN from waves,
  • 106 kN, 244 kN, and 282 kN from wind,
for bulk carriers A, B, and C, respectively.
The thrust force (Tp) and the power (P) supplied to the propellers of the DP system were calculated using Equations (5) and (6). A power-to-thrust efficiency ratio (Cp) of 0.14 was assumed for the power calculations. Then, using Equations (3) and (17), the fuel consumption and the resulting CO2 emissions of bulk carriers A, B, and C were calculated for 250 days a year during concrete loading operations in the CCZ.
Detailed values of fuel consumption resulting from the operation of the DP system to counteract the drift force caused by waves and wind on bulk carriers A, B, and C over 250 days a year during concrete loading in the CCZ can be found in Table A7, Table A8, Table A9 and Table A10 of Appendix B. The values of fuel consumption, which depend solely on Hs, are presented collectively in Figure 13. This Figure illustrates that fuel consumption increases with wave height up to 3 m, reaching a maximum of 600 t, 1250 t, and 2100 t for bulk carriers A, B, and C, respectively. Beyond this level, fuel consumption decreases despite the increase in drift forces. This reduction in fuel consumption occurs because the amount of time the vessel spends in these conditions decreases as the wave height increases, given the lower occurrence of such large waves.
Meanwhile, Figure 14 shows the cumulative CO2 emissions from fuel combustion on ships due to wind and wave action as a function of wave height. For wave heights of up to 4 m, the cumulative CO2 emissions are calculated to be 4900 t, 10,000 t, and 16,500 t for bulk carriers A, B, and C, respectively. These numbers demonstrate that significantly lower amounts of CO2 will be emitted when using type A bulk carriers for the transshipment of extracted concrete. Utilizing type B and C bulk carriers for this mission will result in an increase in CO2 emissions by 100% and 240%, respectively.

3.3. CO2 Emissions during Transport

The number of shipments was calculated using Equation (21), based on the predetermined total mass of PMNs to be transported and the cargo capacities of bulk carriers A, B, and C. The results indicated that there were 264 shipments for bulk carrier A, 63 for B, and 27 for C. Subsequently, utilizing the established routes shown in Figure 1, the distances that bulk carriers A, B, and C need to traverse to reach their designated destination ports were calculated based on Equations (20) and (21). The varying distances, according to the planned routes, are displayed in Table 7, which shows that type A bulk carriers are required to cover the longest distance, while type C has the shortest distance to traverse. These discrepancies increase with route length.
In accordance with a report by the International Maritime Organization [31], the threshold values for the energy efficiency design index (EEDI) were calculated using Equation (19), resulting in 9.48 for bulk carrier A, 4.78 for B, and 3.21 for C. This analysis shows that the largest bulk carriers have the lowest CO2 emissions. Comparatively, type A and B bulk carriers emit nearly 200% and 50% more CO2 than those of type C, respectively.
Based on these data, the total CO2 emissions of bulk carriers A, B, and C, which transport the predetermined mass of PMNs to ports via designated routes, were calculated. The results, graphically presented in Figure 15, indicate that type C bulk carriers will emit the least amount of CO2. Transporting the same amount of extracted PMNs to the farthest destination port using type C bulk carriers will result in CO2 emissions of 998,768 tones, compared to 502,449 tones for type B and as much as 329,338 tones for type A. Therefore, transporting PMNs from the CC zone to the unloading port using type A and B bulk carriers will cause an increase in CO2 emissions by 203% and 53%, respectively, when compared to using type C bulk carriers.

3.4. Assessment of CO2 Emissions from Bulk Carriers during Nodule Loading and Maritime Transport

Based on the CO2 emissions of A, B, and C bulk carriers during the loading of PMNs in the CCZ and during transportation, the total CO2 emissions for the entire mission were calculated. The calculations were conducted under the assumption that the loading in the CCZ is carried out with wave heights up to 3 m and wind speeds up to 12 m/s. The results are presented in Figure 16. On route number 1, which is the shortest, type A, B, and C bulk carriers emit 81,111, 46,058, and 36,902 tons of CO2, respectively. On route number 6, which is the longest, A, B, and C bulk carriers emit 1,002,238, 509,449, and 340,638 tons of CO2 together, respectively. Hence, it follows that, regardless of route length, the smallest bulk carriers emit the largest amount of CO2, and the largest bulk carriers emit the least. Therefore, the amount of CO2 emitted decreases with an increase in bulk carrier size. Moreover, the further the distance to the unloading port of PMNs, the greater the differences in CO2 emissions in favor of large bulk carriers. On the shortest route, type A and B bulk carriers emit 120% and 25% more CO2 than type C bulk carriers, respectively. On the longest route, type A and B bulk carriers emit 194% and 50% more CO2 than type C bulk carriers, respectively.
On the shortest route, the differences in CO2 emissions between small and large bulk carriers are relatively minimal. Figure 17 shows the cumulative CO2 emissions during the extraction and transportation of nodules to the nearest destination port, depending on the wave height (and corresponding wind speed) prevailing during the loading of nodules in the CC zone.
For the smallest bulk carriers, the impact of wind and waves on CO2 emissions is minimal. This is due to the small size of these ships and the need for them to cover longer distances compared to larger ships, which significantly increases CO2 emissions during the voyage. With wave heights up to 1 m and wind speeds up to 7 m/s, these ships emit 77,850 thousand tons of CO2, and with wave heights of 5 m, they emit 83,300 thousand tons of CO2, indicating an increase of only about 7%.
Environmental conditions have a greater impact on medium and large bulk carriers, as the influence of waves and wind on these ships is more significant, increasing CO2 emissions, although the relative emissions during transport are lower. Loading nodules in wave heights up to 5 m and wind speeds up to 19 m/s results in annual CO2 emissions of 51,000 and 45,000 tons for bulk carriers type B and C, respectively. In contrast, loading nodules in wave heights up to 1 m and wind speeds up to 7 m/s results in annual CO2 emissions of 39,500 and 26,100 tons for bulk carriers type B and C, respectively. When compared to emissions at 5 m wave heights, this represents an increase of approximately 29% for bulk carrier type B and as much as 72% for bulk carrier type C. For wave heights up to 2 m and wind speeds up to 10 m/s, the annual emissions are 42,000 and 30,000 tons of CO2 for bulk carriers type B and C, respectively, which corresponds to an increase of approximately 21% for bulk carrier type B and 50% for bulk carrier type C when compared to emissions at 5 m wave heights.
Loading all types of ships in worse weather conditions, i.e., wave heights greater than 5 m and wind speeds above 19 m/s, does not significantly increase CO2 emissions, as such conditions are relatively rare and the duration of the trips in these conditions is short.
Assuming that loading will be conducted in wave conditions up to 1–2 m due to technical issues related to ensuring safe loading, as indicated by [20,21], on the shortest route and on the largest bulk carriers, a significant impact of waves and wind on CO2 emissions can be expected. However, on longer routes, this impact is not substantial.

4. Discussion

The analysis presented in the article indicates that the smallest type A bulk carriers, although characterized by a lower drift force, generate higher CO2 emissions when compared to larger vessels. This is associated with their lower transport efficiency and the greater number of trips required to deliver nodules to the destination port, resulting in higher fuel consumption and, consequently, higher CO2 emissions. In contrast, the largest type C bulk carriers, despite the greater drift force caused by marine environmental effects during loading, emit less CO2 overall, thanks to their higher transport efficiency. These results suggest that selecting larger vessels for the loading and transporting of nodule cargo can positively contribute to reducing the environmental impact of the mission by reducing carbon dioxide emissions.
These studies were conducted with the assumption that all types of ships can conduct loading operations efficiently, even in challenging weather conditions. However, research conducted by Kacprzak [20] has shown that the loading of PMNs on small bulk carriers can only be carried out under limited weather conditions. This means that the use of a fleet of small bulk carriers for the loading and transport of PMNs may lead to an extension of the time required to complete these operations due to the lower number of days in a year that are suitable for loading. These findings suggest that using large bulk carriers for these tasks may be a more efficient solution than using small ones.
Figure 15, Figure 16 and Figure 17 show that the overall CO2 emission related to the transport mission primarily depends on the emissions generated during sea transport. This clearly demonstrates that, although smaller ships may emit less CO2 during the loading operations themselves, their lower transport efficiency, and the need to make a significantly larger number of voyages and thus cover a longer distance, means that when considering the entire venture, their environmental impact is negatively greater than that of large bulk carriers. Combined with the fact that small ships burn relatively more fuel per unit of transported cargo than large ones, this leads to the conclusion that the use of a fleet of large ships could be key from the point of view of ecological efficiency. This discovery may have significant implications for ecological policy in maritime transport. It shows that efforts to limit CO2 emissions in maritime transport should not only include the search for new technologies and fuels, but also the optimization of fleet size and operation planning. Using large bulk carriers for the transportation of all types of cargo could significantly reduce global CO2 emissions from maritime transport, presenting both a challenge and an opportunity for the industry.
Scaling the research findings to encompass general maritime transport clearly shows the substantial potential for reducing CO2 emissions. This is particularly relevant given the feasibility of the annual transportation of up to 11 billion tons of cargo via sea, as reported in UNCTAD [68]. For example, Figure 18 depicts the hypothetical CO2 emissions generated by ships carrying a volume of cargo along a typical route spanning an average of 5000 nautical miles in relation to the ship deadweight capacity. The illustration reveals that utilizing only a fleet of the smallest bulk carriers, each with a deadweight of 1000 tons, would culminate in an annual CO2 emission totaling 1,960,843,563 tons. Conversely, only employing a fleet composed of the largest bulk carriers, each boasting a deadweight of 320,000 tons, would significantly reduce the annual CO2 emission to merely 125,166,010 tons. Therefore, if we assume hypothetically that all the global cargo transported annually was only carried by the smallest bulk carriers, this will increase the emissions by 1,835,677,553 tons of CO2 annually. This roughly equates to the annual burning of 1.142 billion passenger cars with an average consumption of 6 L of diesel per 100 km, assuming that they each travel 10,000 km annually. This example highlights the need for further research and analysis to better understand and utilize the potential for the efficient use of large bulk carriers in ecological maritime transport. The analysis presented in the article potentially represents a significant advancement in understanding the complexities of CO2 emissions in maritime transport. It emphasizes that the ecological efficiency of this sector depends not only on the employed technologies, which primarily involve alternative fuels [69,70,71], but also on operational strategies and planning. Optimizing the size of ships could bring significant benefits to the environment, which should be considered in global efforts towards sustainable development and carbon footprint reduction. However, verifying this approach requires conducting more detailed research in this area, which may be explored in subsequent studies by the authors.

5. Conclusions

The study confirmed the set objectives, providing significant insights into the impact of ship size on CO2 emissions. It introduces new methods for assessing energy efficiency and the use of ANNs to estimate drifting force. The key innovations of the article include the following:
  • the impact of bulk carrier size on CO2 emissions, indicating the possibility of reducing the carbon footprint through the optimal selection of transport units,
  • a new method to assess energy efficiency that takes environmental conditions into account, which can significantly contribute to the optimization of fuel consumption and the reduction of CO2 emissions,
  • the application of ANNs to estimate drifting force, offering new perspectives for ship design and operation, and highlighting the potential of artificial intelligence in the analysis of maritime operations.
The study demonstrated that large bulk carriers, despite the greater drift force acting on them, emit less CO2 than smaller bulk carriers, both when considering the loading at sea of PMNs and their transport to the destination port. This indicates that preferring larger ships in global maritime trade can reduce environmental impacts due to their higher transport efficiency and lower CO2 emissions.
The proposed method to estimate fuel consumption and CO2 emissions, considering environmental conditions with the DP system, could lead to a new energy efficiency indicator. This indicator would aid in assessing and optimizing maritime operations from an environmental perspective, but also requires further research.
Using neural networks (ANNs) to estimate lateral drifting force in regular waves has proven effective, demonstrating its potential in maritime analysis. The developed ANN accurately predicts drifting force based on ship design parameters. However, practical use is limited by the narrow dimension ratio ranges of the training models. Expanding the training dataset with more diverse ship shapes could improve the method’s universality, with further research required.

Author Contributions

All authors contributed equally to this editorial. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All data supporting the reported results are included in the attachments of this article. No new data were created during this study.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

ANNArtificial neural network
AWLateral plane of the vessel above the plane of draught
BBreadth of the ship
BFGSBroyden–Fletcher–Goldfarb–Shanno algorithm
CBBlock coefficient
CCZClarion-Clipperton Zone
CIICarbon intensity indicator
CO2Carbon dioxide
DDisplacement
DPDynamic positioning
DWTDeadweight tonnage
EEDIEnergy efficiency design index
EEOIEnergy efficiency operational indicator
HsSignificant wave height
IMOInternational Maritime Organization
LLength between perpendiculars of the ship
MSEMean squared error
PCCPearson correlation coefficient
PMNsPolymetallic nodules
RMSERoot mean squared error
RRMSERelative root mean square error
SCGScaled conjugate gradient
SGDStochastic gradient descent
TDraught of the ship
TzAverage zero-up-crossing wave period
λWavelength of a regular wave

Appendix A

Table A1. Normalization coefficients for the inputs and output of the ANN.
Table A1. Normalization coefficients for the inputs and output of the ANN.
ShiftScale
L/B.−30.66666666666665.55555555555554
B/T−8.709677419354853.22580645161291
L/T−12.51219512195120.81300813008130
CB−27.000000000000133.3333333333334
l/L−1.538461538461540.76923076923076
Output01.49253731343284
Table A2. Neuron weights and biases in the hidden layer.
Table A2. Neuron weights and biases in the hidden layer.
Neuron in the Hidden Layer
Input123456
L/B0.3012850.3006980.1938480.4277970.160716−0.070965
B/T0.7704410.3707260.7491370.0941640.785890−0.726824
L/T0.7949740.5266351.0353770.8671951.0195680.842477
CB2.5448731.5735342.6298562.7222122.506663−0.238512
l/L3.1111882.3916263.6416225.7050014.0867349.063809
bias0.5584910.4005010.554879−1.2159700.480272−4.862085
Table A3. Neuron weights and biases in the output layer.
Table A3. Neuron weights and biases in the output layer.
Output
Neuron 10.703693
Neuron 2−0.501220
Neuron 31.014894
Neuron 4−1.718017
Neuron 50.080677
Neuron 6−6.121465
bias1.356850

Appendix B

Table A4. Estimated drift force values from beam seas for bulk carrier type A using ANN.
Table A4. Estimated drift force values from beam seas for bulk carrier type A using ANN.
Hs (m)Tz (s)
<44 to 55 to 66 to 77 to 88 to 99 to 1010 to 1111 to 1212 to 13>13
7 to 825663303278719601301861581401284206153
6 to 71965252921341501996659445307217158117
5 to 6144318581568110373248532722616011686
4 to 51002129010897665083362271571118060
3 to 4642826697490325215145100715238
2 to 33614643922761831218256402921
1 to 216020617412381543625181310
0 to 140524431201396432
Table A5. Estimated drift force values from beam seas for bulk carrier type B using ANN.
Table A5. Estimated drift force values from beam seas for bulk carrier type B using ANN.
Hs (m)Tz (s)
<44 to 55 to 66 to 77 to 88 to 99 to 1010 to 1111 to 1212 to 13>13
7 to 81232329140433565271719591393996722532398
6 to 7943251930952729208015001067763553407305
5 to 669318512274200515281102784560406299224
4 to 54811285157913931061765544389282208156
3 to 43088231011891679490348249181133100
2 to 31734635685013822751961401027556
1 to 2772062532231701228762453325
0 to 119516356423122161186
Table A6. Estimated drift force values from beam seas for bulk carrier type C using ANN.
Table A6. Estimated drift force values from beam seas for bulk carrier type C using ANN.
Hs (m)Tz (s)
<44 to 55 to 66 to 77 to 88 to 99 to 1010 to 1111 to 1212 to 13>13
7 to 8888335252565689505540403072229317091283974
6 to 768025664024435538713093235217561309982746
5 to 64991885295732002844227317281290961722548
4 to 5347130920532222197515781200896668501380
3 to 42228381314142212641010768573427321243
2 to 3125471739800711568432322240180137
1 to 2552093293563162531921431078061
0 to 11452828979634836272015
Table A7. Values of fuel consumption FC [t] resulting from the operation of the DP system to counteract the drift force caused by wave action on bulk carrier A.
Table A7. Values of fuel consumption FC [t] resulting from the operation of the DP system to counteract the drift force caused by wave action on bulk carrier A.
Hs (m)Tz (s)
<44 to 55 to 66 to 77 to 88 to 99 to 1010 to 1111 to 1212 to 13>13
7 to 800000100000
6 to 700002332100
5 to 6000391194210
4 to 50031635372310410
3 to 4011463103844215410
2 to 30348147165963610200
1 to 201071116762871000
0 to 103751000000
Table A8. Values of fuel consumption FC [t] resulting from the operation of the DP system to counteract the drift force caused by wave action on bulk carrier B.
Table A8. Values of fuel consumption FC [t] resulting from the operation of the DP system to counteract the drift force caused by wave action on bulk carrier B.
Hs (m)Tz (s)
<44 to 55 to 66 to 77 to 88 to 99 to 1010 to 1111 to 1212 to 13>13
7 to 800000211100
6 to 700004664210
5 to 6000518262111520
4 to 50043074855626921
3 to 40120115215191101381131
2 to 303702683452198625610
1 to 201010321015963173100
0 to 1031093100000
Table A9. Values of fuel consumption FC [t] resulting from the operation of the DP system to counteract the drift force caused by wave action on bulk carrier C.
Table A9. Values of fuel consumption FC [t] resulting from the operation of the DP system to counteract the drift force caused by wave action on bulk carrier C.
Hs (m)Tz (s)
<44 to 55 to 66 to 77 to 88 to 99 to 1010 to 1111 to 1212 to 13>13
7 to 800000332110
6 to 70000713149421
5 to 60008345346251141
4 to 500548137175123592162
3 to 40126183401394223862661
2 to 30391427642452190561330
1 to 2010133335296131378100
0 to 10313146100000
Table A10. Values of fuel consumption (FC) [t] resulting from the operation of the DP system to counteract the drift force caused by wind on bulk carriers A, B, and C.
Table A10. Values of fuel consumption (FC) [t] resulting from the operation of the DP system to counteract the drift force caused by wind on bulk carriers A, B, and C.
Hs (m)vW (m/s)FC (t)
Bulk Carrier ABulk Carrier BBulk Carrier C
7 to 8056128148
6 to 7088202234
5 to 6069159184
4 to 50130297344
3 to 40105241278
2 to 3095218252
1 to 20104239276
0 to 1049112130

References

  1. UNCTAD. Review of Maritime Transport (UNCTAD/RMT/2021). United Nations Conference on Trade and Development. 2021. Available online: https://unctad.org/publication/review-maritime-transport-2021 (accessed on 1 February 2024).
  2. Fu, X.; Lang, J.; Wang, X.; Li, Y.; Lang, J.; Zhou, Y.; Guo, X. The impacts of ship emissions on ozone in eastern China. Sci. Total Environ. 2023, 903, 166252. [Google Scholar] [CrossRef] [PubMed]
  3. Nunes, R.A.O.; Alvim-Ferraz, M.C.M.; Martins, F.G.; Sousa, S.I.V. Environmental and social valuation of shipping emissions on four ports of Portugal. J. Environ. Manag. 2019, 235, 62–69. [Google Scholar] [CrossRef] [PubMed]
  4. Moreno-Gutiérrez, J.; Pájaro-Velázquez, E.; Amado-Sánchez, Y.; Rodríguez-Moreno, R.; Calderay-Cayetano, F.; Durán-Grados, V. Comparative analysis between different methods for calculating on-board ship’s emissions and energy consumption based on operational data. Sci. Total Environ. 2019, 650, 575–584. [Google Scholar] [CrossRef] [PubMed]
  5. Hoang, A.T.; Foley, A.M.; Nižetić, S.; Huang, Z.; Ong, H.C.; Ölçer, A.I.; Pham, V.V.; Nguyen, X.P. Energy-related approach for reduction of CO2 emissions: A critical strategy on the port-to-ship pathway. J. Clean. Prod. 2022, 355, 131772. [Google Scholar] [CrossRef]
  6. Fan, A.; Li, Y.; Liu, H.; Yang, L.; Tian, Z.; Li, Y.; Vladimir, N. Development trend and hotspot analysis of ship energy management. J. Clean. Prod. 2023, 389, 135899. [Google Scholar] [CrossRef]
  7. Barone, G.; Buonomano, A.; Del Papa, G.; Maka, R.; Palombo, A. How to achieve energy efficiency and sustainability of large ships: A new tool to optimize the operation of on-board diesel generators. Energy 2023, 282, 128288. [Google Scholar] [CrossRef]
  8. Nyanya, M.N.; Vu, H.B.; Schönborn, A.; Ölçer, A.I. Wind and solar assisted ship propulsion optimisation and its application to a bulk carrier. Sustain. Energy Technol. Assess. 2021, 47, 101397. [Google Scholar] [CrossRef]
  9. Ytreberg, E.; Astrom, S.; Fridell, E. Valuating environmental impacts from ship emissions—The marine perspective. J. Environ. Manag. 2021, 2021, 111958. [Google Scholar] [CrossRef] [PubMed]
  10. Amudha, K.; Bhattacharya, S.; Sharma, R.; Gopkumar, K.; Kumar, D.; Ramadass, G. Influence of flow area zone and vertical lift motion of polymetallic nodules in hydraulic collecting. Ocean. Eng. 2024, 294, 116745. [Google Scholar] [CrossRef]
  11. Sha, F.; Xi, M.; Chen, X.; Liu, X.; Niu, H.; Zuo, Y. A recent review on multi-physics coupling between deep-sea mining equipment and marine sediment. Ocean. Eng. 2023, 276, 114229. [Google Scholar] [CrossRef]
  12. Cunningham, A. Assessing the feasibility of deep-seabed mining of polymetallic nodules in the Area of seabed and ocean floor beyond the limits of national jurisdiction, as a method of alleviating supply-side issues for cobalt to US markets. Miner. Econ. 2024, 37, 207–226. [Google Scholar] [CrossRef]
  13. Gollner, S.; Haeckel, M.; Janssen, F.; Lefaible, N.; Molari, M.; Papadopoulou, S.; Reichart, G.J.; Trabucho Alexandre, J.; Vink, A.; Vanreusel, A. Restoration experiments in polymetallic nodule areas. Integr. Environ. Assess. Manag. 2022, 18, 682–696. [Google Scholar] [CrossRef] [PubMed]
  14. Bonifácio, P.; Kaiser, S.; Washburn, T.W.; Smith, C.R.; Vink, A.; Arbizu, P.M. Biodiversity of the Clarion-Clipperton Fracture Zone: A worm perspective. Mar. Biodivers. 2024, 54, 5. [Google Scholar] [CrossRef]
  15. Uhlenkott, K.; Meyn, K.; Vink, A.; Martínez Arbizu, P. A review of megafauna diversity and abundance in an exploration area for polymetallic nodules in the eastern part of the Clarion Clipperton Fracture Zone (North East Pacific), and implications for potential future deep-sea mining in this area. Mar. Biodivers. 2023, 53, 22. [Google Scholar] [CrossRef]
  16. Stratmann, T. Role of polymetallic-nodule dependent fauna on carbon cycling in the eastern Clarion-Clipperton Fracture Zone (Pacific). Front. Mar. Sci. 2023, 10, 1151442. [Google Scholar] [CrossRef]
  17. Li, B.; Jia, Y.; Fan, Z.; Li, K.; Shi, X. Impact of the Mining Process on the Near-Seabed Environment of a Polymetallic Nodule Area: A Field Simulation Experiment in a Western Pacific Area. Sensors 2023, 23, 8110. [Google Scholar] [CrossRef] [PubMed]
  18. Katona, S.; Paulikas, D.; Stone, G.S. Ethical opportunities in deep-sea collection of polymetallic nodules from the Clarion-Clipperton Zone. Integr. Environ. Assess. Manag. 2022, 18, 634–654. [Google Scholar] [CrossRef]
  19. Shobayo, P.; van Hassel, E.; Vanelslander, T. Logistical Assessment of Deep-Sea Polymetallic Nodules Transport from an Offshore to an Onshore Location Using a Multiobjective Optimization Approach. Sustainability 2023, 15, 11317. [Google Scholar] [CrossRef]
  20. Kacprzak, P. Assessment of cargo handling operation efficiency in the CCZ for standard bulk carriers in the view of significant amplitudes of roll as a limiting criterion. Sci. J. Marit. Univ. Szczec. 2023, 74, 55–64. [Google Scholar] [CrossRef]
  21. Kacprzak, P. An analysis of shear forces, bending moments and roll motion during a nodule loading simulation for a ship at sea in the Clarion—Clipperton Zone. Sci. J. Marit. Univ. Szczec. 2021, 65, 9–20. [Google Scholar] [CrossRef]
  22. Marques, C.H.; Belchiora, C.R.P.; Capracea, J.-D. Optimising the engine-propeller matching for a liquefied natural gas carrier T under rough weather. Appl. Energy 2018, 232, 187–196. [Google Scholar] [CrossRef]
  23. Jurdziński, M. Processes of a freely drifting vessel. TransNav 2020, 14, 687–693. [Google Scholar] [CrossRef]
  24. Bialystocki, N.; Konovessis, D. On the estimation of ship’s fuel consumption and speed curve: A statistical approach. J. Ocean. Eng. Sci. 2016, 1, 157–166. [Google Scholar] [CrossRef]
  25. Bal Beşikçi, E.; Arslan, O.; Turan, O.; Ölçer, A. An ANN based decision support system for energy efficient ship operations. Comput. Oper. Res. 2016, 66, 393–401. [Google Scholar] [CrossRef]
  26. Petersen, J.P.; Jacobsen, D.J.; Winther, O. Statistical modelling for ship propulsion efficiency. J. Mar. Sci. Technol. 2011, 17, 30–39. [Google Scholar] [CrossRef]
  27. Cepowski, T.; Drozd, A. Measurement-based relationships between container ship operating parameters and fuel consumption. Appl. Energy 2023, 347, 121315. [Google Scholar] [CrossRef]
  28. Panagakos, G.; de Pessôa, T.S.; Dessypris, N.; Barfod, M.B.; Psaraftis, H.N. Monitoring the carbon footprint of dry bulk shipping in the EU: An early assessment of the MRV regulation. Sustainability 2019, 11, 5133. [Google Scholar] [CrossRef]
  29. Kanberoğlu, B.; Kökkülünk, G. Assessment of CO2 emissions for a bulk carrier fleet. J. Clean. Prod. 2021, 283, 124590. [Google Scholar] [CrossRef]
  30. IMO. Resolution MEPC.1/Circ.684 Guidelines for Voluntary Use of The Ship Energy Efficiency Operational Indicator (EEOI); IMO: London, UK, 2009. [Google Scholar]
  31. IMO. Resolution Res MEPC.364(79) Guidelines on the Method of Calculation of the Attained Energy Efficiency Design Index (EEDI) for New Ships; IMO: London, UK, 2022. [Google Scholar]
  32. IMO. Resolution MEPC.346(78) Guidelines for the Development of a Ship Energy Efficiency. Management Plan (SEEMP); IMO: London, UK, 2022. [Google Scholar]
  33. IMO. Resolution MEPC.336(76) Guidelines on Operational Carbon Intensity Indicators and the Calculation Methods (CII GUIDELINES, G1); IMO: London, UK, 2021. [Google Scholar]
  34. Chen, X.; Lv, S.; Shang, W.; Wu, H.; Xian, J.; Song, C. Ship energy consumption analysis and carbon emission exploitation via spatial-temporal maritime data. Appl. Energy 2024, 360, 122886. [Google Scholar] [CrossRef]
  35. Lee, S.-S. Analysis of the effects of EEDI and EEXI implementation on CO2 emissions reduction in ships. Ocean. Eng. 2024, 295, 116877. [Google Scholar] [CrossRef]
  36. Van Laar, G. Sustainable Transport of Polymetallic Nodules [TU Delft]. 2021. Available online: https://repository.tudelft.nl/islandora/object/uuid%3Aac42edba-a4ca-49ec-bb99-1c1b0b203f20 (accessed on 3 February 2024).
  37. Lipton, I.T.; Nimmo, M.J. NI 43-101 Technical Report TOML Clarion Clipperton Zone Project, Pacific Ocean (Issue July). 2016. Available online: https://int.nyt.com/data/documenttools/2021-03-metals-company-technical-report-on-toml-mining-zone-plan/2d5350243bade994/full.pdf (accessed on 8 February 2024).
  38. Dudziak, J. Teoria Okrętu (II); Fundacja Promocji Przemysłu Okrętowego i Gospodarki Morskiej: Gdańsk, Poland, 2008; ISBN 978-83-60584-09-5. [Google Scholar]
  39. Kadomatsu, K.; Inoue, Y.; Takarada, N. On the Required Minimum Output of Main Propulsion Engine for Large Fat Ship with Considering Manoeuverability in Rough Seas. J. Soc. Nav. Arch. Jpn. 1990, 1990, 171–182. [Google Scholar] [CrossRef] [PubMed]
  40. Sprenger, F.; Maron, A.; Delefortrie, G.; Van Zwijnsvoorde, T.; Cura-Hochbaum, A.; Lengwinat, A.; Papanikolaou, A. Experimental Studies on Seakeeping and Manoeuvrability in Adverse Conditions. J. Ship Res. 2017, 61, 131–152. [Google Scholar] [CrossRef]
  41. Shigunov, V.; El Moctar, O.; Papanikolaou, A.; Pottho, R.; Liu, S. International benchmark study on numerical simulation methods for prediction of manoeuvrability of ships in waves. Ocean Eng. 2018, 165, 365–385. [Google Scholar] [CrossRef]
  42. Yasukawa, H.; Hirata, N.; Matsumoto, A.; Kuroiwa, R.; Mizokami, S. Evaluations of wave-induced steady forces and turning motion of a full hull ship in waves. J. Mar. Sci. Technol. 2018, 24, 1–15. [Google Scholar] [CrossRef]
  43. Ueno, M.; Nimura, T.; Miyazaki, H.; Nonaka, K. Steady Wave Forces and Moment Acting on Ships in Manoeuvring Motion in Short Waves. J. Soc. Nav. Arch. Jpn. 2000, 2000, 163–172. [Google Scholar] [CrossRef] [PubMed]
  44. Ueno, M.; Nimura, T.; Miyazaki, H.; Nonaka, K.; Haraguchi, T. Model Experiment on Steady Wave Forces and Moment Acting on a Ship a Rest. J. Kansai Soc. Nav. Archit. 2001, 235, 69–77. [Google Scholar]
  45. Bortnowska, M. The Determination of Power Output of the Tracking Control System of a Mining Ship Zeszyty Naukowe NR 10 Akademii Morskiej W Szczecinie. Sci. J. Marit. Univ. Szczec. 2006, 10, 97–106. [Google Scholar]
  46. ITTC. ITTC—Recommended Procedures and Guidelines-Seakeeping Experiments. In Proceedings of the 27th ITTC Seakeeping Committee, Copenhagen, Denmark, 31 August–5 September 2014; pp. 1–22. [Google Scholar]
  47. Maruo, H. The drift of a body floating on waves. J. Ship Res. 1960, 4, 1–10. [Google Scholar]
  48. Newman, J.N. The drift force and moment on ships in waves. J. Ship Res. 1967, 11, 51–60. [Google Scholar] [CrossRef]
  49. Salvesen, N. Second-Order Steady State Forces and Moments on Surface Ships in Oblique Regular Waves. In International Symposium on Dynamics of Marine Vehicles and Structures in Waves. Mech. Eng. 1974, 22, 225–241. [Google Scholar]
  50. Gerritsma, J.; Beukelman, W. Analysis of the resistance increase in waves of a fast cargo ship 12. Int. Shipbuild. Prog. 1972, 19, 285–293. [Google Scholar] [CrossRef]
  51. Kashiwagi, M. Added Resistance, Wave-Induced Steady Sway Force and Yaw Moment on an Advancing Ship. Ship Technol. Res. (Schistech.) 1992, 39, 3–16. [Google Scholar]
  52. Boese, P. Eine einfache Methode zur Berechnung der Widerstandserhöhung eines Schiffes im Seegang. Ship Technol. Res. 1970, 258, 17. [Google Scholar]
  53. Faltinsen, O.M.; Minsaas, K.J.; Liapis, N.; Skjordal, S.O. Prediction of Resistance and Propulsion of a Ship in a Seaway. In Proceedings of the 13th Symposium on Naval Hydrodynamics, Tokyo, Japan, 6–10 October 1980; pp. 505–529. [Google Scholar]
  54. Papanikolaou, A.; Nowacki, Z. Second-Order Theory of Oscillating Cylinders in a Regular Steep Wave. In Proceedings of the l3th ONR Symp., Tokyo, Japan, 6–10 October 1980; pp. 303–333. [Google Scholar]
  55. Liu, S.; Papanikolaou, A. Prediction of the Side Drift Force of Full Ships Advancing in Waves at Low Speeds. J. Mar. Sci. Eng. 2020, 8, 377. [Google Scholar] [CrossRef]
  56. Cepowski, T. The prediction of ship added resistance at the preliminary design stage by the use of an artificial neural network. Ocean. Eng. 2020, 195, 106657. [Google Scholar] [CrossRef]
  57. Cepowski, T. The use of a set of artificial neural networks to predict added resistance in head waves at the parametric ship design stage. Ocean. Eng. 2023, 281, 114744. [Google Scholar] [CrossRef]
  58. Ahn, Y.; Kim, Y. Data mining in sloshing experiment database and application of neural network for extreme load prediction. Mar. Struct. 2021, 80, 103074. [Google Scholar] [CrossRef]
  59. Dyer, A.S.; Zaengle, D.; Nelson, J.R.; Duran, R.; Wenzlick, M.; Wingo, P.C.; Bauer, J.R.; Rose, K.; Romeo, L. Applied machine learning model comparison: Predicting offshore platform integrity with gradient boosting algorithms and neural networks. Mar. Struct. 2022, 83, 103152. [Google Scholar] [CrossRef]
  60. Citakoglu, H.; Coşkun, Ö. Comparison of hybrid machine learning methods for the prediction of short-term meteorological droughts of Sakarya Meteorological Station in Turkey. Environ. Sci. Pollut. Res. 2022, 29, 75487–75511. [Google Scholar] [CrossRef]
  61. Demir, V.; Citakoglu, H. Forecasting of solar radiation using different machine learning approaches. Neural Comput. Applic. 2023, 35, 887–906. [Google Scholar] [CrossRef]
  62. Premalatha, M.; Naveen, C. Analysis of different combi- nations of meteorological parameters in predicting the horizontal global solar radiation with ann approach: A case study. Renew. Sustain. Energy Rev. 2018, 91, 248–258. [Google Scholar] [CrossRef]
  63. Uncuoglu, E.; Citakoglu, H.; Latifoglu, L.; Bayram, S.; Laman, M.; Ilkentapar, M.; Oner, A.A. Comparison of neural network, Gaussian regression, support vector machine, long short-term memory, multi-gene genetic programming, and M5 Trees methods for solving civil enginee Energies ring problems. Appl. Soft Comput. 2022, 129, 109623. [Google Scholar] [CrossRef]
  64. Broyden, C.G. The convergence of a class of double-rank minimization algorithms 1. General considerations. IMA J. Appl. Math. 1970, 6, 76–90. [Google Scholar] [CrossRef]
  65. Fletcher, R. A new approach to variable metric algorithms. Comput. J. 1970, 13, 317–322. [Google Scholar] [CrossRef]
  66. Goldfarb, D. A family of variable-metric methods derived by variational means. Math. Comput. 1970, 24, 23–26. [Google Scholar] [CrossRef]
  67. Shanno, D.F. Conditioning of quasi-newton methods for function minimization. Math. Comput. 1970, 24, 647–656. [Google Scholar] [CrossRef]
  68. UNCTAD. Review of Maritime Transport (UNCTAD/RMT/2023). United Nations Conference on Trade and Development. 2023. Available online: https://unctad.org/publication/review-maritime-transport-2023 (accessed on 2 March 2024).
  69. Obydenkova, S.V.; Defauw, L.V.E.; Kouris, P.D.; Smeulders, D.M.J.; Boot, M.D.; van der Meer, Y. Environmental and Economic Assessment of a Novel Solvolysis-Based Biorefinery Producing Lignin-Derived Marine Biofuel and Cellulosic Ethanol. Energies 2022, 15, 5007. [Google Scholar] [CrossRef]
  70. Bortuzzo, V.; Bertagna, S.; Bucci, V. Mitigation of CO2 Emissions from Commercial Ships: Evaluation of the Technology Readiness Level of Carbon Capture Systems. Energies 2023, 16, 3646. [Google Scholar] [CrossRef]
  71. Richter, S.; Braun-Unkhoff, M.; Hasselwander, S.; Haas, S. Evaluation of the Applicability of Synthetic Fuels and Their Life Cycle Analyses. Energies 2024, 17, 981. [Google Scholar] [CrossRef]
Figure 1. Location of destination ports.
Figure 1. Location of destination ports.
Energies 17 03383 g001
Figure 2. Drift force coefficients of four ships in regular beam waves at zero speed.
Figure 2. Drift force coefficients of four ships in regular beam waves at zero speed.
Energies 17 03383 g002
Figure 3. General research scheme for the calculation of CO2 emissions during loading.
Figure 3. General research scheme for the calculation of CO2 emissions during loading.
Energies 17 03383 g003
Figure 4. Research scheme for developing an ANN.
Figure 4. Research scheme for developing an ANN.
Energies 17 03383 g004
Figure 5. General research scheme for calculating CO2 emissions during transport.
Figure 5. General research scheme for calculating CO2 emissions during transport.
Energies 17 03383 g005
Figure 6. Values of EEDIc for different types of ships as a function of DWT, assuming a reduction factor X, applicable from 2025.
Figure 6. Values of EEDIc for different types of ships as a function of DWT, assuming a reduction factor X, applicable from 2025.
Energies 17 03383 g006
Figure 7. The general structure of an ANN.
Figure 7. The general structure of an ANN.
Energies 17 03383 g007
Figure 8. Residuals vs. predicted values.
Figure 8. Residuals vs. predicted values.
Energies 17 03383 g008
Figure 9. Scatter plots of predicted vs. measured Cwy values.
Figure 9. Scatter plots of predicted vs. measured Cwy values.
Energies 17 03383 g009
Figure 10. Comparison of ANN predictions with model test results.
Figure 10. Comparison of ANN predictions with model test results.
Energies 17 03383 g010
Figure 11. Comparison of ANN predictions for bulk carriers types A, B, and C with model test results.
Figure 11. Comparison of ANN predictions for bulk carriers types A, B, and C with model test results.
Energies 17 03383 g011
Figure 12. Drift force values from wind for bulk carriers A, B, and C which correspond to wind velocity (vW).
Figure 12. Drift force values from wind for bulk carriers A, B, and C which correspond to wind velocity (vW).
Energies 17 03383 g012
Figure 13. The values of fuel consumption (FC) [t] resulting from the operation of the DP system to counteract the drift force caused by wave and wind action on bulk carriers A, B, and C, corresponding to significant wave height (Hs).
Figure 13. The values of fuel consumption (FC) [t] resulting from the operation of the DP system to counteract the drift force caused by wave and wind action on bulk carriers A, B, and C, corresponding to significant wave height (Hs).
Energies 17 03383 g013
Figure 14. Cumulative CO2 emissions due to combined wave and wind action in relation to significant wave height (Hs).
Figure 14. Cumulative CO2 emissions due to combined wave and wind action in relation to significant wave height (Hs).
Energies 17 03383 g014
Figure 15. Total CO2 emissions of A, B, and C bulk carriers during transportation.
Figure 15. Total CO2 emissions of A, B, and C bulk carriers during transportation.
Energies 17 03383 g015
Figure 16. Total CO2 emissions of A, B, and C bulk carriers during loading and transportation.
Figure 16. Total CO2 emissions of A, B, and C bulk carriers during loading and transportation.
Energies 17 03383 g016
Figure 17. Cumulative CO2 emissions of A, B, and C bulk carriers during loading and transportation on route no. 1, as influenced by environmental conditions represented by significant wave height (Hs).
Figure 17. Cumulative CO2 emissions of A, B, and C bulk carriers during loading and transportation on route no. 1, as influenced by environmental conditions represented by significant wave height (Hs).
Energies 17 03383 g017
Figure 18. Comparative analysis of annual CO2 emissions: smallest vs. largest bulk carriers on a 5000-nautical mile route transporting 11 billion tons of cargo.
Figure 18. Comparative analysis of annual CO2 emissions: smallest vs. largest bulk carriers on a 5000-nautical mile route transporting 11 billion tons of cargo.
Energies 17 03383 g018
Table 1. Technical parameters of the studied bulk carriers, where L—length between perpendiculars, B—breadth, T—draught, CB—block coefficient, D—displacement, DWT—deadweight, and AW—lateral plane of the vessel above the plane of draught.
Table 1. Technical parameters of the studied bulk carriers, where L—length between perpendiculars, B—breadth, T—draught, CB—block coefficient, D—displacement, DWT—deadweight, and AW—lateral plane of the vessel above the plane of draught.
Main ParticularsShip Size
ABC
L [m]103.9185217
B [m]18.224.432.26
T [m]7.05711.0114.02
L/B [-]5.717.586.73
B/T [-]2.582.222.30
L/T [-]14.7216.8015.48
CB [-]0.800.820.85
D [t]11,03641,90085,700
DWT [t]760032,00073,600
AW [m2]74617111977
Table 2. Parameters of the mining ship Hidden Gem.
Table 2. Parameters of the mining ship Hidden Gem.
L [m]228
B [m]42
T [m]12
CB [-]0.82
D [t]96,504
DWT [t]61,042
Table 3. Annual frequency of wave parameter occurrence pwave in the CCZ. Data adapted from [37], where Hs—significant wave height and Tz—average zero-up-crossing wave period.
Table 3. Annual frequency of wave parameter occurrence pwave in the CCZ. Data adapted from [37], where Hs—significant wave height and Tz—average zero-up-crossing wave period.
Hs (m)Tz (s)
<44 to 55 to 66 to 77 to 88 to 99 to 1010 to 1111 to 1212 to 13>13
7 to 8000000.00010.00010.00010.00010.00010
6 to 700000.00020.00050.00070.00060.00040.00020.0001
5 to 60000.00030.00140.00270.00310.00230.00130.00060.0002
4 to 5000.00030.00250.00810.01290.0120.00770.00370.00140.0005
3 to 400.00010.00230.0150.0370.04550.03390.01760.0070.00230.0006
2 to 300.00080.01430.06230.10540.09280.05140.02040.00640.00170.0004
1 to 20.00010.00570.04740.110.10920.06050.02230.00620.00140.00030
0 to 10.00050.00640.01910.01850.00820.00210.00040.0001000
Table 4. Distribution of wind speed. Data adapted from [38], where Hs—the significant wave height, vW—the wind velocity, and pwind—the annual frequency of wind occurrence.
Table 4. Distribution of wind speed. Data adapted from [38], where Hs—the significant wave height, vW—the wind velocity, and pwind—the annual frequency of wind occurrence.
Hs (m)vW (m/s) pwind [-]
17.30.240
210.00.270
312.10.170
415.40.115
519.30.091
622.10.037
724.90.037
827.00.020
Table 5. Design characteristics of ship model subjected to model test on mean drift force, where L—length between perpendiculars, B—breadth, T—draught, and CB—block coefficient.
Table 5. Design characteristics of ship model subjected to model test on mean drift force, where L—length between perpendiculars, B—breadth, T—draught, and CB—block coefficient.
ShipL/B (-)B/T (-)L/T (-)CB (-)References
VLCCb5.523.0116.620.81[43,44]
KVLCC25.522.7915.40.81[40,41]
Bulk carrier5.72.715.390.83[39]
Scb845.522.7915.40.84[42]
Table 6. Model evaluation metrics for the ANN.
Table 6. Model evaluation metrics for the ANN.
PCC train [-]0.98
PCC valid [-]0.99
PCC test [-]0.97
R2 train [-]0.95
R2 valid [-]0.99
R2 test [-]0.95
RMSE train [-]0.05
RMSE valid [-]0.05
RMSE test [-]0.05
RRMSE train [%]19.58
RRMSE valid [%]16.36
RRMSE test [%]14.44
Table 7. Total distance to be covered by A, B, and C bulk carriers.
Table 7. Total distance to be covered by A, B, and C bulk carriers.
Route Number According to Figure 1Route Length from the CC Zone to the Destination Port [NM]Total Distance [NM]
ABC
120441,077,083255,475108,321
224001,264,800300,000127,200
369713,673,506871,325369,442
416,6968,798,5812,086,950884,867
522,51211,863,6662,813,9631,193,120
626,29113,855,4623,286,4001,393,434
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cepowski, T.; Kacprzak, P. Reducing CO2 Emissions through the Strategic Optimization of a Bulk Carrier Fleet for Loading and Transporting Polymetallic Nodules from the Clarion-Clipperton Zone. Energies 2024, 17, 3383. https://doi.org/10.3390/en17143383

AMA Style

Cepowski T, Kacprzak P. Reducing CO2 Emissions through the Strategic Optimization of a Bulk Carrier Fleet for Loading and Transporting Polymetallic Nodules from the Clarion-Clipperton Zone. Energies. 2024; 17(14):3383. https://doi.org/10.3390/en17143383

Chicago/Turabian Style

Cepowski, Tomasz, and Paweł Kacprzak. 2024. "Reducing CO2 Emissions through the Strategic Optimization of a Bulk Carrier Fleet for Loading and Transporting Polymetallic Nodules from the Clarion-Clipperton Zone" Energies 17, no. 14: 3383. https://doi.org/10.3390/en17143383

APA Style

Cepowski, T., & Kacprzak, P. (2024). Reducing CO2 Emissions through the Strategic Optimization of a Bulk Carrier Fleet for Loading and Transporting Polymetallic Nodules from the Clarion-Clipperton Zone. Energies, 17(14), 3383. https://doi.org/10.3390/en17143383

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