An Optimization Model for Tramp Ship Scheduling considering Time Window and Seaport Operation Delay Factors

)e quantity of electrical coal transported through the tramp shipping network is increasing due to the high demands. )is trend has increased the scheduling difficulty combined with the underdevelopment of the private thermal power plant port. )e high coal consumption and low port storage capacity requires the scheduling of the tramp ship to be on a strict time window to ensure the continuous operation of the thermal power plant.)e low port unloading capacity often leads to the port congestion and delay of the unloading operation. )is paper develops a mixed-integer-programming model for the optimization of the tramp ship scheduling to reduce the total operation cost, including the transportation cost and the unloading waiting cost, and the branchand-price algorithm is adopted to solve this large-scale model. )e model and algorithm are tested with historical operation data from the thermal power plant in the southern coastal areas of China. )e optimized scheme significantly reduces the total operation cost by reducing the unloading waiting time and the number of active vessels in certain periods. )e results also demonstrate the algorithm improvement in the aspects of the optimization quality and efficiency comparing with the heuristic solution.


Introduction
e supply of electric coal (coal used to generate electricity) in the southern coastal areas of China follows the basic pattern of transporting coal from north to south [1]. Electric coal is mainly transported through waterways from the main producing regions of northern China to the numerous thermal power plants located in the southern coastal areas [2]. In this vast and complex supply system, large-scale power generation companies (who typically operate approximately 10-20 power plants) play the role of transport organizers. ese companies are responsible for not only organizing the electric coal transport fleet through the international ship charter market but also setting the scheduling and operation plan of the fleet. Specifically, these companies make decisions regarding the electric coal transportation plan, ship routing, and scheduling. In the past, the experience and subjective judgment of decision makers (heads of the transportation departments of the power generation companies) have been the primary basis for formulating the electric coal transportation plan [3]. Due to the limited amount of transport and the small size of the fleet, this practice did not present severe flaws at the time and was able to ensure that transportation plans based on the human experience achieved a delicate balance between efficiency and cost.
However, in recent years, the electric coal transportation system has undergone a disruptive change with further increases in the amount of thermal power generation. To meet the growing electricity demand, since 2010, the Chinese government has stepped up efforts to build thermal power units in the southern coastal areas. In turn, the volume of coal and the size of coal transportation fleets have also increased rapidly. In 2010, the company transported only 30,000 tons of electric coal per day and used only five vessels every day on average [4]. In contrast, in 2016, this company transported approximately 70,000 tons per day and used over 20 vessels every day [4].
With the rapid growth of transportation needs, several characteristics of the near-sea electric coal transportation have drawn the attention of many researchers. e electric coal transportation is strictly based on pre-made plans, which includes the loading and unloading task, fleet management, and fleet scheduling. Because of the high consumption volume and the low storage capacity of the small thermal power plant port, the transportation tasks have a high requirement on the time sensitivity to ensure the continuous supply of the electrical coal for the thermal plant. e fleets consist of different types and number of ships with different transportation capacities and often supported with the renting of the additional vessel for temporary high capacity demands [5,6]. In electric coal transportation systems, the operating capacity of these ports is typically extremely limited because unloading ports are mainly constructed and operated by thermal power plants. ese ports can only provide unloading services for one ship at a time, often resulting in long queues of ships waiting to be unloaded (i.e., the port congestion problem). Lastly, the railway has always been one of the essential modes of transportation in the southern Chinese coastal electric coal transportation system [7]. Since the Chinese railway is the dominant transportation mode, the establishment of the modern integrated transportation hubs has also enabled the land-and-water coordinated transport [8,9]. ree factors should be considered for the formulation of the electric coal transportation plan (referred to as the electric coal ship-scheduling problem, ECSSP). e first factor is the time-window constraints of electric coal transportation. Since the storage capacity of the unloading port is quite limited compared to the electric coal consumption, even the minor delay of the transportation could cause the shutdown of the powerplant, which in turn leads to economic and social loss. us, it is crucial for the shipment to be completed within their time windows. e second factor is the congestion of the unloading port. e congestion could cause port operation delays, which could not be neglected due to the short voyage duration and high operation frequency [10]. Because of the strict time window and the interdependent scheduling for the whole transportation operation, the delay could have a chain reaction effect, which would decrease the effectiveness of the entire transportation network. e third factor is the synergy between water transportation and railway transportation. e complementarity and cooperation of the railway with waterway transportation are necessary guarantees to ensure the stability of electric coal transportation. e abovementioned three factors make the current ECSSP highly complex. Decision makers should not only consider the constraint of the shipping time windows on scheduling and the ship unloading processes but also the reasonable cooperative application of water transportation and railway transportation.
us, it is more complicated than the traditional tramp ship routing and scheduling problem and deserves further study [11]. erefore, how to design a reasonable electric coal transportation plan has become a significant issue for major power generation companies [3].
Based on the unique characteristics of the electric coal transportation system in the southern coastal areas of China, this paper addresses the ECSSP under the background of hard time-window constraints, unloading port congestion, and synergy between water transportation and railway transportation. e main contributions of this paper are as follows: (1) A new tramp ship routing and scheduling optimization method considering port congestion is proposed.
(2) A mixed-integer-programming model is developed to optimize the shipping routing, the corresponding scheduling, and the share rate for water and railway transportation simultaneously (e.g., the beginning time point of the transportation task and the necessary waiting times at unloading ports). (3) A new algorithm based on branch-and-price is designed, which can solve ECSSP efficiently.
is paper is organized as follows. In Section 2, we review the relevant literature, clarify the field of ECSSP research, and analyze the difference and connection between the ECSSP and the corresponding classic problems. In Section 3, a detailed description of ECSSP is provided. A mixed-integer-programming model is developed in Section 4. e algorithm of this model is designed in Section . In Section 6, the model and the algorithm are tested with historical data. e performance of the algorithm is presented, and the practical value of the model and algorithm is discussed based. Conclusions are provided in Section 7.

Literature Review
According to Ronen's definition, the ECSSP belongs to the tramp ship routing and scheduling problem (TRSP) [12]. Systematic reviews of the literature on the TRSP have been provided by Christiansen et al. [11,13], Ronen [12,14], Lane, et al. [15], and Meng, et al. [16]. Christiansen and Fagerholt [17] noted that approximately 40% of the TRSP models use the set partitioning (SP) model framework, which can characterize the TRSP effectively. e core idea of the SP model framework is to translate the TRSP into a classic set partitioning problem by assuming that all the feasible shipscheduling plans are preknown. A feasible prerequisite for this framework is that decision makers must be able to provide an estimation of the feasible ship-scheduling plans set in advance.
For the TRSP, the early representative studies include the following research. Brown et al. [18] conducted a thorough study of the crude oil transportation problem. In this study, carriers must ship crude oil from the Middle East to Europe and North America with hard time-window constraints. For a known fleet, the decision maker must design a crude oil transport scheduling for each vessel. In response to this problem, Brown et al. [18] constructed an integerprogramming model based on the SP framework. Due to the relatively small size of the problem, the feasible shipscheduling plans set can be obtained by enumeration, which means that the model can be solved efficiently. erefore, we refer to this kind of method the enumeration SP algorithm (ESPA). Bausch et al. [19] and Perakis and Bremer [20] also used the ESPA to solve the TRSP in their respective studies. Christiansen and Fagerholt [17] discussed an extension of the TRSP where decision makers must consider the impact of port nonoperating periods and climatic conditions on the ship-scheduling plan and also built a model based on SP and solved it using the ESPA.
For the port congestion problem, while many studies indicate it is a significant and complex issue, only limited solutions are proposed. Jafari [21] studied the causes of the lag and halt in port operation, including document incompleteness, shortage of ship and truck, the unpreparedness of goods owners, and labor issues. ese factors could have significant effects on port operation efficiency. Rahman et al. [22] also study the factors that caused the port operation delay and found that a variety of aspects of ports, ships, cargo owners, and environmental uncertainty all contribute to the port congestion. eir research suggested that the subfactor of foul weather and tide prediction is the most significant factor for port operation delay, which is also the main impact factor for the fleet scheduling problem. For electrical coal transportation, the unloading port for the corresponding thermal power plant is relatively small and only has 1 or 2 berths with low port storage capacity. Once something unexpected happens, subsequent ships can only wait for the available berths for the loading and unloading operations. us, the port congestion could cause a significant delay for the following operation. To solve this complex issue, several studies have attempted to provide effective solutions. For example, Kavirathna et al. [23] proposed the vessel transportation policy for the terminals' cooperative management to increase the overall efficiency of the port. e berth scheduling problem is also studied to achieve the minimized total port operation delay [24,25]. e problem of port operation delay can also be regarded as the unbalanced distribution of logistic resources in time and space. Some scholars adopted the method of state-space-time (SST) network to promote the collaboration of the logistic network and optimize the distribution of logistic resources in time and space. Wang et al. [26] adopted the state-space-time (SST) network to solve the problem of collaborative two-echelon multicenter vehicle routing optimization. By promoting collaboration among service providers in the logistics network, logistic costs are reduced, transportation efficiency is improved, and transportation emissions are reduced [26]. e green logistics route planning problem with ecopackages could also be optimized by identifying the best location strategy with the state-space-time (SST) network [27]. While those methods show great potential, they are not suitable for the Chinese electric coal transportation issues for the relatively small unloading ports with only limited unloading capacity. erefore, from the perspective of fleet management, introducing the operation delay time into the TRSP could be an innovative solution for the port congestion issue.
In recent years, scholars have reached a consensus on the modeling of the TRSP [11], and the standard TRSP model has been proposed. However, because the fleet sizes mentioned by these TRSP issues are continuously expanding and the demands for transport are diversified, it has become increasingly difficult to enumerate all feasible ship-scheduling plans. e traditional ESPA-based enumeration method is no longer applicable, and how to solve the TRSP efficiently has become the research focus. According to the standard TRSP model, Brønmo et al. [28] proposed a multistage local search heuristic algorithm and compared the advantages and disadvantages of the new algorithm through four examples. is study proved that the new algorithm could effectively solve the large-scale TRSP within a feasible timeframe. Kobayashi and Kubo [29] proposed a column generation solution algorithm (CGSPA) for the standard TRSP model, in which the shortest path algorithm considering the constraints of time windows were incorporated. Fagerholt et al. [30] studied a variant of the TRSP in which there is a particular mechanism for addressing a shortage in transport capacity.
is mechanism allows decision makers to address the shortage of fleet capacity by temporarily renting vessels. However, the chartered vessels must complete at least two transport missions before releasing. For this problem, they proposed an integerprogramming model and a tabu search algorithm (denoted as TSPA) based on the work of Korsvik et al. [31]. Norstad et al. [32] considered the factor of ship speed while studying the TRSP and proposed a local heuristic search algorithm (named multistart) to solve this problem. Hemmati et al. [10] proposed a two-phase hybrid heuristic algorithm for the maritime inventory and routing problem that can also be used to solve the TRSP. Wang et al. [33] also adopted the integerprogramming-model-based approach for searching the reasonable solutions for each vehicle in the distribution center network, which provides the best sequential transportation coalition and efficient model performance. Similar studies include the work of Sopot and Gribkovskaia [34], Gansterer et al. [35], and Cuesta et al. [36].
A comprehensive analysis of the existing relevant literature illustrates that the majority of studies assume that the ships in the fleet do not affect each other during the transportation task. However, in the ECSSP proposed in this paper, the operational statuses of each ship are interrelated. Due to the operation delay caused by port congestion, the deviation of the transportation plan of any ship may have unexpected effects on the transport of other ships. erefore, the problem proposed in this paper represents a special type of TRSP, in which the scheduling plans of the vehicles (ships) interfere with and affect each other. Solving the ECSSP could enrich the research on the routing and scheduling optimization method with the consideration of the port congestion. e mixed-integer-programming model proposed in this research considered characteristics such as the shipping routing, related scheduling, and land-and-water coordinated transport and achieved simultaneous optimization. is paper also proposed a branchand-price-based algorithm for the effective solving of ECSSP, which would be beneficial to solve other complex models.

Problem Description
e China Southern Coastal Electric Coal Transportation System (CSTS) studied in this paper is a complex transportation system involving multiple links, multiple constraints, and multiple transportation modes.
e system consists of the electric coal supply port (referred to as the loading port), electric coal unloading port (referred to as the Journal of Advanced Transportation unloading port), water transport subsystem, and rail transport subsystem. Electric coal must be shipped from the loading port to reach a specific unloading port by the waterway or rail subsystem. e CSTS has five unique features that cannot be ignored.
First, while the spatial distribution of loading ports is concentrated, the spatial distribution of the unloading ports is scattered; thus, the water transport network is star shaped. As shown in Figure 1, in the CSTS, the loading ports are densely distributed in a small area (Supply Port Group) in the western region of the Bohai Sea. ese ports typically have sufficient transport capacity with the major coalproducing regions and, thus, have a steady coal supply capacity. In addition, due to the dense distribution of loading ports, the western Bohai Sea region has an extremely high port service capacity and can meet the demand for electric coal loading operations of any volume at almost any time. In contrast, because most of the unloading ports are built near power plants, they are widely clustered in the coastal areas of South China. In terms of the network structure, if we abstractly view the supply port group as a virtual loading port, then the network can be considered as a star-shaped transport network, with the center as the virtual loading port and the peripheral nodes as the unloading ports.
Second, electric coal transportation has strict time and capacity requirements. In practice, due to the limitation of coal storage capacity, thermal power plants must regularly replenish coal inventories and prepare a detailed task list of electric coal purchases and transportation in advance.
e task list provides details of each transport mission, including the type of coal to be transported, the volume of coal, the loading port, the unloading port, and the time window for loading and unloading. e time window is the time duration between the earliest possible loading time and the latest possible uploading time of the ship loading and unloading operation according to the berth operation schedule and ship operation schedule. Because the implementation status of the task list directly determines whether the relevant thermal power plants can operate stably, decision makers often demand that the electric coal transportation system strictly fulfills the listed transportation tasks, even without cost considerations. erefore, the electric coal transport system must be able to provide adequate capacity to ensure that the electric coal is shipped and delivered within the prescribed time window.
ird, electric coal transportation ships have been switched between the status of commissioning and storage. In the CSTS, decision makers usually need to manage a selfowned fleet of 20 to 30 ships of various types. In order to reduce the operating costs of the fleet, decision makers must implement temporary decommissioning of some vessels and choose when to use them again based on the demand for electric coal transportation. erefore, in ECSSP, decision makers must not only consider the ship's scheduling plan but also must plan the operation status (commissioning and decommissioning) of each ship. Because different types of ships have significantly different decommissioning costs and operating costs, how to formulate an appropriate operation plan based on the needs of electric coal transportation has also become an important issue.
Fourth, the transportation pattern of single demand, single mode, and single voyage is adopted.
is transportation pattern is the standard ship-scheduling rules for thermal coal transportation in production and transportation practice with significant advantages. Since 3%-8% cargo loss will occur for the transportation mode of the multiple bulk cargo, the adopted transportation could efficiently reduce the cargo loss for the frequent electric coal transportations in the network. Besides, transporting multiple cargoes in a single voyage or transporting the cargoes to multiple unloading nodes could bring significant difficulties for the statistical management of the company. us, the transportation pattern of single demand, single mode, and single voyage is adopted in the Coastal Electric Coal Transportation System (CSTS). In the CSTS, this pattern is adopted to facilitate the settlement of freight, avoid mutual pollution between different types of coal, and reduce unnecessary coal losses due to transport. In this pattern, only one mode of transportation can be used for certain transportation tasks. e transportation task cannot be split in any way, and a ship can complete at most one transportation task in a single voyage. is pattern forces decision makers to make a careful choice between water transportation and rail transportation when setting up a fleet to match the ship's capacity with transportation demands. How to make efficient use of the transportation resources of the system has become one of the critical issues that decision makers must solve on the premise of ensuring that transportation tasks are completed on time.
Fifth, the factors involved in transport organization are complex. is feature is mainly reflected in the following two interactive relationships: the interaction between rail transportation and waterway transportation. e interaction relationship selected for this paper is because of the following reasons. Even though the waterway transportation has the advantage of lower transportation batch cost, thermal coal transportation should consider other factors rather than a simple economic factor. e supply of thermal coal must guarantee the continuous operation of thermal power plants as the premise. However, the disruption of the coal transportation could cause the shutdown of the thermal power plant due to the low port storage capacity of the thermal coal port and the high coal consumption rate. As discussed above, the electoral coal transportation by waterway could suffer from the port congestion issue, whereas rail transportation can ignore uncertain factors such as the port congestion.
us, it is necessary to consider the combination of water transportation and other transportation modes to complete the thermal coal transportation to ensure that thermal coal can be supplied to thermal power plants in time. In the CSTS, the water transportation subsystem is designed, managed, and operated by the power generation group (i.e., the decision maker). Decision makers will form a fleet in advance according to the electricity production plan and then formulate and implement a dispatch plan for electric coal ships based on the transportation task list and the operational capabilities of each port. e railway transportation subsystem is operated by the government department of railway transportation. e decision makers obtain the necessary transportation services through purchase. e advantages of the waterway transportation subsystem are large transportation capacity and low unit transportation cost; however, the transportation organization is very complicated due to the impact of the port congestion.
e advantage of the railway transportation subsystem is that the transportation capacity is abundant, and the transportation process is reliable. e decision maker does not need to consider the transportation scheduling problem, nor is it affected by the factors of port congestion, but the transportation cost is high. erefore, how to make full use of the complementary advantages of the two transportation modes to improve transportation efficiency and reduce transportation costs is another crucial issue that decision makers must consider.
e interactions between ship-scheduling plans amplify the impact of port congestion. In the traditional TRSP problem, the scheduling process of a single ship is relatively independent, and the operation of each ship does not affect the others. However, in the electric coal ship-scheduling problem, the port congestion has a chain reaction effect on the operation processes of all ships in the plan period, which causes difficulty in the electrical coal supply for the subsequent thermal power plant. For example, when ship A encounters congestion in port i, this congestion may delay its arrival at subsequent port j, which in turn affects ships B and C that use port j as the unloading port. e impact of congestion in a certain port, with the unloading port as a medium, creates a continuous and complex chain reaction causing more significant operation delays throughout the fleet and planning period and ultimately disrupts the scheduling of every ship within the system. Under this mechanism, the improper scheduling of any one vessel may disrupt the operation of the entire system through the amplification of a chain reaction. erefore, the scheduling problem of ECSSP with the consideration of operation delay is far more complicated than that of the traditional TRSP problem. Decision makers must not only focus on the conventional decision-making problems, such as the matching of ships and transportation tasks and the order of ships to complete transportation tasks, but also accurately consider when each task should begin, when it is completed, and the waiting period at each port, thus reducing unnecessary economic losses.
Based on the abovementioned facts, we present the following description of the ECSSP proposed in this paper. Given are (1) the planning period, (2) the task list of electric Journal of Advanced Transportation coal transportation (including transportation volume, loading and unloading port, and loading and unloading time window), and (3) the information regarding renting ships (e.g., daily rent, earliest available time, dead weight, and average speed). Assumptions: there is only one hub-type loading port in the system, and the loading and unloading capacity of the port is sufficient for the transportation needs, and it is the starting point for all of the electric coal transportation tasks; there are several unloading ports in the system, each port has several berths, and berthing and unloading operations are arranged in accordance with the principle of first-come first-serve; electric coal transportation tasks can use either waterway or railway transportation; electric coal transportation adopts the pattern of a single-demand, single-mode, single-batch mode; all ships are activated at the beginning of the planning period, and the decision makers can choose the opportunity to switch (enable or store) according to their needs. However, since this switching requires a long preparation time and complicated procedures, it is assumed that the operating status of each ship can only be switched at most two times during the observation period. Decision making: to achieve the goal of the lowest total cost of the transportation system, the following decisions are required: choosing the transportation method for each transportation task; assigning transportation tasks to each ship; determining the task-execution sequence, the best start and end time; and switching the ship's operating status between active and inactive.

Water Voyage and Land Voyage.
To facilitate the reader's understanding, we first introduce the concepts of water voyage and land voyage. Water voyage refers to the procedure by which a ship departs from a loading port, arrives at a designated unloading port, and then returns to the loading port (assuming that the ship must return to the loading port to be on standby after transportation). e water voyage can be further divided into delivery subvoyage and return subvoyage. e former refers to the process of the ship's departure from the loading port to the unloading port, and the latter refers to the process of returning the ship from the unloading port to the loading port. Land voyage refers to the process by which electric coal is transported by the railway system.
Both water voyage and land voyage are derived from the electric coal transportation task. Assume that an electric coal transportation task refers to electric coal being shipped from the port a in period t 1 -t 2 and arriving at the port b in period t 3 -t 4 , where the transportation volume is d tons. en, the water voyage derived from this task is the process of the loading of d tons of electric coal at the port a in period t 1 -t 2 , then starting the delivery, arriving at the port b in period t 3t 4 , and finally, returning to port a after unloading. e land voyage for this task is d tons of coal are loaded from the port a in period t 1 -t 2 and arrive at the port b in period t 3 -t 4 . To simplify the problem, we set the following terms for the land voyage: (1) for any transportation task, the electric coal transportation can always be delivered on time by land voyage and that land voyage is always feasible; however, the cost is higher; (2) because of the single-demand, singlemode, single-voyage mode of transportation, water and land voyages derived from the same transportation are mutually exclusive, and decision makers can choose only one way to fulfill the transportation task. From now on, we shall refer to the water and land voyages derived from the same transportation task as a complementary voyage. For the ease of reading, the sets, parameters, and variables that appear multiple times in this article are listed in Table 1.
It should be noted that this article additionally defines two virtual water voyages, which are, respectively, denoted as 0 and e. Voyage 0 means that the ship starts to carry out water transportation tasks from an inactive status, while voyage e means that the ship no longer performs any water transportation tasks and enters the storage state. For example, if ship A performs voyage 0 at time t 1 and voyage e at time t 2 , it means that ship A switched from the storage state to the operating state at time t 1 and from the operating state to the storage state at t 2 . For voyage 0 and voyage e, the values of t ls i , t us i , T l i , and T u i are all 0, and the value of t le i and t ue i is the length of the observation period.

Feasible Single-Vessel Scheduling
Plan. e concept of a single-vessel scheduling plan can be defined based on the concept of a water voyage. e single-ship-scheduling plan is the work plan and the schedule of the electric coal transportation ship, which determines the operation timing, operation order, and specific water voyage nodes. Let vector ((α vp ), (y ip ), (t s ip ), (T wl ip ), (T wu ip )) represent a known singleship-scheduling plan p. Whether the design of a single-shipscheduling plan p is reasonable depends mainly on whether the key parameters satisfy the following constraints: Constraint group 1: the constraints of the task selection and the sequencing of the water voyage. Equations (1)-(4) provide the first set of constraints that ensure that the single-ship-scheduling plan is feasible, which mainly describes the conditions that a feasible plan p should meet in the selection and sequencing of the water voyage. x ijp is a binary variable; if the ship in the single-ship-scheduling plan p starts water transportation task j immediately after it finishes task i, it takes a value of 1; otherwise, it is 0. Equations (1) and (2) ensure that the designs for voyage 0 and voyage e are satisfied, which means the plan p starts at virtual node 0 and ends at the virtual node e. Equation (3) is used to ensure that any water voyage i has presequence and postsequence voyages if it is to be carried out. e algebraic relationship between y ip and x ijp is defined in equation (4). Journal of Advanced Transportation Constraint group 2: the constraints on the time point when the water voyage is undertaken. e relationship between the three variables t s ip , T wl ip , and T wu ip and the key time nodes of the water transportation task is shown in Figure 2. For example, in the plan p, if the first water voyage executed is i, and then, the t s ip is the starting time of the water voyage i. Since the execution order of the task is depicted by x ijp , the values of t s ip and x ijp are tightly related. In addition, to ensure that the single-ship-scheduling plan is feasible, the values of t s ip , T wl ip , and T wu ip should also meet the loading and unloading time-window constraints related to the water transportation task. I e set of known water voyages, each element of which is denoted as i; denote 0 ∈ I, e ∈ I(0 and e as virtual nodes, the meaning will be explained in the following section) G e set of unloading ports, each element is denoted as g V All ships that can be rented, each element is denoted as v P e set of single-ship-scheduling plan, each element is denoted as p Q e set of fleet-scheduling plan, each element is denoted as q, Q⊆P Parameters: Based on the abovementioned analysis, the second set of constraints that ensure the feasibility of the single-shipscheduling plan is shown in equations (5)- (11), where M represents a maximal positive number. e binary variable z ipt represents whether the corresponding ship of the plan p is performing water voyage i and at unloading status at the time point t (dayt); if so, then, z ipt takes a value of 1; otherwise, it is 0. Equations (5) and (6) depict the algebraic relationship between t s ip and x ijp ensuring the start time of a water voyage is the end time of its preceding voyage, where t 0 v represents the initial time for the ship v to be operational. Equations (7) and (8) ensure that the settings of T wl ip , T wu ip , and t s ip satisfy the timewindow constraints of the relevant water voyage. Equation (9) requires that the value of z ipt is 0; when y ip � 0, namely, if the plan p does not implement the water voyage, the ship used in the plan shall not be in the unloading state of the water voyage at any time. Equations (10) and (11) determine the logical relation between T wl ip , T wu ip , and z ipt at time t i . Equation (12) makes sure that the time for ship docked in the port is not less than the port unloading time used in the plan p. Constraint group 3: constraints on ship selection considering the single-vessel scheduling. Equations (13)- (15) are the third group of constraints, which mainly describe the conditions that the feasible plan p should meet in choosing the vessels. Among them, w v indicates the transportation capacity of the ship v; the variable y ipv indicates whether the plan p uses the ship v to perform the water voyage i; if it is, it is taken as 1; otherwise, it is taken as 0; d i represents the task volume of the water voyage i. Equation (13) is used to ensure that plan p uses only one vessel. Equation (14) is used to ensure that the capacity of the vessel is larger than the transportation demand. Equation (15) is the mathematical relationship among y ipv , y ip , and α vp ; that is, when y ip and α vp are both equal to 1, y ipv is taken as 1; otherwise, it is 0. It should be noted that, in order to reduce the complexity of solving the model and construct a linear programming optimization model, this article linearizes equation (15) to obtain equation (15a). e apostrophe in the equation label indicates it is the linearization version of the equation without the apostrophe. In the following, we define a plan p that satisfies equation (1)-(15) as a feasible single-ship-scheduling plan. e complete set of such feasible plans is denoted as P. (15a)

Feasible Fleet Scheduling Plan
As mentioned at the end of Sections 1 and 2, there will be an interactive relationship between different singleship feasible plans when they are implemented. e fourth set of constraints will be constructed for the impact between different single-ship feasible solutions. Constraint group 4: constraints on the interactive relationship between multiple single-ship plans. In the present study, we define a subset Q(Q⊆P) as the fleet scheduling plan. e elements of the subset Q are a single-ship plan and should satisfy the following two conditions. (1) the single-ship single-task constraint: in Q, a ship can only carry out a single-ship-scheduling plan once during the planning period. (2) e exclusivity constraint of port operation: the number of unloading ships should not exceed the number of berths designed in the port at the same time.
Based on the abovementioned analysis, we propose the feasible conditions for the subset Q shown in equations (16) and (17). Among them, the variable x p represents whether the feasible single-ship-scheduling plan p ∈ P is selected in the subset Q; if it is, x p is taken as 1; otherwise, x p is 0. e indicator variable β ig indicates whether the port is an unloading port in the waterway voyage i. If it is, it is taken as 1; otherwise, it is taken as 0. b g represents the number of berths in the ports g. Equation (16) ensures that each ship can only be operated once with a single-ship-scheduling plan, which means parts of the fleets are allowed to not participate in the plan. Equations (16a) and (16b) are the linearization versions of equation (16) for the same constraint. Equation (17) ensures that the number of unloading ships does not exceed the number of berths designed in the port at the same time. Equations (17a) and (17b) are the linearization versions of equation (17) for the same constraint.

Electric Coal Fleet Transportation Plan Optimization Model
Based on the abovementioned constraints, the scheduling process of waterway transportation has been fully characterized. e last set of constraints is introduced for the constraints of the choice of waterway and rail transportation. Constraint group 5: constraints on the choice of waterway and rail transportation. According to the assumptions, any water transportation voyage can be replaced by a complementary railway land voyage, so the decision maker can choose either of them. e fifth set of constraints contains only equation (18). Equations (18a) and (18b) are the linearization versions of equation (18) for the same constraint. It shows that if any feasible single-vessel solution for the waterway voyage i is not selected for set Q, (y ip x p � 0), the task of corresponding electric coal transportation must be undertaken by the complementary railway land plan (c i must be 1).
Based on the definitions of constraint 1-5, we can provide a complete mathematical model (denoted as ECSSM). e model is also improved from the common SPP model. Let c op p represent the operating costs of the plan p; c a i represents the cost of the land voyage associated with the water voyage i; c r v represents the daily rent of ship v(which equals to the difference between the daily operating cost when the ship v is activated and the daily management cost when ship v is sealed). Equation (19) is the objective function of the model to minimize the total system transportation cost, which equals to the sum of the cost of water voyage and railway land. Equation (20) s.t. equation (1)～ (18).
e ECSSM is a large-scale linear integerprogramming model. e ECSSM takes into full consideration the synergetic transportation of sea and rail, the time-window constraints of transportation tasks, and the congestion in the port. It can adequately describe the characteristics of the electrical coal transportation in the southern coastal areas of China. However, the ECSSM still have some limitation for the wide application of costal waterway transportations. With the consideration of the electric coal transportation's operational practice in the Chinese coast area, the transportation pattern of single vessel, single mission, and single voyage is adopted for the tramp ship schedule. Although this assumption is suitable for the scope of this study, it also limits the application of the ECSSM for the applications with other characteristics, such as multiple-voyage or multiple types of cargoes. Furthermore, the ECSSM also adopts the assumption of fixed vessel speed, which is a common assumption for the waterway transportation with the consideration of the economical speed. However, since the thermal power plant could not risk the shutdown, the tramp ship would consider changing the vessel speed to deliver the electric coal in time for extreme situations. erefore, adopting the assumption of fixed speed ignores the other possible measure to mitigate the damage caused by the operation delay. Although the ECSSM has such limitations; these limitations only become significant in the other cases with distinct differences as well as the extreme application cases. For the purpose of studying the electric coal transportation in the Chinese coastal area, the ECSSM could accurately describe the characteristics of transportation and suitable for further study.
e crux of the model's solvability lies in the size of the single-ship-scheduling plan set P. However, the plan set P involves the setting and combination of many factors such as waterway voyage selection, waterway voyage sequencing, waterway operation time node, and waiting time. In reality, the size of the plan set P is extremely large, and the enumeration algorithm cannot meet the needs of practical application.
us, we must design a targeted solution to solve the ECSSM efficiently.

Algorithm Design
e ECSSM is a large-scale mixed-integer optimization model. At present, there are three methods for solving such models. One is the traditional calculation method based on branch and bound. is method has low efficiency and is powerless for many complicated problems. e second is the heuristic algorithm based on the evolutionary framework. e efficiency of this kind of algorithm is relatively high, but it cannot ensure that the global optimal solution is obtained, and the algorithm has poor robustness and unstable calculation results. e third is the branch-and-price algorithm (B&P) which has been widely adopted in recent years. Because of the high computational efficiency of the branch-and-price algorithm and the theoretically accurate solution of large-scale models, it has become an ideal tool for solving large-scale mixed-integer optimization. erefore, this paper adopts the branch-and-price algorithm to solve the abovementioned model. e core of the branch-and-price algorithm is column generation. e idea of column generation evolved from the original Dantzig-Wolfe decomposition idea. It decomposes the original model into two models, called the master model (MM) and submodel (SM). e decision variables in the main model are only a subset of all variables in the original model. e role of the submodel is to continuously introduce new variables into the main model until the optimal solution of the original model is obtained. e objective function of the submodel is the reduce cost of the introduced variable. e flow of the branch-and-price algorithm can be briefly described as follows: (1) Construct an initial solution and solve the main model. In theory, any feasible solution satisfying all constraints can be taken as the initial solution. In addition, by solving the main model, the dual price of each constraint in the model is also obtained. (2) Solve the submodel and obtain the reduce cost of the variable to be introduced, that is, the objective function of the submodel. If the reduce cost is negative, the variable (called column in the main model) is added to the main model. (3) Solve the main model and submodel repeatedly until the reduce cost of the submodel is nonnegative. In theory, the optimal solution of the main model at this time is the optimal solution of the original model. e abovementioned process is shown in Figure 3.

Mathematical Expression of the Main Model.
e mathematical expression of the main model (denoted as MM) is shown in equations (21)- (25). Among them, P represents the currently known feasible set of single-ship plan; the rest of the symbols are the same as before. e structure of this model is similar to that of the ECSSM; the difference is that the MM uses P to replace the original full set P of feasible solutions of the single-ship plan. Since P is known, the variables of this model are x p and c i . Since P⊆P, the optimal solution of the MM should theoretically be the upper bound of the ECSSM. e scale of P is much smaller than that ofP, so we use Gurobi to solve the model.

Mathematical Expressions of Submodels.
e submodel (denoted as SM) is as shown in equations (26)-(41). P represents the new feasible single-ship-scheduling plan (i.e., a new column) to be generated. π (23) v , π (24) lt , and π (25) i represent the shadow price from equations (23)-(25) of the MM, and the rest of the symbols are the same as before. Equation (26) is the objective function of the SM, which means whether the single-ship feasible plan has a further probability for improvement. If the value of equation (26) is negative, P needs to be further optimized. Equations (27)-(41) are the constraint conditions for the SM. ey are the corresponding equations of equations (1)- (15). e difference is that the P is used to replace the p in original equations, which means the main characterization object of the SM is P. ECSSM-SM: In the submodel, each waterway voyage i is abstracted as a transportation node and the executing order of the voyage is regarded as the transportation between the node. Virtual nodes 0 and e are added additionally. en, the subproblem in this section will be converted to an elemental shortest path problem with resource constraint (ESPPRC). erefore, it can be solved using the common multidimensional labeling method. In this ESPPRC, the total cost of the whole route equals to the total costs of all nodes (water voyages) of the route.

5.2.
Design of Multidimensional Labeling. Multidimensional labeling requires labeling for each transportation node. e label L i used in this article contains a set of operational data for a v-type shuttle tanker from the port to the node pos(L v i ). e labels in this chapter include three different types of information related to the problem attributes.
where t j represents the operation time for the water voyage j. It should be noted that the decision variables of the original problem include the waiting time in the water voyage, so there will be multiple possibilities for the time and cost of completing the water voyage i and the water voyage j in sequence. In this chapter, the waiting time is measured in days, and the discrete processing is conducted for various situations, corresponding to the different transportation routes between the two nodes.
After defining the update rule, the domination rule is defined. For the same ship type v, when the resident node of the labels L v 1 and L v 2 are the same FPSO, namely, , if the following conditions are met, the label L v 2 is dominated by the label L v 1 . e specific conditions are (1) less time consumption t(L v i ) ≤ t(L v j ); (2) less operation cost c(L v i ) ≤ c(L v j ).

Numerical Experiments
6.1. Experimental Design. At present, there are only two feasible algorithms for the ECSSP other than the CGSPA; one is the ESPA, and the other is a heuristic algorithm based on experience (denoted by PHA). Practice rule-based heuristic algorithm is the heuristic algorithm based on the operation rules derived from the practical transportation operation. e practical rules enable the PHA to effectively optimize the transportation task volume and maintain the advantage in enhancing computing efficiency by reducing computational time [37]. e heuristic algorithm is proved to be effective in solving the complex transportation network optimization problem involving multiple network nodes and multiple vehicles [37]. us, PHA is a suitable algorithm for the computing efficiency comparison.
In the first phase of our experiment, we demonstrate the advantages of the CGSPA by analyzing the performance of the CGSPA, ESPA, and PHA in terms of computational accuracy and efficiency. In the second stage of this experiment, we perform an optimization calculation for a real historical instance and analyze the results.
We use the transport demand and fleet data of a Chinese thermal power generation group in May and June 2015 (totaling 60 days) as data. Please refer to Table 2 for further details. e ports involved in the experiment are shown in . We assume that 30 vessels can be rented, including three vessels of 20,000 tons, ten vessels of 20,000 to 60,000 tons, and 17 vessels of 60,000 to 90,000 tons. Please refer to Table 3 for further details regarding each ship (e.g., daily rent, available time, and detailed transportation capacity). In the first stage of the experiment, based on the abovementioned data, we randomly generate test data for the comparative calculation. In the second phase of the experiment, we directly use the abovementioned data for the optimization calculation.

Algorithm Efficiency.
e calculation process of the ESPA is relatively simple. First, we obtain the set P by enumeration, and then, we use Gurobi to solve the ECSSM model. In theory, the ESPA algorithm can solve the ECSSP accurately and is an ideal reference to evaluate the accuracy of the CGSPA. e calculation process of the PHA is relatively complex. is algorithm is the primary method used by power generation companies to solve the ECSSP. A summary of the PHA calculations is provided in Table 4, where Step 2.1 is the heart of the algorithm. Step 2.1 is used to estimate whether the ship v′ can complete the task i k on time. For details, please refer to Table 5.
To compare the performance of the CGSPA, ESPA, and PHA, we generate 15 sets of test data based on the data given in Section 6.1. For convenience, we characterize these sets by the number of unloading ports (P), the number of vessels that can be rented (S), the length of the planning period (D), and the number of tasks (T). For example, if a dataset has five unloading ports, four ships, a planning period of 30 days, and ten tasks, this set is named P5S4D30T10. We use C++ to program all the algorithms and use a PC which has an Intel I5-6500 CPU and 16 GB DDR4-2138 memory to run them. e integer program solver used by the CGSPA and ESPA is Gurobi 7.5. Table 6 shows the calculation results and efficiency comparison of the three algorithms for each test. CDPT represents the total computation time of the CGSPA for solving the submodel, CGUT represents the total computation time of the algorithm calling Gurobi to solve the main model, and CGT represents the computation time of the entire calculations. ENT represents the time consumed by the ESPA to enumerate set P, EGUTrepresents the total time that the ESPA spends using Gurobi to solve the ECSSM, and EGT represents the total time that the ESPA takes to complete the calculations. %CGSPA− ESPA and % CGSPA-PHA present the gaps corresponding to the relative difference between the solution values of the CGSPA and ESPA and between the CGSPA and PHA, respectively. OM represents that the relevant calculation is aborted because the memory capacity of the computer is exceeded.
As shown in Table 6, the CGSPA is far superior to the ESPA in terms of computational efficiency. For example, in the small-sized example P4S7T10D30, the CGSPA obtained the optimal solution in 54.71 s, whereas the ESPA took 5.02 s. In the mid-size example P11S17T30D30, because of the size of set P, the ESPA did not obtain an optimization result in an acceptable time (3 h), whereas the CGSPA completed the calculation in only 176.1 s. In terms of accuracy, the CGSPA is far superior to PHA. For example, in the mid-sized example P9S13T20D40, the result of the CGSPA was 17.6% more accurate than the PHA result, whereas in the larger example P11S17T30D40, the CGSPA's result was 19.2% more accurate than the PHA result. e abovementioned results illustrate that compared with the existing algorithms, the CGSPA has notable advantages in terms of computation time and calculation accuracy.
In addition, to test the robustness of the CGSPA, we perform multiple calculations by adjusting the rentable ship set. e results are shown in Table 7, where NUM of P′ represents the total number of P′ (the number of the columns) created by the CGSPA during the calculation. e number of rentable vessels is not positively related to the total time spent by the CGSPA. For example, P6S9T20D40 (7 ships) has a shorter computation time than P6S7T15D40 (9 ships). Further analysis indicates that the efficiency of the submodel solution and the total amount of P′ generated are the key factors determining the efficiency of the CGSPA. e results shown in Table 7 illustrate that the CGSPA has better robustness and can still solve the problem after specific parameters are changed.

Scenario Analysis.
In the second phase of the experiment, we use the model and algorithm proposed in this paper to solve a real historical instance and compare the optimal solution with the actual electric coal transportation plan used by power generation enterprises.
is issue considers all the rentable vessels given in Table 2, all the tasks given in Table 3, and all the ports given are in Section 6.1 with a 70-day planning horizon. e relevant rental fee of August and September in 2016 is set as the average relevant ship rental rate. e railway rates are taken from the China Railway Corporation announced at the same time as the tariffs.
Comparisons of the key parameters of the electric coal transportation optimization plan and the actual transportation plan are shown in Table 8 and Figures 4-6 , where OPR represents the optimized plan, RPR represents the actual transportation plan, NUM represents the number of vessels used, GRT represents the total renting duration of the vessels in days, GWT represents the total waiting time of the vessels in days, RTU represents the ratio of GRT to GWT,  Table 4: A summary of the PHA calculation process.

PHA calculation process
Step 0: according to the lower bound of the unloading time window, sort the water voyages and obtain an ordered set I � i k , where i k ∈ I, k represents the number of the water voyage in I, and for any k, td e i k ≤ td e i k+1 Step 1: select the ship v′ with the highest shipping capacity among the currently available ship set V and obtain the currently available time point t′ of the ship Step 2: if I ≠ ϕ and V ≠ ϕ, then let k � 1; else, go to step 3 Step 2.1: if the transportation volume of the task i k is no less than 60% of the capacity of v′, the ship v′ may complete the task as required (see Table 5); then, delete i k in I, update the ship available time point t′, and record that ship v′ performs task i k ; else, let k � k + 1 Step 2.2: if k > |I|, or t′ > |T|, then remove v′ from ship set V and return step 2.1 Step 3: the calculation is completed Table 5: Procedure for determining whether vessel v can process mission i at time t.
Step 1: if the transportation loading of i k is larger than the ship capacity of vessel typev, then go to step false Step 2: if t ≤ t ls i k , thent wl i k � t ls i k − t, t � t ls i k (calculating the waiting time at the loading port) Step 3: if t > t le i k , then go to step false Step 4: if t + t l i k + t r1 i k > t ue i k , then go to step false Step 5: if t′| i k � ∅, then go to step false Step true: take t′ � min t′

Return true
Step false: t � t − t wl i k , return false (Return value: true, possible; false, and impossible) AW represents the average volume of transportation tasks undertaken by the vessels of the relevant type, RW represents the average utilization rate of the relevant vessels (i.e., the average ratio of the volume of completed traffic to the capacity of the relevant vessels), and % OPR− RPR represents the relative standard deviation of the optimized plan and actual plan in total transportation costs. As shown in Table 6, the performance of the optimized plan in terms of capacity allocation is highly similar to that of the actual plan. Both use all the available ships and reach a high level of utilization (both exceeding 90%). ere was not a significant waste of capacity in either plan. However, the optimization plan is significantly better than the actual plan in terms of the time of use of the ships, with the effect being particularly evident for large electric coal transportation ships. Specifically, in the practical plan, the waiting time for large ships in the port accounts for more than 21% of the total renting duration, compared with only 0.01% in the optimized plan.
To further analyze the advantages of the optimized plan in reducing the waiting time of the vessels, Figure 5 compares the number of waiting ships on each observation point (day) between the optimized plan and the actual plan. e waiting time of ships during the planning period is significantly reduced in the optimized plan, with a daily average number of waiting ships of only 0.30 compared to 2.98 in the actual plan. e abovementioned analysis    illustrates that the optimized plan can effectively reduce the waiting time of the vessels at the unloading port while ensuring that the transportation tasks are completed as required, thus reducing the fleet operation and rental costs. Figures 6 and 7 compare the differences between the optimized and actual plans in terms of fleet size, the time point to start renting a ship, and the time point to stop renting it. As shown in the figure, during phase A of the observation period (from day 0 to day 25), the overall size of the fleet used by both programs was controlled at a scale of 20-25 ships. e options for the two plans at the point of chartering are also similar. However, significant differences were observed during phase B of the observation period (day 25 to day 70). Due to the rationalization of the transportation assignment by the optimized plan, many vessels were allowed to be released on day 25, whereas in the actual plan, releasing was not started until day 45, when a large number of vessels were suspended. As a result, the average fleet size of the optimized plan in phase B was only 14.30, compared to that of 19.65 in the actual plan. e smaller average fleet size combined with the lesser average vessel chartering time makes the optimized plan far superior to the actual plan in terms of electric coal transportation costs, with an approximate 22.60% reduction (% OPR-RPR).
Neither the actual plan nor the optimized plan has applied rail transportation. In other words, water transportation is highly competitive in the actual cases covered in this paper. To determine the advantages of the competitiveness of water transportation, we use the model and algorithm proposed in this paper to study the impact of the floating rent of vessels on the competitiveness of water transportation.
e results of the analysis are shown in Figure 8. Railway transportation begins to be employed when the daily rent increases by 50%. e railway transportation share rate increases to 3.82% and 10.52% when the daily rent increases by 300% and 450%, respectively. is result illustrates that, in the cases discussed in this experiment, the water transportation system is highly competitive and will not fundamentally change the condition that the electric coal transportation is dominated by water transportation even though the rents of ships have risen substantially.

Conclusions
In this paper, a new mathematic model for the ECSSP is developed with the consideration of the synergetic transportation of sea and rail, the time-window constraints of transportation tasks, and the congestion in the port. e model can be used to optimize fleet scheduling planning while considering the multiple transportation plan constraints. A precise algorithm is designed based on the architecture of the branch-and-price algorithm to effectively solve relatively large-scale optimization. e reliability of the model and algorithm is demonstrated through algorithm analysis and a practical case in the southern coastal areas of China. Results indicate that the proposed model and algorithm can effectively solve the ECSSP for the increasing     transformation need in the coastal areas of Southern China. e total transportation cost is significantly reduced by reducing the total waiting time and the active vessels throughout the transportation task.
ere are still issues that need further study. For example, in the transportation process, we assume that the transport system adopts the pattern of a single-vessel, singlemission, single-voyage system. Although this mode of transportation is a convention for the transportation of coal in China, it may lead to wasted shipping capacity in many cases. erefore, how to optimize electric coal transportation under the condition of allowing multiple ships to participate in a set of voyage tasks is an interesting and challenging issue. Meanwhile, the ship's speed is assumed to be fixed, whereas a ship can adjust the speed in practice. Further studies should investigate the transportation of coal considering the variable speed.
Data Availability e transport demand and fleet data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare no conflicts of interest.