Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

Stochastic optimization model for ship inspection planning under uncertainty in maritime transportation


  • Maritime transportation plays a significant role in international trade and global supply chains. Ship navigation safety is the foundation of operating maritime business smoothly. Recently, more and more attention has been paid to marine environmental protection. To enhance maritime safety and reduce pollution in the marine environment, various regulations and conventions are proposed by international organizations and local governments. One of the most efficient ways of ensuring that the related requirements are complied with by ships is ship inspection by port state control (PSC). In the procedure of ship inspection, a critical issue for the port state is how to select ships of higher risk for inspection and how to optimally allocate the limited inspection resources to these ships. In this study, we adopt prediction and optimization approaches to address the above issues. We first predict the number of ship deficiencies based on a k nearest neighbor (kNN) model. Then, we propose three optimization models which aim for a trade-off between the reward for detected deficiencies and the human resource cost of ship inspection. Specifically, we first follow the predict-then-optimize framework and develop a deterministic optimization model. We also establish two stochastic optimization models where the distribution of ship deficiency number is estimated by the predictive prescription method and the global prescriptive analysis method, respectively. Furthermore, we conduct a case study using inspection data at the Hong Kong port to compare the performances of the three optimization models, from which we conclude that the predictive prescription model is more efficient and effective for this problem.

    Citation: Ran Yan, Ying Yang, Yuquan Du. Stochastic optimization model for ship inspection planning under uncertainty in maritime transportation[J]. Electronic Research Archive, 2023, 31(1): 103-122. doi: 10.3934/era.2023006

    Related Papers:

    [1] Ran Yan, Shuaian Wang . Ship detention prediction using anomaly detection in port state control: model and explanation. Electronic Research Archive, 2022, 30(10): 3679-3691. doi: 10.3934/era.2022188
    [2] Simon Tian, Xinyi Zhu . Data analytics in transport: Does Simpson's paradox exist in rule of ship selection for port state control?. Electronic Research Archive, 2023, 31(1): 251-272. doi: 10.3934/era.2023013
    [3] Xi Luo, Ran Yan, Shuaian Wang, Lu Zhen . A fair evaluation of the potential of machine learning in maritime transportation. Electronic Research Archive, 2023, 31(8): 4753-4772. doi: 10.3934/era.2023243
    [4] Hongyu Zhang, Yiwei Wu, Lu Zhen, Yong Jin, Shuaian Wang . Optimization problems in liquefied natural gas transport and storage for multimodal transport companies. Electronic Research Archive, 2024, 32(8): 4828-4844. doi: 10.3934/era.2024221
    [5] Shuaian Wang, Xuecheng Tian, Ran Yan, Yannick Liu . A deficiency of prescriptive analytics—No perfect predicted value or predicted distribution exists. Electronic Research Archive, 2022, 30(10): 3586-3594. doi: 10.3934/era.2022183
    [6] Liling Huang, Yong Tan, Jinzhu Ye, Xu Guan . Coordinated location-allocation of cruise ship emergency supplies under public health emergencies. Electronic Research Archive, 2023, 31(4): 1804-1821. doi: 10.3934/era.2023093
    [7] Xuecheng Tian, Ran Yan, Shuaian Wang, Yannick Liu, Lu Zhen . Tutorial on prescriptive analytics for logistics: What to predict and how to predict. Electronic Research Archive, 2023, 31(4): 2265-2285. doi: 10.3934/era.2023116
    [8] Yiwei Wu, Yao Lu, Shuaian Wang, Lu Zhen . New challenges in fleet deployment considering EU oil sanctions. Electronic Research Archive, 2023, 31(8): 4507-4529. doi: 10.3934/era.2023230
    [9] Yiwei Wu, Yadan Huang, H Wang, Lu Zhen . Nonlinear programming for fleet deployment, voyage planning and speed optimization in sustainable liner shipping. Electronic Research Archive, 2023, 31(1): 147-168. doi: 10.3934/era.2023008
    [10] Xu Xin, Xiaoli Wang, Tao Zhang, Haichao Chen, Qian Guo, Shaorui Zhou . Liner alliance shipping network design model with shippers' choice inertia and empty container relocation. Electronic Research Archive, 2023, 31(9): 5509-5540. doi: 10.3934/era.2023280
  • Maritime transportation plays a significant role in international trade and global supply chains. Ship navigation safety is the foundation of operating maritime business smoothly. Recently, more and more attention has been paid to marine environmental protection. To enhance maritime safety and reduce pollution in the marine environment, various regulations and conventions are proposed by international organizations and local governments. One of the most efficient ways of ensuring that the related requirements are complied with by ships is ship inspection by port state control (PSC). In the procedure of ship inspection, a critical issue for the port state is how to select ships of higher risk for inspection and how to optimally allocate the limited inspection resources to these ships. In this study, we adopt prediction and optimization approaches to address the above issues. We first predict the number of ship deficiencies based on a k nearest neighbor (kNN) model. Then, we propose three optimization models which aim for a trade-off between the reward for detected deficiencies and the human resource cost of ship inspection. Specifically, we first follow the predict-then-optimize framework and develop a deterministic optimization model. We also establish two stochastic optimization models where the distribution of ship deficiency number is estimated by the predictive prescription method and the global prescriptive analysis method, respectively. Furthermore, we conduct a case study using inspection data at the Hong Kong port to compare the performances of the three optimization models, from which we conclude that the predictive prescription model is more efficient and effective for this problem.



    Maritime transportation is the backbone of globalized trade and the manufacturing supply chain [1,2,3,4], as more than 80% of world merchandise trade by volume is carried by sea, as reported by the United Nations Conference on Trade and Development [5]. Large shipping volumes also render substantial and diverse navigation, safety and marine environmental protection concerns to the high seas, territorial water, and ports and terminals [6,7,8]. According to the International Maritime Organization (IMO), the best way of improving safety at sea is developing international regulations that should be followed by all shipping nations [9,10]. In recent years, environmental sustainability issues in shipping and port activities have become a major policy concern all over the world [11,12,13]. Several mandatory and voluntary regulations are imposed by the IMO [14] as well as local governments [8,15,16]. If a ship's conditions, crew or operations are substantially below the standards required by these regulations, the ship is deemed a substandard ship [17]. A ship's flag state is the nationality or tribe a ship belongs to and under whose laws the ship is plying in international waters [18,19]. Therefore, the ship's flag state is considered the first line of defense against substandard shipping. However, the nature of international shipping leads to the fact that ships may not frequently visit ports that belong to their flag states thus heavily weakening the power of ship flag states in catching substandard ships. As a complement of flag state control, port state control (PSC) is implemented on a regional level, where the port authorities in a certain region are rendered the ability to inspect foreign ships visiting their ports to ensure that they will not cause danger and pollution to their waters.

    Before carrying out one day's PSC inspections at a port, the first thing for the port authority to do is to select those ships with the highest risks for inspection among all the foreign visiting ships. For example, at the Hong Kong port, within the Tokyo Memorandum of Understanding (MoU), a ship risk profile (SRP) under the new inspection regime (NIR) is adopted for high-risk ship selection [20]. Then, each of the selected ships will be inspected by one available ship inspector at the port state, who is also called a PSC officer (PSCO). A selected ship will only be inspected once during this port visit, and one PSCO can inspect several ships per working day. During the onboard inspection, a condition found that does not comply with relevant conventions or regulations is recorded as a deficiency. If one or more very severe deficiencies are found, an intervention action called detention can be taken by the port state to prevent the ship from sailing and causing potential danger and pollution to their waters. As port inspection resources (mainly refering to the available PSCOs) are quite limited, not all visiting ships can be inspected during the normal working period per working day (which is usually 8 hours) of PSCOs [21]. Overtime pay should be made to a PSCO if he/she cannot finish the workload assigned within the normal working period. Meanwhile, ship noncompliance found in inspections can reduce risks and pollution, and thus can be seen as a reward for PSC inspection, where such noncompliance can be represented by ship deficiencies detected. Therefore, the port state needs to make a trade-off between whether to inspect foreign visiting ships and the workload assigned to available PSCOs, so as to maximize the total inspection benefits (i.e., the total number of deficiencies identified) while minimizing the PSCOs' overtime pay.

    There have been some studies on the ship risk prediction problem by adopting empirical analyses [22] and machine learning methods [23,24], which treat this problem as a classification task about whether a ship should be inspected. Recently, some pioneer studies consider combining the ship risk prediction results with PSCO assignment optimization models so as to rationalize the decision of PSC inspection [25]. However, these studies cannot meet the real world ship inspection decision requirements described above for two main reasons: First, they try to optimize ship selection and PSCO assignment decisions with fixed inspection resources without considering the flexible workload assignment and the corresponding inspection cost. Second, most studies follow a predict-then-optimize framework, where the predicted parameters in the following optimization model are treated in a deterministic manner, and the optimization results largely depend on the quality of the prediction models. This study tries to bridge these gaps by integrating prediction with optimization to address the above decision-making problem while considering the trade-off between more detected deficiencies and less overtime pay.

    A k nearest neighbor (kNN) model using real PSC inspection records at the Hong Kong port is first adopted to predict the number of ship deficiencies. Then, three optimization models are established to maximize the total reward from deficiencies identified while minimizing PSCOs' total overtime pay. To be more specific, in the first model, we develop a deterministic optimization model based on the predicted mean of each unknown parameter (i.e., ship deficiency). In the other two models, unknown parameters in the optimization models are modeled as random variables, where the predictive prescription method and the global prescriptive analysis method are adopted to estimate the distributions of the unknown parameters, respectively. We use the Hong Kong port as a case study to construct the kNN model and evaluate the three optimization models. The results show that a kNN model with k=3 can well predict the number of ship deficiencies. Moreover, the optimization model with the predictive prescription method performs best among the three models, and it is suitable for problems where the objective function has a nonlinear term of the unknown parameters.

    According to a review of academic literature conducted by [26], existing studies on PSC can be roughly divided into four categories: factors influencing PSC inspection results, ship selection scheme for inspection, effects of PSC inspection, and suggestions to improve PSC inspection. This research lies in the area of ship selection scheme for inspection, and thus we first review the existing literature within this category.

    Historical inspection records are the main data source used to predict ship risk in the near future or ship inspection performance with the current PSC inspection practice. Based on such historical inspection records as well as external data sources on ship specifications (e.g., ship characteristics and information on ship management parties), empirical analysis is conducted for ship condition prediction. For example, ship risk dimensions are quantified by considering the combination of past incidents and detention information to target high-risk ships. Five combined methods were proposed to realize the combined ship classification methods, which are expected to have predictive power for incidents [22]. Furthermore, the authors quantified risk areas, including collisions, drift groundings, powered grounding, engine-related failures, fire and explosion, hull-related failures, pollution and loss of life, on different ship types such as oil tanker, dry bulk, chemical tanker, general cargo and container ships. The risk areas can be used to guide the assignment of ship inspection priority in PSC inspections [27]. In addition to empirical analyses, machine learning models are also popular for ship risk prediction. Widely used machine learning models include support vector machine, which can be found in [24], [23] and [28] for ship detention prediction; Bayesian network, which can be found in [29] for ship detention probability prediction of bulk carriers and [30] for ship deficiency number prediction; and tree-based model. Multi-target regression random forests were developed to predict ship deficiency conditions in [25]. An XGBoost model incorporating shipping domain knowledge was proposed in [31] for ship deficiency number prediction. A balanced random forest model was used in [32] for ship detention prediction, considering the fact that the un-detained ships highly outnumber detained ships in PSC.

    In more recent years, an emerging research stream in PSC is to rationally allocate inspection resources at a port state for ship inspection, and the inspection resources here mainly refer to PSCOs. One typical example is [25], which aimed to match the expertise of PSCOs with the ships assigned to them, so as to maximize the inspection benefits represented by the total number of deficiencies that can be identified. To achieve the objective, three ship condition prediction models were developed, and the prediction results were combined with PSCO assignment models. Prediction targets in the first prediction model were natural, the numbers of deficiencies under each deficiency category, and the training target was to minimize the prediction error. Prediction targets of the second prediction model were parameters of the decision variables in the following PSCO assignment model. Prediction targets of the third prediction model were also the number of deficiencies under each deficiency category, but the training target was motivated by the structure of the PSCO assignment model to optimize the decision performance. [31] developed an XGBoost model for ship deficiency prediction by considering shipping domain knowledge regarding ship flag performance, recognized organization performance and company performance along with several PSCO scheduling models considering inspection templates and un-dominated inspection templates to maximize the inspection benefit evaluated by the total number of deficiencies. Both studies mentioned above aimed to allocate limited inspection resources at a port authority to maximize inspection benefits. There is another working paper aiming to design inspection resource allocation schemes for ship maintenance for ship operators while considering three types of costs: inspection cost, repair cost and risk cost. A smart predict-then-optimize (SPO) method using an ensemble of SPO trees was developed to incorporate ship maintenance decisions into the prediction of ship detention probability [33].

    The above review shows that when evaluating inspection benefits for port authorities, the benefits only consider ship deficiency number detected, while the costs of conducting an inspection, e.g., the total payment or overtime pay to PSCOs, have not been considered yet. In practice, there is a trade-off between inspecting more ships and reducing inspection costs, and thus taking both factors into account can no doubt increase model applicability. To bridge this gap, optimization models considering more realistic factors, including the practical inspection period onboard, PSCOs' overtime pay, and reward for identifying a ship deficiency, are proposed in this study. Furthermore, in all the ship inspection resource allocation models discussed so far, the optimization models are treated in a deterministic manner where the unknown parameters are the predicted means conditional on the auxiliary data, i.e., historical inspection records and ship specifications. However, due to the fact that the data volume is limited, there might be some missing features, and there are inevitable errors in the data collected. Machine learning models developed for ship condition prediction are not perfect, indicating that the predicted conditional mean is not equal to the actual conditional mean value. To address this issue, this study adopts a kNN model for ship deficiency number prediction and considers data uncertainties by turning the original optimization model into a predictive prescription model based on kNN and a global prescriptive analysis model considering training error based on kNN.

    According to the working process of the PSC authorities, on the morning of each working day, there is a set of foreign visiting ships (denoted by S) coming to the port state that are candidates to be inspected. There is a set of PSCOs (denoted by P) who will be assigned to inspect ships with higher risk. Every detected deficiency on a ship will yield a reward (denoted by c). Practically, the inspection time for a ship sS is considered to be strongly related to its deficiency number ns, because whenever a deficiency, especially a detainable deficiency, is discovered, the PSCO has to take a certain amount of time to provide convincing evidence, as otherwise, it may be appealed by the ship captain or shipping company [34]. Furthermore, the PSCO needs to communicate with the ship captain and crew after detecting a deficiency. Moreover, other categories of deficiencies related to the detected one might also be checked. Therefore, we assume the ship inspection time is linear in its number of deficiencies and can be written as a+bns, where a and b are empirical values gained from historical data. Notably, we can also adopt a more complex function of ship inspection time and deficiency number by taking factors such as ship type and size into account. However, since these factors are known characters and are the same for one ship in all optimization models, they would have similar influences on our models and would not affect the relative performances of the models. So, we do not consider them to quantify ship inspection time in this study for simplicity. The total time initially assigned to each PSCO pP for ship inspection on each working day is tp. We note that if a PSCO spends more time than tp on ship inspection work, there will be overtime pay. Moreover, we assume that the overtime pay is a piecewise linear function of overtime working hours, i.e., 300 HKD/h is required for overtime within 1 hour, 400 HKD/h is for 1–2 hours, and 600 HKD/h is for more than 2 hours. Our model can also accommodate other forms of overtime payment function by modifying the corresponding constraints. The main purpose of the PSC authorities is to maximize the inspection benefits, which consist of the reward for detected deficiencies and the negative inspection costs. Since these two objectives can be converted to the same measurement, i.e., cost, we adopt a parameter λ to transform them into a single objective which aims to minimize the total costs by the weighting method. To summarize, our decision problem is to rationally assign the available PSCOs to the selected ships so as to achieve a trade-off between the total reward for identified deficiencies and the total payment for overtime work of PSCOs.

    The detailed mathematical notations used in our problem are provided below.

    Known parameters:

    S: the set of candidate ships for inspection, indexed by s.

    P: the set of PSCOs to be assigned for ship inspection, indexed by p.

    tp: the total time initially assigned to PSCO p for ship inspection time on each working day.

    a: the basic inspection time of a ship.

    b: the inspection time increment with each deficiency of a ship.

    c: the reward for each detected deficiency.

    λ: a known parameter used to normalize the objective, which lies in (0,1).

    Unknown parameters:

    ns: the number of deficiencies of ship s.

    Decision variables:

    zps: a binary variable that equals 1 if PSCO p is assigned to inspect ship s, or 0 otherwise.

    xp: a continuous variable that indicates the total working time of PSCO p.

    yp: a continuous variable that indicates the overtime pay of PSCO p.

    The model can be written as follows:

    minx,y,zλcpPsSnszps+(1λ)pPyp (3.1a)
    s.t.pPzps1,sS, (3.1b)
    xp=sS(a+bns)zps,pP, (3.1c)
    yp300(xptp),pP, (3.1d)
    yp400(xptp1)+300,pP, (3.1e)
    yp600(xptp2)+300+400,pP, (3.1f)
    xp,yp0,pP, (3.1g)
    zps{0,1},pP,sS. (3.1h)

    Objective (3.1a) maximizes the inspection benefit, which consists of the reward for identifying the (estimated) deficiencies of the ships inspected and the negative overtime pay to PSCOs. Constraint (3.1b) guarantees that each ship is inspected by no more than one PSCO. Constraint (3.1c) presents the total working time of each PSCO. Constraints (3.1d)–(3.1f) specify the overtime pay of each PSCO. Constraints (3.1g)–(3.1h) define the domains of variables.

    In our prediction and optimization approaches, a prediction model is first developed to predict the key unknown parameters (i.e., the number of deficiencies of each ship) in the optimization model by leveraging auxiliary data (i.e., PSC inspection records published by MoUs). Based on the prediction model, three optimization models are then constructed to generate decisions. We note that the overtime pay in the objective is nonlinear with the ship deficiency numbers, which are the unknown parameters predicted by the prediction model.

    A couple of machine learning approaches, such as kNN, decision trees, and random forests, can be used in this context to estimate the distribution of the deficiency numbers of ships [25,31,32]. In this study, we adopt a kNN model for ship deficiency number prediction because it is one of the most intuitive and easy-to-implement machine learning methods for prediction tasks, and it can be easily understood and applied to practical operations management planning in the maritime industry. Moreover, kNN is suitable to generate neighborhood points and to predict the distribution of deficiency numbers in the predictive prescriptive methods [35,36], which will be discussed later. Therefore, we develop a kNN model to predict the number of deficiencies for each ship based on historical records and ship-related features. The data used to predict the deficiency number are introduced in Section 4.1.1. Then, the detailed kNN model construction process is presented in Section 4.1.2.

    The data set used in this study is derived from the Asia Pacific Computerized Information System (APCIS)* supported by the Tokyo MoU and the World Register of Ships (WRS) database. Specifically, APCIS provides detailed inspection records at the Hong Kong port and PSC-related information of the inspected ships within the Tokyo MoU, while WRS provides comprehensive ship-related features: ship construction, engine, dimensions, registration, ownership, fixtures, class, etc. We searched the most frequently used features in current literature from WRS and combined it with the APCIS using ship IMO number to form a unified data set, which contains a total of 4,404 inspection records at the Hong Kong port from 2015 to 2019 with the corresponding ship specification of each record. Based on this data set, 14 features (denoted by w), which are considered to have high correlations with ship deficiency numbers in current literature [25,26,31,32], were used as the input of the kNN model, and the output is the deficiency number of each ship denoted by n. The selected 14 features include ship-related features (e.g., ship age, gross tonnage (GT), length, depth, beam, and type), PSC-related features (e.g., flag performance, recognized organization (RO) performance and company performance over the past three years reported by the Tokyo MoU) and ship historical inspection features (e.g., period since the last inspection, deficiency number of the last inspection, total detentions in all previous inspections, flag change times and casualties in last five years). Then, we randomize the whole data set and divide it into the training set, validation set and test set, containing 60%, 20% and 20% of all records i.e., (2643,881, and 881 records), respectively. Since there are situations where ship deficiencies are extremely large (e.g., 63, 77 and 87), 0.5% records in the training data with the largest ship deficiencies were filtered in order to avoid the influence of outliers on the model.

    *https://apcis.tmou.org/

    http://www.tokyo-mou.org/inspections_detentions/psc_database.php

    To be specific, we list the detailed interpretation and statistical information of the 14 features on the whole data set in Table 1. Since there are a few incomplete records, such as ships with undefined flag performance, RO performance or company performance, empty ship length or depth values, and no last inspection record, we fill the missing values by using the statistical value (e.g., mean, median, and mode) of the corresponding feature in the training set, as is shown in the column Fill-in value for missing values. Moreover, as required by the Tokyo MoU, from the best to the worst, the states of ship flag performance are "White, Grey" and "Black, " and the states of ship RO and company performances are "High, " "Medium, " "Low, " and "Very low, " respectively.

    Table 1.  Feature explanation, encoding method, and descriptive statistics.
    Feature name Meaning Encoding Fill-in value for missing values Mean value1 Max value1 Min value1
    Ship age (years) The time interval between a ship's keel laid date and the current inspection date. No encoding. / 11.42 48.00 0
    GT (100 cubic feet) The internal volume of a ship. No encoding. / 43,388.43 217,612.00 299.00
    Length (meters) The length of a ship. No encoding. 211.19 (mean)2 212.59 400.00 40.75
    Depth (meters) The vertical length of a ship from the top of the keel to the underside of the upper deck. No encoding. 17.5 (mean) 17.52 38.00 3.30
    Beam (meters) The width of a ship. No encoding. / 31.49 63.10 7.80
    Type Ship type searched from the database and according to the annual report on PSC from [37]. One-hot encoding: Bulk carrier: 1 for bulk carrier and 0, otherwise; Container ship: 1 for container ship and 0, otherwise; General cargo/multipurpose: 1 for general cargo/multipurpose and 0, otherwise; Passenger ship: 1 for passenger ship and 0 otherwise; Tanker: 1 for tanker and 0 otherwise; Other: 1 for other ship types and 0 otherwise. / / / /
    Flag performance The flag state performance of a ship calculated by the historical performance of the ships under a flag over the past three years reported in [38]. Label encoding: 'White' 33, 'Grey' 2, 'Black' 1. White (mode) / / /
    RO performance The recognized organization performance of a ship determined by the inspection and detention history of its ships over the last three years [39]. Label encoding: 'High' 4, 'Medium' 3, 'Low' 2, 'Very Low' 1. High (mode) / / /
    Company performance The company performance of a ship based on the Tokyo MoU database considering its ships' inspection performance in a running 36-month period [20]. Label encoding: 'High' 4, 'Medium' 3, 'Low' 2, 'Very Low' 1. Medium (mode) / / /
    The time since the last inspection (months) The time interval between the last inspection date and the current inspection date of a ship. No encoding. 5 (median) 8.56 178 0
    The number of deficiencies in the last inspection The number of deficiencies of a ship identified in the last initial PSC inspection within the Tokyo MoU. No encoding. 0 (mode) 3.28 55 0
    The number of total detentions The total detention times of a ship in all PSC inspections by all PSC authorities. No encoding. / 0.70 19 0
    The number of flag changes The total flag changing times of a ship since its keel laid date. No encoding. / 0.73 8 0
    Casualties in the last five years Whether a ship has encountered casualties in the last five years. One-hot encoding: 1 for any casualty that occurs in the last five years and 0, otherwise. / / / /
    1 : These three columns represent the average, maximal, and minimal values of different features of ships in the entire data set, respectively.
    2 : This indicates that the empty ship length is filled by 211.19, which is the mean value of this feature among all the ships whose length value is not empty in the training set.
    3 : This indicates that state 'White' is encoded to be 3.

     | Show Table
    DownLoad: CSV

    kNN is a powerful non-parametric supervised learning method which was first proposed by [40]. The only hyperparameter in kNN is k, which is the cardinality of the neighborhood considered by the algorithm. As one of the most commonly used machine learning methods, it is widely applied to classification and regression tasks. In this article, we adopt kNN to predict the deficiency number of a ship, which is calculated by the average values of its k nearest neighbors. Moreover, Euclidean distance is used as the distance metric to find the k nearest neighbors.

    Specifically, the training data set contains N inspection records and is denoted by (w1,n1),(w2,n2),,(wN,nN), where wi, i1,,N, is the input feature vector of a ship, and ni is the output value, which is its deficiency number recorded in data record i. The working process of kNN is relatively simple, because the training phase is just to memorize the training set. Then, given a new test sample with feature vector w0, its predicted deficiency number ˆn0 is obtained by looking for the k training samples wi that are the closest to w0 under the defined distance metric (i.e., Euclidean distance) and then averaging their associated targets ni. Mathematically, the prediction function of the target of a new sample w0 using kNN can be presented as follows:

    ˆn0=iNk(w0)nik, (4.1)

    where Nk(w0)={i=1,,N:Nj=1I[w0wiw0wj]k} is the neighborhood set of k data points that are the closest to w0, and N is the number of inspection records in the training set. In addition to being used for ship deficiency number prediction, kNN is also used in prescriptive methods to derive weights in the objective function, which will be discussed later.

    The prediction model is validated using two common metrics for regression problems in machine learning: mean squared error (MSE) and mean absolute error (MAE). Given a total of J samples in the data set, the real output nj and the predicted output ˆnj for a sample j1,,J, the definitions of MSE and MAE are as follows:

    MSE=1JJj=1(ˆnjnj)2, (4.2)
    MAE=1JJj=1|ˆnjnj|. (4.3)

    In M1, we treat each unknown parameter ns in the optimization model as a deterministic value (i.e., the predicted value ˆns). To be specific, we assume that the distribution of ns can be written as

    Pr(ns=ˆns)=1. (4.4)

    Therefore, the optimization model can be reformulated as

    [M1]zDeterministicargminx,y,zλcpPsSˆnszps+(1λ)pPyp, (4.5)
    s.t.xp=sS(a+bˆns)zps,pP, (4.6)
    (3.1b),(3.1d)(3.1h). (4.7)

    In M2, the unknown parameters are modeled as random variables. We adopt the predictive prescription method proposed by [35] to capture the uncertainty, where weight functions of unknown parameters are generated in a data-driven manner based on kNN. First, the distribution of deficiency number ns for ship s with feature vector ws can be written as

    Pr(ns=ni)={1k,ifiNk(ws),0,ifi{1,,N}Nk(ws). (4.8)

    Then, the vector of ship deficiency number for all ships denoted by n would have k|S| possible values, which is an exponential function of S. Thus, we adopt the sample average approximation (SAA) approach and randomly select a set of scenarios U from the Cartesian product of the neighborhood set of every ship {Nk(w0)××Nk(ws)} to approximate the distribution of n. Mathematically, the distribution of n can be written as follows:

    Pr(n=nu)={1|U|,ifuU,0,otherwise, (4.9)

    where the value of n under scenario uU is nu. The decision problem of interest can be reformulated as the following stochastic optimization problem with the above distribution (4.9) of its parameters:

    [M2]zwSAAargminx,y,zλcpPsSˆnszps+1λ|U|uUpPyup (4.10a)
    s.t.pPzps1,sS, (4.10b)
    xup=sS(a+bnus)zps,pP,uU, (4.10c)
    yup300(xuptp),pP,uU, (4.10d)
    yup400(xuptp1)+300,pP,uU, (4.10e)
    yup600(xuptp2)+300+400,pP,uU, (4.10f)
    xup,yup0,pP,uU, (4.10g)
    zps{0,1},pP,sS. (4.10h)

    Note that we employ the predicted deficiency number of each ship ˆns in the first term in the objective function since it provides an exact approximation, as is proved below.

    Proof:

    Since the first term of the objective is linear, we have

    λc|U|uUpPsSnuszps=λcpPsS(uUnus|U|zps).

    If the number of selected scenarios is large enough, and set U is the Cartesian product of the neighborhood set of every ship, then

    λcpPsS(uUnus|U|zps)=λcpPsS(u{Nk(w0)××Nk(ws)}nusk|S|zps)=λcpPsS(iNk(ws)ns,ikzps)=λcpPsSˆnszps.

    Therefore, pPsSˆnszps is an exact approximation of the first term in the objective.

    A global method to estimate the distributions of the unknown parameters relies on the prediction error of each record in the training set. To be specific, the deficiency number of each ship can be first predicted by kNN and then revised by the prediction error in the training set. Therefore, given the auxiliary data ws and its predicted output ˆns, the distribution of ns taking a total of N possible values can be approximated by

    Pr(ns=ˆns+ei)=1N,i{1,,N}. (4.11)

    Then, the total number of scenarios for n would be N|S|. Similar to (4.4), we randomly select a small subset U of all possible scenarios. We denote eu as the error between the predicted and actual parameter values under scenario uU. Thus, the distribution function of n can be written as follows:

    Pr(n=ˆn+eu)={1|U|,ifuU,0,otherwise. (4.12)

    Based on this distribution function, the stochastic optimization model can be written as

    [M3]zGlobal,SAAerrargminx,y,zλcpPsSˆnszps+1λ|U|uUpPyup, (4.13)

    subject to constraints (4.10b)–(4.10h).

    In this section, we use the Hong Kong port as a case study to construct the kNN model and evaluate the three optimization models.

    For the kNN model, the only hyperparameter used to control the learning process is k. Theoretically, when the value of k is smaller, the training error will decrease, but the model is more complex, and the generalization ability might be weakened. Thus, k is empirically set in the range [1,8] with an interval of 1, from which the optimal k is chosen by considering the two metrics (i.e., MSE and MAE) on the validation set.

    The values of other parameters in the optimization models are set according to the real situation at the Hong Kong port and are presented in Table 2. We assume that 20 foreign ships come to the port in one day, and 4 PSCOs are available for inspection. The 20 ships that are candidates for inspection are randomly selected from the test set, and the number of deficiencies of them is first predicted by the kNN model. The basic inspection time for a ship is set as 1 hour, and the inspection time will increase by 0.25 hours whenever there is one more deficiency detected. Moreover, each PSCO is assumed to be the same, and the initially assigned time to PSCO p for inspecting the ships tp is equally set to 8 hours. In addition, in order to get the reward for each detected deficiency, we first calculate a PSCO's hourly wage by dividing the average monthly salary of a PSCO (i.e., 85,870 HKD) by the working hours in a month (i.e., 176 hours under the assumption that a PSCO works 22 days a month except weekends and works 8 hours a day), and the result is 488 HKD/h. Then, we can get the average time spent for inspecting one deficiency as (1+0.25×4.8)/4.8=0.46 h, which is calculated by the average deficiency number of one inspection and the linear function of the inspection time. Thus, the cost of inspecting each deficiency is 488×0.46= 225 HKD, and the reward for each detected deficiency (i.e., c) is set to twice of it as 450 HKD. Furthermore, we set the values of |U| as 300,500, and 1000, respectively, so as to explore the impacts of the selected scenarios in M2 and M3. To balance the reward for deficiency detection and PSCOs' overtime pay, parameter λ is set to 0.3, 0.5 and 0.8, respectively.

    https://www.mardep.gov.hk/theme/maritime_industry/en/marine_officer.html

    Table 2.  Settings of parameters related to ship inspection at the Hong Kong Port.
    Parameter Value
    |S| 20
    |P| 4
    tp 8 h
    a 1
    b 0.25
    the average monthly salary of a PSCO 85,870 HKD
    the average deficiency number of one inspection1 4.8
    c 450
    λ 0.3, 0.5, 0.8
    |U| 100,200,300
    1 : This is calculated on the whole data set.

     | Show Table
    DownLoad: CSV

    The experiments are performed on a computer with an Intel Core i5 CPU and 8 GB memory under the Mac operating system. The models are implemented in the Python programming language using Gurobi 9.1.2 as the solver.

    Following the procedure introduced in Section 4.1.2, we developed kNN models with k varying from 1 to 8 and recorded the MSE and MAE on the validation set. The results are shown in Table 3. We can find the optimal k is 3, as it leads to the minimal MSE (34.65) and MAE (3.35). The computation times of different k values are similar and are below 0.4 seconds.

    Table 3.  Prediction performances of kNN with different k values on the validation set.
    k MSE MAE Computation time (s)
    1 42.37 3.71 0.23
    2 35.41 3.42 0.35
    3 34.65 3.35 0.25
    4 35.33 3.40 0.21
    5 36.50 3.46 0.21
    6 36.03 3.47 0.23
    7 36.09 3.46 0.24
    8 35.75 3.45 0.22

     | Show Table
    DownLoad: CSV

    For comparison, we also develop a support vector regression (SVR) model, which is suitable for nonlinear regression problems and is extended from a support vector machine for classification problems [41]. The kernel function is the most significant factor that affects the performance of the SVR. Therefore, we tune the kernel functions by choosing one from radial basis function ('rbf'), polynomial ('poly') and sigmoid, respectively, which are popular kernel functions in SVM [42]. The other hyperparameters are set as the default values in sklearn for consistency. The SVR models are tested on the same data sets as kNN, and the prediction performances on the training and validation set are shown in Table 4. We can observe that 'rbf' performs best among the three kernel functions, with lower MSE and MAE on both training and validation sets. Therefore, we choose 'rbf' as the kernel function of SVR. Moreover, we can also observe that the MSE and MAE of the SVR model with kernel function 'rbf' on the validation set are still larger than those of the kNN model with k=3. In addition, the computation times of both the kNN and SVR models are small.

    Table 4.  Prediction performances of SVR with different kernel functions on the training and validation sets.
    Kernel functions Training set Validation set
    MSE MAE Computation time (s) MSE MAE
    'rbf' 24.83 2.96 0.27 41.64 3.43
    'poly' 30.46 3.27 0.29 49.35 3.63
    'sigmoid' 9975.59 71.88 0.56 10402.97 75.27

     | Show Table
    DownLoad: CSV

    Given the chosen k=3 of kNN and kernel function 'rbf' of SVR, we divided the whole test set into 10 groups and ran each group to record its MSE and MAE. The minimum, maximum, mean and variance of MSE and MAE on the 10 groups of test sets are shown in Table 5. It can be seen that the mean values of MSE and MAE of kNN on the test sets are 30.42 and 3.32, respectively, which are slightly less than those on the validation set. By comparing the results of the two models, we can find that all measurements of kNN are better than SVR with respect to MSE. For MAE, the mean and minimal values of kNN are slightly larger than those of SVR, while the maximal value and standard derivation of kNN are smaller than those of SVR, which indicates that kNN is more stable for all test sets. To summarize, the kNN model with k=3 is effective and stable in predicting the deficiency number of a ship. We note that more complex models such as tree-based models can be built to predict the ship deficiency number and may perform better on the data sets by carefully adjusting model hyperparameters. However, these models with various intricate hyperparameters are not as straightforward as the kNN model and cannot be easily understood by the port officials. Consequently, they may be reluctant to use these models.

    Table 5.  Prediction performances of kNN and SVR on the test set.
    Model Metric Min Max Mean Standard deviation
    kNN MSE 17.60 48.61 30.42 10.98
    MAE 2.77 3.72 3.32 0.29
    SVR MSE 9.47 52.75 32.96 16.05
    MAE 2.11 3.96 3.11 0.49

     | Show Table
    DownLoad: CSV

    Based on the kNN prediction model, we compared the performances of the three optimization models introduced in Section 4.2. We first randomly selected 10 groups of instances, each of which consists of 20 records from the test set. Then, we tested the three optimization models on these instances in turn and set the maximal computation time limit for each instance as 600 seconds. In addition, to avoid an unreasonable number (i.e., negative or too large of a number) of ship deficiencies, we set a range [0,ˉns+3×σ(ns)] for the possible value of ns in M3, where ˉns and σ(ns) are the mean and standard deviation values of the ship deficiency number in the training set, respectively. Specifically, if the calculated number of ship deficiencies is larger than the upper bound, which is a statistical measurement commonly used when the data set is large enough, we set it as the upper bound; similarly, if it is smaller than zero, we set it as zero, as ship deficiency number cannot be negative. The comparison results are shown in Table 6. Notably, the optimal decisions can be made if the deficiency number of each ship is predicted accurately in advance (denoted by "perfect-forecast policy"). However, this is an ideal situation that cannot be achieved in practice. We use the objective of the perfect-forecast policy as a lower bound in theory to compare the performances of the three proposed optimization models.

    Table 6.  Comparison of computational performances of M1, M2 and M3.
    Performance λ M1 |U|=100 |U|=200 |U|=300
    M2 M3 M2 M3 M2 M3
    Average total computation time (S) 0.3 0.05 19.30 155.60 176.50 514.89 323.99 565.70
    0.5 0.03 41.11 422.34 401.19 600.00 520.17 600.00
    0.8 0.01 39.90 427.12 382.68 563.32 600.00 600.00
    Standard deviation of computation time (S) 0.3 0.02 13.78 175.17 189.52 151.16 222.65 77.00
    0.5 0.02 36.32 196.32 192.48 0.00 141.01 0.00
    0.8 0.01 33.57 248.10 186.83 121.67 0.00 0.00
    Average number of detected deficiencies 0.3 83.30 82.60 81.80 83.20 80.90 83.80 80.70
    0.5 95.90 96.40 91.10 96.40 92.40 96.40 91.10
    0.8 101.10 101.10 101.10 101.10 101.10 101.10 101.10
    Average overtime pay 0.3 2,560.00 2,252.50 2,507.50 2,115.00 2,362.50 2,480.00 2,262.50
    0.5 5,037.50 4,995.00 3,990.00 5,122.50 4,367.50 5,322.50 4,357.50
    0.8 6,525.00 6,362.50 6,410.00 6,325.00 6,302.50 6,560.00 6,615.00
    Average absolute objective value1 0.3 9,453.50 9,574.25 9,287.75 9,751.50 9,267.75 9,577.00 9,310.75
    0.5 19,058.75 19,192.50 18,502.50 19,128.75 18,606.25 19,028.75 18,318.75
    0.8 35,091.00 35,123.50 35,114.00 35,131.00 35,135.50 35,084.00 35,073.00
    Average gap from the perfect-forecast policy 0.3 17.92% 16.80% 19.56% 15.46% 19.59% 16.84% 19.36%
    0.5 6.98% 6.34% 9.11% 6.69% 8.63% 7.07% 10.02%
    0.8 1.04% 0.94% 0.94% 0.90% 0.89% 1.03% 1.11%
    1 : To avoid the minus sign for simplicity, this row shows the average absolute values of the objectives, where larger values are better.

     | Show Table
    DownLoad: CSV

    From Table 6, we can observe that much less time is required to solve M1 compared to that of M2 and M3 in all cases. The computation time of M2 is also much smaller than M3 under the same parameter settings. Moreover, the average computation times of both M2 and M3 would increase greatly as expected as the number of selected scenarios |U| increases. Most instances in M2 and M3 would reach the time limit (i.e., 600 seconds), when |U| increased to 200 and 300, respectively. The standard deviation of computation time for different models also vary a lot. While the standard deviation of computation time in M1 is quite small, it is much larger in M2 and M3. This is because more time is required to solve M2 and M3, and small changes in parameters can greatly affect the computation time. When the number of selected scenarios |U| increases, there is also a rapid upward trend in the standard deviation of computation time in M2 and M3. Particularly, in M2 and M3, when the value of |U| increases to a certain level that the solution times of all instances reach the time limit, the standard deviation of computation time is about zero. In addition, there is no obvious relationship between the value of λ and the computation time.

    Taking the perfect-forecast policy as a benchmark, we can further compare the performances of the three models under different parameter settings as shown in the last row of Table 6. When the value of |U| is fixed, and the value of λ increases, the performances of the three models improve, as the gaps between the objectives and the optimal policies drop from nearly 20% to about 1%. This verifies our intuition that the larger the coefficient of the linear term of deficiency number is, the better the results the models can provide. However, there is no distinct downtrend of the gaps with the increasing value of |U| as expected. Two factors may contribute to this result: First, although the value of |U| increases, the sampling scenarios are still very small compared to the possible values of ship deficiency numbers, which would reach k|S| or N|S|; second, too many sampling scenarios would increase the computational difficulty of the model and may lead to a sub-optimal solution within the limited time. In addition, M2 outperforms M1 in most cases (except when λ=0.5 and |U|=300) and is also much better than M3 (except when λ=0.8 and |U|=200). Especially, when the value of λ is 0.8, the average gap between the perfect-forecast policy and M2 is slightly smaller than that of M1. However, the average gap between the perfect-forecast policy and M2 is much smaller than that of M1 when the value of λ decreases to 0.3, as the proportion of the nonlinear term of uncertain parameters is larger. Meanwhile, M3 always fails to provide a satisfactory result compared to M1 and M2, which may be because of the limited sampling number compared to the extremely large number of possible scenarios in M3.

    Moreover, λ is a parameter used to achieve a trade-off between the reward for detected deficiency number and the cost of overtime pay. When λ increases from 0.3 to 0.8, the total number of detected deficiencies shows an increasing trend, as suggested by the optimal solution, since the weight of reward for each detected deficiency becomes larger. Then, a longer time is also needed to detect such deficiencies, and thus the overtime pay of PSCOs for ship inspection is also increased. Overall, there is a downward trend of the objective value, i.e., the total ship inspection cost, as λ increases, because the weight of reward for detecting the ship deficiencies increases.

    To summarize, the average model computation times and their standard deviation increase in the order of M1, M2 and M3. In addition, the model computation times of M2 and M3 are sensitive to the value of |U|, because it increases the number of constraints in the optimization model. In almost all cases, M2 outperforms M1 and M3 and has smaller gaps between the objective achieved and the perfect-forecast policy. Moreover, the gaps between M2/M3 and the perfect-forecast policy could decrease if there is a larger number of selected scenarios and no computation time limitation. We can therefore conclude that M2 is the most effective, stable and flexible among the three models. This reveals the benefits of a predictive prescription approach toward a model where the objective has a nonlinear term of the unknown parameters (the deficiency numbers of ships).

    Maritime safety and marine environment protection are two main topics in maritime study that have been gaining wide attention over the world for many years, which can be seen from the increasing number of and more and more stringent regulations and conventions implemented. PSC inspection is an effective safety net to catch substandard ships that do not comply with the requirements in the regulations and conventions. To improve the efficiency of PSC inspection, higher-risk ship selection and PSCO resources assignment are addressed as essential problems. Whereas previous studies model this problem in a deterministic manner to maximize the reward for detected ship deficiencies based on the prediction model of ship deficiency number, this study adopts a more comprehensive objective, which considers not only the benefits of identified deficiencies, but also the cost of ship inspection. Based on this target, a deterministic model (M1), and two stochastic models (i.e., M2, based on the predictive prescription method, and M3, based on the global prescriptive analysis method) were established using a kNN prediction model as the predictor of ship deficiency number. Extensive numerical experiments were conducted using inspection data at the Hong Kong port to evaluate the performances of the three models and study the impacts of significant parameters.

    To be specific, we tested the performance of the kNN model with different values of k, from which the optimal k=3 was chosen by minimizing two metrics for regression (i.e., MSE and MAE) on the validation set. Based on this kNN model, we compared the average performances of the three optimization models with ten groups of randomly generated instances from the test sets. The results show that M2 stands out by the minimal average gap from the perfect-forecast policy and the relatively small average total computation time. Meanwhile, the results also show that M2 is particularly suitable for the problem, where the objective is nonlinear with the predicted uncertain parameters. By contrast, M1, which is widely used in current studies, bears larger optimality gaps than M2 in almost all cases, although it has the least average total computation time among the three models. In addition, M3 cannot perform well in those cases, which may be due to the use of all training data errors in estimating the distribution of the number of ship deficiencies, the issue of over-fitting to the predictions and the limit of computation time. However, there is still an improvement possibility for M3 if we design more suitable filter methods for the errors used to estimate the distribution of unknown parameters. Moreover, by conducting sensitivity analysis on important parameters, we can observe that there is a falling trend in the average gap between the perfect-forecast policy as the value of λ increases in all of the three models, and there is an increase in the average total computation time for M2 and M3 when the value of |U| becomes larger.

    The main contributions of this paper lie in the following aspects:

    1) Taking an important policy in the maritime industry into account, we propose prediction and optimization approaches which consider the trade-off between detecting more deficiencies and reducing overtime workload. This helps port states to identify high-risk ships effectively and assign the PSCOs more flexibly. Therefore, the main objectives of PSC to eliminate substandard shipping and safeguard the sea can be enhanced.

    2) Theoretically, we implement both a predictive prescription model and a global prescriptive analysis model based on kNN for ship deficiency number prediction and inspection resource assignment. To the best of our knowledge, this study is the pioneer in integrating prediction with optimization in the ship inspection application and can be easily understood and adopted by port authorities for practical ship inspection planning.

    3) In the numerical tests, we conduct extensive numerical experiments using inspection data at the Hong Kong port to compare M1, based on the traditional deterministic model, with M2, based on the predictive prescription method, and M3, based on the global prescriptive analysis method. The experiment reveals the effectiveness of a predictive prescription approach towards ship inspection planning where the reward for ship deficiencies detected and cost of PSCOs' overtime pay are taken into account.

    However, there are also some limitations of our work. First, to achieve a trade-off between detecting more deficiencies and reducing inspection cost, we only use a weighting method to unify these two objectives without developing other models. Second, the performance of the global prescriptive analysis model M3 cannot perform well as expected, since the predicted distribution function of ship deficiency number is revised by the prediction errors on the whole training set. For example, if the prediction error of one training sample is too large, revising the prediction result with this error may result in an unreasonable predicted ship deficiency number. Therefore, future research could be conducted from the following aspects: From the modeling perspective, we could adopt a dual-objective optimization model or a two-stage optimization model to achieve a trade-off between the total reward for deficiencies identified and the total overtime pay; from the methodological point of view, more scenario selection strategies could be adopted to predict the distribution function of M3 and thus improve the effectiveness of the model.

    The authors declare there are no conflicts of interest.



    [1] C. Li, Y. Xie, G. Wang, X. Zeng, H. Jing, Lateral stability regulation of intelligent electric vehicle based on model predictive control, J. Intell. Connected Veh., 4 (2021), 104–114. https://doi.org/10.1108/JICV-03-2021-0005 doi: 10.1108/JICV-03-2021-0005
    [2] S. Wang, D. Zhuge, L. Zhen, C. Y. Lee, Liner shipping service planning under sulfur emission regulations, Transp. Sci., 55 (2021), 491–509. https://doi.org/10.1287/trsc.2020.1010 doi: 10.1287/trsc.2020.1010
    [3] L. Zhen, Y. Hu, S. Wang, G. Laporte, Y. Wu, Fleet deployment and demand fulfillment for container shipping liners, Transp. Res. Part B Methodol., 120 (2019), 15–32. https://doi.org/10.1016/j.trb.2018.11.011 doi: 10.1016/j.trb.2018.11.011
    [4] K. Wang, S. Wang, L. Zhen, X. Qu, Cruise service planning considering berth availability and decreasing marginal profit, Transp. Res. Part B Methodol., 95 (2017), 1–18. https://doi.org/10.1016/j.trb.2016.10.020 doi: 10.1016/j.trb.2016.10.020
    [5] S. N. Sirimanne, J. Hoffman, W. Juan, R. Asariotis, M. Assaf, G. Ayala, et al., Review of Maritime Transport 2019, UNCTAD, New York, USA, 2019.
    [6] N. Lyu, Y. Wang, C. Wu, L. Peng, A. F. Thomas, Using naturalistic driving data to identify driving style based on longitudinal driving operation conditions, J. Intell. Connected Veh., 5 (2021), 17–35. https://doi.org/10.1108/JICV-07-2021-0008 doi: 10.1108/JICV-07-2021-0008
    [7] W. Yi, L. Zhen, Y. Jin, Stackelberg game analysis of government subsidy on sustainable off-site construction and low-carbon logistics, Cleaner Logist. Supply Chain, 2 (2021), 100013. https://doi.org/10.1016/j.clscn.2021.100013 doi: 10.1016/j.clscn.2021.100013
    [8] D. Huang, S. Wang, A two-stage stochastic programming model of coordinated electric bus charging scheduling for a hybrid charging scheme, Multimodal Transp., 1 (2022), 100006. https://doi.org/10.1016/j.multra.2022.100006 doi: 10.1016/j.multra.2022.100006
    [9] IMO, Maritime Safety, 2019. Available from: https://www.imo.org/en/OurWork/Safety/Pages/default.aspx.
    [10] L. Wu, Y. Adulyasak, J. F. Cordeau, S. Wang, Vessel service planning in seaports, Oper. Res., 70 (2022), 2032–2053. https://doi.org/10.1287/opre.2021.2228 doi: 10.1287/opre.2021.2228
    [11] L. Zhen, Y. Wu, S. Wang, G. Laporte, Green technology adoption for fleet deployment in a shipping network, Transp. Res. Part B Methodol., 139 (2020), 388–410. https://doi.org/10.1016/j.trb.2020.06.004 doi: 10.1016/j.trb.2020.06.004
    [12] S. Wang, R. Yan, A global method from predictive to prescriptive analytics considering prediction error for "predict, then optimize" with an example of low-carbon logistics, Cleaner Logist. Supply Chain, 4 (2022), 100062. https://doi.org/10.1016/j.clscn.2022.100062 doi: 10.1016/j.clscn.2022.100062
    [13] J. Qi, S. Wang, H. Psaraftis, Bi-level optimization model applications in managing air emissions from ships: A review, Commun. Transp. Res., 1 (2021), 100020. https://doi.org/10.1016/j.commtr.2021.100020 doi: 10.1016/j.commtr.2021.100020
    [14] IMO, Initial IMO GHG Strategy, 2019. Available from: https://www.imo.org/en/MediaCentre/HotTopics/Pages/Reducing-greenhouse-gas-emissions-from-ships.aspx.
    [15] R. Yan, S. Wang, Integrating prediction with optimization: Models and applications in transportation management, Multimodal Transp., 1 (2022), 100018. https://doi.org/10.1016/j.multra.2022.100018 doi: 10.1016/j.multra.2022.100018
    [16] L. Zhen, Y. Wu, S. Wang, Y. Hu, W. Yi, Capacitated closed-loop supply chain network design under uncertainty, Adv. Eng. Inf., 38 (2018), 306–315. https://doi.org/10.1016/j.aei.2018.07.007 doi: 10.1016/j.aei.2018.07.007
    [17] IMO, Resolution A. 1155(32) Adopted on 15 December 2021 (Agenda items 12 and 14) Procedures for port state control, 2021, 2022. Available from: https://www.register-iri.com/wp-content/uploads/A.115532.pdf.
    [18] S. Wang, L. Zhen, D. Zhuge, Dynamic programming algorithms for selection of waste disposal ports in cruise shipping, Transp. Res. Part B Methodol., 108 (2018), 235–248. https://doi.org/10.1016/j.trb.2017.12.016 doi: 10.1016/j.trb.2017.12.016
    [19] W. Yi, S. Wu, L. Zhen, G. Chawynski, Bi-level programming subsidy design for promoting sustainable prefabricated product logistics, Cleaner Logist. Supply Chain, 1 (2021), 100005. https://doi.org/10.1016/j.clscn.2021.100005 doi: 10.1016/j.clscn.2021.100005
    [20] Tokyo MoU, Information sheet of the new inspection regime (NIR), 2014. Available from: http://www.tokyo-mou.org/doc/NIR-information%20sheet-r.pdf.
    [21] C. Chai, Z. Zhou, W. Yin, H. David, S. Zhang, Evaluating the moderating effect of in-vehicle warning information on mental workload and collision avoidance performance, J. Intell. Connected Veh., 5 (2022), 49–62. https://doi.org/10.1108/JICV-03-2021-0003 doi: 10.1108/JICV-03-2021-0003
    [22] C. Heij, S. Knapp, Shipping inspections, detentions, and incidents: an empirical analysis of risk dimensions, Marit. Policy Manage., 46 (2019), 866–883. https://doi.org/10.1080/03088839.2019.1647362 doi: 10.1080/03088839.2019.1647362
    [23] R. Xu, Q. Lu, K. Li, W. Li, Web mining for improving risk assessment in port state control inspection, in 2007 International Conference on Natural Language Processing and Knowledge Engineering, IEEE, (2007), 427–434. https://doi.org/10.1109/NLPKE.2007.4368066
    [24] R. Xu, Q. Lu, W. Li, K. Li, H. S. Zheng, A risk assessment system for improving port state control inspection, in 2007 International Conference on Machine Learning and Cybernetics, IEEE, 2 (2007), 818–823. https://doi.org/10.1109/ICMLC.2007.4370255
    [25] R. Yan, S. Wang, K. Fagerholt, A semi-"smart predict then optimize"(semi-SPO) method for efficient ship inspection, Transp. Res. Part B Methodol., 142 (2020), 100–125. https://doi.org/10.1016/j.trb.2020.09.014 doi: 10.1016/j.trb.2020.09.014
    [26] R. Yan, S. Wang, Ship inspection by port state control—review of current research, in Smart Transportation Systems 2019, Springer, (2019), 233–241. https://doi.org/10.1007/978-981-13-8683-1_24
    [27] S. Knapp, C. Heij, Improved strategies for the maritime industry to target vessels for inspection and to select inspection priority areas, Safety, 6 (2020), 18. https://www.mdpi.com/2313-576X/6/2/18
    [28] S. Wu, X. Chen, C. Shi, J. Fu, Y. Yan, S. Wang, Ship detention prediction via feature selection scheme and support vector machine (svm), Marit. Policy Manage., 49 (2022), 140–153. https://doi.org/10.1080/03088839.2021.1875141 doi: 10.1080/03088839.2021.1875141
    [29] Z. Yang, Z. Yang, J. Yin, Realising advanced risk-based port state control inspection using data-driven bayesian networks, Transp. Res. Part A Policy Pract., 110 (2018), 38–56. https://doi.org/10.1016/j.tra.2018.01.033 doi: 10.1016/j.tra.2018.01.033
    [30] S. Wang, R. Yan, X. Qu, Development of a non-parametric classifier: Effective identification, algorithm, and applications in port state control for maritime transportation, Transp. Res. Part B Methodol., 128 (2019), 129–157. https://doi.org/10.1016/j.trb.2019.07.017 doi: 10.1016/j.trb.2019.07.017
    [31] R. Yan, S. Wang, J. Cao, D. Sun, Shipping domain knowledge informed prediction and optimization in port state control, Transp. Res. Part B Methodol., 149 (2021), 52–78. https://doi.org/10.1016/j.trb.2021.05.003 doi: 10.1016/j.trb.2021.05.003
    [32] R. Yan, S. Wang, C. Peng, An artificial intelligence model considering data imbalance for ship selection in port state control based on detention probabilities, J. Comput. Sci., 48 (2021), 101257. https://doi.org/10.1016/j.jocs.2020.101257 doi: 10.1016/j.jocs.2020.101257
    [33] X. Tian, R. Yan, Y. Liu, S. Wang, A smart predict-then-optimize method for targeted and cost-effective maritime transportation, Transp. Res. Part B Methodol., 142 (2020), 100–125. https://doi.org/10.1016/j.trb.2020.09.014 doi: 10.1016/j.trb.2020.09.014
    [34] IMO, Procedures for Port State Control, 2019. Available from: https://www.imo.org/en/OurWork/IIIS/Pages/Port%20State%20Control.aspx.
    [35] D. Bertsimas, N. Kallus, From predictive to prescriptive analytics, Manage. Sci., 66 (2020), 1025–1044. https://doi.org/10.1287/mnsc.2018.3253 doi: 10.1287/mnsc.2018.3253
    [36] L. Galli, T. Levato, F. Schoen, L. Tigli, Prescriptive analytics for inventory management in health care, J. Oper. Res. Soc., 72 (2021), 2211–2224. https://doi.org/10.1080/01605682.2020.1776167 doi: 10.1080/01605682.2020.1776167
    [37] Tokyo MoU, Annual report 2016 on port state control in the Asia-Pacific region, 2017. Available from: http://www.tokyo-mou.org/doc/ANN16.pdf.
    [38] Tokyo MoU, Black–Grey–White lists, 2017. Available from: http://www.tokyo-mou.org/doc/Flag%20performance%20list%202017.pdf.
    [39] Paris MoU, Criteria for responsibility assessment of recognized organizations (RO), 2013. Available from: https://www.parismou.org/sites/default/files/RO%20responsibility%20rev11_0.pdf.
    [40] E. Fix, J. L. Hodges, Discriminatory analysis: nonparametric discrimination, consistency properties, Int. Stat. Rev., 57 (1989), 238–247. https://doi.org/10.2307/1403797 doi: 10.2307/1403797
    [41] A. J. Smola, B. Schölkopf, A tutorial on support vector regression, Stat. Comput., 14 (2004), 199–222. https://doi.org/10.1023/B:STCO.0000035301.49549.88 doi: 10.1023/B:STCO.0000035301.49549.88
    [42] H. Bhavsar, M. H. Panchal, A review on support vector machine for data classification, Int. J. Adv. Res. Comput. Eng. Technol., 1 (2012), 185–189.
  • This article has been cited by:

    1. Ran Yan, Shuaian Wang, Lu Zhen, Shuo Jiang, Classification and regression in prescriptive analytics: Development of hybrid models and an example of ship inspection by port state control, 2024, 163, 03050548, 106517, 10.1016/j.cor.2023.106517
    2. Ewa Chodakowska, Darius Bazaras, Edgar Sokolovskij, Veslav Kuranovic, Leonas Ustinovichius, TRANSPORT RISKS IN THE SUPPLY CHAINS – POST COVID-19 CHALLENGES, 2024, 25, 1611-1699, 211, 10.3846/jbem.2024.21110
    3. Hongyan Dui, Huiting Xu, Haohao Zhou, Reliability analysis and resilience measure of complex systems in shock events, 2023, 31, 2688-1594, 6657, 10.3934/era.2023336
    4. Xu Zhan, Yang Yong, Wang Xiao, Phased mission reliability analysis of unmanned ship systems, 2023, 31, 2688-1594, 6425, 10.3934/era.2023325
    5. Xizi Qiao, Ying Yang, King-Wah Pang, Yong Jin, Shuaian Wang, Ship Selection and Inspection Scheduling in Inland Waterway Transport, 2024, 12, 2227-7390, 2327, 10.3390/math12152327
    6. Ying Yang, Ran Yan, Shuaian Wang, Prescriptive analytics models for vessel inspection planning in maritime transportation, 2024, 190, 03608352, 110012, 10.1016/j.cie.2024.110012
    7. Mahdi Jahangard, Ying Xie, Yuanjun Feng, Leveraging machine learning and optimization models for enhanced seaport efficiency, 2025, 1479-2931, 10.1057/s41278-024-00309-w
    8. Qi Hong, Mo Jia, Xuecheng Tian, Zhiyuan Liu, Shuaian Wang, A Surrogate Piecewise Linear Loss Function for Contextual Stochastic Linear Programs in Transport, 2025, 13, 2227-7390, 2033, 10.3390/math13122033
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2221) PDF downloads(130) Cited by(8)

Figures and Tables

Tables(6)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog