Routing a Fleet of Drones from a Base Station for Emission Detection of Moving Ships by Genetic Algorithm

: A fleet of drones is considered in the routing problems with an offshore drone base station, considering the simultaneous movements of drones and ships. A model, entitled meeting model, between a drone and a moving ship is devised, and an extended model is developed based on the vehicle routing problem model. A genetic algorithm based on a sequential insert heuristic (SIH) is designed to solve the model as a holistic framework with two strategies to determine the sequential assignments of ships to drones, namely, the DroneByDrone, and ShipByShip strategies. The proposed models and solution algorithms are demonstrated and verified by experiments. Numerical studies show that the DroneByDrone strategy can overperform the ShipByShip strategy regarding traveling distances. In addition, when considering the simultaneous movement of the ship and drone, improving the drone flying speeds can reduce the flying time of drones rather than optimizing the ship’s moving speed. The managerial implications and possible extensions are discussed based on modeling and experimental studies.


Introduction
As an essential component of logistics transportation, maritime transportation plays a crucial role in connecting global markets and promoting international trade.Global maritime cargo trade is expected to nearly triple between 2020 and 2050 [1].As a result, the growth in shipping demand is causing shipping pollution concerns.According to the findings of the International Maritime Organization, the shipping industry's carbon dioxide emissions constitute approximately 2.2% of worldwide total emissions.Ship fuel combustion also produces nitrogen oxides, sulfur oxides, and other pollutants.The oxides and pollutants negatively impact the environment and human health [2].The pollutants emitted from ships at high altitudes can spread to an area with air currents, causing more severe impacts on the ground and atmosphere, especially in Port City.Therefore, the International Maritime Organization (IMO) institutes Emission Control Areas (ECAs) globally to regulate maritime emissions.
ECA aims to significantly reduce ship pollutant emissions and promote a green and sustainable shipping industry.Within ECAs, maritime regulations mandate the utilization of fuels with a sulfur content that does not exceed a specified threshold.Nonetheless, the absence of a cost-benefit for low sulfur oil compared to heavy oil leads numerous shipping enterprises and vessel owners to favor heavy oil for marine fuel.Consequently, the monitoring of ship emissions in ECAs becomes critically essential.The traditional detection method is to take samples on board, send them to a third party for inspection, or use a fast detector for detection.However, the operation has low detection efficiency and is complex.Hence, innovative emission detection and regulation methods are among the most essential ways to control shipping pollution.Among them, onboard drone-based sniffing technology can realize fast and remote online monitoring and can be used as a feasible and reliable means of pollution monitoring.However, drones' initial acquisition costs and ongoing maintenance are relatively high.As a result, reasonable scheduling and optimizing the routes for drones should be developed to ensure cost-effectiveness.
This study contributes to the related studies in the following three aspects.First, we consider the simultaneous movements of drones and ships according to their respective moving directions and speeds and propose a meeting model to optimize the catch-up time during drone ship emissions.Second, since a fleet of drones' detection scenario for ships is much like the classical VRP, we formulated the routing problem based on VRP models.However, the customer points visited by the vehicles in the VRP are stationary while the ships are moving.The general VRP models cannot be fully adapted to this problem.We extend the VRP models by proposing an extended model considering "moving customers".Finally, we develop a sequential insert heuristic (SIH) and SIH-based bi-stage solution algorithm to compute the drones' and ships' meeting times and positions and then create an SIH-based genetic algorithm to calculate the visiting order of ships by the drone for emissions detection.
The rest of the sections are organized as follows.Section 2 reviews the related studies on ECA, ship emission detection, and drone routing problems.Section 3 introduces the problem and investigates two routing modes, the TSP-based and VRP-based modes.Section 4 formulates the mathematical models.Section 5 presents a genetic algorithm based on a sequential insert heuristic algorithm.Then, we conduct a series of numerical experiments to examine the proposed models and algorithms and discuss the managerial implications in Section 6.We conclude this study in Section 7.

Related Studies
This study featured the drone routing problem inspired by drone applications in emission detection in emission control areas (ECAs).

Emission Control Area
An ECA regulates the release of harmful gases from ships, aiming to diminish the emissions of sulfur dioxide, nitrogen oxides, and particulate matter regarding the effects of ship emissions on human health and the environment [3].Data from the IMO indicate that the implementation of ECAs in Europe has resulted in a 60% decrease in sulfur dioxide emissions from ships [4].Implementing ECA has driven the development of new technologies, such as exhaust gas scrubbers and liquefied natural gas (LNG), which can effectively control ship emissions.
Table 1 summarizes the research on ECA in terms of the research questions, methods, and corresponding regions or ports.ECA is a universal emission control solution; so, most governments and ports have taken it as a policy-making and administration tool.The research in Table 1 investigated various features, including ECA formation, impact assessment, scheduling optimization based on ECA, and stakeholder behaviors.Four methods are used in the research: data-driven [5], mathematical analysis [6], optimization [7], and evaluation methods [8].In Table 1, emission control is significant, and various technologies were applied.Drones should emerge as new devices for ECA management, challenging drone scheduling, and routing problems.

Study
Research Problem Method Region/Port [9] Design China's domestic ECA BLP China [6] The accomplishment of the ECA program development A China [10] Assess noxious gases of ship operations in a potential ECA D Incheon Port [11] The spatial relationship between the net benefit and the size of ECA Simulation China [7] Locate ECA to minimize the impacts of Sulfur emissions on human health MILP China, Africa [12] Dual environmental effects of ECA and reduced speed zones MILP Pacific Ocean [5] Investigate the heterogeneity of the ECA policy effectiveness D Yangtze River Delta [8] Examine whether ECA regulations impact the efficiency of ports DEA European Union, North America [13] Impact of Sulfur ECA on port state control's inspection outcome A - [4] Influence of sulfur emission control areas on particulate matter emission A China [14] Optimize ship route and speed with ECA simultaneously A - [15] Global impacts of ECA A Asia, Europe [16] Decision on liner shipping considering ECA regulations MINP North America [17] Analyze the ship's response strategy to ECA MINP North America

Ship Emission Detection
Ship emission detection technologies play a critical role in reducing the impact of ship emissions on the environment and public health [3].Ship emission detection technologies are still developing, and future research should focus on improving the technology's accuracy, efficiency, and cost-effectiveness [18].Drones equipped with sniffing technology have been tested and applied to detect atmospheric pollution on ships.Using drones for ship emission detection can achieve remote, online, and real-time detection.It is more efficient, flexible, and convenient than conventional detection methods.
Table 2 summarizes some groundbreaking research on ship emission detection in three areas: the research problem, research method, and detection method.The research problem includes ship emission detection methods, detection accuracy optimization, and detection effectiveness evaluation.Research methods related to ship detection issues include neural networks [18], simulation analysis [19], machine learning [20], optimization [21], and utilizing data such as Automatic Identification System (AIS) data [22], satellite images, and videos.Detection methods involve remote sensing, onboard measurement [23], and imaging pattern recognition [24].This study applied a fleet of drones to ship emission detection and investigated the scheduling and routing problem of drones equipped with specialized emission detection devices.

Study
Research Problem Detection Solution [18] Individual ship pollution or ship plume evaluation Onboard NN [20] A big data analytics framework for ship performance monitoring Onboard ML [25] Apply an energy consumption approach to estimate container ship emissions Onboard EA [22] A method evaluating emissions in detailed Onboard EA [23] Precise identification of carbon emission hotspot areas and ship types in ports Onboard D [19] Monitor the SO2 concentration in ship exhausts Imaging S [26] Estimate SO2, CO2, and NO concentrations of exhaust gas by imaging methods Imaging NN [27] Investigating the impacts of satellite fire observation accuracy on the top-down nitrogen oxide emission estimation in northeastern Asia.RS D [21] Develop a marine diesel SOX online detection system RS H [28] Remotely detect sulfur contents in fuel oil used by ships in ECA RS S [24] Detect NO2 pollution plumes from ships with the TROPOMI/S5P satellite sensors.RS D

Drone Routing Problem
Several complex stages are involved during ship emission detection, including law enforcement personnel boarding, document verification, fuel sampling, and actual detection [29].Drones offer numerous benefits, including their high mobility, precise information gathering, and prompt data transmission, making them widely applicable in various fields [30].Drones can detect and monitor ship emissions in the emission control area, constrained by energy consumption and cost [31].Most importantly, drones used for ship emission detection must handle the simultaneous movements of ships, which is a distinct feature in the literature and challenges the existing routing models and solution algorithms.
Most of the literature formulates the drone routing problem as a mixed-integer linear programming (MILP) to determine the task assignment and the visited orders for each drone.Table 3 lists the research on drone routing problems under various scenarios.The VRPs are NP-hard; so, general MILP solvers can only solve small-scale tasks.Hence, to solve medium-and large-scale tasks, researchers usually develop heuristics [32] (e.g., greedy, neighborhood search), a metaphor approach [33,34] (e.g., genetic algorithm, Ant colony optimization), and a decomposition-based algorithm [35] (e.g., Branch-and-Cut, Branch-and-Price).
In terms of application, the delivery and distribution (mainly including emergency [36], last-mile delivery [33], and monitoring [37]) can inspire the model and algorithm developments.We can observe distinct features of the drone scheduling problems: endurance and collaborative delivery with trucks [38].This study formulated the drone routing problem based on VRP modes and their concepts.In ship emission detection, the tasks performed by a fleet of drones and a drone may visit more than one ship.
Applying drones to ship emission detection has the characteristics of a basic VRP.Specifically, a drone flies from a drone base station to the ship, whose position can be determined by the AIS receivers connected with the station [22].Namely, the ship to be detected corresponds to the customer in a typical VRP.As a difference, the ship is not a stationary object but moves with time.These distinct features challenge the existing models and solution algorithms studied in Table 3.To model the proposed problem, we review the related VRP with drones to provide a fruitful methodological reference.

Problem Statement
In this section, we consider a routing problem with a single drone base station and a fleet of drones for ship emission detection (Section 2.1) in an ECA (Section 2.2).As studied in Section 2, the drone routing problem for ship emission detection is featured by the dynamic movements of drones and ships in an ECA.Therefore, it is a complicated extension of the vehicle routing problem, and it challenges the modeling and algorithm development due to the new feature.
The ships' sailing direction and speed are known, and the drones' speed is significantly higher than the ships' sailing speed.The drones fly toward the ships' direction of travel and meet them on the routes.The objective is to minimize the total flying time of the drones for efficient, speedy, and cost-effective ship emission detection.We analyzed a VRP-based mode for the routing problem of a fleet of drones.
The VRP-based mode employs a fleet of drones, and each drone visits a sub-set of ships, as depicted in Figure 1.Based on factors such as the number of ships, their locations, the movement directions, and speeds, it decides the routes for each drone to detect the ships.In the VRP-based mode in Figure 1, two drones start from the base station, and ship 4 is detected before ship 5 and the detection task is completed shortly after their journey begins.Many factors are involved in ship pollution emission detection by drones: drone cost, endurance distance, drones, fly distance, detection efficiency, ships, efficiency, fly times, and emergency capability.The number of drones in the fleet will impact these factors.Changing the number of drones in the base station significantly affects the waiting time and detection sequence of ships in ECA.For example, more drones in the fleet may increase detection efficiency while incurring higher drone investment costs.

Modeling
Using drones for ship emission detection is an innovative and effective attempt to improve traditional emission detection.To better leverage the utility of drones, we propose a drone routing problem for ship emission detection.Considering the synchronized motion of drones and ships when ships arrive in the ECA, we propose a meeting mode between drones and ships and integrate the motion constraints of drones and ships into the classic VRP, constructing an extended model for this problem.

Meeting Model between a Drone and a Moving Ship
As studied above, when the drone navigates to the ship's predetermined "present" position, the ship will have advanced to a new location along its trajectory, as defined by its current and intended positions, as demonstrated in Figures 2. Therefore, it is advantageous to consider the movements of both the drone and ship simultaneously.We calculate the rendezvous point, considering their respective trajectories, and then direct the drone to this estimated position.
For a ship  situated at coordinates (  ,   ), its movement is directed toward a target location (  ,   ) at a speed   .Meanwhile, a drone base station k is positioned at ( ̅  ,  ̅  ).If station k deploys a drone at speed  ̅  to approach ship , the drone must travel for a duration  to intersect with the moving ship at its target trajectory.
In Figure 2, the point of intersection between the drone and the ship is depicted.The ship moves from point  (  ,   ) to point  (   ,    ).Two circles can be drawn: one centered at ( ̅  ,  ̅  ) for the drone with a radius of  ̅  ,, and another for the ship centered at (  ,   ) with a radius of   .The meeting point is where these two circles intersect on the line ab.This intersection is calculated using model (1), where    ,    are specified in equations (2) and (3).The [Mmeet] model can be solved using Equations ( 4)- (10).A feasible meeting time is then ascertained by Equation ( 11), disregarding any negative values.The meeting coordinates (  ,   ) are determined by the meeting time t, as calculated in (1).
The meeting position (, ) and flying time   (the valid value in (10)) can be computed for a drone at the base station  and a ship  by the above meeting model, denoted as follows: ,   ,   ←  [meet] ( ̅  ,  ̅  ,   ,   ,    ,    ,  ̅  ,   ).

Baseline VRP Model
The VRP represents a crucial and fundamental problem in combinatorial optimization, gathering extensive attention from numerous researchers [47,48].In a VRP, the objective is to optimize the routes for a fleet of vehicles to minimize the total travel costs, considering a set of designated customer locations.Using a drone fleet for ship emission detection parallels a VRP scenario, where each drone embarks from a base station and sequentially visits a series of ships, ensuring a single visit per ship.Therefore, we will refer to the VRP models and establish the baseline model.
The baseline model encompasses three distinct sets: a set of ships denoted as  = {1,2, ⋯ }, indexed by ; a set of drones denoted as  = {1,2, ⋯ }, indexed by ; and a combined set of ships and drones denoted as  + = {1,2, ⋯ }, also indexed by .The variable   represents the time taken for a drone to fly from ship  to ship .D denotes the total number of drones available at the station.The binary decision variable   ∈ {0, 1} is set to 1 if drone k is assigned to visit ship  immediately after ship , and 0 otherwise.Similarly, the binary decision variable   ∈ {0,1}, is 1 if drone  is scheduled to visit ship , and 0 otherwise.Additionally, the binary decision variable   signifies the sequential order in which ship  is monitored by drone .Utilizing these data and variables, the baseline model is designed to optimize the allocation of ships to drones.min ∑     ,, (11) Subject to + 1 ≤   + (1 −   ), ∀ ∈  + ,  ∈ ,  ∈  (16)   ,   ∈ {0,1}, ∀,  ∈  + ,  ∈  (17)   > 0,  ∈  + ,  ∈  (18) In the baseline model, the objective function (11) represents the minimum flight distance of drones; Constraint (12) ensures that each ship must be assigned to an AGV for detection.In Constraint (13), the assignment variable (  ) is presented by the sequencing variable (  ).Namely, ship  is detected by drone  if ship  is assigned to drone .Constraint (14) ensures that the in-degree and out-degree of visiting the ship are either one or zero.In Constraint (15), the number of ships detected by each drone is no more than the specific value.Constraint ( 16) sequences the drone operations, namely, the orders of ships the drones visited.The constraint avoids loops within a route by the order variable of a drone () visiting two ships (, ).Here,  = ||.When   = 1, the constraint is   + 1 ≤   , namely, the orders of visiting the ships , .Constraint (17) and Constraint (18) constrain the ranges of the variables.

Extended Model
In the baseline model,   represents the traveling time (distance or other costs and metrics) from a position  to another position .In other words,   is constant because the positions of ,  are determined and not changed in the model.However, in the ship emission detection scenario, ,  are two moving ships, and the value of   depends on the time of the drone visiting the ships.As studied in Section 3, in the ship emission detection scenario,   represents the drone's flying time from a ship  to another ship .A drone visits  just after visiting .A route from the drone base station  can be denoted by   =< ,  1 ,  2 , ⋯ ,   ,  > , where  represents the drone and a virtual depot, and  1 ,  2 , ⋯ ,   represents the visiting sequence of ships.Apparently,     +1 depends on the sequence < ,  1 , ⋯ ,   ,  +1 ].Formally, we introduce a function   (  , , ) to present the dependence,   , (  ,   ) ←   (  , , ).

Solution Methods
We develop an SIH-based genetic algorithm to solve this problem.A sequential insert heuristics algorithm is devised to compute the meeting positions and times for a sequence of ships to be visited.The genetic algorithm is used for the iterative solution.

Sequential Insert Heuristics (SIH)
Based on the meeting model, a sequential insert heuristics algorithm is devised to compute the meeting positions and times for a sequence of ships to be visited.First, initialize the meeting times  and the meeting positions (, ).The following information is available: the present and target positions and speed of the ship; the position of the drone base station and the drone's speed; and a ship sequence to be visited by the drone.Then, by setting the iteration variable , Set  =  1 to compute the first meeting time    The detailed process is shown in Algorithm 1.The algorithm takes a sequence of ships as inputs to generate the times and positions meeting with the drone .In the given six steps,  [meet] (⋅) is invoked to compute the meeting time and position for the ship in the given ship sequence.Then, calculate the meeting time and position with the next ship in this sequence.By  SIH (•) , the meeting times  and positions (, ) between a drone and a sequence of ships can be computed.This is a prerequisite for computing the drone's and target ship's total detection time in the following sections.The time complexity of Algorithm 1 is (), where  represents the scale of the number of ships.

SIH-Based Bi-Stage Solution Algorithm
About the VRP, we have formulated the baseline model to determine the sequence in which a drone visits ships, assuming that the ships are stationary and awaiting detection.This model is then expanded, considering the simultaneous movements of both drones and ships.Subsequently, a two-stage algorithm is developed for this purpose.The first stage involves deriving an initial solution by solving the baseline model, which does not consider the concurrent movements of drones and ships.The second stage employs Sequential Insert Heuristics-Algorithm 1 (SIH)-to translate the ship sequences from the initial solution into actual drone routes, this time factoring in the simultaneous movements of the drones and ships.
The SIH-based bi-stage solution algorithm (Algorithm 2) first solves the VRP to obtain the visit sequence  of drones to ships (Step 1).Then, for each drone visit sequence in turn, the algorithm applies the SIH algorithm- SIH (⋅)-to determine the meeting times and positions of each drone visiting the ships (Step 2).In Algorithm 2, the second stage can improve the first stage's solution, while its optimality is constrained by the first stage, which does not consider the simultaneous movements of drones and ships.

SIH-Based Genetic Algorithm
As a population-based optimization algorithm, the genetic algorithm (GA) is widely used in routing and scheduling problems [49].During the solving process, the GA first encodes a chromosome string, imitating the biological evolution process, and performs selection, crossover, and mutation to generate chromosomes that are different from and superior to their parents.The new chromosomes are then subjected to the same operations of selection, crossover, and mutation, and this iteration continues until the optimal solution is found in the population.This study encoded ships as a sequence in the chromosome genes, representing their numeric order.The decoding process utilizes a sequential insertion algorithm (Algorithm 1) to translate these sequences into drone detection routes, thereby establishing the meeting positions and times.
The characteristics of the problem are integrated into the GA encoding and decoding techniques.As outlined in Algorithm 2, the solution is encoded in the form of a sequence of ships.The decoder then utilizes this sequence to calculate the specific time and position for the drone's detection of each ship.The GA devised in this study can be denoted by  GA with three sets of algorithmic hyperparameters () and two problem-related parameters (ℎ, ).Executing  GA can produce the meeting time and position for each ship, the assigned drone, and the drones' visiting sequences of ships.

The Chromosome Presentation
In this study, the drones used are homogenous, and investment-related decisions predetermine the number.The chromosome consists of a sequence of ship identifiers.For example, a chromosome encoded as integers for 12 ships is as follows: 9 → 8 → 12 → 10 → 1 → 4 → 6 → 5 → 2 → 3 → 11.The sequence is a permutation of the ships' identifiers, 1,2, ⋯ ,12.The sequence indicates the priorities of ships to be considered for the assignments of ships to drones and the ships' emission detection orders.The decoding algorithm (Section 5.3.3) will determine the assignments of ships to drones and each drone's visiting orders to the ships.

The Decoding Algorithm
The decoding algorithm transforms genetic information generated by GAs into solutions for the given problem.This study uses a real-number encoding scheme.Given a chromosome (a sequence of ships), the ships to be inspected are divided into subsequences for drones using either the drone-by-drone sequence division strategy (DroneByDrone) or the ship-by-ship sequence division strategy (ShipByShip).Then, Algorithm 2 (SIH) can decode each subsequence into a feasible solution with determined meeting times and positions for drones and ships.
In summary, the decoding scheme consists of two stages.First, the ships are assigned to drones based on the priorities determined by the chromosome.This stage is entitled "sequence to subsequences" in Section 5.3.3.Second, the ships assigned to a drone constitute a subsequence of ships, which can be used to compute the meeting times and positions by Algorithm 2. In Algorithm 3, the two steps are formulated.The first stage (Step 2 of Algorithm 3) elucidated above is denoted by (, ).In the second stage (Step 3 of Algorithm 3), the meeting times and positions are computed by   (⋅) for each route generated by the first stage.Step 2 Divide the ship sequence into sub-sequences (see Algorithms 4 and 5)

Sequence-to-Subsequence Division Strategies
Sequence-to-subsequence division strategies divide a sequence into subsets under specific principles and rules.Section 5.3.3introduces two methods, namely, DroneByDrone and ShipByShip.
The DroneByDrone strategy involves computing the average number of ships that each drone needs to inspect, given a known sequence of drones performing ship pollution inspections and the number of drones available.This average value, denoted by , is used to group about  ships as a batch and assign them to the same drone.Thus, the algorithmic hyperparameters are given as follows:  = [  ,   ,   ,   ,  = DroneByDrone].Algorithm 4 gives an implementation of the DroneByDrone strategy.Step 2.
Set the sub-sequence for each drone base station.

Return {𝑅 𝑘 |∀𝑘}
The ShipByShip strategy is an alternative approach to the sequence-to-subsequence division strategies.Under the ShipByShip strategy, ships are sequentially assigned to their corresponding drones based on a regular arithmetic progression of the detected ship order.Algorithm 5 implements the ShipByShip strategy with  = [  ,   ,   ,   ,  = ShipByShip].Step 2.

Numerical Experiments
This study analyzes the optimization problem of routing a fleet of drones for ship pollution emission detection in ECAs.To validate the effectiveness of the devised models (Section 4) and the GA and its invoking sub-algorithms (Section 5), experimental studies were conducted using Python programming language on a computer with 11th Gen Intel(R) Core (TM) i5-11300H configuration (Intel Corporation, Santa Clara, CA, USA).

Dataset Generation
In this study, the drone's speed is fixed at  0 = 25 m per second.Ships within the ECA navigate at speeds uniformly distributed between 10 and 20 knots, roughly equivalent to   ~ [10,20] meters per second.Here, [a, b] represents a uniform distribution with the lower and upper limits , , respectively.The ECA is considered to be located at the Yangtze River estuary, spanning an area of 20 km by 10 km.To abstract from specific real-world locations, this virtual area represents a generic ECA, emphasizing methodological generalizability.The drone's endurance range varies depending on battery technology, with different types having different capabilities.It is assumed that the drone's endurance is sufficient for this study, which may influence the setup of ships and the virtual area.
Datasets are generated with the following specifications: The drone base station is positioned at coordinates (0, 0).For each ship, the present and target positions are generated within the 20 × 10 (km × km) range.The position farther from the drone base station is designated as the present position, while the closest is the target.The pair of present and target positions determine the ship's direction of movement.The dataset is named in the format "KkNnVvXxYY", where K indicates the number of drones at the base station, N the number of ships, V is the drone's flying speed, X is the horizontal range, and Y the vertical range.For instance, "K2N5V25X20Y10" denotes an ECA spanning 20 × 10 km 2 with two drones at the base station tasked to detect emissions from five ships, and the drones traveling at a speed of 25 m/s.

Algorithmic Parameter Tuning
The hyperparameters of GA include the number of population iterations (), the population size (), the crossover rate (), and the mutation probability ().The parametertuning experiments used the dataset "K3N30X20Y10V25".The numerical values of the algorithm parameters are determined through parameter-tuning experiments, where each parameter is adjusted and, others are fixed, and the values that produce the best results are selected.
Figure 3 depicts the optimal values of the GA for various combinations of crossover and mutation probabilities.These two parameters jointly influence the computational efficiency and performance of the GA.The results reveal that a crossover rate of 0.8 and a mutation probability of 0.7 produce the best optimization outcomes.The gray area in Figure 4 shows the parameter combination with the highest objective values.
As shown in Figure 4, the mean and distribution of the detection time for multiple drones and ships were investigated for the iterations of the GA.As indicated by the results, we can choose the settings of 250 iterations.
Table 4 shows the impact of population size on the optimal value of the GA.The analysis indicates that the smallest objective function value of the optimal solution is obtained when the population size is 45.As the population size increases, the objective function value remains almost constant while the computation time increases.
These results demonstrate the vital role of parameter tuning and iteration adjustment in achieving optimal performance for the GA.Based on the results, the following settings are obtained: (, , , ) = (250, 45, 0.7, 0.8).

Demonstration of the Models and Algorithms
The dataset "K2N10V10X20Y10" demonstrates the effectiveness of the proposed models and the GAs by using four solutions, as presented in Figure 5.For the demonstration purpose, the present positions of the ships are all set to (0,   ) and   is distributed from 5 to 10 uniformly.Moreover, the target positions of the ships share the same values of   .Notably, the devised model and algorithms do not require the ship movement tracks to be paralleled horizontally.So, the ships should move horizontally.
In Figure 5a, the ships are presumed to be stationary, not moving, and just waiting for emission detection.Therefore, the meeting positions are still located at (0,   ).In Figure 6b, the bi-stage SIH is used.In the first stage, the ships are presumed to be stationary, and the baseline model is used to obtain the assignments of ships to the two drones and the sequence of ships for each drone.So, in Figure 5b, the assignments and sequences of ships to each drone are the same as in Figure 5a.However, the second stage algorithm [SIH] changes the meeting positions and times.In Figure 5b, the movements of ships can be observed.In Figures 5c,d, the ShipByShip and DroneByDrone-based GAs are executed to generate demonstrative results.In these two figures, we can observe no differences.However, when the algorithms are repeated 50 times, the average objectives will incur minor differences, as presented in Table 5.The two decoding strategies will also incur discrepancies in convergence speeds and quality, as shown in Figure 6.The DroneByDrone strategy can achieve better solutions while lowering convergency performance more than the ShipByShip strategy.

Comparisons of the Division Strategies
The impacts of dataset sizes on the division strategies are investigated in Figure 7.The solution methods, [GA + DroneByDrone] and [GA + ShipByShip], are executed for the datasets with different numbers of drones (  ) and ships ( || ),  = 2,3, and || = 10, 15, ⋯ , 50.With the same number of ships to be detected, the total time for detecting ship pollutants is significantly shorter when using three drones compared to that when using two drones.However, under the same drone quantity condition, adopting the DroneByDrone strategy results in shorter drone flight times than when using the ShipBy-Ship strategy.As described in Section 5.3.3, the DroneByDrone strategy prioritizes scheduling drone visits to ships based on the number of ships a drone can access, resulting in a sequence of nearby visits.In contrast, the ShipByShip strategy assigns each ship to a corresponding drone.Compared with the DroneByDrone strategy, the ShipByShip strategy

Sensitivities of Ships and Drones' Speeds
The speeds of both drones  0 and ships   are contingent upon technological advancements and investments.Increasing the speeds of ships and drones results in higher energy expenditures and augmented emissions.Additionally, technological limitations can restrict the maximum achievable speeds of drones.Accelerating drone speeds might necessitate innovation and further investment.In the subsequent analysis, we explore how variations in the speeds of drones and ships affect the total operational time.
Table 6 presents the sensitivity of the ships' speeds (  ) by varying them −5% and 5% and the sensitivity of the drones' speeds (  ) by varying them −5% and 5%.From the information in Table 6, we can see that the increase or decrease in the shipping speed does not show an apparent correlation regularity on the effect of the objective function, while the rise in the drone flight speed can effectively decrease the value of the objective function; the decline in the drone flight speed will make the value of the objective function grow larger.As presented in Figure 7, improving the flight speed of drones is probably more effective for optimizing the detection time of multi-drone ship pollution emissions.As presented in Figure 8, it is evident that the drones' speed significantly impacts the flight time required to complete ship-pollution-monitoring tasks.The effect of drone speed on flight time was analyzed for two and three drones, with task numbers N ranging from 10 to 50.The results show that increasing drone speed leads to a noticeable decrease in flight time, whereas decreasing drone speed leads to a significant increase in flight time.Figure 9 reveals that under the same conditions of drone quantity,   (1 − 5%),   , and   (1 + 5%) almost overlap.This indicates that changes in ship speed do not significantly impact the flight time required for drones to complete ship-pollution-monitoring tasks.

Discussions and Managerial Implications
Using drones for ship pollution monitoring in ECA is an innovative way to monitor ship emissions at sea, and it offers higher monitoring efficiency than traditional monitoring methods.The establishment of ECA has set clear targets and regulations for ship emission reduction in specific areas, but practical implementation of these policies requires efficient regulatory methods and execution capabilities.To achieve effective regulation, it is necessary to optimize the deployment of drones for ship monitoring tasks.In this study, we formulate the problem of efficient and reasonable ship emission detection as a multidrone routing optimization problem.
We draw the following conclusions as managerial implications based on the experimental results.
(1) The optimization of drone flight routing for ship emission detection is a dynamic problem involving the simultaneous movement of drones and ships.The proposed meeting model can appropriately characterize the dynamic problem.The model is an extension and theoretical increment to the static VRP problem.In addition, the meeting model complements the study of operations research theory on dynamic problems.In practice, compared to single drones, multi-drone and ship meeting models based on the VRP model can overcome the difficulties of low detection efficiency and long flight distances.Experiments have verified their advantages in applications.(2) Regarding detection efficiency, the DroneByDrone strategy has an advantage over the ShipByShip strategy in terms of performance.The reason is that with drones traveling at much higher speeds than ships, the DroneByDrone strategy prioritizes the drone detection sequence, thereby reducing the detection time.In addition, in the meeting problem, fast-moving objects should be used as the entry point for algorithm design.Therefore, the DroneByDrone strategy should be prioritized to handle the drone routing problem.However, the meeting model is nonlinear, and it is not easy to find an exact solution.Thus, linearization methods need further research to save costs.(3) Our research focuses on the scheduling problem of drones at the operational level, with the understanding that coordination mechanisms and information systems are vital determinants of operational performance.First, an effective drone scheduling system must dynamically track the positions of ships using data acquired from AIS receivers.Second, drone base stations must exchange scheduling information and integrate their schedules into a centralized, top-level scheduling system.This integration facilitates the global deployment of drones.Third, local ECA authorities should collaborate in sharing the schedules for ship emission detection and the resulting data, ensuring comprehensive monitoring and management.

Conclusions
The adverse effects of ship emissions on the ecological environment and human health in coastal areas have prompted a need for innovative ship emission detection methods.In this regard, we investigated the optimization problem of drone routing in ship emission detection.Firstly, effective detection methods are necessary for the ECA to achieve good ship emission monitoring and control effects, thereby reducing the impact of ship emissions on people's health in the coastal areas.Traditional detection methods suffer from low efficiency and cumbersome procedures.Using drones for ship emission detection is an innovative and effective attempt to improve traditional emission detection.To better leverage the utility of drones, we propose a multi-drone routing optimization problem for ship emission detection.Secondly, considering the synchronized motion of drones and ships when ships arrive in the ECA, we propose a meeting mode between drones and ships and integrate the motion constraints of drones and ships into the classic VRP, constructing an extended model for this problem.Thirdly, we develop an SIH-based genetic algorithm to solve this problem.The sequence insertion algorithm determines the meeting time and location of drones and ships, and the genetic algorithm is used for the iterative solution.Finally, taking the Yangtze River Delta ECA as an example, we develop a dataset generation method and conduct experiments to demonstrate and study the proposed model, algorithm, and sensitivity analysis.
Our study proposes an innovative and practical approach to ship emission detection using drones and develops a multi-drone routing optimization model and an SIH-based genetic algorithm for solving it.The experimental results demonstrate the effectiveness and feasibility of the proposed method, which can be applied to other ECA or environmental monitoring tasks.The baseline Model (Section 4.2) can obtain routes while not considering dynamic simultaneous movements of ships and drones; it, coupled with SIH, can constitute the bi-stage solution algorithm.Further, a GA can use SIH as a decoding algorithm with two strategies to allocate ships to drone routes, DroneByDrone, and ShipByShip.The bi-stage and the GAs (with DroneByDrone or ShipByShip) can solve the drone routing problem, while the bi-stage is not competitive compared with the latter two methods.
In future research, this study can be expanded in several ways: (1) The base station for drone flights is currently fixed and singular.Future research could explore the possibility of having multiple base stations in an ECA and optimizing the selection of base station locations in conjunction with drone detection routing to minimize drone detection routing and achieve better management and economic value.Additionally, the genetic algorithm based on sequence insertion used in this study is effective but cannot guarantee optimality.Future research could develop more efficient algorithms to solve this problem.(2) Implementing machine learning and artificial intelligence techniques could be explored to improve the accuracy and efficiency of the drone detection system.This could involve training models to identify patterns in drone flight routing and predict potential breaches of the ECA.Additionally, the use of advanced sensors and imaging technologies could be investigated to enhance the system's detection capabilities.(3) The potential for integrating the drone detection system with other security systems could be explored.For example, the system could be integrated with perimeter intrusion detection or video surveillance systems to provide a comprehensive security solution.This would require careful consideration of data integration and sharing protocols to ensure seamless integration and optimal performance.Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.Written informed consent has been obtained from the patient(s) to publish this paper.

Figure 1 .
Figure 1.A conceptual diagram of the VRP-based mode.

Figure 2 .
Figure 2. A meeting model for a drone and directional mobile ship.

Algorithm 3
Decoding algorithm Input : A chromosome of a sequence of ships (, ): an operator to divide  into || sub-sequences for the stations   ,   ,    ,    ,   , ∀: the present and target positions and speed of ship   ̅  ,  ̅  ,  ̅  : the position of the drone base station and the drone's speed Output : the total flying time [  ,   ,   |∀]: The meeting times and positions of the drones Steps Step 1. Initialize the total flying time. = 0.

Figure 3 .
Figure 3. Objectives of various mutation and crossover probabilities.

Figure 4 .
Figure 4. Distributions of objectives for different iterations.

Figure 5 .
Figure 5. Demonstrations of four solution methods.

Figure 6 .
Figure 6.Trace plots of the GAs with DroneByDrone and ShipByShip.
detection of ships that are further apart from each other, and as ships and drones move synchronously, the flight time during the drone's visit to a ship increases when the ship being visited is farther away.

Figure 7 .
Figure 7. Comparisons of the DroneByDrone strategy and ShipByShip strategy.

Figure 9 .
Figure 9.A diagram of ships' speed sensitivity analysis.

Table 1 .
Pioneering studies on ECA.

Table 2 .
Pioneering studies on Ship emission detection.

Table 3 .
Pioneering studies on drone routing problems.

Table 4 .
Average objective values and computing times for various population sizes.

Table 5 .
The objective values for four different solution approaches.
Funding:Institutional Review Board Statement: Not applicable.