Next Article in Journal
The Digital Economy and Agricultural Modernization in China: Measurement, Mechanisms, and Implications
Previous Article in Journal
Determining the Interaction of Social Sustainability with the Physical Environment: A Case Study of Adana, Turkey
Previous Article in Special Issue
Assessing the Belt and Road Initiative’s Impact: A Multi-Regression Model Based on Economic Interaction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Path Optimization for Collaborative UAVs and Mothership Monitoring of Air Pollution from Port Vessels

School of Shipping Economics and Management, Dalian Maritime University, Dalian 116026, China
*
Author to whom correspondence should be addressed.
Sustainability 2024, 16(12), 4948; https://doi.org/10.3390/su16124948
Submission received: 26 April 2024 / Revised: 21 May 2024 / Accepted: 4 June 2024 / Published: 9 June 2024

Abstract

:
The seriousness of vessel air pollution has forced the International Maritime Organization (IMO) to introduce a series of relevant laws and regulations. This paper proposes a monitoring scheme based on the synergistic operation of motherships and UAVs. This scheme innovatively adopts a harbor sea patrol vessel or the other official vessel (mothership) as the mobile power supply base for UAVs and realizes efficient and accurate monitoring of vessel air pollution in the pre-monitored area at sea by carrying multiple UAVs. The focus of this paper is on the path optimization problem for multi-UAV collaboration with mothership (MUCWM) monitoring, where the objective is to minimize the total monitoring time for MUCWM. The following three main aspects are studied in this paper: (1) multi-UAV monitoring path optimization; (2) the collaboration mechanism between the mothership and multiple UAVs; and (3) mothership traveling path optimization. In order to effectively solve the above problems, this thesis constructs a path optimization model for multi-UAV collaborative mothership monitoring of air pollution from vessels in port waters; solves the model using the improved adaptive differential evolution (IADE) algorithm; and verifies the effectiveness of the model and the algorithm by using the position data in the Automatic Identification System (AIS) of vessels in Ningbo Zhoushan Port. Through the performance comparison and sensitivity analysis of the algorithm, it is confirmed that the algorithm can effectively solve the path planning problem of the collaborative operation between the mothership and multiple UAVs. The research results in this paper not only help to reduce the air pollution level of harbor vessels and improve the efficiency of sea cruising but also play an important supporting role in the enforcement of relevant emission regulations.

1. Introduction

In recent years, international shipping has rapidly expanded as a low-cost way to transport bulk commodities [1], accounting for over 90% of global trade [2]. According to the United Nations’ 2023 Review of Maritime Transport report, maritime trade is expected to grow by 2.4% in the post-pandemic era [3]. However, as the shipping industry flourishes, the environmental issues associated with marine fuels are becoming more apparent [1]. Vessel air pollution is now recognized as the third primary source of pollution, after vehicle exhaust and industrial emissions. [4]. Vessels emit significant pollutants such as carbon dioxide (CO2), nitrogen oxides (NOX), sulfur oxides (SOX), and particulate matter (PM) while navigating. Global data indicate that SOX emissions from vessels make up 10–15% of total human activity emissions [5]. These pollutants not only lead to ocean acidification [1,6] but also cause various respiratory diseases in humans. Studies show that vessel pollution has become a major source of air pollution in port cities, particularly from container vessels powered by high-sulfur fuels exceeding 3.50% sulfur content, which severely pollute nearby waters and air [7].
Serious environmental and health concerns have compelled the International Maritime Organization (IMO) to take decisive actions. Since 2015, the IMO has established four Emission Control Areas (ECAs) in the Baltic Sea, the North Sea, North America, and the U.S. Caribbean. The organization intends to establish additional Emission Control Areas (ECAs) in various regions to minimize pollution from vessels [8]. The International Maritime Organization (IMO) now mandates that vessels entering these ECAs utilize fuel containing less than 0.1% m/m sulfur content [9]. In 2020, the IMO’s Marine Environment Protection Committee (MEPC) set a global limit on sulfur content in marine fuels at 0.50% m/m [10]. In the same year, the Maritime Safety Administration under China’s Ministry of Transportation declared the “2020 Global Marine Fuel Sulfur Restriction Implementation Plan” along with the “Guidelines for the Supervision and Management of Air Pollutant Emissions from Ships”. According to these guidelines, the sulfur content of fuel used by vessels in inland waterway emission control areas must be no more than 0.1% m/m. Additionally, the sulfur content in the fuel used by ships operating within coastal air pollution control areas and on international journeys in non-emission control areas should not surpass 0.5% m/m [11]. Despite these regulations, the extra costs of purchasing gas purification devices or using low-sulfur fuels lead some shipping companies to violate these rules to cut costs. Methods include using cheaper, high-sulfur fuels, shutting down pollution control equipment in ECAs, and falsifying vessel logs and fuel consumption records [1]. Historical data indicate that non-compliance rates can reach as high as 12.3% [11].
There are two main traditional methods for monitoring air pollution from vessels. The first is manual boarding, which includes taking oil samples, checking oil change records in the logbook, and using a rapid sulfur detector for marine fuel oil. The second method involves real-time tracking of air quality using land-based environmental monitoring stations. However, both methods have drawbacks. They are inefficient, susceptible to environmental influences, and offer limited monitoring coverage. This is particularly problematic during emergencies such as epidemics, when manual boarding and testing become almost impossible.
UAVs are characterized by high flexibility, no geographical restrictions, and high operability [12]. As UAV technology continues to develop, its use in the shipping industry is becoming increasingly widespread. UAVs are now used for port infrastructure security inspections, maritime oil spill monitoring, ship piloting, maritime smuggling detection, and hazardous chemical inspections at ports [13,14,15,16,17]. Additionally, UAVs have started to be employed in vessel air pollution monitoring. For instance, on 20 May 2020, a DJI M210 UAV, named “Rhinoceros V2” and equipped with a system for detecting air pollution, carried out exhaust gas measurements close to the Ganjiang Marine Office of the Yangzhou Maritime Bureau in China [18]. On 5 January 2021, Zhangjiagang Maritime Bureau used a UAV-mounted gas detector to detect a vessel with excessive sulfur emissions [19]; on 26 June 2022, Changshu Maritime Bureau used a UAV to monitor air pollution from vessels and inspected and found a vessel with an exhaust sulfur value of 3.0% m/m or more [20]; on 20 October 2022, the Port of Rotterdam Authority attempted to apply the new Aera3 UAV for monitoring air pollution emissions from vessels [21]. These instances show that UAVs are a feasible and promising tool for monitoring air pollution from vessels [7].
However, the above-mentioned studies are still nascent, with several critical issues needing resolution:
(1)
The energy constraints of UAVs [22].
(2)
The limitation of the size of the sea area that can be monitored [8].
(3)
The limitation on the number of vessels that can be monitored [8].
To further optimize the research on monitoring air pollution from vessels using UAVs, it is essential to thoroughly explore and address these issues. This paper proposes a monitoring scheme involving a mothership and multiple UAVs working in concert. In this case, the mothership could be a port sea patrol ship or the other official ship. We have developed a path optimization model for MUCWM monitoring of air pollution from vessels in port waters, aiming to minimize the total monitoring time.
The core focus of this paper is the path planning problem for UAVs collaborating with a mothership to monitor air pollution from vessels. This involves three main aspects: (1) optimizing the monitoring paths for multiple UAVs, (2) developing a synergy mechanism between the mothership and multiple UAVs, and (3) optimizing the travel path of the mothership. We have created a path planning model for MUCWM monitoring designed to minimize the total time needed to complete the monitoring task.
The purpose of this study is to explore how the mothership and UAVs can monitor air pollution from port vessels in a collaborative manner. Based on the self-constructed path optimization model for the joint operation of the mothership and UAVs to monitor the air pollution from port vessels and solved by using the improved adaptive differential evolution (IADE) algorithm, the optimization of the monitoring paths of the UAVs, the collaborative work between the mothership and the UAVs, and the optimization of the travel paths of the mothership are implemented. In the case study, the effectiveness of the proposed model and algorithm is verified by using the vessel position data from the Automatic Identification System (AIS) of Ningbo Zhoushan Port. The experimental results show that our algorithm can effectively plan the traveling paths of the mothership and UAVs in collaborative operations, which is of great theoretical significance and application value for the study of UAVs monitoring air pollution from vessels in port.
The contributions of this study are significant: (1) We propose a novel path planning model for UAV collaborative mothership monitoring of air pollution from vessels. This model optimizes paths for both UAVs and the mothership working together. (2) We employ the IADE algorithm to solve the model. We have verified the superiority of this algorithm through performance comparisons and sensitivity analysis. (3) We conducted numerical experiments using vessels in the sea area around Zhoushan Port in Ningbo, China. We gathered position data for 100 ships from the Automatic Identification System (AIS) to validate the feasibility of our algorithm.
The innovations of this paper are as follows: (1) Model innovation: This study proposes a path optimization model for multi-UAV collaborative with mothership (MUCWM) surveillance. The model also optimizes the collaborative monitoring routes of the mothership and multiple UAVs, and this method can effectively reduce the overall operation time and improve monitoring efficiency. Therefore, the solution strategy of path planning in this paper can provide a useful reference for other scholars engaged in related research. (2) Algorithm innovation: In this paper, we propose an improved adaptive differential evolution (IADE) algorithm based on the collaborative operation of the mothership and multiple UAVs, which successfully solves the path planning problem of UAVs collaborating with the mothership to monitor the air pollution of port vessels. After comparative analysis, the algorithm demonstrates strong search capability and convergence speed, which significantly improves the efficiency and performance of the problem, and thus better meets the requirements of practical applications. The algorithm is especially suitable for harbor environments, where fast problem-solving is required.
The remaining sections of this paper are structured as follows: Section 2 outlines the related literature. The methodology for establishing the MUCWM model and implementing the algorithm is detailed in Section 3 and Section 4, respectively. Section 5 presents the experiments and the results, with the conclusion following in Section 6.

2. Literature Review

2.1. Vehicle-Mounted Drone Scheduling

The vehicle-mounted drone scheduling problem is typically modeled as a variant of the Traveler Problem (TSP) or the Vehicle Path Problem (VRP), with the choice depending on the number of vehicles involved. The Traveling Salesman Problem with Drones (TSP-D) addresses scenarios where one vehicle and one or more drones are collaboratively optimized for operation. In contrast, the Vehicle Routing Problem with Drones (VRP-D) applies to cases involving multiple vehicles and multiple drones, all optimized together. The primary optimization objectives include minimizing the total time, reducing the joint routing costs, and decreasing the customer’s waiting time.
In the problem presented in this paper, the mothership can be viewed as a truck for joint distribution with multiple UAVs, and the decision objective is also to minimize the total time for mothership and multi-UAV monitoring, which can be viewed as a variant form of the TSP-D scientific problem.

2.1.1. Vehicle-Mounted Drone Scheduling (TSP-D)

The TSP-D draws inspiration from Murray and Chu [23], who introduced two related problems: the Flight Side Kick Traveling Salesman Problem (FSTSP) and the Parallel Drone Scheduling Traveling Salesman Problem (PDSTSP). These aim to minimize the total delivery time for a single vehicle and drone. They proposed two Mixed Integer Programming (MIP) formulations and two simple heuristics, tested on scenarios with up to 10 customers. Agatz et al. [24] explored a similar problem called TSP-D, which differs from FSTSP by allowing the return to previously visited nodes and waiting to receive drones at the launch site. It also considers the drones’ maximum flight distance and uses an Integer Linear Program (ILP) to choose the best sequence of operations to minimize completion time.
Bouman et al. [25] developed a dynamic programming-based exact approach to TSP-D, capable of solving much larger instances. Ha et al. [26] introduced a Minimum-Cost TSP-D, aimed at reducing the total transportation cost through two algorithms: a local search-based TSP-LS and a Greedy Random Adaptive Search Procedure (GRASP). Yurek and Ozmutlu [27] devised an iterative optimization algorithm for TSP-D that breaks down the problem into finding truck routes and then optimizing these routes using a mixed-integer linear program. Tu et al. [28] extended the TSP-D to include multiple UAVs traveling with the truck, creating a new variant called “TSP-mD”.

2.1.2. Vehicle-Mounted Drone Scheduling (VRP-D)

The VRP-D extends the classic Vehicle Routing Problem (VRP) by allowing both trucks and drones to independently serve customers, creating a many-to-many relationship between these vehicles. Wang et al. [29] first proposed the VRP-D, which focuses on delivering goods using multiple vehicles and drones. The objective is to minimize the total task duration, assigning multiple customers to each drone per dispatch. The model tests several worst-case scenarios to explore the potential time savings from integrating trucks with drones compared to using trucks alone. Poikonen et al. [30] extended this work by considering the battery life of drones and expanding the range of the worst-case scenarios examined. Kitjacharoenchai et al. [31] developed a model that coordinates trucks and multiple drones to synchronize package delivery, aiming to minimize the time it takes for both vehicle types to reach the warehouse. Dayarian et al. [32] introduced a “drone resupply” model, termed the “Vehicle Route Problem for Drone Resupply” (VRPDR). This model involves a collaboration between a group of drones and vehicles to deliver online orders from a fulfillment center to a home. Karak and Abdelghany [33] introduced the “Hybrid Vehicle Drone Routing Problem” (HVDRP), in which several drones are deployed from a central mother ship to execute numerous pickups and deliveries at the same time. Poikonen and Golden [34] recently developed the “k-Multi-visit Drone Routing Problem” (k-MVDRP), which considers tandem connections between trucks and k drones, allowing drones to deliver one or more packages to customers.
Another significant research direction in the VRP-D is the Vehicle Routing Problem with Time Window and Drones (VRPTWD). Guerriero et al. [35] proposed a VRPTWD model that considers a soft time window and focuses on customer satisfaction. Phan and Suzuki [36] investigated a variant of this problem that includes simultaneous receiving and sending constraints with a multi-objective optimization approach. Schermer et al. [37], Wang et al. [38], and Sacramento et al. [39] have also proposed variant forms of the VRP-D.

2.2. Research Status of UAV Electric Energy Mobile Supply

Typical approaches to tackle the UAV charging problem involve the use of stationary charging stations and mobile charging vehicles (MCVs) [40]. Recently, significant research has focused on mobile replenishment strategies for UAVs [41]. Yu et al. [42] utilized unmanned ground vehicles (UGVs) to replenish UAVs. Their goal was to allow energy-restricted UAVs to visit a series of stations in the least amount of time, introducing an autonomous UAV charging model grounded in the generalized traveling salesman problem (GTSP). They categorized the problem into three types: Multiple Stationary Charging Stations (MSCS), Single Mobile Charging Stations (SMCS), and Multiple Mobile Charging Stations (MMCS). Inspired by aerial refueling, Zhu et al. [43] introduced a novel concept for aerial charging of UAVs on missions through wireless power transmission. They categorized UAVs into charging UAVs (CUAVs) and mission-performing UAVs (MUAVs) and developed a Deep Reinforcement Learning (DRL)-based scheduling algorithm to optimize the flight paths of both UAV types. Zhou et al. [44] showcased a collaborative network involving UAVs and ground vehicles. In this arrangement, a sub-network of UAVs in the air assists a sub-network of ground vehicles via air-to-air and ground interactions. The ground vehicles function as energy supervisors and chargers for the UAVs. This network structure has been effective in areas like pollution monitoring, disaster response, and information distribution. Luo et al. [45] investigated a two-echelon cooperative routing issue involving ground vehicles (GVs) and their onboard unmanned aerial vehicles (UAVs). In this configuration, the GVs function as mobile support units for the UAVs. They formulated a novel 0–1 integer planning model that accounts for the spatio-temporal cooperative constraints affecting the routes of both GVs and UAVs. To address the model, they introduced two heuristic methods: the first creates a comprehensive route covering all targets, while the second devises a travel path for the GVs and synchronizes the flight paths of the UAVs.

2.3. A Study of Path Planning for Monitoring Air Pollution from Vessels with UAVs

The use of unmanned aerial vehicles (UAVs) to monitor air pollution from vessels has attracted considerable academic interest. Scholars have conducted research in this area: Xia et al. [8] solved the path planning problem for UAV monitoring of air pollution from vessels in emission control areas (ECAs). They proposed a mixed-integer linear programming model based on time-expanded networks and developed a solution method using Lagrangian relaxation techniques. Shen et al. [46] addressed the path planning problem of multiple UAVs for collaborative monitoring of marine air pollution in harbors. They constructed a dynamic multi-objective path planning model and proposed a UAV path planning algorithm for dynamic environments. Yuan et al. [47] developed an enhanced tracking algorithm for UAV gas sensor systems to monitor marine vessel emissions. Sun et al. [1] improved the ACO algorithm to solve the UAV path planning problem for monitoring vessel air pollution. They introduced a hierarchical pheromone update strategy and a partition-based pheromone management mechanism to enhance the traditional ACO algorithm. Luo et al. [48] also developed an improved ACO algorithm for scenarios involving a large number of vessels and UAV base stations. Table 1 summarizes the studies on UAV path planning for monitoring air pollution from vessels. Sun et al. [7] introduced a dynamic scheduling strategy based on reinforcement learning (RL) that takes into account the fluctuations in the speed and direction of ships at sea and thus significantly reduces the flight speed of the UAV.
Currently, most of the research on UAV monitoring of air pollution from vessels focuses on the shore-based replenishment method. This method limits the monitoring range of the UAV, and in this mode, the UAV needs to return to shore frequently to recharge, consuming a large amount of electrical energy for the round trip. This reduces the efficiency of energy utilization and increases the time required for the UAV to complete the monitoring task. To address this current problem, Shen et al.’s study [12] explored the use of a ship as a mobile supply station for UAVs. However, the problems not addressed in this study include the following: (1) the possibility of releasing the UAVs prematurely; (2) the charging time of the UAV is not considered; (3) the travel route of the loaded ship is determined at the initial stage, and the paths of the loaded ship and the UAVs are not planned at the same time; and (4) the simulated annealing artificial swarm algorithm, which is designed to be influenced by the initial solution, is highly affected by the initial solution.

2.4. Literature Summary

Currently, research on using UAVs to monitor vessel air pollution is still in the early stages. Although there have been several related studies, challenges such as UAV electrical energy constraints remain. This paper utilizes the mothership as the mobile electric energy supply base for UAVs to address the issue of monitoring vessel air pollution. This approach has relevance to research on vehicle-mounted drone scheduling and UAV mobile resupply. Section 2.1 indicates that the vehicle-mounted UAV scheduling problem typically represents a variant of the Traveling Salesman Problem with Drones (TSP-D) or the Vehicle Routing Problem with Drones (VRP-D). Section 2.2 indicates that most existing studies on mobile supply bases for UAV electrical energy are based on shore, with only a few attempts to utilize a loading ship as a mobile supply base for UAVs. Optimization was carried out to address the possible problems in the only article proposing the study of mobile resupply bases, Shen et al. [12]. In the model of this paper, (1) the release and recovery operations of the UAVs are not synchronized, but are carried out independently according to the actual monitoring needs; (2) the UAV charging time is taken into account for the modeling; (3) the collaborative monitoring routes of the mothership and multiple UAVs are optimized simultaneously, which can minimize the overall operation time; (4) the improved adaptive differential evolution (IADE) algorithm is applied to this problem scenario, the search performance and convergence ability of the algorithm are enhanced by the improved strategy, and finally, the effectiveness of the IADE algorithm is verified by an algorithm performance comparison and sensitivity analysis. The monitoring scheme in this paper provides a more ideal solution for multi-UAV collaborative monitoring of air pollution from vessels in a harbor and improves the monitoring efficiency of air pollution from vessels in harbor waters. Table 1 provides a summary of research on UAV monitoring of air pollution from vessels.

3. Construction of MUCWM Model

In this section, we established a path optimization model for multiple UAVs to collaborate with a mothership to monitor vessel air pollution in the port.

3.1. Problem Description

The workflow of the monitoring of air pollution from vessels in the harbor implemented by the mothership in collaboration with multiple shipborne UAVs is shown in Figure 1, where M denotes the mothership, v1-v10 denote the monitored ships, u1–u3 denote the UAVs, c1–c9 denote the candidate supply bases, the blue solid line in the figure denotes the traveling path of the mothership, the solid lines of other colors denote the sailing path of the UAVs, and the dashed lines denote the paths of the UAVs flying to the supply points. First, the mothership carries the UAVs from the harbor shore to cruise along the predetermined monitoring route. The shipborne UAVs take off from the mothership, fly to the location of the vessel to be monitored in the harbor waters, monitor the concentration of its emissions in the vessel’s plume, and fly to the location of the next vessel immediately after the monitoring until the energy is not enough to maintain the subsequent monitoring tasks; then they return to the mothership to replenish the electric power and continue to perform the monitoring tasks when all of the monitoring tasks are completed. After completing the monitoring tasks of all the vessels, the UAVs will be recovered and carried away from the harbor area by the mothership and return to the departure point of the harbor.
The main decision problem for the above monitoring effort is how to plan the cruise paths of the mothership and the multiple shipborne UAVs in such a way as to minimize the time of the total monitoring mission. Issues to be considered in this problem include the following:
(1)
How to allocate and optimize the UAV monitoring sequence, that is, the collection of vessels monitored by UAVs and the monitoring sequence.
(2)
How to optimize the driving path of the mothership. This problem involves the selection of the UAV energy supply location.
(3)
How to realize the cooperation between the mothership and UAV. The goal to be achieved in this problem is to achieve collaboration between multiple UAVs and the mothership to minimize the total monitoring time.

3.2. Model Assumptions

In the real port environment, it is very complicated to use UAVs to monitor the air pollution of vessels. This paper mainly considers non-environmental factors and the path optimization problem of vessel air pollution monitoring under relatively ideal conditions. Therefore, the following hypotheses are put forward:
(1)
Both the mothership and the UAV can move freely in R 2 (Euclidean plane) (that is, within the port sea area).
(2)
Both the mothership and the UAV sail in a straight line.
(3)
The UAV travels at maximum speed.
(4)
The charging time of each UAV and the time for monitoring the vessel to be monitored are fixed.
(5)
Each UAV can monitor multiple vessels, and each vessel can only be monitored once.
(6)
The position of the vessel to be monitored is considered unchanged.
(7)
The speed of the UAV > the speed of the mothership.
(8)
Ignore the flight distance of the UAV during the take-off and landing phases.
(9)
The mothership’s fuel and the power to replenish the UAV are unlimited.

3.3. Symbol Description

The following lists the definitions and expressions of parameters and variables used in constructing a multi-UAV collaborative mothership monitoring port vessel air pollution path optimization model:
The auxiliary variables are shown in Table 2.

3.4. MUCWM Path Optimization Model

This paper chooses to construct the objective function based on the travel time of the mothership. The total time to complete the monitoring task consists of three parts:
(1) The time spent by the mothership while traveling. (2) The time spent by the mothership waiting for the UAV. (3) The time spent by the mothership charging the UAV. Based on this, a path optimization model for UAVs to collaborate with the mothership to monitor air pollution from vessels in ports is constructed:
Obj:
min T t o t a l = c i , c j E C d c i , c j × z c i , c j + d l M , S t a S M + u i U c j C t ( u i , c j ) t u i , c j M + u i U v j , v k E V y v j , v k u i × T U
The muti-UAV path optimization model constraints are as follows:
u i U v j , v k E V x v j , v k u i = N U ,   i = 1,2 , , N U ,   j , k = 1,2 , , N V
u i U v j V x S t a , v j u i = 1 ,   i = 1,2 , ,   N U ,   j = 1,2 , ,   N V
u i U v j V x v j , l u i u i = 1 ,   i = 1,2 , ,   N U ,   j = 1,2 , ,   N V
u i U x v j , v k u i = u i U x v k , v j u i , i = 1,2 , , N U ,   j , k = 1,2 , ,   N V ,   j k ,   v j , v k E V
u i U w ( u i , v j ) = 1 ,   j = 1,2 , , N V
The muti-UAV collaborative with mothership optimization model constraints are as follows:
0 < d v j , c n × y v j , v k u i × z c m , c n < T m a x L u i , v j / S u N u i , v j V × T V   and   c m = p u i , v j , u i U ,   i = 1,2 , , N U ,   v j , v k E V ,   j , k = 1,2 , , N C , c m , c n E C ,   m , n = 1,2 , , N V
z c j , c k × t ( u i , c j ) t u i , c j M > 0 , u i U ,   i = 1,2 , , N U ,   c j , c k E C ,   j , k = 1,2 , , N C
z c j , c k × t ( M , c k ) t ( u i , c j ) × S U d c j , c k > 0 , u i U ,   i = 1,2 , , N U ,   c j , c k E C ,   j , k = 1,2 , , N C
x v j , v k u i 0,1 ,   u i U ,   i = 1,2 , , N U ,   v j , v k E V ,   j , k = 1,2 , , N C
y v j , v k u i = 1   0 < T m a x L u i , v j / S u N u i , v j V × T V d v j , v k × x v j , v k u i < d v k , p u i , v j 0 T m a x L u i , v j / S u N u i , v j V × T V d v j , v k × x v j , v k u i d v k , p u i , v j , v j , v k E V ,   j , k = 1,2 , , N C ,   u i U ,   i = 1,2 , , N U
z c i , c j 0,1 ,   c i , c j E C ,   i , j = 1,2 , , N C
w ( u i , v j ) 0,1 ,   u i U ,   i = 1,2 , , N U ,   v j V , j = 1,2 , , N V
Constraint (2) means that there are N V UAVs in total that start from the port’s initial point at the same time to monitor the air pollution of vessels in the port sea area. Constraint (3) is the starting route constraint of each UAV, which means that each UAV can only start from the starting point of the port once. Constraint (4) is the constraint that each UAV passes through the last supply point, which means that each UAV can only reach the last supply node once. Constraints (3) and (4) simultaneously constrain each UAV to have only one flight path. Constraint (5) is used to ensure the continuity of the flight path of each UAV; that is, each monitored vessel node has the same degree of entry and exit for its monitoring UAV. Constraint (6) is used to ensure that each monitored vessel can be monitored by a UAV and can only be monitored once. Constraint (7) is the optional supply point constraint and the service sequence constraint provided by the mothership for UAV charging. Constraint (8) is the time window constraint for the mothership to arrive at the supply point m, which is used to ensure that the mothership must arrive at the target supply point to replenish power for the UAV before the UAV’s power is consumed. Constraint (9) is used to determine whether the mothership can reach the target supply node at its maximum speed within the remaining endurance of the UAV. Constraints (10)–(13) are 0,1 variable constraints, indicating the value range of the decision variable.

4. Model Solution Based on Improved Adaptive Differential Evolution Algorithm

In this paper, an improved DE algorithm is utilized to solve the model. In order to make the DE algorithm more applicable to the problem scenarios in this paper, adaptive and arithmetic improvements are made on the basis of the inherent advantages of the original DE algorithm to increase its performance and ability to solve the problems in this paper, and the improved algorithm is referred to as the improved adaptive differential evolution (IADE) algorithm.

4.1. Adaptive Parameter Settings

(1)
Adopt the lifespan mechanism and extinction mechanism for population NP.
This mechanism regulates population size and maintains diversity. The lifespan mechanism operates as follows: Initially, the age of all individuals is set to zero, aligning with the number of iterations. In each generation, when a new individual is created, its lifespan is assigned based on the value of its fitness function. An individual is removed from the population once its age exceeds its predetermined lifespan. The formula for calculating lifespan is as follows:
L F m i n + η × f ( x i ) f m i n f a v g f m i n f a v g f ( x i ) 1 2 × ( L F m i n + L F m a x ) + η × f ( x i ) f a v g f m a x f a v g f a v g < f ( x i )
Among them, f a v g represents the average value of the fitness function value of the current population and f m i n and f m a x represent the minimum value and maximum value of the fitness function in the current population, respectively. The fitness function value of the algorithm in this article is the minimum value of the total time consumption of the monitoring task. That is, f = T t o t a l . L F m a z represents the maximum value that the life span can reach, that is, the maximum number of iterations, and L F m i n represents the minimum value that the life span can reach. At the same time, η = 1 2 × ( L F m a x L F m i n ) . The extinction mechanism strategy is as follows: When the population fitness function value does not change for p × N P times in a row, it means that the population in the current iteration may fall into a local optimum. At this time, the extinction mechanism strategy can be used to increase the diversity of the population. The parameter p is the extinction coefficient. Generally, the larger the extinction coefficient, the higher the probability of an individual being eliminated or reset; the smaller the extinction coefficient, the higher the probability of an individual being retained. This article sets the extinction coefficient p = 0.3 to ensure the timely generation of new individuals. After an individual becomes extinct, a new population of individuals needs to be generated. Here, elite individual replication is used to generate a new population of individuals. The specific formula is as follows:
n u m = ζ   ·   ( N P m a x N P G )
ζ is the replication rate. This article sets the replication rate ζ = 0.9 to ensure the majority of elite individuals. N P m a x is the maximum value of the population size, and N P G is the current population size.
When the lifespan mechanism comes into play, the diversity of the population may decline sharply, which increases the risk of the population falling into a local extreme. Once this situation occurs, the extinction mechanism will be activated to enhance the diversity of the population by eliminating some individuals and using the replication strategy of elite individuals. This mechanism helps prevent premature convergence of the population, thereby improving the search efficiency and convergence of the algorithm. Therefore, the lifespan mechanism and the extinction mechanism collaborate with each other to realize the adaptive adjustment of N P at different stages.
(2)
Introducing an adaptive parameter adjustment strategy of Beta distribution for scaling factor F and crossover probability C R .
Beta distribution is a common continuous probability distribution in the field of probability statistics, and its value range is limited to the [0, 1] interval. The Beta distribution is often used to describe the probability distribution characteristics of random variables, such as describing the prior distribution of probability, proportion, or parameters. In many fields, such as machine learning and mathematical statistics, Beta distribution plays an important role. The specific form of the Beta distribution is given by its probability density function, which describes the probability distribution characteristics of the Beta distribution in detail. The probability density function of the Beta distribution is as follows:
f ( x | a , b ) = x α 1 ( 1 x ) β 1 B ( a , b )
where x is the value of the random variable, α and β are the two shape parameters of the Beta distribution, and B a , b is the Beta function: B ( a , b ) = 0 1 x α 1 ( 1 x ) β 1 d x . When α = β = 1 , the Beta distribution degenerates into a uniform distribution. When α and β are large, the Beta distribution presents a bell-shaped shape concentrated in the middle. When α and β are small, the Beta distribution presents a long-tail shape.
This article introduces a method of dynamically adjusting F and C R ; that is, the dynamic adaptive adjustment method of Beta distribution is introduced to improve the performance of the algorithm. Specifically, the Beta distribution is used to generate random numbers, and then the generated random numbers are mapped to the value ranges of F and C R . Through this dynamic adjustment method, the DE algorithm is able to select appropriate values of F and C R at different stages. In the initial state, we set α 1 = β 1 = 1 so that the probability distribution function presents a (0, 1) uniform distribution. During the population evolution process, the parameters F and C R change based on the following Beta probability distribution formula:
f ( x | α G , β G ) = x α G 1 ( 1 x ) β G 1 B ( α G , β G )
where α G is the number of individuals with better performance in the current F (or C R ) (such as successful evolution or better fitness function value), and β G is the number of individuals with poorer performance in the current F (or C R ) (such as failure to evolve successfully or making the fitness function value better, which does not make the fitness function value more optimal). We substitute the number of individuals with better performance α G and the number of individuals with poor performance β G in the G generation population into Equation (16), respectively, to calculate the values of C R and F in generation G + 1 , and make constraint boundaries for F and C R ( F 0 ,   2 , C R [ 0 ,   1 ] ), so that F and C R can change adaptively according to the iteration results, effectively improving the convergence speed of the algorithm and the diversity of the population. The calculation process of adaptive F and C R is shown in Figure 2.

4.2. Hybrid Mutation Improvement Strategy

The traditional differential evolution algorithm has limitations in mutation improvement strategies, which limits its versatility. To address this problem, this paper proposes a mutation strategy based on “DE-BB” (differential evolution bare bones) to optimize mutation operations. This strategy is an innovative combination of the backbone particle swarm algorithm and the differential evolution algorithm. Experimental results show that this strategy effectively utilizes the advantages of the backbone particle swarm optimization algorithm in the local depth search, thereby improving the convergence accuracy of the differential evolution algorithm. In the strategy, we incorporate elements of the classic differential mutation and combine them with the adaptive improvement strategy of the F-value to generate individual trial vectors. The specific mutation improvement strategy formula is as follows:
u i , j , G = p i , j , G + r a n d × x r 1 , j , G x r 2 , j , G i f   r a n d ( 0 ,   1 ) C R x i , b e s t , G + F × x r 1 , j , G x r 2 , j , G + F × x r 3 , j , G x r 4 , j , G o t h e r w i s e
In the formula, r 1 ~ r 4 are random integers that are different from each other and randomly selected from the population set 1,2 , , N P , and F is the scaling factor, where p i , j , G = r j × p b e s t i , j , G + ( 1 r j ) × N P b e s t j , G , where r j is also a random decimal in the interval [0, 1]. p b e s t i , j , G is the historical local optimal value of the individual, and N P b e s t j , G is the optimal individual value of the current population.

4.3. Hybrid Cross Improvement Strategy

In traditional differential evolution algorithms, the binomial crossover method (binomial) is often used to cross the contribution vector v generated by the mutation operation with the current target vector x to generate a test vector u . Subsequently, the algorithm competes between the trial vector u and the target direction x through a greedy selection mechanism, and only the solution vector with better fitness is retained to enter the next round of iteration, while the worse solution vector is eliminated. However, when the contribution vector v is already very close to the global optimal position, the crossover operation may destroy its complete solution structure. If the test vector u can fully absorb the information of vector v , it will help increase the diversity of potential solutions and reduce the risk of the algorithm falling into a local optimum. To this end, this article quotes the crossover strategy u u i , t of the multi-selection mechanism and combines it with the adaptive improvement strategy of C R to obtain the crossover improvement strategy formula as follows:
u u i , j , G + 1 = u i , j , G i f   r a n d θ   a n d   r a n d C R v i , j , G o t h e r w i s e
where θ represents the mixed crossover probability threshold; u i , j , G represents the value of the j dimension of individual i , calculated by the traditional binomial crossover method. θ represents the choice probability of performing a binomial crossover or a fully absorbing crossover. According to the experiment of the proposer of this strategy, when the value of parameter θ is set in the range of [0.4, 0.6], the average convergence speed and convergence accuracy of the algorithm are better than those in other intervals. Therefore, this article sets the parameter θ to 0.5.

4.4. IADE Algorithm Solution Steps

4.4.1. Codec Strategy

The encoding method used in this article is natural number encoding. Each individual is composed of two chromosomes. The natural number x k in the first chromosome represents the serial number of the UAV, that is, x 1 , x 2 , , x k , , x N V , where N V represents the number of monitored vessels; the natural number y k in the second chromosome represents the serial number of the monitored vessel; that is, y 1 , y 2 , , y k , , y N V also represents the number of monitored vessels. Among them, the upper and lower correspondence between chromosomes represents the vessel to be monitored by the UAV. The chromosome encoding strategy is shown in Figure 3. The UAV with serial number 1 corresponds to the serial number of the monitored vessel {3, 8, 1, 10}, which means that the order of the vessels monitored by UAV No. 1 is {3→8→1→10}.
Decoding requires designing an appropriate decoding solution based on the characteristics of the problem and converting the encoding after the algorithm is executed into the solution to the problem. The specific steps are as follows:
Step 1: According to the order in which each UAV monitors the vessel, consider the UAV endurance time T e n d u . If, before the UAV goes to the location of the next monitored vessel, T e n d u < d v n o w ,   v n e x t + d v n e x t , c p r e + T V , that is, the UAV cannot return to the supply point after monitoring the next vessel, it is necessary to insert a new supply point gene after the gene of the current vessel position, as shown in Figure 4. T e n d u represents the remaining range of the UAV, d ( x , y ) represents the distance between the two positions x and y , v n o w represents the current location of the vessel being monitored by the UAV, and s e q n e x t represents the position of the next vessel that needs to be monitored, and c p r e represents the location of the previous charging point of the UAV; T V represents the monitoring time.
Step 2: Select a charging point that meets the conditions for the UAV based on the remaining power of the UAV. If T e n d u > d v n o w ,   C is satisfied, add the current charging point to the gene sequence and add the feasible charging points to the UAV monitoring sequence that is added to the UAV feasible monitoring sequence set r u , as shown in Figure 4. Finally, a recovery and supply point need to be added at the end of the genetic sequence of each UAV. After the UAV completes its monitoring mission, the mothership is used to recover it.
Step 3: According to the time window constraint in u _ r o u t e , that is, t c n o w > t + d ( c n o w , c p r e ) , select the charging point driving path that meets the mothership speed constraint, as shown in ①→②→⑤ in Figure 5 and the ①→④→⑤ sequence. t c represents the maximum time point when the UAV reaches supply point c and hovers above it, t represents the current time point, and S M represents the speed of the mothership.
Step 4: Calculate the total duration of the monitoring task based on the charging point driving path of the mothership, return the fitness function value, and end.
Table 3 gives the solution form of MUCWM path optimization.

4.4.2. Algorithm Process

The process of implementing the IADE algorithm is as follows:
Step 1: Initialize parameters, including population size N P , extinction coefficient p , elite individual replication rate ζ , individual component value range b o u n d s , mutation scaling factor F , shape parameters α 1 , β 1 , crossover probability C R , mixed crossover probability threshold θ , the maximum number of iterations G m a x .
Step 2: Generate an initial population, with each individual containing 10 D parameters, and the value range of the parameters is b o u n d s .
Step 3: Decode the initial population, multiply the first dimension of the individual by N U , and perform a rounding operation to determine the set of vessels monitored by the UAV. N U represents the number of UAVs and corresponds to the second-dimension parameters of the individual. Attach a set of sequences 1,2 , , N V to the position, sort the parameters from small to large, retain the attached serial number in the sorting, and use this sorting as the order of the UAV monitoring vessel collection; N V means the number of vessels to be monitored.
Step 4: Let G = 1 represent the first iteration.
Step 5: Let i = 1 represent the first individual.
Step 6: Mutation operation: randomly select 4 individuals x r 1 , x r 2 , x r 3 , and x r 4 in the population, and generate mutated individuals according to the mixed mutation improvement strategy of Formula (19) v i , G .
Step 7: Crossover operation: according to the hybrid crossover improvement strategy of Formula (19), generate new test individuals u i , G .
Step 8: Greedy selection: calculate the fitness function value (minimum value of monitoring task completion time); if the fitness function value f ( u i , G ) < f x i , G 1 , update the individual of the previous generation x i , G 1 = u i , G , otherwise the individual of the previous generation is retained.
Step 9: If i < N P , let i = i + 1 and return to step 5, otherwise go to step 10.
Step 10: Adaptive parameter adjustment. According to the fitness function value, use Formulas (15) and (16) to calculate the adjusted population; use Formula (18) to calculate the adjusted mutation scaling factor F and crossover probability C R .
Step 11: If G < G m a x , G = G + 1 , go to step 4, otherwise go to step 11.
Step 12: Output the historical optimal solution f and the algorithm ends.
The pseudo-code of the improved adaptive differential evolution algorithm is shown in Algorithm 1.
Algorithm 1: Improved adaptive differential evolution algorithm
1 Begin
2 Initialize related parameters: N P , p , ζ , b o u n d s , F , α 1 , β 1 , C R , θ , G m a x
3 Generate initial population
4 Decoding
5 for  G in  G m a x
6      for  i  in  N P
7             v i , G = The optimal individual and four random individuals generate mutant individuals through Formula (17)
8            Boundary constraint check after mutation
9           u i , G = Generate new test individuals through crossover using Formula (18)
10            if  f ( u i , G ) < f ( x i , G 1 )  then
11                 x i , G 1 = u i , G
12            end if
13      end for
14       N P = Calculate the adjusted population by Formulas (14) and (15)
15       F , C R = scaling factor and crossover probability adaptively adjusted through Formula (16)
16 end for
17 Output the optimal solution f , the optimal individual N P b e s t , the algorithm execution time and the algorithm convergence chart
18 Finish
The flow chart of the algorithm is shown in Figure 6.

4.5. Algorithm Design for MUCWM Path Optimization Problem

This section mainly describes the basic ideas for solving the MUCWM path optimization problem in three aspects: optimization of UAV monitoring paths, collaboration between motherships and UAVs, and optimization of mothership monitoring paths.

4.5.1. UAV Monitoring Path Optimization

First, the mothership carries multiple UAVs and sets off from the port starting point S t a at the same time. The UAVs will take off from the mothership at the beginning to perform monitoring tasks. At this time, monitoring needs to be planned for multiple UAVs. To achieve this goal, two problems need to be solved: (1) How to determine the collection of UAVs to monitor vessels; (2) how to determine the order of UAVs to monitor vessels.
Ignoring the power constraints of the UAV, that is, treating the UAV power as infinite, we randomly assign vessel collections and monitoring sequences to the UAV. The specific strategy is as follows: First, add an extra dimension 1,2 , , N V sequence to the test individual u i , G . Secondly, after multiplying the values of the first dimensional elements of the test individual u i , G by N U 1 , respectively, since the value range of the individual elements belongs to 0,1 , and after rounding off the elements and adding 1, the value range of the first dimension element value of u i , G will become [ 1 ,   N U ] , q , where N U is the number of UAVs, and then we sort the second-dimensional elements of u i , G from large to small. During the sorting process, the corresponding positions of the additional sequences remain unchanged. Through the above two strategies, the random allocation of the UAV monitoring vessel collection and monitoring sequence is achieved. The random allocation algorithm for UAV monitoring is shown in Algorithm 2.
Algorithm 2: Random allocation algorithm for UAV monitoring
1   u i , G ( 1 ) = np.round ( N U 1 )
2   s e q = np.arange ( N V ) . reshape ( 1 , N V ) Generate additional sequences
3   u i , G = The trial individual formed by appending s e q to the third dimension of u i , G
4Count different values in u i , G 1
5   u n i q u e   = np.unique u i , G ( 1 )
6  for k  in enumerate ( u n i q u e )
7         u i , G = Split u n i q u e _ v a l u e s according to the value of unique_values
8  end for
9  Sequential allocation strategy for UAV monitoring of vessels:
10 for  u i , G ( 2 ) ( m )  in  u i , G ( 2 )
11         u i , G = Sort according to the value of u i , G 2 m
12  end for
Figure 7 shows a random allocation strategy in a scenario with 2 UAVs and 5 monitored vessels. The sequence of UAV monitoring vessels is divided into the following: u 1 : v 2 v 4 v 1 , u 2 : v 5 v 3 .

4.5.2. Collaboration between the Mothership and UAVs

Since, in the MUCWM path optimization problem, the mothership plays the role of mobile resupplying for UAVs, this means that the mothership needs to meet the UAVs at the supply point and replenish power for the UAVs in time. Therefore, the mothership and the UAVs need to face the problem of human–machine collaboration, which will be realized based on the selection of UAV power supply points.
First, divide the monitored sea area into equal-area grids, and use the midpoints of all grids as the set of candidate power supply points for UAVs: C = 1,2 , , N C , where N C is the divided grid quantity. Based on the value of the individual u i , G after the random allocation strategy is executed, and based on the maximum endurance time of the UAV T m a x , it is judged whether it is necessary to insert a power supply point for the UAV between the two vessels to be monitored during the monitoring process. Take the UAV u 1 in Figure 8 as an example to describe the supply point selection strategy, taking into account the UAV power constraints during the UAV monitoring of vessels.
In the initial case, set f l a g = 1 as the flag that shows that the UAV is fully charged and starts to perform monitoring tasks, and loop through each UAV u i ,   u i U , i 1,2 , 3 , , N U in sequence; then select the supply points it needs to pass during its voyage. Let T c o s t = 0 ; T c o s t represents the time consumed by the current UAV to perform the task, loop through the vessel nodes v j , v j V , j 1,2 , 3 , , N V that the UAV needs to monitor, and judge whether it requires the insertion of feasible supply points between two consecutive monitored vessels. The specific strategies are as follows:
(1) Determine whether the vessel node is the first node monitored by the UAV when it departs from the port, that is, whether it meets f l a g = 1 and j = = 1 at the same time. If the conditions are met, it is judged that the UAV flies to the first node. The vessel node to be monitored will return to the port S t a after monitoring it, that is, T m a x will be greater than or equal to 2 d ( S t a , v j ) / S U + T V , where S U represents the UAV flight speed; if the conditions are met, the UAV will fly to the location of v j and monitor it, as shown in Figure 8, the process of UAV u 1 flying to v 2 . In the algorithm, it is necessary to record the remaining endurance time of the UAV T e n d u , the UAV monitoring sequence r u i and T c o s t and set the flag to 0, indicating that the next vessel node is not the first vessel in the UAV monitoring sequence. It is worth mentioning that since the UAV may have multiple candidate power supply points that meet the conditions in subsequent monitoring tasks, the latter node will be calculated from the previous supply point. Therefore, the solution space is an unmanned power supply point. The first monitored vessel node in the UAV monitoring path is a tree structure array with the root as the root. Therefore, the current vessel node v j is assigned to r u 1 as the root node of the UAV monitoring path r u 1 .
(2) Determine whether the next vessel node is the last vessel node in the UAV monitoring sequence. If not, we need to determine whether the UAV can return to the UAV after flying to the next node and monitoring the vessel. Then we determine the last power supply point of the machine, that is, determine whether T e n d u is greater than or equal to [ d ( v j , v j + 1 ) + d ( v j + 1 , c p r e ) ] / S U + T V ; if it can, the UAV will fly to the next supply point. In Figure 8, the UAV u 1 flies from the location of v 4 to the location of v 1 for monitoring. The reason for using the above strategy to determine whether the UAV needs to replenish electrical energy is that it must ensure that the UAV has a candidate supply point where it can land. Even if the UAV reaches the next vessel node and completes the monitoring task, the UAV’s remaining power is not satisfied. The recharge point with electric energy constraint can also be returned to the previous recharge point for charging. This strategy ensures that the UAV can land and charge at any point in time and will not run out of power and have nowhere to land during the monitoring mission. If it is not satisfied, it means that the remaining range of the UAV is not enough to support its subsequent monitoring tasks. At this time, all electric energy supply points c k , c k C , k 1,2 , 3 , , N C are traversed, and the unmanned machine selects feasible power supply points. If T e n d u is greater than or equal to d ( v j + 1 , c k ) , then c k is added to the solution of the tree structure, and branches may appear, as shown in Figure 8. The area in Q 1 , the supply point in the area, represents the position of the supply point that can be reached from the location of the vessel v 2 to be monitored under the constraints of UAV endurance. Then we update T e n d u and T c o s t . At this time, T c o s t = T c o s t + T u , where T u represents the charging time of the UAV, and we let flag = 1. Finally, the time consumption T e n d u from each time the UAV passes through the power supply point is added to the time window set at t w i n d o w in order to solve the path planning of the mothership. In Figure 8, the UAV u 1 is traveling from the location of the vessel v 2 to the vessel v 4 . Between the locations, the c 1 supply point needs to be inserted. At this time, the next node is the first vessel node monitored by the UAV after it is fully charged and takes off from the mothership; that is, it satisfies flag = = 1 and j > 1. Then it is judged whether the UAV can return to the location of the current supply point after flying to the first vessel node to be monitored and monitoring it, that is, T m a x 2 d ( c n o w , v j ) / S U + T V ; if satisfied, then we update T e n d u and T c o s t , and set the flag to 0. As shown in Figure 8, the UAV u 1 flies from c 1 to the vessel node v 4 to be monitored.
(3) If the next node is the last node in the UAV monitoring sequence, select an electric energy supply point for the UAV that satisfies the remaining power constraint (recovery point: the UAV will be recovered by the mothership at this node, and the mothership returns to the port together and ends the monitoring task), then c k is added to the solution of the tree structure, branches may appear, and then T e n d u , T c o s t , and t w i n d o w are updated. At this time, T c o s t = T c o s t + d ( v j , c k ) / S U .
The algorithm for selecting UAV power supply points is shown in Algorithm 3 (the first node time monitored when the UAV takes off from the port).
Algorithm 3: The algorithm for selecting UAV power supply points
1  Begin
2  flag = 1
3  for i in  N U
4   T c o s t = 0
5      for j in  N V
6                    The first node monitored when the UAV takes off from the port
7          if flag = =1 and j = =1
8              if  T m a x 2 d ( S t a , v j ) / S U + T V
9                   T e n d u = E U 2 d ( s t a r t , v j ) / s u t D
10                  T c o s t = T c o s t + T e n d u
11                  f l a g = 0
12                  r u 1 v j Add to solution structure
13             end if
14         end if
15       end for
16  end for
17  Output  r u , T e n d u , t w i n d o w
18  Finish

4.5.3. Mothership Traveling Path Optimization

The results in Section 4.5.2 are used as parameters and input into the optimization algorithm of the mothership monitoring path, which is illustrated in Figure 9. The optimization strategy of the mothership monitoring path is as follows:
(1) First, select one of the feasible monitoring paths of each UAV for calculation. Sort r u in time order according to t w i n d o w to get s e q M , which is the sequence in which the mothership passes through each supply point. As shown in Figure 9, the navigation path of the mothership M is S t a c 4 c 2 c 1 S t a . Then assign the supply point corresponding to the sorted t w i n d o w to the set W . Traverse the elements in W, and use the speed constraint of the mothership to determine whether the mothership can reach the corresponding supply point to charge the UAV before the UAV’s power is exhausted, that is, whether d ( c k , c k + 1 ) / S M is less than t w i n d o w ( c k ) + T e n d u ; if satisfied, add c k to the mothership’s travel path r M , and update T c o s t , as shown in Figure 9, each time the mothership moves from the current supply point to the next supply point, such as c 4 c 2 .
(2) Since the mothership and the UAV can wait for each other at the supply point, it is necessary to determine the order in which the mothership and the UAV arrive at the supply point and then update the time. Specifically, it needs to be judged whether t c k M is greater than t w i n d o w ( c k ) , where t c k M represents the time when the mothership arrives at the supply point. If it is satisfied, it means that the mothership arrived after the UAV arrived at the supply point c k . At this time, T c o s t needs to be updated, T c o s t = T c o s t + t c k M t w i n d o w c k , and the UAV waits for the mothership’s time to be added to the time consumption.
(3) After executing a set of data, assign T c o s t + d ( l m , S t a ) / S M to f . f represents the completion time of the monitoring task. l m represents the mothership passing the last supply point. This formula represents the mothership passing through the last supply point. The time when the vessel-borne UAV returns to the port is added to T c o s t and then assigned to f . After all the data have been traversed, the smallest f is selected as the optimal solution T t o t a l of the algorithm. Finally, output r M and T t o t a l .
The pseudocode of the mothership traveling path optimization algorithm is shown in Algorithm 4.
Algorithm 4: The Mothership Traveling Path Optimization Algorithm
1  Begin
2  Input  r u , T e n d u , t w i n d o w
3  for i in  r u 1
4      for j in  r u 2
5      …
6          for n in  r u N U
7               s e q M = Sort r u in time order according to t w i n d o w
8              W = corresponding supply point sequence after sorting
9               T c o s t = 0
10            for k in W
11                if  T c o s t   + d ( c k , c k + 1 ) / S M < t w i n d o w ( c k ) +   T e n d u
12                      r M c k
13                      T c o s t = T c o s t +   T c o s t / S M + T V
14                end if
15             end for
16              T t o t a l T c o s t + d ( l m , S t a ) / S M
17       end for
18       
19       end for
20  end for
21   m i n T t o t a l = m i n ( f )
22  Output  r M ,   m i n T t o t a l
23  Finish

5. Experiment and Analysis

In this section, in order to prove the effectiveness of the improved adaptive differential evolution algorithm, this article uses Python 3.12.0 to implement the solution of the algorithm. This experiment was run on a computer with an Intel (R) Core (TM) i5-7200U CPU at 2.50 GHz to 2.70 GHz and a memory of 8 GB.

5.1. Data Preparation and Parameter Settings

This paper conducted a numerical experiment on the vessels entering and leaving the port in a certain port area of Ningbo Zhoushan Port. At 17:40 on 3 April 2023, 100 vessels in the above port were obtained from the automatic identification system (AIS). The position data of vessels sailing in the sea area are partly shown in Table 4. The table gives the vessel serial number, vessel name, vessel type, vessel location (including longitude and latitude), and the update time of the current position. In order to more intuitively display the position of the vessel in the port sea area, a real-time satellite image of the selected sea area is given, as shown in Figure 10. The area within the dotted line is the monitored port sea area.
Before parameter setting, we consider some more realistic factors, such as whether the sea state and weather environment have any influence on the path optimization problem of air pollution monitoring of port vessels by shipboard UAVs.
First of all, the performance and stability of the current DJI UAV are very high, and as an example, the wind resistance of the DJI Warp M300 RTK model UAV can reach 12 m/s, i.e., level 6. According to the survey, in some open sea areas, such as the central part of the ocean or some particularly windy sea areas (such as the roaring westerly wind belt in the South Atlantic Ocean), the average wind speed may be between 7 m/s (about 25 km/h) and 10 m/s (about 36 km/h). Wind speeds are somewhat lower in harbor-based waters, so the monitoring mission is largely unaffected by winds at sea. Secondly, wave heights average a few centimeters to 20 m and can only reach 30 m in very extreme natural conditions [49]. However, when monitoring vessels, UAVs need to fly into the vessel’s plume to do so, and the average vessel’s air height (i.e., the height of a vessel in the ocean between sea level and the very top of the vessel) can vary greatly depending on the type and use of the vessel, and the following are the approximate ranges of air heights for some common types of vessels: container vessels: larger container vessels, usually in the range of 40 to 60 m; bulk carriers: medium-sized bulk carriers may be between 20 and 30 m; tankers: roughly between 30 and 50 m; cruise vessels: usually over 50 m; and medium-sized passenger vessels: these are generally between 20 and 30 m. The vessel’s plume is at a position even higher than the highest point of the vessel, so under normal circumstances, the plume height is higher than the waves, and therefore, the UAV monitoring process will not be affected by the waves as long as it maintains the flight height at the plume position.
The default values of the parameters involved in the experiment are shown in Table 5.
The description of the main parameters in Table 5 is as follows:
Number of UAVs: in the scenario of the vessel exhaust gas measurement trial carried out by the Hanjiang Maritime Safety Department of the Yangzhou Maritime Safety Administration, the number of UAVs was established at three.
Number of equal area divisions in the sea area (number of candidate supply points): In the sea area, it is not advisable to select too many candidate supply points, because too many supply points will greatly increase the running time of the algorithm; at the same time, it is not advisable to select too few candidate supply points, because there may be no solution if there are too few supply points. Therefore, this article divides the sea area into nine equal areas; that is, there are nine candidate supply points.
Port starting point coordinates: set according to the coordinates of a central pier in the port.
Maximum UAV flight speed and maximum UAV endurance time: Based on the UAV produced by DJI that can be used in water conservancy: M300 RTK flight parameters: 55 min maximum flight time, 23 m/s maximum flight speed, and 12 m/s wind resistance. Taking into account the impact of the environment and special circumstances, the maximum flight speed of the UAV is set to 18 m/s, and the maximum endurance time of the UAV is set to 45 min.
Average monitoring time of each vessel: based on previous monitoring experience, the time for UAV monitoring of each vessel is set to 3 min [50].
UAV charging time: In order to save the time required for the total monitoring task, the battery is directly replaced on the mothership to replenish the UAV’s power. Therefore, the battery replacement time is set to 2 min.
The maximum sailing speed of the mothership: Calculated based on the previous average speed of China Maritime Safety Administration port patrol vessels of 20.13 knots. After converting the unit into m/min, it is 621 m/min.
Table 6 shows the parameter settings in the improved adaptive differential evolution algorithm.

5.2. Solution to MUCWM Path Optimization Algorithm Based on IADE Algorithm

According to the actual geographical location of the port, this article scales the monitored sea area to the x o y plane rectangular coordinate system of 200 × 200 to run the algorithm of this article. Following scaling, the initial positions of the 100 vessels to be monitored, one mothership, and three UAVs were determined, represented by (x, y) coordinates. These data are partially displayed in Table 7. Initially, the three UAVs and the mothership share the same location. Figure 11 visually illustrates the initial locations of the vessels to be monitored, the mothership, and the UAVs.
The data are substituted into the improved adaptive differential evolution algorithm to solve the problem. Figure 12 shows the MUCWM path optimization results obtained in the above scenario. Different colors represent different UAV flight paths, and the dotted line is when the UAV goes to the supply point for charging. Since the figure involves many line segments, the path of the mothership will be displayed in the table. Table 8 shows the specific results of path optimization for the mothership and three UAVs.

5.3. Analysis of Results

5.3.1. Algorithm Performance Comparison

In order to verify the accuracy of the algorithm in small-scale calculation examples, the results of IADE solving this problem in small-scale cases are compared with the Gurobi running results. We set the problem scales, N U = 1 ,   2 ,   3 ,   N V = { 5 ,   10 ,   20 ,   50 } , to conduct 12 independent experiments for each scale, and we recorded the GAP between the algorithm and Gurobi in each experiment, where GAP = (f(IADE) − f(Gurobi))/f(Gurobi) × 100%. It is stipulated that if Gurobi runs for more than 30 min and still does not obtain a result, it is considered that Gurobi cannot obtain the solution to the calculation example. The IADE algorithm records the results and running time of each generation and selects the running time of the generation with the smallest difference from the Gurobi result. The comparison results are shown in Table 9.
It can be seen from the table that, in the independent experiment with N U = 3 ,   N U = 20 , the GAP value between IADE and the accurate result is less than 5%, and the solution speed has been greatly improved. Although Gurobi can obtain accurate results, it takes a long time to solve, so the comparison with the Gurobi results proves the effectiveness of this algorithm in solving the collaborative routing optimization problem of UAVs and motherships, and it has practical application significance in real-life scenarios.
When solving the model, in order to verify the superiority of the improved adaptive differential evolution algorithm compared to the traditional differential evolution algorithm, this paper compares the parameters of the adaptive improved differential evolution algorithm with the basic differential evolution algorithm. The group size and the number of iterations are set to N P = 100 and G m a x = 300, respectively. The resulting operation iteration diagram is shown in Figure 13.
From Figure 14, it can be seen that the DE algorithm, after the adaptive improvement strategy, shows a trend of being better than the traditional DE algorithm in terms of convergence speed and optimal solution performance.
On the basis of the adaptive parameter differential evolution algorithm, after adding the hybrid mutation improvement strategy and the hybrid crossover improvement strategy, the test example is also run when N P = 100 and G m a x = 300, and the resulting convergence diagram is shown in Figure 14. It can be seen from the figure that although the convergence speed has dropped slightly, the improved algorithm has strong search capabilities and performs well in manifesting the optimal solution. Therefore, the overall solution performance of the improved adaptive differential evolution algorithm is higher than that of the adaptive differential evolution algorithm.
We took the test case data from Shen et al. [12], as shown in Figure 15, and adjusted the parameters to the values in the article, as shown in Table 10. Since the vessel distribution is relatively concentrated, the number of candidate supply points can be increased to improve the efficiency of the algorithm. Since Shen et al. [12] did not consider the UAV charging time, the UAV charging time was set to 0. Bringing the adjusted data into the algorithm of this article, the result is 116.76982 min, which is smaller than the optimal solution of 126.69466 min in the article. The optimization performance of the improved adaptive differential evolution algorithm is further demonstrated. Figure 16 shows the MUCWM path optimization results obtained by using the algorithm of this article to execute the test case data in Shen et al. [12].

5.3.2. Sensitivity Analysis

Figure 17a illustrates that with an increasing number of UAVs, the overall time for MUCWM monitoring decreases, although the rate of reduction in total monitoring time gradually diminishes. Thus, beyond a certain number of UAVs, further increases will not significantly decrease the total duration of the MUCWM monitoring task. Figure 17b reveals that as the number of UAVs grows, the total execution time of the algorithm consistently increases, suggesting a roughly positive correlation between the number of UAVs and the algorithm’s running time.
Figure 18a demonstrates that as the number of grid divisions rises, there is a decrease in the total MUCWM monitoring time, particularly noticeable when the grid count increases from 1 to 9. However, with further increases in grid divisions, the rate of reduction in the total monitoring task time starts to taper off. Hence, beyond a certain number of grid divisions, the advantages of further increases in reducing the overall time of the MUCWM monitoring task become marginal. Figure 18b shows that as the number of grid divisions grows, the total execution time of the algorithm also increases, with the rate of increase accelerating, suggesting that the number of grid divisions significantly affects the algorithm’s running time.
Figure 19 indicates that extending the maximum endurance time of the UAV has a minimal effect on the total duration of the MUCWM monitoring task. This minimal impact is largely due to the introduction of the mothership serving as a mobile power supply for the UAVs, which significantly mitigates the power limitation issues of the UAV. Figure 20 additionally reveals that, while increasing the sailing speed of the mothership reduces the overall monitoring time of MUCWM, this reduction tends to plateau. Thus, after the mothership’s speed reaches a certain level, further increases in speed yield diminishing returns in reducing the total monitoring time. Figure 21 illustrates that, with varying numbers of UAVs, when the number of vessels to be monitored rises from 30 to 150, the algorithm’s running time increases by 2487.962 s. Notably, when the count of vessels monitored in the port surpasses 120, there is a significant spike in the algorithm’s running time. Consequently, the algorithm exhibits improved computational efficiency when the number of vessels monitored remains below 120.
In particular, we discussed the impact on the total monitoring time in the case of leakage of 5–10% of the monitored vessels, and 10 sets of experiments were conducted for each value of leakage, and the results obtained after averaging are shown in Figure 22. From the figure, it can be observed that the total monitoring time tends to decrease as the number of missed vessels increases, and we find that the impact on the total monitoring task is more significant when it is reduced to 10%. This provides a compromise for the port; the port can tradeoff between the time of the total monitoring task and the coverage of the monitored vessels and can choose to give up a small portion of the vessels to obtain a shorter monitoring time, and this analysis provides a certain reference value for the monitoring program of the vessels in the port’s sea area.

6. Conclusions

This paper presents an innovative monitoring solution for vessel air pollution, utilizing a mothership and multiple UAVs to conduct efficient and comprehensive monitoring in maritime pre-monitoring areas. This significantly enhances monitoring efficiency and extends the monitoring scope. The core focus of this research is the path optimization problem for multi-UAV collaborative mothership (MUCWM) monitoring, which encompasses three critical issues: (1) optimizing the monitoring paths of multiple UAVs to ensure efficient and accurate coverage of the target area; (2) developing a collaboration mechanism between the mothership and multiple UAVs to synchronize their monitoring efforts; and (3) optimizing the mothership’s traveling path to maximize overall monitoring efficiency.
To address these challenges, the paper develops a path optimization model for multi-UAV collaborative mothership monitoring of vessel air pollution in port waters. The model is solved using an improved adaptive differential evolution (IADE) algorithm. Utilizing vessel position data from the Automatic Identification System (AIS) of Ningbo Zhoushan Port, we demonstrate the model and algorithm’s effectiveness through an example. Performance comparisons and sensitivity analyses of the algorithm highlight its significant advantages in solving the path planning problems of collaborative operation between the mothership and multiple UAVs.
However, the research also acknowledges certain limitations: the mission optimization model lacks practical constraints such as UAV flight altitude, sea wind speed, and air control. Furthermore, the model designates static supply points as the midpoints of a divided grid network. Future research will explore dynamically adjusting the supply points based on navigation conditions to achieve a better optimal solution.

Author Contributions

Conceptualization, J.S. and L.S.; methodology, J.S.; software, D.Y.; validation, L.S., J.S. and D.Y.; formal analysis, J.S.; data curation, J.S.; writing—original draft preparation, L.S.; writing—review and editing, L.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Liaoning Province Social Science Planning Fund, grant number 2022-ZSK078, the National Natural Science Foundation of China, grant number 72171032, the National Natural Science Foundation of China, grant number 71702019, and the Humanities and Social Science Foundation of the Ministry of Education, grant number 21YJAZH070.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are available in the paper.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Sun, Z.H.; Luo, X.; Wu, E.Q.; Zuo, T.Y.; Tang, Z.R.; Zhuang, Z. Monitoring Scheduling of Drones for Emission Control Areas: An Ant Colony-Based Approach. IEEE Trans. Intell. Transp. Syst. 2021, 23, 11699–11709. [Google Scholar] [CrossRef]
  2. Sirimanne, S.N.; Hoffman, J.; Juan, W.; Asariotis, R. Review of Maritime Transport 2019; United Nations: New York, NY, USA, 2019; Volume 9. [Google Scholar]
  3. United Nations. Review of Maritime Transport 2023: United Nations Conference on Trade and Development. 2023. Available online: https://unctad.org/publication/review-maritime-transport-2023 (accessed on 27 September 2023).
  4. Liang, G.Y. Neglected pollution sources. Environment 2016, 6, 17–19. [Google Scholar]
  5. Transport & Environment, Air Pollution from Ships. 2018. Available online: https://www.transportenvironment.org/articles/air-pollution-ships (accessed on 20 January 2018).
  6. Li, L.; Gao, S.; Yang, W. The enforcement of ECA regulations: Inspection strategy for on-board fuel sampling. J. Comb. Optim. 2022, 44, 2551–2576. [Google Scholar] [CrossRef]
  7. Sun, P.Z.; Luo, X.; Zuo, T.; Bao, Y.; Sun, Y.N.; Law, R.; Wu, E.Q. Emission Monitoring Dispatching of Drones Under Vessel Speed Fluctuation. IEEE Trans. Intell. Transp. Syst. 2022, 23, 21833–21847. [Google Scholar] [CrossRef]
  8. Xia, J.; Wang, K.; Wang, S. Drone scheduling to monitor vessels in emission control areas. Transp. Res. Part B Methodol. 2019, 119, 174–196. [Google Scholar] [CrossRef]
  9. IMO Low Sulfur Fuel Emission Control Partial State Requirements. 2018. Available online: https://www.sohu.com/a/284957167_120015135 (accessed on 27 December 2018).
  10. China MSA. The Latest Global Emission Control Areas and Low-Sulfur Oil Areas Data Collation and Sink. 2022. Available online: http://www.360doc.com/content/22/0823/07/77972649_1044930083.shtml (accessed on 23 August 2022).
  11. OECD. Reducing Sulphur Emissions from Ships. 2018. Available online: https://www.oecd-ilibrary.org/transport/reducing-sulphur-emissions-from-ships_5jlwvz8mqq9s-en (accessed on 20 January 2018).
  12. Shen, L.; Hou, Y.; Yang, Q.; Lv, M.; Dong, J.X.; Yang, Z.; Li, D. Synergistic path planning for ship-deployed multiple UAVs to monitor vessel pollution in ports. Transp. Res. Part D Transp. Environ. 2022, 110, 103415. [Google Scholar] [CrossRef]
  13. Stein, M. Conducting safety inspections of container gantry cranes using Unmanned Aerial Vehicles. In Dynamics in Logistic, Proceedings of the 6th International Conference LDIC 2018, Bremen, Germany, 20–22 February 2018; Springer International Publishing: Cham, Switzerland, 2018; pp. 154–161. [Google Scholar]
  14. Urbahs, A.; Zavtkevics, V. Oil pollution monitoring of sea aquatorium features with using unmanned aerial vehicles. In Proceedings of the 17th International Conference Transport Means, Kaunas, Lithuania, 23–24 October 2014; pp. 75–78. [Google Scholar]
  15. Zhang, S.; Xin, M.; Wang, X.; Zhang, M. Anchor-free network with guided attention for ship detection in aerial imagery. J. Appl. Remote Sens. 2021, 15, 024511. [Google Scholar] [CrossRef]
  16. Primeau, N.; Abielmona, R.; Falcon, R.; Petriu, E. Maritime smuggling detection and mitigation using risk-aware hybrid robotic sensor networks. In Proceedings of the 2017 IEEE Conference on Cognitive and Computational Aspects of Situation Management (CogSIMA), Savannah, GA, USA, 27–31 March 2017; pp. 1–7. [Google Scholar]
  17. Kim, J.S.; Lee, M.J.; Nam, H.; Do, S.; Lee, J.H.; Park, M.K.; Park, B.H. Indoor and Outdoor Tests for a Chemi-capacitance Carbon Nanotube Sensor Installed on a Quadrotor Unmanned Aerial Vehicle for Dimethyl Methylphosphonate Detection and Mapping. ACS Omega 2021, 6, 16159–16164. [Google Scholar] [CrossRef] [PubMed]
  18. Soarabolity. How does UAV Estimate the Sulfur Content of Ship Fuel by Monitoring Ship Exhaust [EB/OL]. Available online: https://www.sohu.com/a/399172599_100226347 (accessed on 2 June 2020).
  19. Flyable technology Soarability. Maritime Administration Explores Implementation of Drone Ship Exhaust Monitoring. 2022. Available online: https://baijiahao.baidu.com/s?id=1738314327938212810&wfr=spider&for=pc (accessed on 14 July 2022).
  20. Jiangsu Maritime Bureau of the People’s Republic of China, News Events. 2022. Available online: https://www.js.msa.gov.cn/art/2022/7/11/art_47_1356157.html (accessed on 11 July 2022).
  21. China seShipping. Port of Rotterdam Is Trialling New Drones to Monitor Ship Emissions. 2022. Available online: https://www.chineseshipping.com.cn/cninfo/News/202210/t20221021_1370484.shtml (accessed on 21 October 2022).
  22. Explanation of the Endurance of Multi-Rotor Drones and Fixed-Wing Drones. Military/Avionics. 2023. Available online: https://m.elecfans.com/article/2238023.html (accessed on 7 September 2023).
  23. Murray, C.C.; Chu, A.G. The flying sidekick traveling salesman problem: Optimization of drone-assisted parcel delivery. Transp. Res. Part C Emerg. Technol. 2015, 54, 86–109. [Google Scholar] [CrossRef]
  24. Agatz, N.; Bouman, P.; Schmidt, M. Optimization approaches for the traveling salesman problem with drone. Transp. Sci. 2018, 52, 965–981. [Google Scholar] [CrossRef]
  25. Bouman, P.; Agatz, N.; Schmidt, M. Dynamic programming approaches for the traveling salesman problem with drone. Networks 2018, 72, 528–542. [Google Scholar] [CrossRef]
  26. Ha, Q.M.; Deville, Y.; Pham, Q.D.; Hà, M.H. On the min-cost traveling salesman problem with drone. Transp. Res. Part C Emerg. Technol. 2018, 86, 597–621. [Google Scholar] [CrossRef]
  27. Yurek, E.E.; Ozmutlu, H.C. A decomposition-based iterative optimization algorithm for traveling salesman problem with drone. Transp. Res. Part C Emerg. Technol. 2018, 91, 249–262. [Google Scholar] [CrossRef]
  28. Tu, P.A.; Dat, N.T.; Dung, P.Q. Traveling salesman problem with multiple drones. In Proceedings of the 9th International Symposium on Information and Communication Technology, Da Nang, Vietnam, 6–7 December 2018; pp. 46–53. [Google Scholar]
  29. Wang, X.; Poikonen, S.; Golden, B. The vehicle routing problem with drones: Several worst-case results. Optim. Lett. 2017, 11, 679–697. [Google Scholar] [CrossRef]
  30. Poikonen, S.; Wang, X.; Golden, B. The vehicle routing problem with drones: Extended models and connections. Networks 2017, 70, 34–43. [Google Scholar] [CrossRef]
  31. Kitjacharoenchai, P.; Lee, S. Vehicle routing problem with drones for last mile delivery. Procedia Manuf. 2019, 39, 314–324. [Google Scholar] [CrossRef]
  32. Dayarian, I.; Savelsbergh, M.; Clarke, J.P. Same-day delivery with drone resupply. Transp. Sci. 2020, 54, 229–249. [Google Scholar] [CrossRef]
  33. Karak, A.; Abdelghany, K. The hybrid vehicle-drone routing problem for pick-up and delivery services. Transp. Res. Part C Emerg. Technol. 2019, 102, 427–449. [Google Scholar] [CrossRef]
  34. Poikonen, S.; Golden, B. Multi-visit drone routing problem. Comput. Oper. Res. 2020, 113, 104802. [Google Scholar] [CrossRef]
  35. Guerriero, F.; Surace, R.; Loscri, V.; Natalizio, E. A multi-objective approach for unmanned aerial vehicle routing problem with soft time windows constraints. Appl. Math. Model. 2014, 38, 839–852. [Google Scholar] [CrossRef]
  36. Phan, D.H.; Suzuki, J. Evolutionary multiobjective optimization for the pickup and delivery problem with time windows and demands. Mob. Netw. Appl. 2016, 21, 175–190. [Google Scholar] [CrossRef]
  37. Schermer, D.; Moeini, M.; Wendt, O. A matheuristic for the vehicle routing problem with drones and its variants. Transp. Res. Part C Emerg. Technol. 2019, 106, 166–204. [Google Scholar] [CrossRef]
  38. Wang, Z.; Sheu, J.B. Vehicle routing problem with drones. Transp. Res. Part B Methodol. 2019, 122, 350–364. [Google Scholar] [CrossRef]
  39. Sacramento, D.; Pisinger, D.; Ropke, S. An adaptive large neighborhood search metaheuristic for the vehicle routing problem with drones. Transp. Res. Part C Emerg. Technol. 2019, 102, 289–315. [Google Scholar] [CrossRef]
  40. Zajac, S.; Huber, S. Objectives and methods in multi-objective routing problems: A survey and classification scheme. Eur. J. Oper. Res. 2021, 290, 1–25. [Google Scholar] [CrossRef]
  41. Qin, W.; Shi, Z.; Li, W.; Li, K.; Zhang, T.; Wang, R. Multiobjective routing optimization of mobile charging vehicles for UAV power supply guarantees. Comput. Ind. Eng. 2021, 162, 107714. [Google Scholar] [CrossRef]
  42. Yu, K.; Budhiraja, A.K.; Tokekar, P. Algorithms for routing of unmanned aerial vehicles with mobile recharging stations and for package delivery. arXiv 2017, arXiv:1704.00079. [Google Scholar]
  43. Zhu, K.; Yang, J.; Zhang, Y.; Nie, J.; Lim, W.Y.B.; Zhang, H.; Xiong, Z. Aerial refueling: Scheduling wireless energy charging for UAV enabled data collection. IEEE Trans. Green Commun. Netw. 2022, 6, 1494–1510. [Google Scholar] [CrossRef]
  44. Zhou, Y.; Cheng, N.; Lu, N.; Shen, X.S. Multi-UAV-aided networks: Aerial-ground cooperative vehicular networking architecture. IEEE Veh. Technol. Mag. 2015, 10, 36–44. [Google Scholar] [CrossRef]
  45. Luo, Z.; Liu, Z.; Shi, J. A two-echelon cooperated routing problem for a ground vehicle and its carried unmanned aerial vehicle. Sensors 2017, 17, 1144. [Google Scholar] [CrossRef]
  46. Shen, L.; Wang, Y.; Liu, K.; Yang, Z.; Shi, X.; Yang, X.; Jing, K. Synergistic path planning of multi-UAVs for air pollution detection of ships in ports. Transp. Res. Part E Logist. Transp. Rev. 2020, 144, 102128. [Google Scholar] [CrossRef]
  47. Yuan, H.; Xiao, C.; Wang, Y.; Peng, X.; Wen, Y.; Li, Q. Maritime vessel emission monitoring by an UAV gas sensor system. Ocean. Eng. 2020, 218, 108206. [Google Scholar] [CrossRef]
  48. Luo, X.; Sun, Z.H.; Qiu, S. Ant colony system based drone scheduling for ship emission monitoring. In Proceedings of the 2021 IEEE Congress on Evolutionary Computation (CEC), Kraków, Poland, 28 June–1 July 2021; pp. 241–247. [Google Scholar]
  49. Drones for Ship Exhaust Monitoring, DroneNet. Available online: https://www.163.com/dy/article/IEOJKCVR05149OCK.html (accessed on 16 September 2023).
  50. Sea Wave. Key Laboratory of Coastal Disaster and Protection (Hohai University), Ministry of Education. 2019. Available online: https://coast.hhu.edu.cn/2020/0512/c2585a203485/page.htm (accessed on 11 November 2019).
Figure 1. Vessel air pollution monitoring route map.
Figure 1. Vessel air pollution monitoring route map.
Sustainability 16 04948 g001
Figure 2. Calculation process of adaptive F and C R .
Figure 2. Calculation process of adaptive F and C R .
Sustainability 16 04948 g002
Figure 3. Chromosome coding strategy.
Figure 3. Chromosome coding strategy.
Sustainability 16 04948 g003
Figure 4. Replenishment point gene insertion strategy.
Figure 4. Replenishment point gene insertion strategy.
Sustainability 16 04948 g004
Figure 5. Viable paths for UAV monitoring.
Figure 5. Viable paths for UAV monitoring.
Sustainability 16 04948 g005
Figure 6. Flowchart of improved adaptive differential evolution algorithm.
Figure 6. Flowchart of improved adaptive differential evolution algorithm.
Sustainability 16 04948 g006
Figure 7. A random assignment strategy for monitoring vessels by UAVs.
Figure 7. A random assignment strategy for monitoring vessels by UAVs.
Sustainability 16 04948 g007
Figure 8. Synergistic strategy for UAVs and the mothership.
Figure 8. Synergistic strategy for UAVs and the mothership.
Sustainability 16 04948 g008
Figure 9. Monitoring path optimization strategy for the mothership.
Figure 9. Monitoring path optimization strategy for the mothership.
Sustainability 16 04948 g009
Figure 10. Schematic diagram of the location of the test vessels in Ningbo Zhoushan Port.
Figure 10. Schematic diagram of the location of the test vessels in Ningbo Zhoushan Port.
Sustainability 16 04948 g010
Figure 11. Schematic of the vessels’ positions after zooming.
Figure 11. Schematic of the vessels’ positions after zooming.
Sustainability 16 04948 g011
Figure 12. Schematic diagram of MUCWM path optimization results.
Figure 12. Schematic diagram of MUCWM path optimization results.
Sustainability 16 04948 g012
Figure 13. Adaptive strategy iteration graph.
Figure 13. Adaptive strategy iteration graph.
Sustainability 16 04948 g013
Figure 14. Iterative diagram of operator improvement strategy.
Figure 14. Iterative diagram of operator improvement strategy.
Sustainability 16 04948 g014
Figure 15. Test data from Shen et al. [12].
Figure 15. Test data from Shen et al. [12].
Sustainability 16 04948 g015
Figure 16. MUCWM path optimization results for running test data.
Figure 16. MUCWM path optimization results for running test data.
Sustainability 16 04948 g016
Figure 17. The relationship between the number of UAVs and the optimal solution and the running time of the algorithm: (a) the relationship between the total time of monitoring tasks and the number of UAVs; (b) the relationship between the algorithm running time and the number of UAVs.
Figure 17. The relationship between the number of UAVs and the optimal solution and the running time of the algorithm: (a) the relationship between the total time of monitoring tasks and the number of UAVs; (b) the relationship between the algorithm running time and the number of UAVs.
Sustainability 16 04948 g017
Figure 18. The relationship between the number of grids and the optimal solution and the running time of the algorithm: (a) the relationship between the total time of monitoring tasks and the number of grids; (b) the relationship between the algorithm running time and the number of grids.
Figure 18. The relationship between the number of grids and the optimal solution and the running time of the algorithm: (a) the relationship between the total time of monitoring tasks and the number of grids; (b) the relationship between the algorithm running time and the number of grids.
Sustainability 16 04948 g018
Figure 19. The relationship between the maximum endurance of UAVs and the total time of monitoring tasks.
Figure 19. The relationship between the maximum endurance of UAVs and the total time of monitoring tasks.
Sustainability 16 04948 g019
Figure 20. The relationship between the maximum speed of the mothership and the total time of monitoring tasks.
Figure 20. The relationship between the maximum speed of the mothership and the total time of monitoring tasks.
Sustainability 16 04948 g020
Figure 21. The relationship between the number of vessels to be monitored and the algorithm running time.
Figure 21. The relationship between the number of vessels to be monitored and the algorithm running time.
Sustainability 16 04948 g021
Figure 22. Relationship between the proportion of vessels to be monitored that are missed and the total time of the monitoring tasks.
Figure 22. Relationship between the proportion of vessels to be monitored that are missed and the total time of the monitoring tasks.
Sustainability 16 04948 g022
Table 1. Summary of research on UAV monitoring of air pollution from vessels.
Table 1. Summary of research on UAV monitoring of air pollution from vessels.
LiteratureCollaboration ObjectElectricity Supply MethodUAV RoutingCarrier RoutingCollaborative Approach
Xia et al. [8]NShore-basedYNN
Shen et al. [22]NShore-basedYNN
Yuan et al. [23]NNNNN
Sun et al. [1]NShore-basedYNN
Luo et al. [24]NShore-basedYNN
Sun et al. [7]NNYNN
Shen et al. [12]ShipMobile supplyYYStep by step
This studyMothershipMobile supplyYYSynchronize
Y: Yes; N: No
Table 2. Model symbol description.
Table 2. Model symbol description.
SymbolMeaning
N C Number of supply points (equal area grids).
N U Number of UAVs.
N V Number of vessels to be monitored.
C C = c 1 , c 2 ,   , c C is a set of N C supply points, c i is the i t h supply point, i   { 1,2 , 3 ,   ,   N C } , c i C .
U U = u 1 , u 2 ,   , u N U is a set of N U UAVs, u i is the i t h UAV, i   { 1,2 , 3 ,   ,   N U } , u i U .
V V = { v 1 , v 2 ,   , v N V } is the set of N V vessels to be monitored, v i is the i t h vessel to be monitored, i   { 1,2 , 3 ,   ,   N V } , v i   V .
S t a Initial point.
d m i , n j The distance between m i node and n j node, where m i and n j can be the location of the vessel to be monitored or the location of the supply point, m i , n j { V C } .
v i , v j E V v i , v j represents the arc connecting the vessel nodes v i , v j to be monitored, v i , v j V ; E V is the set of arcs between the vessels to be monitored.
c i , c j E C c i , c j represents the arc connecting supply points c i , c j , c i , c j C ; E C is the set of arcs between supply points.
T U Charging time of UAV.
T m a x Maximum endurance of UAV.
T V Time for the UAV monitoring of the vessel.
t ( u i , c j ) The time when UAV u i arrives at supply point c j , u i U , c j C .
t u i , c j M The time when the mothership arrives at the supply point c j to serve the UAV u i , u i U , c j C .
l M The last supply node passes by mothership, l M C .
l u i The last supply node passes by UAV u , l u i C ,   u i U .
S M The speed of the mothership.
S U The speed of the UAV.
p ( u i , v j ) The previous supply point when UAV u i travels to vessel node v j , u i U , v j V .
L ( u i , q j ) The distance traveled by the UAV u i when it is fully charged and travels to q j , u i U , where q j can be the location of the supply point or the location of the vessel q j { V C } .
N ( u i , q j ) V The number of vessels monitored when UAV u i is fully charged and travels to q j , u i U , where q j can be the location of the supply point or the location of the vessel q j V C .
x v j , v k u i = 1 0 A value of 1 means that the UAV u i passes through the arc v j , v k between vessel nodes, 0 means no.
y v j , v k u i = 1   0 A value of 1 means that the UAV u i needs to insert a supply point between vessel nodes v j , v k , 0 means no.
z c i , c j = 1 0 A value of 1 means that the mothership passes through the arc c i , c j between supply point nodes, 0 means no.
w ( u i , v j ) = 1 0 A value of 1 means that the vessel node v j is monitored by the UAV u i , 0 means no.
Table 3. Solution form of MUCWM path planning.
Table 3. Solution form of MUCWM path planning.
UAV monitoring sequence u 1 v 3 v 8 c 1 v 1 v 10
u 2 v 6 v 2 c 2 v 7 c 4 v 9
u 3 v 5 v 4
Mothership sailing sequence M c 2 c 1 c 4
Minimum time consumption of monitoring tasks m i n T t o t a l
Table 4. Data from selected test vessels.
Table 4. Data from selected test vessels.
Serial NumberTypeNameLongitudeLatitudeUpdate Time
1cruise vesselMINFUDINGYU02378122.21036729.8765517:38:30
2freighterZHESHENGYU07203122.36769329.87685717:39:23
3other typesSHUNAN118122.16776529.87763517:39:24
4passenger vesselZHELINGYU96006122.16863329.87823317:38:30
5cruise vesselMINFUDINGYU04373122.24871829.87940317:22:29
6other typesZHELINGYUYUN20079122.31313329.88215317:39:19
7cruise vesselZHEXIANGYUYUN01008122.31312529.88227517:34:58
8passenger vesselZHEPUYU81003122.31258729.88262317:14:21
9other typesMINFUDINGYU09892122.3008829.88293717:24:20
10cruise vesselZHESANYUYUN80399122.30063329.88319817:39:20
11freighterXIN HAI 125122.30081729.883217:33:44
12other types02387-21 51%122.37322329.88388817:11:32
13other typesFUDINGYU07855-6 90%122.30591529.8844716:51:14
14cruise vesselZHEPUYUYUN98788122.28733529.88652217:37:51
15other types07853-25 80%122.21089829.88774516:57:28
16port suppliesRUN ZHOU 7122.18638829.88834517:33:42
17passenger vesselMA ZHI DU 2122.37259229.88837717:39:06
18freighterSUQIYU02669122.37233729.88843717:34:28
19port suppliesPU MAO GONG 4122.3723229.88859317:24:03
20other typesSUQIYU02669-33%122.26059529.89108817:30:34
Table 5. Parameter settings involved in the test.
Table 5. Parameter settings involved in the test.
ParameterValue
Number of UAVs3
Number of equal-area divisions within the sea area (number of candidate supply points)9
Harbor initial point coordinates(153.8975, 172.3765)
UAV maximum flight speed18 m/s
Maximum endurance of UAV45 min
Average monitoring time per vessel3 min
UAV charging time2 min
Mothership maximum sailing speed621 m/min
Number of vessels to be monitored100 vessels
Table 6. Parameterization of improved adaptive differential evolutionary algorithms.
Table 6. Parameterization of improved adaptive differential evolutionary algorithms.
ParameterDefinitionValue
NPNumber of populations100
G_maxThe maximum number of iterations300
FVariation scaling factorAdaptive (initially 1)
CRCrossover probabilityAdaptive (initially 1)
boundsElement upper and lower bounds[0, 1]
Table 7. Partially scaled vessels’ position data.
Table 7. Partially scaled vessels’ position data.
VesselXY
v 1 46.998504593025540.0728363324858847
v 2 133.204317830522033.016281062558292
v 3 181.79624765949187.069408740377959
v 4 167.3523756864698613.041988003445184
v 5 91.4837201703954416.482433590408455
v 90 146.9947346595266178.70608397601762
v 91 148.79173630574027179.34875749784328
v 92 154.04953692653433182.60497000858064
v 93 6.501878683514789188.69751499570967
v 94 20.46319916557688188.9203084832821
Harbor initial point153.8975172.3765
Table 8. Specific results of MUCWM path optimization.
Table 8. Specific results of MUCWM path optimization.
UAV monitoring sequence u 1 Sta→v80→v68→v60→v75→c4→v22→v66→v36→v84→v17→c6→v46→v74→v39→v28→v61→c5→v24→v42→v12→v81→c9→v63→v85→v96→v41→v14→c3→v21→v55→v73→c5→v99→v34→v100→v65→v3→c3→v13→v58→v87→c4
u 2 Sta→v93→v83→v57→v26→v23→v27→c5→v88→v77→v50→v76→v72→c5→v47→v32→v67→v40→c6→v35→v7→v78→v59→v48→c5→v97→v62→v91→v82→v95→v1→c3→v10→v30→v70→v86→c8→v44→v56→v54→c8→v89→v5→c7
u 3 Sta→v49→v51→v11→v20→c5→v45→v2→v53→c2→v4→v6→v71→v90→c7→v64→v69→v92→v29→c3→v18→v43→v25→v52→c4→v9→v94→v15→v38→c4→v33→v19→v31→v37→c8→v79→v98→v16→v8→c7
Mothership sailing sequence M Sta→c4→c5→c6→c2→c7→c6→c9→c3→c5→c3→c4→c3→c5→c8→c7→c4→Sta
Minimum time consumption of monitoring tasks112.25619 min
Table 9. Experimental comparison of the results of IADE and Gurobi.
Table 9. Experimental comparison of the results of IADE and Gurobi.
Problem SizeOptimal Solution (min) CPU Running Time (s)
This ArticleGurobiGAPThis ArticleGurobi
N U = 1 N V = 5 83.2665783.2665701.524.98
N V = 10 93.8387593.8387508.9321.92
N V = 15 100.64764100.64764036.54118.34
N V = 20 115.90341115.903410114.06523.03
N U = 2 N V = 5 69.3349869.3349803.5827.05
N V = 10 80.9303480.93034012.9599.07
N V = 15 99.0637197.436521.67%40.65256.83
N V = 20 105.99024103.952771.96%136.98976.64
N U = 3 N V = 5 66.9834766.9834706.9887.54
N V = 10 75.5469174.969640.77%18.97386.83
N V = 15 84.9644784.863481.19%49.541605.75
N V = 20 93.71631--145.86-
Table 10. Test data parameter settings.
Table 10. Test data parameter settings.
ParameterValue
Number of UAVs3
Number of equal-area divisions within the sea area (number of candidate supply points)25
Harbor initial point coordinates(0, 125)
UAV maximum flight speed1000 m/min
Maximum endurance of UAV35 min
Average monitoring time per vessel2 min
UAV charging time26 m/min
Mothership maximum sailing speed90 vessels
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shen, L.; Sun, J.; Yang, D. Research on Path Optimization for Collaborative UAVs and Mothership Monitoring of Air Pollution from Port Vessels. Sustainability 2024, 16, 4948. https://doi.org/10.3390/su16124948

AMA Style

Shen L, Sun J, Yang D. Research on Path Optimization for Collaborative UAVs and Mothership Monitoring of Air Pollution from Port Vessels. Sustainability. 2024; 16(12):4948. https://doi.org/10.3390/su16124948

Chicago/Turabian Style

Shen, Lixin, Jie Sun, and Dong Yang. 2024. "Research on Path Optimization for Collaborative UAVs and Mothership Monitoring of Air Pollution from Port Vessels" Sustainability 16, no. 12: 4948. https://doi.org/10.3390/su16124948

APA Style

Shen, L., Sun, J., & Yang, D. (2024). Research on Path Optimization for Collaborative UAVs and Mothership Monitoring of Air Pollution from Port Vessels. Sustainability, 16(12), 4948. https://doi.org/10.3390/su16124948

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop