1. Introduction
The shipping industry makes a significant contribution to global trade, with over 80% of goods being transported by sea annually. Although carbon emissions from shipping constitute only a small portion of global carbon emissions, the share of shipping emissions in global anthropogenic emissions increased from 2.76% in 2012 to 2.89% in 2018 [
1]. According to the United Nations Conference on Trade and Development’s (UNCTAD) “Review of Maritime Transport 2023,” this proportion has further risen to 3% [
2]. In response, the International Maritime Organization (IMO) set forth short-term and long-term decarbonization targets in 2023: a 40% reduction in carbon intensity by 2030 and near-zero emissions by around 2050 [
3]. The short-term targets include measures such as the Energy Efficiency Existing Ship Index (EEXI), Carbon Intensity Indicator (CII), and Ship Energy Efficiency Management Plan (SEEMP) [
4]. At the same time, as the global economy enters a downturn and the shipping market continues to decline with decreasing freight rates, fuel costs account for more than 50% of the daily operational costs of ships [
5]. Currently, emission reduction measures for existing operational ships primarily include ship intelligence, the application of energy-saving devices, route optimization, and speed optimization. Among these, speed optimization [
6] has a particularly significant impact on energy savings and emission reductions for existing ships, simultaneously lowering operational costs and addressing the increasingly stringent emission targets set by the IMO.
In recent years, many scholars have conducted studies on ship speed reduction. Degiuli [
7] studied speed reduction in Panamax container ships in specific areas of the Mediterranean Sea, and the author indicated that a 13.6% reduction in speed could save up to 31% of fuel consumption and greenhouse gas emissions in all surveyed areas. Farkas et al. [
8] conducted studies on Panamax container ships on fixed routes, analyzing speed reductions of 10% to 40% annually and monthly. They found that an annual speed reduction of 40% could reduce fuel consumption per nautical mile by up to 29.2% on the North Atlantic route.
Currently, there are three main types of ship fuel consumption models: white-box models, black-box models, and grey-box models [
9]. The white-box model is based on the principles of ship propulsion. Tillig et al. [
10] established a white-box model for the main engine’s fuel consumption based on the principle of ship–engine–propeller matching. This model is very intuitive, allowing a direct functional relationship between the main engine’s fuel consumption and the ship’s speed [
11]. Researchers use the principle of ship–engine–propeller matching to build a white-box model using simulation platforms. Tran et al. [
12] developed a ship energy consumption model directly using the SIMULINK R2023a platform. However, in actual ship operations, the sea conditions are very complex, and the white-box model often cannot adequately reflect the impact of the external environment on the main engine. Therefore, the accuracy of the white-box model is less ideal under high wind and wave conditions [
13]. The black-box model is a purely data-driven fuel consumption model, mainly divided into statistical models and machine learning models [
14]. Unlike the white-box model, which is entirely based on physical principles, the black-box model cannot directly show the functional relationship between fuel consumption and ship speed. Ma et al. [
15] used the main engine’s speed as the optimization variable, optimizing for fuel consumption and sailing time. They employed the K-means clustering algorithm to classify sea conditions based on wind speed, wind direction, and current speed and used the multi-objective particle swarm optimization (MOPSO) algorithm to optimize fuel consumption across different segments. Due to the lack of a clear functional relationship, the accuracy of the black-box model is highly dependent on the settings of model parameters, which significantly affect the results [
16].
The grey-box model combines physical principles with data analysis [
17], offering the advantages of both white-box and black-box models. Grey-box models can be structured in series or parallel configurations. Cai [
18] used big data from ships to develop an efficient and reliable fuel consumption prediction model, exploring the impact of data diversity, quality, and quantity on black-box and grey-box models. They demonstrated that grey-box models could achieve higher accuracy with less data compared to black-box models. Ma et al. [
19] optimized ship speed and routes using the NSGA2 algorithm, targeting shipping costs and emissions as optimization objectives. These studies demonstrate that it is feasible to optimize ship energy consumption by controlling the main engine’s speed through the establishment of grey-box models for ship fuel consumption based on different sea conditions. Current research on ship energy optimization based on ship–engine–propeller matching mostly focuses on overall speed reduction for routes or speed optimization for specific segments. However, the frequency of high wind and wave conditions on major routes is increasing in some cases. For example, during winter, it is nearly impossible to completely avoid high wind and waves when navigating across the Pacific and Atlantic Oceans in mid to high latitudes. Successive cyclones or cyclone groups can cause severe sea conditions spanning thousands of nautical miles, making avoidance impractical. Therefore, this study establishes a ship energy consumption model that considers the increased impact of high wind and wave conditions on fuel consumption by using the main engine’s speed as a variable. It provides an operational solution for actual severe sea conditions.
In this study, a white-box model for ship fuel consumption was established on the SIMULINK simulation platform based on the principle of ship–engine–propeller matching. On the MATLAB platform, ship noon reports (including the daily average main engine power, daily average speed, daily average fuel consumption, daily average ship speed, and 24 h sailing distance), AIS data (including heading, ship speed, and other navigational information), onboard sensor data, weather forecasts, and data from the fifth generation of the ECMWF atmospheric reanalysis dataset (ERA5) [
20] were integrated. The K-Medoids clustering algorithm was used to cluster the collected sea condition data. Using the Random Forest Regression (RFR) algorithm, ship fuel consumption and speed were predicted. The non-dominated sorting genetic algorithm II (NSGA-II) [
21] was optimized, incorporating variable weighting, discrete optimization, and the ideal point method. For different sea conditions, the weights were adjusted accordingly. Using the main engine’s speed as the optimization variable and fuel consumption per nautical mile and ship speed as optimization objectives, a black-box model for the main engine’s fuel consumption was established. The black-box model was combined with the white-box model in series to develop a grey-box model for optimizing ship fuel consumption under high wind and wave conditions.
3. Improvement in the Non-Dominated Sorting Genetic Algorithm
The NSGA-2 algorithm [
19], a commonly used multi-objective optimization algorithm, can find better solutions and convergence distributions in most problems. It is suitable for various engineering applications, providing multiple solutions for conflicting objective functions, and has good applicability in solving speed optimization problems [
17]. Based on the black-box model for fuel consumption established in the previous section, the results can form a functional relationship between the main engine speed as the optimization variable and the fuel consumption per nautical mile and speed as optimization objectives. The equation for fuel consumption per nautical mile under a specific sea condition is shown in Equation (5):
where
is the speed under the corresponding sea condition, and
is the main engine’s speed under the corresponding sea condition.
The optimization objective functions shown in Equations (6) and (7) can use the main engine’s fuel consumption
and sailing speed
as optimization objectives during the actual navigation process, with speed
as the optimization variable.
The objective function for bi-objective optimization is as follows:
Since there may be multiple solutions with the same non-dominance level, the initial NSGA-2 algorithm can yield multiple feasible solutions, such as solutions that favor reducing fuel consumption per nautical mile, increasing speed, or balancing both. These three types of solutions may each have multiple instances, leading to different results in each calculation, causing instability and non-uniqueness in the solutions, and making it difficult to determine the optimal speed for each sea condition. To address the shortcomings of the NSGA-2 algorithm in this study, the algorithm was optimized by constraining the main engine speed to determine the optimal speed for the main engine under various severe sea conditions. The specific optimization steps are as follows:
- 1.
Calculate the maximum and minimum main engine speeds for the current sea condition’s allowable critical safe speed under the main engine’s propulsion power using the method proposed by Aertssen [
28] and the minimum safe propulsion power provided by IMO [
29], and then normalize these values.
In Equation (9),
represents the ship’s voluntarily reduced speed,
denotes the ship’s design service speed, and ‘M’ and ‘N’ are determined by the encounter angle and wind speed, respectively.
In Equation (10), DWT represents the ship’s deadweight tonnage, and ‘a’ and ‘b’ are two coefficients recommended in the guidelines.
- 2.
Apply the variable weighting method to increment and from 0 to 1 at intervals of 0.01 in the multi-objective optimization problem, using the linear weighted sum method to transform the multi-objective speed optimization problem into a single-objective speed optimization problem. Restrict the weights of and under different sea conditions: prioritize fuel consumption per nautical mile optimization ( > ) for sea conditions 1 to 3 and prioritize navigation safety (<) for sea conditions 4 to 6.
- 3.
Discretize the main engine speed in 1 r/min increments and introduce auxiliary decision variables according to the target requirements, further transforming the weighted single-objective nonlinear optimization problem into a 0–1 mixed-integer linear programming problem.
- 4.
Solve the problem multiple times under different weights to obtain the Pareto optimal solution set.
- 5.
Use the ideal point method to balance the decisions from the obtained Pareto optimal solution set to determine a unique optimal solution and ultimately identify the best optimization speed. The ideal point principle is shown in Equation (11), where the distance of each point in the Pareto optimal solution set to the ideal point with the minimum means of the two objectives is calculated, and the point with the shortest distance is chosen as the final optimal solution.
In Equation (11), is the distance between the Pareto frontier point and the ideal point; is the objective function value for objective 1 in the Pareto optimal solution set; is the minimum value of objective 1; is the objective function value for objective 2 in the Pareto optimal solution set; and is the minimum value of objective 2.
4. Results
The previous section introduced the steps for improving the NSGA-2 algorithm. The key to optimization lies in adjusting the objective weights and under different sea conditions, which significantly impacts the ship’s speed and fuel consumption. Below, sea conditions 1 and 6 are used as examples to study the effect of different weight coefficients on the ship’s speed and fuel consumption per nautical mile.
Sea condition 1 is characterized by an average wave height of 2.48 m, an average wind speed of 6.92 m/s, and a relative wind direction of 320°. For a Handymax bulk carrier, navigation under this condition is relatively safe, so more emphasis is placed on fuel consumption per nautical mile. The value of
is set from 0.5 to 1 in intervals of 0.01, corresponding to a value range of 0.5 to 0 for
. As shown in
Figure 8, the maximum allowable speed for the ship under sea condition 1 is 13.5 knots, and the minimum speed is 10 knots. The speed increases as the fuel consumption weight
decreases and the speed weight
increases. In
Figure 8, it can be seen that fuel consumption under sea condition 1 initially increases with an increase in
, reaching a maximum at
= 0.65, and then decreases, reaching a minimum near
= 0.85. Combining
Figure 8 and
Figure 9, it is evident that simply reducing the speed does not necessarily decrease the main engine’s fuel consumption. The relationship between speed and fuel consumption lies between a quadratic and cubic function, closer to a cubic function, and is still consistent with the classical approximate relationship between fuel consumption and speed.
Under sea condition 6, with an average wave height of 5.3 m, an average wind speed of 13.67 m/s, and a relative wind direction of 45°, navigational safety is paramount. If head or oblique waves are encountered on the route, the relative speed between the waves and the ship is high, causing significant impact on the hull. In severe cases, this can result in large rolling motions, significant deck wetness, bottom slamming, and propeller racing. Therefore, the speed must be reasonably reduced while considering the highest possible speed to quickly exit the severe wave area. In this case, the value of
ranges from 0.5 to 0, and
ranges from 0.5 to 1. From
Figure 10 and
Figure 11, it can be seen that under a wave height of 5 m, fuel consumption decreases with an increase in
, reaching a minimum at
= 0.88, and then increases. However, considering both speed and fuel consumption, it can be observed that when
ranges from 0.5 to 0.7, corresponding to speeds of 5.5 to 6.5 knots, the relationship between fuel consumption and speed is not the classic approximate cubic function but rather an approximate linear relationship. When
ranges from 0.7 to 1, it becomes more like a quadratic function relationship.
By determining the optimal objective weights,
and
, under different sea conditions, the optimal speed that balances both the navigation time and main engine’s fuel consumption was found, addressing the shortcomings of the NSGA-2 algorithm. Finally, the weights
and
for sea conditions 1 to 6 are shown in
Table 8.
5. Discussion
By comparing the pre-optimization ship and engine data in
Table 9 with the post-optimization data in
Table 10, the following can be observed: Firstly, from the perspective of RPM, the optimized speeds for all sea conditions are concentrated in the range of 96–114 r/m compared to the pre-optimization range of 89–110 r/m, showing a more unified and stable characteristic post-optimization. Taking sea condition 1 as an example, the engine speed before optimization was 89 rpm, which is 77% of the rated speed. After optimization, it increased to 96 rpm, which is 83% of the rated speed, enabling the engine to operate at its optimal condition for this sea state with an optimization rate of 10%. The optimal engine speeds for sea conditions 1 to 3 are 83% of the rated speed; for sea condition 4, it is 86%; and for the more severe sea conditions, 5 and 6, the optimal engine speeds are 97% and 98% of the rated speed, respectively. This concentrated and unified speed setting helps improve the operational efficiency and reliability of the main engine.
Secondly, in terms of speed, the optimized data show a significant increase in speed under most sea conditions. For example, the speed for sea condition 1 increased from 11 knots pre-optimization to 12.4 knots post-optimization, representing an improvement of 10.9%, and the speed for sea condition 2 increased from 10.5 knots to 11.8 knots, representing an improvement of 12.3%. Even under low-speed conditions (such as sea condition 6), the speed increased from 6 knots to 7 knots, representing an improvement of 16.7%. These improvements not only help reduce the navigation time but also reduce the operational costs of the ship. The most significant change is in fuel consumption, which significantly decreased under all sea conditions post-optimization. The most notable optimization is seen in sea condition 6, where fuel consumption decreased from 137.41 kg/nm pre-optimization to 107.41 kg/nm post-optimization, representing a reduction of 21.9%. The lowest fuel consumption in sea condition 2 also decreased from 85.80 kg/nm to 74.80 kg/nm, representing a 12.8% improvement, as shown in
Figure 12. This indicates that by optimizing the main engine’s speed, the engine can operate at the current optimal condition, significantly improving fuel efficiency and reducing fuel consumption per nautical mile, effectively lowering operational costs and environmental pollution. Optimizing the main engine’s speed to optimize the ship’s speed significantly improved fuel efficiency and navigation time. The optimized data indicate that without reducing the speed, fuel consumption was significantly reduced, and the navigation time was shortened. This optimization method is of great importance for enhancing the economic efficiency and environmental performance of ships.
6. Conclusions
This study investigated the impact of the main engine’s speed variations on the fuel consumption of bulk carriers under severe sea conditions from the perspective of ship–engine–propeller matching using a 55,000-ton Handymax bulk carrier as the research subject. A grey-box fuel consumption model was used, with the white-box model being based on the principle of ship–engine–propeller matching. The wave-added resistance was determined using the method developed by Liu and Papanikolaou, and the wind resistance was calculated using the wind resistance coefficients provided by ITTC. The main engine parameters were determined using a combination of navigation data and bench tests, and the propeller model was determined using the regression coefficient method. Wind and wave data, as well as navigation data, were obtained from ship sensors, AIS, the ERA5 database, and ship noon reports. The K-Medoids clustering algorithm was used to classify severe sea conditions. The Random Forest Regression algorithm was used to construct the black-box energy consumption model with the main engine’s speed as the variable. The fuel consumption per nautical mile and speed were the optimization objectives for the NSGA-2 algorithm. Finally, the white-box and black-box models were combined to establish the grey-box fuel consumption model. The following conclusions were drawn:
Under severe sea conditions, reasonably adjusting the engine speed to ensure the engine operates at its optimal condition can reduce fuel consumption per nautical mile. However, excessively reducing the engine speed can lead to increased fuel consumption and decreased ship speed.
Under severe sea conditions, increasing the main engine’s speed while ensuring navigational safety and without reducing speed can significantly reduce fuel consumption per nautical mile. In sea conditions with a wave height of 5.30 m and a wind speed of 13.67 m/s, increasing the speed from 6 knots to 7 knots not only increased the speed by 16.7% but also reduced fuel consumption per nautical mile by 21.9%.
Under severe sea conditions, the classic cubic relationship between speed and fuel consumption may not hold. In conditions with a wave height of 5.30 m and a wind speed of 13.67 m/s, the relationship between speed and fuel consumption was a combination of an approximate quadratic and linear relationship.
This study investigated the optimization of ship engine energy consumption in rough sea conditions, demonstrating the importance of adjusting the engine speed according to sea conditions. By altering the engine speed to operate at its optimal condition, the per nautical mile fuel consumption in adverse sea conditions can be significantly reduced, as well as the sailing time through these areas. Additionally, a preliminary study on the relationship between speed and fuel consumption in rough sea conditions was conducted. Future research will collect more navigation and sea condition data, further classify sea conditions, and consider the operational optimization of auxiliary engines to achieve overall energy consumption optimization for ships. Further studies will also investigate the relationship between engine fuel consumption and sea conditions in rough seas. Moreover, fuel consumption optimization for other major ship types besides bulk carriers will be conducted to ensure the model’s applicability to the three main ship types.