一种基于AIS数据的船舶航行特征及偏好分析方法与流程

本发明属于船舶交通领域，涉及一种基于ais数据的船舶航行特征及偏好分析方法。
背景技术：
：船载的船舶自动识别系统(automaticidentificationsystem，ais)可以每隔几秒钟向岸上基站或卫星发送船舶的ais数据，ais数据包括三种信息类型，即静态信息、动态信息和航行相关信息。静态信息包括船名、船舶类型、船舶尺度、货物种类等；动态信息包括时间、位置、船舶速度和船舶航向等；航行相关信息包括吃水、目的地等。航行相关信息由船上人员手动更新，数据可靠度不高，在ais数据分析中，一般不考虑航行相关信息。船舶ais数据的动态信息与静态信息是分离的，因此对于船舶ais数据分析方法，大多通过分析船舶ais数据的动态信息，为后续研究奠定基础，如船舶碰撞风险分析、船舶航迹预测。未见对于船舶ais数据的动态信息和静态信息相结合的船舶航行特征及偏好分析方法。技术实现要素：为了实现上述目的，本发明提供一种基于ais数据的船舶航行特征及偏好分析方法。通过对船舶ais数据中代表船舶自身属性的静态信息与代表船舶航行状态的动态信息相结合进行数据分析，明了直观地分析了拥有不同静态信息的船舶在航道中的航行特征及航行偏好，确定船舶的航行特征可以为航道、锚地通过能力仿真动态模拟提供输入参数，从而为航道规划、船舶组织调度提供定量参考。此外，了解船舶的航行偏好可以为海事管理提供参考，如对偏离航行偏好的船舶给予重点关注，从而降低航行安全风险。本发明的技术方案为：一种基于ais数据的船舶航行特征及偏好分析方法，包括以下步骤：步骤1，统计船舶ais数据在航道中任意位置划定长度200m、宽度为航道宽度(若航道外有船舶航行，则宽度为航道宽度基础上再向外延伸0.5倍航道宽度)的矩形区域作为船舶ais数据统计区，若同一艘船舶单次航程在统计区内有多次ais数据，保留任意一次ais数据即可，统计得到一个时间范围内(如一季、一年)通过ais数据统计区的所有船舶的ais数据。步骤2，确定船舶ais数据中的静态信息和动态信息在静态信息中，因货物种类、船名可能发生改变，因此静态信息只考虑船舶类型和船长。在动态信息中，因航道多为直水道，不考虑船舶航向，因此动态信息只考虑时间、船舶位置和船速。步骤3，确定不同船舶类型的船速特征ais数据中的船舶类型未对货船进一步划分，如集装箱船、散货船等，因此船舶类型分为货船、油船和客船。根据货船ais数据中的时间数据和船速数据，确定货船在日间(6:00～18:00)的船速最大值和最小值，将最值船速区间等分为多个速度区间(速度区间数量根据最值船速区间大小确定，无固定标准)，统计得到货船在日间(6:00～18:00)的船速分布图并拟合得到船速的正态分布概率密度曲线(正态性检验p值大于0.5)。货船在夜间(18:00～6:00)的船速正态分布概率密度曲线确定方法与日间(6:00～18:00)相同。油船和客船的船速特征分析方法与货船相同。不同船舶类型的船舶在不同时间段的船速正态分布概率密度曲线可以为航道、锚地通过能力仿真动态模拟提供输入参数，从而为航道规划、船舶组织调度研究提供定量参考。步骤4，确定不同船舶类型的船舶在航道中航行位置偏好如步骤3，将船舶类型分为货船、油船和客船。利用k-means方法对货船ais数据中的位置数据进行聚类分析，以平行于航道轴线的平行线作为聚类中心，聚类中心数量从1开始不断增加，计算不同聚类中心数量条件下任一聚类的畸变函数，根据“肘部法则”确定聚类中心数量，进而得到各聚类中心在航道内位置及各聚类货船数量在货船总量中所占的比例。油船和客船位置聚类分析方法与货船相同。油船、客船均对航行安全要求较高，因此明确不同类型船舶在航道中的位置偏好有助于海事部门指挥船舶航行，降低船舶碰撞概率。步骤5，确定各类型船舶中不同船长船舶的船速偏好如步骤3，将船舶类型分为货船、油船和客船。根据货船ais数据中的船长数据，将船长均匀分为多个类别(具体类别数量根据航道中船舶长度及分析需求综合确定)；再根据货船ais数据的船速数据，将船速均匀分为多个等级(具体等级数量根据航道中船舶速度及分析需求综合确定)，统计特定船长类别下处于不同船速等级的货船数及特定船长类别的货船总数，并计算特定船长类别下处于不同船速等级的货船数占特定船长类别下货船总数的比例。油船和客船的不同船长船舶的船速分析方法与货船相同。不同船舶类型、不同船长的船舶制动距离有所区别，因此明确各类型船舶中不同船长船舶的船速偏好有助于海事部门根据船舶类型和船长重点关注不满足相应船速偏好的船舶，避免船速过快或过慢，在保障航道通航效率的同时降低船舶碰撞风险。本发明的有益效果：本发明所述的基于ais数据的船舶航行特征及偏好分析方法，将船舶ais数据包含的静态信息和动态信息相结合，针对船舶在航道中的航行作业过程，探究不同类型船舶在不同时段的船速特征以及不同类型、不同船长船舶的航行偏好，确定了航道、锚地通过能力仿真动态模拟研究所需参数，加强了船舶航行中海事管理的针对性并提高效率。附图说明图1为船舶定线制航行示意图，船舶双向通航。图2为船舶航行速度分布图；其中，(a)为日间货船航行速度分布图；(b)为夜间货船航行速度分布图；(c)为日间油船航行速度分布图；(d)为夜间油船航行速度分布图。图3为对船舶位置分析结果图；其中，(a)为货船位置分析结果图；(b)为油船位置分析结果图。图4为船长和船速关系图；其中，(a)为货船船长与船速关系图；(b)为油船船长与船速关系图。具体实施方式下面结合附图和技术方案，进一步说明本发明的具体实施方式。所需要的数据从ais船舶自动识别系统中提取，数据来源于2018年某港区全部进出港船舶。一种基于ais数据的船舶航行特征及偏好分析方法，步骤如下：如图1所示，定线制航路满足船舶双向通行，船舶分别由东侧入口进入向西航行和由西侧入口进入向东航行，此说明以2018年自定线制航路东侧入口进入的船舶为分析对象。步骤1，定线制航路分隔带禁止船舶航行，在定线制航路东侧入口处划定长度200m、宽度为在航道宽度基础上向北延伸0.5倍航道宽度作为船舶ais数据统计区，保留船舶单次航程在数据统计区内任意一次ais数据即可，去掉多余数据，统计得到2018年通过数据统计区的所有船舶ais数据。步骤2，保留船舶ais数据中的船舶类型、船长、时间、船舶位置和船速，去掉船名、货物种类和船舶航向等其他信息。步骤3，此港区内无客运码头，航道中无客船行驶，因此将船舶类型分为货船和油船两类。以货船为分析对象，提取货船ais数据中时间处于日间(6:00～18:00)时的船速数据，确定这些船速数据的最大值和最小值，将最大值和最小值间的船速区间等分为20个速度区间，统计处于各速度区间的货船数并拟合得到货船在日间(6:00～18:00)行驶时船速的正态分布概率密度曲线，如图2中(a)所示。采用同样的方法得到货船在夜间(18:00～6:00)行驶时船速的正态分布概率密度曲线，如图2中(b)所示。以油船为分析对象，采用与货船相同的数据分析方法，分别得到油船在日间(6:00～18:00)和夜间(18:00～6:00)行驶的船速正态分布概率密度曲线，如图2中(c)和(d)所示。以上得到的正态分布概率密度曲线方程可以作为航道、锚地通过能力仿真动态模拟所需输入参数，仿真结果为航道规划、船舶组织调度提供定量参考。步骤4，将船舶类型分为货船和油船。以货船为分析对象，利用k-means方法对货船ais数据中的位置数据进行聚类分析，以平行于航道轴线的平行线作为聚类中心，聚类中心数量从1开始不断增加，计算不同聚类中心数量条件下任一聚类的畸变函数，根据“肘部法则”确定聚类中心数量为5，得到5个聚类中心在航道内位置及各聚类货船数量在货船总数中所占比例。其中，类别1船舶属于航道内接近分隔带的船舶，类别2船舶和类别3船舶属于航道中心位置的船舶，类别4船舶属于航道边缘的船舶，类别5船舶属于航道外的船舶。如图3中(a)所示，类别1船舶和类别2船舶的数量最多，分布也最密集；越靠近航道边缘，船舶分布越稀疏，并且在航道外有一定数量的船舶分布。以油船为分析对象，采用与货船相同的船舶位置数据分析方法，确定聚类中心数量为5，得到5个聚类中心在航道内位置及各聚类油船数量在油船总数中所占比例。各类别船舶在航道中所处位置与货船相同，如图3中(b)所示，油船位置分布趋势与货船相同。海事部门根据以上分析结果，指挥油船沿航道边缘行驶，降低船舶碰撞概率。步骤5，将船舶类型分为货船和油船。以货船为分析对象，根据货船ais数据中的船长数据，确定这些船长数据的最大值和最小值，将最值船长区间等分为5个类别；根据货船ais数据中的船速数据，确定这些船速数据的最大值和最小值，将最值船速区间等分为5个等级，具体数值见表1。统计特定船长类别下处于各船速等级的货船数及特定船长类别下的货船总数，并计算特定船长类别下处于各船速等级的货船数占特定船长类别下货船总数的比例，如图4中(a)所示，对于船长为类别1～类别3的货船，当船速位于船速等级1或船速等级5时，海事部门应给予重点关注。以油船为分析对象，船长等级对应的船长范围以及船速等级对应的船速范围确定方法同货船，具体数值见表2。统计特定船长类别下处于各船速等级的油船数及特定船长类别下的油船总数，并计算特定船长类别下处于各船速等级的油船数占特定船长类别下油船总数的比例，如图4中(b)所示，对于船长为类别3、类别4的油船，当船速位于船速等级1或船速等级4或船速等级5时，海事部门应给予重点关注。表1船长类别对应的船长范围以及速度等级对应的速度范围(货船)船长类别船长范围(m)速度等级速度范围(m/s)类别10<l≤113.2等级12.88<v≤6.46类别2113.2<l≤226.4等级26.46<v≤10.02类别3226.4<l≤339.6等级310.02<v≤13.58类别4339.6<l≤452.8等级413.58<v≤17.14类别5452.8<l≤566.0等级517.14<v≤20.7表2船长类别对应的船长范围以及速度等级对应的速度范围(油船)本发明提出的上述描述是为了便于该研究领域的普通研究人员能理解和运用本方法进行计算分析。熟悉本领域的人员可以对这些分析进行修改，并把在此说明的数据分析方法应用到其他例如不同类型船舶对航线的选择偏好中。本领域研究人员根据本发明的揭示，不脱离本发明范畴所做出的改进和修改都应该在本发明的保护范围之内。当前第1页12