Optimal Ship Deployment and Sailing Speed under Alternative Fuels

: The European Union (EU) has implemented a sub-quota of 2% for renewable marine fuels to be utilized by vessels operating within its jurisdiction, effective starting from 2034. This progressive policy signiﬁes a signiﬁcant leap towards reducing carbon emissions and promoting sustainable development. However, it also presents notable challenges for shipping companies, particularly in terms of fuel costs. In order to support shipping companies in devising optimal strategies within the framework of this new policy, this study proposes a mixed-integer linear programming model. This model aims to determine the optimal decisions for fuel choice, sailing speed and the number of vessels on various routes. Furthermore, we showcase the adaptability of our model in response to ﬂuctuations in fuel prices, relevant vessel costs, and the total ﬂeet size of vessels. Through its innovative insights, this research provides invaluable guidance for optimal decision-making processes within shipping companies operating under the new EU policy, enabling them to minimize their total costs effectively.


Introduction
The global shipping industry plays a vital role in ensuring the movement of goods and commodities across the world [1][2][3].However, traditional fuels used in this sector, primarily fossil fuels, have raised significant concerns due to their adverse environmental impacts, including carbon emissions and air pollution [4].As societies worldwide strive to combat climate change and transition towards sustainable practices, the exploration and integration of renewable fuels have emerged as a pivotal solution for the maritime industry [5][6][7].Renewable fuels, also known as biofuels or alternative fuels, are derived from organic matter such as plants, algae, and waste materials.Unlike traditional fossil fuels, these sources are considered sustainable as they can be replenished through natural processes, reducing dependency on finite resources [8].The integration of renewable fuels into the maritime sector holds tremendous potential to greatly mitigate the environmental footprint of shipping operations while promoting sustainable development.Furthermore, the utilization of renewable energy sources aligns with international efforts to achieve the objectives of the Paris Agreement and the International Maritime Organization's (IMO) related emission reduction targets [9,10].Understanding the significance of renewable fuels in maritime transportation is significant for guiding policy decisions and driving sustainable practices within the shipping industry.
According to [11], the EU has established a sub-quota of 2% for renewable marine fuels to be utilized by vessels operating within its jurisdiction, effective from the year 2034.This means that starting from 2034, at least 2% of the fuel used by vessels during voyages within the EU must be derived from renewable sources.Note that for these voyages linking the EU and non-EU areas, half of the fuel consumption generated is attributed to the EU.Undoubtedly, this exerts a crucial impact on the optimal decision-making of shipping companies.
In this paper, we put forth mathematical frameworks aimed at facilitating the attainment of optimal decision-making in the pursuit of cost minimization for shipping companies, while simultaneously ensuring compliance with the renewable fuel policy recently introduced by the EU.To the utmost extent of our understanding, our study represents a groundbreaking inquiry, uniquely incorporating this freshly established EU policy.Specifically, our investigation addresses the ensuing research inquiries: 1. What are the optimal fuel choices within different types of areas that result in the minimal overall expenditure while adhering to the renewable fuel target stipulated by EU policy? 2. What are the optimal sailing speeds, taking into account both traditional and renewable fuels, within the EU and non-EU areas, to minimize the aggregate expenses incurred by shipping companies in accordance with the recently proposed EU renewable fuel regulations?3. What is the optimal number of deployed vessels in diverse shipping routes, including considerations for the chartering in or chartering out ships, that will culminate in the lowest total costs while simultaneously meeting the EU's 2% renewable fuel target?4. How do the optimal fuel choices within the EU and non-EU areas, sailing speeds associated with different fuel types, as well as the optimal number of vessels to be equipped in each shipping route, including the consideration of vessels chartering in or chartering out, vary in response to fluctuations in fuel prices, relevant vessel costs, and the total fleet size of vessels?
To tackle the four research inquiries mentioned above, we initially introduce a sophisticated mixed integer nonlinear optimization (MINLP) model characterized by its intricate nature and challenging problem-solving complexity.Subsequently, we employ advanced mathematical techniques to convert the MINLP model into a mixed-integer linear programming (MILP) formulation.This transformation enables the utilization of readily available optimization solvers for solving the MILP model.Lastly, we conduct a series of comprehensive experiments and thorough sensitivity analyses to assess the model's performance in response to various parameter variations.

Literature Review
We review two streams of literature closely related to our study: (i) the renewable fuel in shipping; (ii) the optimal decisions in shipping.

The Renewable Fuel in Shipping
In the level of policy, many countries and organizations have shed light on the transformative potential of renewable fuels and their role in shaping the future of sustainable shipping [12].In detail, at the IMO, there are ongoing deliberations regarding the implementation of Market-Based Measures (MBMs) to enhance the economic desirability of lowand zero-carbon fuels compared to fossil fuels.The proposed MBMs encompass a range of approaches, including the imposition of global levies on marine fuels, the establishment of an Emissions Trading System (ETS), and the exploration of other hybrid mechanisms [13].Moreover, the FuelEU Maritime Directive exemplifies a joint endeavor to decrease the greenhouse gas (GHG) of energy employed aboard maritime vessels, elucidating a profound ambition of achieving a 75% reduction by 2050 [14].This multipronged objective can be actualized by means of actively advocating and extensively embracing renewable and low-carbon fuels [14].In addition, a comprehensive overhaul of the European Energy Taxation Directive (EU ETD) looms on the horizon, enacting thresholds for the minimum taxation rates on bunker fuel.Noteworthy is the fact that fossil fuels bear the brun of the highest minimum tax rate, valiantly standing at EUR 10.75/GJ, while renewable fuels are subject to the lowest rate, valiantly striking at EUR 0.15/GJ [15].In addition, the revised renewable energy directive (RED II) propounds an exacting target of a minimum 13% reduction in GHG intensity within the transport sector by 2030 [15].Moreover, it establishes sub-targets, tailored specifically to advance biofuels and non-biological renewable fuels [16].
In the aspect of renewable fuel utilization, various studies have explored the potential of different alternatives.Hydrogen fuel, due to its notable efficiency and environmental advantages, stands as a paramount sustainable fuel option.In practical applications, methanol is commonly employed to generate hydrogen.The alcohol-hydrogen fuel mixed by a series of high-temperature catalytic reactions can be used in the shipping industry [17].Moreover, biodiesel emerges as a sustainable energy source that exhibits exceptional biodegradability and possesses low toxicity, making it an excellent alternative to fossil fuels across various sectors [18].Additionally, a comprehensive study [19] focused on evaluating the environmental and economic aspects of using Liquefied Natural Gas (LNG) as a ship fuel.The results from this extensive life-cycle analysis and cost assessment demonstrated that LNG generates lower greenhouse gas emissions, up to 28% less than heavy fuel oil, while producing slightly higher nitrogen oxide emissions.From the prospective of technology, LNG, Liquefied Petroleum Gas (LPG), and methanol are considered more mature technologies, along with biodiesel, hydrogen, and ammonia fuels, which exhibit greater potential for future development [10].However, the present utilization of a diverse array of low-carbon fuels presents certain drawbacks that impede their immediate substitution for conventional fossil fuels.Addressing this challenge entails the establishment of multi-period energy planning, which not only facilitates adaptation to demand fluctuations and evolving emission restrictions at distinct stages, but also aids in energy projection and future investment strategizing [20].What is more, from an economic standpoint, the dual fuel propulsion system stands out as the most viable and cost-efficient alternative for container vessels in the present.By ingeniously alternating between conventional fossil fuels and renewable energy sources, this system optimally adheres to forthcoming regulatory requirements and standards [21].

The Optimal Decisions in Shipping
In the realm of maritime transportation, the pursuit of attaining the optimal trajectory, sailing speed, fuel consumption, and other pertinent optimal decisions constitutes the focus in cost minimization.Within this context, the exploration of optimal decisions has been undertaken to curtail operational expenses [22][23][24][25].
Several studies have been conducted to optimize sailing speed and fuel consumption for the purpose of reducing shipping costs.For the relationship between fuel consumption and sailing speed, the prevailing consensus posits that there exists a cubic relationship between fuel consumption and sailing speed [26].Laporte et al. [27] conducted a comprehensive study on speed optimization problems within the liner transportation network, taking into account time window constraints.Arijit et al. [28] and other researchers extensively investigated the implementation of a slow streaming policy as a means to minimize fuel consumption and associated costs.The adoption of slower speeds has proven effective in reducing fuel consumption; however, it can also result in delivery delays.Moreover, reducing the sailing speed necessitates deploying more vessels for the service, thereby augmenting the operational expenses of the ships [29].In addressing this trade-off, He et al. [30] proposed a speed optimization problem where a set of speeds for each segment of a given route is determined to minimize costs while considering time windows and speed limits for each segment.Aydin et al. [31] conducted an assessment of speed optimization in liner shipping, incorporating stochastic port times into their model.Their objective was to minimize total fuel consumption while upholding schedule reliability.Li et al. [32] proposed an innovative approach that combines the optimization of sailing routes and speeds, taking into account the interrelation between them as well as environmental factors.
In the context of model formulation involving multiple optimal decision variables, Lu et al. [33] introduce a dual-objective optimization model for ship speed, with a focus on the influence of the ECA control area, berthing costs, and AMP systems.By employing the multi-objective PSO (Particle Swarm Optimization) algorithm, the study aimed to identify the Pareto solution set that simultaneously minimizes ship operating costs and carbon emissions.Subsequently, the TOPSIS (Technique for Order Preference by Similarity to Ideal Solution) algorithm was applied to evaluate and extract the optimal compromise solution from the Pareto set.Mahsa et al. [34] formulated a bi-objective programming model that addresses the joint optimization of ship scheduling and sailing speeds for a specific shipping service.This model takes into account factors such as non-identical stream flow speeds.Iris et al. [35] focused their research on the well-known berth allocation problem (BAP).They proposed a novel mathematical formulation that extends the traditional BAP to encompass multiple ports within a shipping network.Furthermore, their model implementation demonstrates that precise speed discretization can lead to significantly improved economic and environmental outcomes.Dulebenets et al. [36] introduced an MINLP model aimed at minimizing the overall cost of liner shipping routes.Meng et al. [37] evaluated an optimization model for fuel consumption and ship speed in liner transport, accounting for deviations from the scheduled speed.Lee et al. [38] investigated a dynamic planning model to determine fuel consumption under uncertain fuel costs.The trade-off between vessel quantity and speed was further elucidated in [39][40][41], emphasizing the significance of selecting the optimal number of ships for liner shipping.Sheng [40] investigated optimal vessel speed and fleet size considerations for services operating within emission control areas (ECAs).
In general, prior research efforts have predominantly focused on optimizing decisions related to sailing speed, fuel consumption, deployed vessels across different routes, and fleet size.These studies aim to provide more effective decision support for shipping companies by jointly considering multiple optimal decision variables.However, limited attention has been paid to the fuel selection and corresponding sailing speeds for different legs of the journey.Our research endeavors to address this gap by integrating the latest renewable fuel policies implemented by the EU with existing research findings.In this paper, we propose an MILP model, which aims to determine the optimal fuel selection, speed decisions, and the corresponding deployment quantity of vessels, simultaneously achieving cost minimization and promoting environmental conservation and sustainable development.

Research Contributions
1. Theoretical contributions.The present research addresses a significant research gap by focusing on the optimal selection of fuel, sailing speed, and the number of ships under the newly proposed EU policy.Notably, the existing literature has overlooked this specific aspect.To the best of our knowledge, this study represents the pioneering effort in establishing mathematical models aimed at minimizing the overall costs incurred by shipping companies while considering the implications of the new EU policy.The proposed approach employs an MILP model to determine the optimal decisions for shipping companies.By conducting rigorous experiments and sensitivity analyses, this study yields specific solutions while evaluating the impacts of various parameters.2. Practical contributions.This research contributes valuable practical insights into the development of optimal strategies for shipping companies to effectively minimize costs and ensure compliance with the new EU emissions policy.The obtained results possess practical implications for fostering sustainable growth within the shipping industry, facilitating its alignment with environmental regulations.Notably, the proposed mathematical model serves as a decision-making tool, providing shipping companies with a framework to navigate the challenges presented by the new EU emissions policy effectively.
The rest of the paper is organized as follows.Section 2 describes the research problem in detail and develops the mathematical model.Section 3 proposes solution methods for addressing the initial proposed model.Section 4 conducts experiments and sensitivity analysis.Finally, conclusions are drawn in Section 5.
The main notations used in this study are summarized in Table 1.The cost of chartering in a vessel K The total fleet size of vessels T The service frequency of each shipping route µ 1 The traditional fuel price per tonne µ 2 The renewable fuel price per tonne b m The total berthing time at all ports on liner shipping route m, m ∈ M L m i The total length of the shipping route m within the i type of areas The minimum sailing speed v max The maximum sailing speed J The integer used to discretize sailing speed, j = 0, 1, ...., J Function f (v 3 ) Fuel consumption rate at the sailing speed of v Decision variables x m The number of vessels to be deployed on each liner shipping route m, m ∈ M x out The number of chartering out vessels x in The number of chartering in vessels l tm i The sailing length in the i type of areas on each liner shipping route m with traditional fuel, i = 0, 1, 2, m ∈ M l rm i The sailing length in the i type of areas on each liner shipping route m with renewable fuel, i = 0, 1, 2, m ∈ M z tmj i Binary decision variable that equals 1 if vessels sail using traditional fuel on liner shipping route m in the i type of areas with speed v j and 0 otherwise, Binary decision variable that equals 1 if vessels sail using renewable fuel on liner shipping route m in the i type of areas with speed v j and 0 otherwise, i = 0, 1, 2, m ∈ M, j = 0, 1, ...., J

Problem Description and Model Development
In this study, we consider vessel company decisions on the sailing speeds within different areas and the choice of fuel, i.e., the choice of traditional fuel and renewable fuel.Our primary objective is to assist vessel companies in achieving cost savings while simultaneously meeting the EU's 2% renewable target requirement.
We consider a shipping network with multiple routes; M denotes the set of shipping routes, m ∈ M. The liner shipping route m ∈ M calls at several ports in Europe and Asia.The vessel company serves these routes by employing vessels equipped with dualfuel engines, which allows them the operation on either traditional fuel or renewable fuel (though the latter option incurs higher costs).There are three distinct types of voyages in the considered container liner service network: (1) voyages within non-EU areas; (2) voyages linking the EU and non-EU areas; (3) voyages within the EU areas.We use v tm 0 , v tm 1 , and v tm 2 to denote the sailing speed using traditional fuel in voyages within non-EU areas, linking the EU and non-EU areas, and within the EU areas, respectively.And v rm 0 , v rm 1 , and v rm 2 denote the sailing speed using renewable fuel in voyages within non-EU areas, linking the EU and non-EU areas, and within the EU areas, respectively.Referring to [42], the fuel consumption and sailing speed have a cubic relationship.Therefore, we use av 3 to denote vessel fuel consumption (note that v could be v tm 0 , v tm 1 , v tm 2 , v rm 0 , v rm 1 , and v rm 2 ).We use l tm 0 , l tm 1 , and l tm 2 to denote the total length of voyages using traditional fuel within non-EU areas, linking the EU and non-EU areas, and within the EU areas, respectively.And l rm 0 , l rm 1 , and l rm 2 denote the total length of voyages using renewable fuel within non-EU areas, linking the EU and non-EU areas, and within the EU areas, respectively.
The vessel company possesses a total fleet size of K vessels.Additionally, the company has the option to charter out any surplus vessels, generating additional income represented by c out .On the other hand, if the existing fleet of K vessels is insufficient to serve the shipping network, the company can charter in additional vessels at a fee denoted by c in .These chartering options provide the company with flexibility in managing its vessel resources and optimizing its operations based on the specific demands and requirements of the shipping network.We have c in > c out .The operating cost of vessels on each liner shipping route m is c m .
In order to minimize the overall total costs, the vessel company must make strategic decisions regarding the following key factors: 1.The number of vessels to be deployed on each liner shipping route m, denoted by x m ; 2. The number of chartering in or chartering out vessels, denoted by x in and x out , respectively; 3. The sailing speed on each liner shipping route m with different types of fuel, i.e., v tm 0 , v tm 1 , v tm 2 , v rm 0 , v rm 1 , and v rm 2 ; 4. The sailing length on each liner shipping route m with different types of fuel, i.e., l tm i and l rm i , i = 0, 1, 2. Additionally, the vessel company needs to meet the following constraints: 1.The EU's 2% renewable target requirement; 2. The service frequency of each shipping route m, denoted by T m .That is, the time for finishing a single trip should be within T m days; 3. The restriction of the number of vessels; 4. The restriction of vessel sailing speed.
The optimization model is formulated as follows: Objective function (1) consists of four parts.Firstly, ∑ m∈M c m x m calculates the total operation cost of vessels on m different liner shipping routes.Secondly, c in x in represents the cost of chartering in additional vessels when the existing fleet of K vessels is insufficient to serve the shipping network; on the contrary, the third part, c out x out , is the additional revenue of chartering out surplus vessels.As Objective function (1) calculates the minimum overall total costs, we subtract this additional income.Fourthly, ∑ m∈M ∑ 2 i=0 (µ 1 ) represents the total fuel costs, including traditional fuel cost and renewable fuel cost, where µ 1 and µ 2 denote the traditional fuel and renewable fuel price, respectively.Constraint (2) meets the requirement of at least 2% of the renewable fuel utilized by vessels during voyages within the EU area starting from 2034.Constraints (3) restrict the service frequency of each shipping route m.Constraint (4) restricts the total number of vessels on m shipping routes, only involving the existing vessels and additional vessels that charter in or charter out, which is subtracted.Constraints (5) regulate the total voyage in every area on each shipping route m, consisting of the length of voyage using renewable fuel and traditional fuel.Constraints (6) and Constraints (7) give the domain of v tm i and v rm i , respectively, indicating the maximum and minimum of v rm i and v rm i .Constraint (8), Constraint (9) and Constraint (10) regulate the number of deployed vessels in route m, and the number of additional vessels should be positive integers.Constraints (11) and Constraints (12) restrict the length of each leg in different types of areas on m routes, and they should be positive.For parameters, i = 0 indicates areas within the EU; i = 1 represents the areas linking the EU areas and non-areas; i = 2 denotes areas within non-EU.µ 1 (USD/tonne) and µ 2 (USD/tonne) denote the fuel price of traditional fuel and renewable fuel.According to practice, µ 1 < µ 2 .b m represents the total berthing time at all ports.L m i denotes the total length of the shipping route m of different areas.The decision variables include v tm i , v rm i , l tm i , and l rm i .

Solution Methods
Model (M1) is hard to solve due to the operation of multiplying and dividing the decision variables.We next develop methods to address the nonlinear terms and transform (M1) into an MILP model, which improves computational efficiency.
We discretize sailing speed v by 0.01 knot.We define and we set j = 0, 1, ..., J. Therefore, the sailing speed v can be discretized to v j : We further adopt binary decision variables to indicate which discretized sailing speed is chosen on each shipping route m with different types of fuel.To be more specific, z tmj i , i = 0, 1, 2, m ∈ M, and j = 0, ..., J denotes which sailing speed v j is chosen on shipping route m on area i using traditional fuel; z rmj i , i = 0, 1, 2, m ∈ M, and j = 0, ..., J denotes which sailing speed v j is chosen on shipping route m on area i using renewable fuel.Therefore, we use the new binary decision variables z tmj i and z rmj i to replace v tm i and v rm i , and (M1) can be transformed into the following (M2).
subject to 0.5 In Model (M2), we discretize sailing speed v to facilitate our research, improving computational efficiency.Specifically, we utilize v tmj i and v rmj i to represent the sailing speeds for different types of voyages using traditional fuel and renewable fuel, respectively.Furthermore, we set the value of W to L m i .Broadly speaking, Model (M2) encompasses three categories of decision variables.The first category involves integer decision variables x m , x in , and x out , which denote the number of vessels to be deployed on each liner shipping route m, as well as the quantities of vessels to be chartered in or out.The second category features binary decision variables z tmj i and z rmj i , comprising a total of 6mJ binary decision variables.When z tmj i = 1, the corresponding sailing speed v tmj i is selected for voyages utilizing traditional fuel, resulting in the adoption of traditional fuel consumption.On the other hand, when z rmj i = 1, the corresponding sailing speed v rmj i is chosen for voyages using renewable fuel.The renewable fuel and traditional fuel consumption are determined by expression av 3 .Lastly, the third category consists of continuous variables l tm i and l rm i , totaling 6m continuous variables.By transforming Model (M1) into Model (M2), we can enhance computational efficiency.
Model (M2) is still hard to solve because of terms l tm i z tmj i and l rm i z rmj i in Objective function (14) and Constraints (16).As mentioned earlier, l tm i , l rm i , z tmj i and z rmj i are decision variables, so the multiplications of l tm i and z tmj i , l rm i and z rmj i lead to the formation of nonlinear terms within the objective function and constraint conditions, thereby transforming the model into an MINLP problem, which is inherently challenging to solve.As a nonlinear problem can encompass multiple feasible regions or sets of similar values for the decision variables that satisfy all constraints, each feasible region may contain multiple "peaks" (in maximization problems) or "valleys" (in minimization problems).Determining which peak is the tallest or which valley is the deepest lacks a general approach.Additionally, there can exist spurious peaks or valleys referred to as "saddle points."Due to these possibilities, nonlinear optimization methods offer limited guarantees in terms of identifying the "true" optimal solution.As a nonlinear problem can encompass multiple feasible regions or sets of similar values for the decision variables that satisfy all constraints, each feasible region may contain multiple "peaks" (in maximization problems) or "valleys" (in minimization problems).Determining which peak is the tallest or which valley is the deepest lacks a general approach.Additionally, there can exist spurious peaks or valleys referred to as "saddle points".Due to these possibilities, nonlinear optimization methods offer limited guarantees in terms of identifying the "true" optimal solution.We define and The following inequalities hold: where Q m i represents an exceedingly large value that bounds the upper limit of l tm i and l rm i .We set Q m i to be equal to L m i .Typically, a smaller value of Q m i is preferred (while still ensuring the correctness of the model), as it usually leads to shorter computational time compared to a model with larger Q m i .Regarding Constraints (30), since the binary variable z tmj i can take values of either 0 or 1, when z tmj i equals 1, Constraints (30) can be expressed as l tmj i ≤ L m i , m ∈ M, j = 0, . . ., J, i = 0, 1, 2. This condition is evidently valid since l tm i is bounded by the total voyage length L m i .Conversely, when z tmj i equals 0, Constraints (30) become 0 ≤ l tm i , m ∈ M, j = 0, . . ., J, i = 0, 1, 2, which is also true.
Concerning Constraints (31), given the aforementioned definition, η , m ∈ M, j = 0, . . ., J, i = 0, 1, 2. Consequently, Q m i becomes the upper bound for l tm i , equating to the total leg length L m i .Regarding Constraints (32), we can discuss the case where z tmj i is equal to 1 or 0 separately.When z tmj i is equal to 1, Constraints (32) can be transformed into l tmj i ≥ l tm i , m ∈ M, j = 0, . . ., J, i = 0, 1, 2. It is evident that this inequality holds.On the other hand, when z tmj i is equal to 0, the inequality becomes 0 ≥ l tm i − Q m i , m ∈ M, j = 0, . . ., J, i = 0, 1, 2, implying that the minimum value of Q m i is l tmj i , consistent with Constraints (31).The validity of Constraints (33) is unquestionable, as both l tmj i and z tmj i are greater than or equal to 0. Similarly, Constraints (34)-(37) hold for the same reasons.
subject to 0.5 With the help of Constraints ( 30)-( 37) and discretization, the original optimization model is transformed into an MILP programming model, which can be solved by the off-the-shelf optimization solvers, such as CPLEX and Gurobi.To validate our models, we design a shipping network involving four shipping routes for the experiment, setting the parameters, e.g., c m , c in , c out , µ 1 , and µ 2 according to practice.

Experiment Settings
The experiments were run on a laptop computer equipped with 2.60 GHz of Intel Core i7 CPU and 16 GB of RAM, and Model (M3) was solved by Gurobi Optimizer 10.0.2 via Python API.

Selected Shipping Routes
We select four routes from Asia to northern Europe 1 to test the performance of Model (M3).Details are shown in Table 2.

Route ID
Port Rotation (City)

. Parameter Settings
We first set the values of parameters for drawing the basic results, and we conduct sensitivity analysis to examine the impacts of these parameters.
1.The operation cost c m .Referring to [43], we first set c m = USD 180,000 per week for a 5000-TEU (Twenty-foot Equivalent Unit) container ship.2. The fee of chartering in a vessel c in .Referring to [44], we set c in = USD 120,000 per week.3. The fee of chartering out a vessel c out .Referring to [44], we set c out = USD 100,000 per week.4. The traditional fuel price µ 1 .Referring to [45], we set µ 1 to be an average value of 600 (USD /tonne). 5.The renewable fuel price µ 2 .Referring to [45], we set µ 2 to be an average value of 1000 (USD /tonne).6.A company's total fleet size K. Referring to [46], we set K to be an average value of 60. 7. Referring to [43], We set v max = 18 knots and v min = 13 knots.8. Referring to [42], we set f (v 3 ) = 0.00043 × v 3 , a = 0.00043.

Basic Results
Based on the routes presented in Table 2 and the parameter settings, we conducted numerical experiments and obtained the results in Table 3.As outlined in Section 2, the decisions regarding vessel sailing speeds and fuel choices play a significant role.Therefore, we analyzed the sailing speeds with traditional and renewable fuels, as well as voyages in different types of areas.Specifically, during the course of the experiments, as we mentioned earlier, Model (M1) entails nonlinear terms and divisions between variables, rendering it a nonlinear model that poses significant challenges for solving.However, by discretizing the sailing speed and introducing variables z The optimal value of the objective function for Model (M3) is denoted as "OBJ," and distances are measured in nautical miles (nm).From Table 3, we observe that the number of deployed vessels is comparable across the four routes.However, Route 1 has slightly more ships (13) compared to the other routes, primarily due to its longer total voyage distance among the four routes.Additionally, the existing fleet of vessels is sufficient, leading to the chartering out of 11 ships.In general, vessels use traditional fuel in all types of areas, including those within the EU, those connecting the EU and non-EU regions, and non-EU areas.However, renewable fuel consumption is concentrated within the EU areas.As mentioned earlier, vessel companies must consider the EU's 2% renewable fuel target requirement.For voyages linking the EU and non-EU areas, half of the fuel consumption is attributed to the EU, while for voyages within the EU areas, all fuel consumption is attributed to the EU.To meet the 2% renewable fuel requirement at a lower cost, vessel companies prioritize using renewable fuel predominantly within the EU areas.This approach minimizes fuel consumption while still meeting the renewable target requirement, considering that renewable fuel prices are higher than traditional fuel prices.Additionally, a vessel exhibits a slower speed when using sustainable energy compared to its speed when utilizing traditional fuel.

Sensitivity Analysis
The unit traditional fuel price and renewable fuel price may change as more policies are issued to promote the use of renewable fuel.In addition, in the basic analysis, some important parameters, e.g., the total fleet size of vessels, the operating cost per ship and the revenue of chartering out a vessel, are set to be deterministic.However, these parameters often fluctuate in real life.Therefore, sensitivity analyses on these parameters are conducted to investigate the influences of these parameters on the operation decisions.In the sensitivity analysis, the parameters are divided into three sorts.The first one is the fuel price, including traditional fuel price and renewable fuel price; the second one is the relevant cost of the vessels, involving the operating cost of a vessel, the revenue of chartering out a vessel as well as the cost of chartering in a vessel; the third one is the total fleet size of vessels.

Impact of the Fuel Price
Given the implementation of the EU's 2% renewable fuel policy, the selection and pricing of different fuel types have significantly impacted optimal operational decisions of vessel companies.Consequently, in this section, we specifically explore the effects of traditional and renewable fuels on such decision-making processes.In our initial analysis, the unit price of traditional fuel (µ 1 ) is established at 600 dollars per tonne.In sensitivity analysis, we set µ 1 varying between 500 and 800 dollars per tonne.The computational results are reported in Table 4.
According to the findings presented in Table 4, the objective value exhibits an upward trend as the price of traditional fuel increases.Furthermore, an increased allocation of vessels is observed across the four shipping routes, accompanied by a decrease in sailing speeds as traditional fuel prices soar.Given the cubic relationship between fuel consumption and sailing speed, it follows that higher vessel speeds result in greater fuel consumption.Consequently, in response to the rising unit price of traditional fuel, shipping companies are inclined to reduce sailing speeds to mitigate traditional fuel consumption and thereby achieve cost savings in their operations.Additionally, if the unit price of traditional fuel (µ 1 ) reaches excessively high levels, shipping companies may opt to deploy additional vessels to maintain the desired weekly service frequency.
In basic analysis, the unit price of renewable fuel (µ 2 ) is set at 1000 dollars per tonne.However, as stated in [21], the growing emphasis on green and sustainable development has led to the formulation of favorable policies aimed at promoting the utilization of renewable fuel.This suggests the potential future decline in renewable fuel prices.Thus, we set the range of µ 2 to span from 800 to 1000 dollars per tonne.Computational findings are presented in Table 5.
Regarding renewable fuel prices, similar patterns emerge as with traditional fuel prices, albeit to a lesser extent.Specifically, the total cost increases as the price of renewable fuel rises.However, due to the minimal proportion of renewable fuel in the overall fuel consumption, the impact of changing renewable fuel prices is relatively minor compared to the effects observed with traditional fuel.Due to the minimal impacts of fluctuations in cost of chartering in a vessel (c in ) and the revenue of chartering out a vessel (c out ) (where c in resemble ship operation cost, resulting in an increase in total cost as c in rises, while the total cost decreases with an upturn in c out ), our analysis focuses solely on ship operation cost.
In the aforementioned analysis, the c m predetermined weekly operating cost for a vessel stands at 180,000 dollars.Nevertheless, the value of c m is subject to fluctuation due to the impact of various unpredictable factors and risks [47], such as the global outbreak of COVID-19 in 2020, which reportedly sparked an escalation in ship operating costs [48], or the gradual impact of technological advancements, which may potentially reduce these costs.Consequently, the range of c m is defined as 160,000 dollars to 240,000 dollars, with corresponding results meticulously documented in Table 6.
Analysis of Table 6 reveals a direct correlation between the increase in the operating cost of a vessel (c m ) and the corresponding rise in the objective value.This signifies that as the total operating costs for vessels surge, the overall objective value exhibits an upward trajectory.Furthermore, as the value of c m intensifies, a prudent approach to saving on operating costs is witnessed through a reduction in the number of vessels deployed across the four routes.In the initial analysis, the total fleet size of vessels, denoted as K, is assumed to have an average value of 60.However, it is important to note that the total fleet size varies among different companies and at different stages of development.Consequently, the value of K varies according to the scale of each company.Therefore, it becomes imperative to conduct a sensitivity analysis of the total fleet size of vessels.In this experiment, we consider a range of fleet sizes, from 40 to 70 ships.It is worth mentioning that the vessel type under investigation is a 5000-TEU (Twenty-foot Equivalent Unit) container ship, which typically does not exceed 70 in number.
Analyzing Table 7, we observe that as the total fleet size of vessels increases, the objective value decreases.This is attributed to a shift in the vessel situation of companies, transforming it from having inadequate vessels to having surplus ships.Thus, additional revenue is generated through chartering out surplus vessels.Moreover, it is notable that the allocation of ships to the four shipping routes remains unchanged regardless of variations in the total number of ships.The shipping companies only need to determine the number of ships to be chartered out or chartered in based on the balance between the sum of ships allocated to the four routes and the company's existing total number of ships.

Conclusions
This research investigates the shipping company's optimal strategy regarding fuel selection, sailing speed, ship deployment, and the number of ships chartered in or chartered out, considering the EU's proposed new policy.Initially, we present an innovative MINLP model, subsequently converting it into an MILP model through the application of advanced mathematical techniques.Through rigorous experiments, our proposed model demonstrates its efficacy, leading to the following conclusions: (i) Due to the lower prices compared to those of renewable fuels, traditional fuels are employed in all three types of regions, while renewable fuels are selectively used only in specific segments within the EU area.(ii) In terms of sailing speed, given the cubic relationship between speed and fuel consumption, as well as the price disparities between the two fuel types, ships tend to operate at higher speeds when using traditional fuels, thus reducing operational costs compared to utilizing sustainable fuels.(iii) The number of vessels employed and the mileage of routes are interconnected, with decisions regarding vessel chartered in or chartered out contingent upon striking a balance between the total number of vessels employed to each route and the overall fleet size of vessels.Furthermore, we analyze the influence of variations in fuel prices, relevant costs of vessels, and the total fleet size of vessels on optimal decisions.In general, increases in fuel prices, vessel-related costs, and fleet size result in rising total costs.Specifically, when the prices of both fuel types increase, shipping companies, seeking cost reduction, tend to lower sailing speeds to minimize fuel consumption, while also employing a higher number of vessels to satisfy weekly service frequencies.Moreover, when the operational costs of vessels increase, the number of vessels employed decreases.It is worth noting that changes in the total fleet size of vessels do not affect the number of vessels employed to individual routes but rather impact decisions regarding vessel leasing and chartering.
This investigation offers valuable insights into shipping company strategies under the new policy.Overall, our study contributes to the understanding of how shipping companies can make informed decisions in response to the EU's new policy.By providing efficient solution methods and examining the sensitivity of optimal decisions to various factors, our research offers practical guidance for navigating the challenges and opportunities presented by the policy.
Despite the significant contributions of this study, there exist certain limitations that require further investigation.
Firstly, the proposed model assumes static conditions and does not account for dynamic factors such as changing market trends, weather conditions, or evolving regulations.Incorporating these dynamic elements would be valuable for a more comprehensive analysis.
Secondly, our research focuses primarily on the shipping company's perspective and optimal decision-making.However, future studies could consider the broader EU's policy on the maritime industry as a whole, including the implications for sustainability, environmental protection, and the overall supply chain efficiency.
Additionally, our study assumes perfect information availability and precise parameter estimations.In reality, uncertain data and imperfect information are common challenges.Future research endeavors could explore robust optimization techniques or apply stochastic programming to address these uncertainties and enhance the reliability of decisionmaking processes.
In conclusion, this study provides valuable insights into the optimal strategies for shipping companies under the new EU policy.However, addressing the aforementioned limitations and pursuing further research in the suggested directions will yield a more comprehensive understanding of the complex interplay between fuel selection, sailing speed, and fleet management, contributing to the advancement of sustainable and efficient maritime operations.

0
The sailing speed on liner shipping route m within non-EU areas with traditional fuel, m ∈ M v tm 1 The sailing speed on liner shipping route m within linking the EU and non-EU areas with traditional fuel, m ∈ M v tm 2 The sailing speed on liner shipping route m within EU areas with traditional fuel, m ∈ M v rm 0 The sailing speed on liner shipping route m within non-EU areas with renewable fuel, m ∈ M v rm 1 The sailing speed on liner shipping route m within linking the EU and non-EU areas with renewable fuel, m ∈ M v rm 2 The sailing speed on liner shipping route m within EU areas with renewable fuel, m ∈ M v min transformed the MINLP model (M1) into the MILP model (M3), facilitating the determination of optimal values during the experiments process.Hence, we employ Model (M3) to ascertain the optimal decisions for sailing speed, travel distance, fuel selection, and other relevant factors.

Table 2 .
Summary of shipping routes.

Table 4 .
Impact of the unit price of traditional fuel.

Table 5 .
Impact of the unit price of renewable fuel.

Table 6 .
Impact of the cost of operating a vessel.

Table 7 .
Impact of the total fleet size of vessels.