传统的海上交通调查手段复杂且耗时费力，船舶自动识别系统（Automatic Identification System，AIS）不仅能提高船舶航行效率与安全，而且能提供大量的海上交通特征.AIS的信息采集系统能实时地接收各港口AIS接收点采集到的船载AIS数据.这些数据具有实时性强、数量大、内容丰富等特点，可以供基于电子海图的船舶监控系统使用，也可以作为历史数据对海上交通特征进行分析.
船舶在港口航道（以下简称“航道”）中航行，需要选择合适的航速.若航速过低，则舵效较差，很难有效地控制船首向和船位.若航速过高，则船体的下沉量较大，船体水下的富余水深随之减少，致使船舶所受的阻力增大，舵效变差.由于船舶惯性较大，航速较高时，减速所需距离较大，不利于在避让或紧急情况下采取有效措施.因此，船舶进出港时的航速受到诸多约束，同时也呈现出一定的规律.
大型船舶具有惯性大、余速消除慢等特点，控制船舶进出港航速是保证船舶安全的一种重要方法.本文通过构建AIS信息数据库，经统计分析得到船舶进出港各阶段的航速分布，从而提出航速控制建议，为船舶自动进出港提供理论支撑，为新引航员及船长提供航速参考，有效指导船舶交通管理（Vessel Traffic Services，VTS）部门对速度异常的船舶进行监视，提高VTS的主动服务效率.
1船舶进出港各阶段航速统计实现
1.1船舶进出港航速统计总体流程
船舶进出港各阶段航速统计流程见图1.图1中：ETL（Extract Transform Load）为数据的提取、转换和加载技术；OLAP（OnLine Analytical Processing）为在线联机分析处理技术.
1.2主要计算公式
利用所采集的有效AIS信息，只 需输入船舶在码头的起算点经纬度，即可利用式（1）求得该时刻船舶与码头间距离为船长的几倍，此距离所对应的航速即为AIS信息中该时刻船舶的航速.
DL=1 852.25D（（φstart，λstart），（φend，λend））L （1）
式中：DL为舶舶与码头间距离所等效的船长倍数；D为海图上两经纬度之间的距离，m；L为船长，m；（φstart，λstart）表示码头起算点的纬度和经度；（φend，λend）表示该时刻船舶的纬度和经度.
1.3船舶进出港航速统计具体步骤
步骤1AIS信息采集与预处理.
通过定时地从AIS数据接收服务器中取出所需AIS消息报文文件，运用AIS信息解码模块实现解码和数据入库，并在数据库中预处理使用前的无效数据，使得到的AIS信息基本表更加有效.预处理的主要工作有：删除动态信息表中MMSI为0的记录；删除明显错误的AIS数据，如速度、经纬度超过合理值（如出现负值速度或速度>100 kn，或经度>180°，或纬度>90°等）；根据AIS_NAME表删除超出研究水域范围的记录等.基于导入的原始数据，建立数据库，围绕主题进行数据建模和统计分析.
步骤2定量划分各阶段.
船舶操纵运动不但与船舶运动控制有关，还与船型参数等有关.不同排水量、不同船型参数的船舶的运动性能也不尽相同.根据参考文献[11]中“利用距离码头几倍船长控制船舶速度区间范围，进而完成对船舶系泊操纵”的论述，本文同样采用几倍船长定量划分船舶航行的各个阶段，进而统计得出各阶段船舶航速分布.
步骤3船舶最远距离的确定.
①船舶进港过程为减速操纵，船舵控制航向的能力随着船速的降低而减弱；反之，出港过程为加速操纵，船舵控制航向的能力随着船速的提高而增强.相比较而言，船舶出港操纵相对较容易，故在此以船舶进港操纵过程确定研究距离.
②在船舶由沿海水域驶入港口水域并向停泊位置接近的过程中，由于港内航行需要频繁改变船速，船舶应处于备车航行状态，是否降速取决于距离停泊位置的距离、船舶吨位、操纵性、通航环境以及船型等因素.如大型集装箱船，由于其操纵性能较好，一般在10 n mile左右时进行备车，此时距离泊位大约为50L. 步骤4统计各阶段航速.
①确定所需研究的船舶吨级和港口，确定划分阶段最远距离；
②进行AIS数据预处理，根据①在数据库中搜索出在一定时间空间内满足条件的船舶AIS动态表数据，存入临时表A中；
③对临时表A按MMSI分组查询，每组按时间先后顺序排序，计算前后船舶AIS记录的时间差、经度差、纬度差，选择这些差值小于其设置的最大阈值、航速小于其设置的最大阈值的数据，存入临时表B中；
④根据数学模型，计算表B中实时船舶位置距离码头起算点距离，换算成船长倍数（根据历史船舶交通流，选取船舶主要靠泊位置的中心点为船舶距离泊位的起算点），并将进港船舶数据存入表C中，出港船舶数据存入表D中；
⑤根据研究需要，对表C和D进行相应的深入分析.
步骤5获取拟合曲线.通过上述步骤可以得到该类船舶进出港航速概率点分布（各阶段的航速都处在一定的范围区间内）.为得到最佳的拟合曲线，采用最小二乘法对数据进行回归拟合.因为接收到的船舶AIS信息已经考虑风流综合影响，所以在实际的船舶操纵过程中应保持船舶对地航速服从最佳拟合曲线分布.
2实例研究
2.1数据选取及统计分析结果
厦门港是我国综合运输体系的重要枢纽、集装箱运输干线港、东南国际航运中心和对台航运主要口岸，位居世界集装箱港口前20强，主要船型为1万TEU左右（12万吨级）的集装箱船舶.嵩屿港集装箱码头在厦门所有集装箱码头中最靠近主航道，为便于研究船舶进出港各阶段航速，选择嵩屿港集装箱码头船型1万TEU左右的集装箱船舶为研究对象.选取2013年1月，4月，6—10月的AIS数据.
根据条件滤除错误数据，从数据库中得出符合上述条件的有效船舶的AIS数据，共计船舶44艘次，其中：进港船舶24艘次；出港船舶20艘次.船舶进港和出港速度分布及其最佳拟合曲线见图2.
航迹分布和距离嵩屿港码头约30L处的门限，其中：A为厦门港主航道，B为距离嵩屿港码头约30L处的门限.
由图2可知，当DL<15时，船舶进出港航速的变化服从对数分布，但当DL≈30时，由于厦门港主航道两侧青屿和浅滩的影响，船舶需采取安全航速通过此门限，船舶最终航速服从多项式分布.
为进一
步分析航速的分布特征，通过获取的最佳拟合曲线，分别得到12万吨级集装箱船舶进港和出港时的速度分布，即
Y1=-0.7D5L-0.5D4L-0.000 5D3L-
0.056D2L+1.600 9DL （2）
Y2=-0.7D5L-0.4D4L+0.007D3L-
0.221 7D2L+2.857 9DL （3）
式中：Y1和Y2分别为进港和出港船舶速度分布.
综合本航道的特点，提出对该吨级船舶进出港速度控制的建议，见表1和2.
从表1可知：对于进港船舶，当DL∈[30，40）时，由于门限两侧青屿和浅滩的影响，船舶减速并将速度控制在7～8 kn，在通过门限后，为保证进港时的效率，可适度加车进港；当DL∈[10，20）时，速度达到最大值，控制在12～13 kn，之后开始减速.从表2可知，对于出港船舶，
2.3结果验证
为验证上述数据的准确性，请具有丰富海上实践经验、持有船长适任证书的人员在大型船舶操纵模拟器上按上述建议航速进行进出港模拟.模拟平台为大连海事大学航海技术研究所研制的大型船舶操纵模拟器.该模拟器已取得中国海事局和IMO认证，满足本次模拟的各项要求.模拟结果表明，按照该建议航速操纵船舶进出港，能够获得较好的效果.
3结束语
基于AIS数据，通过运用数据库的统计和分析功能，进而求得船舶进出港的速度分布.对厦门嵩屿港集装箱码头12万吨级集装箱船舶进出港各阶段的航速分布进行统计分析，并根据速度的概率分布，给出该类船舶进出厦门嵩屿港集装箱码头的建议航速.该研究可以为引航员有效引导船舶进出港提供参考，为船舶自动进出港航速提供新思路，并为港口主管部门更好地监管通航环境提供理论依据.此外，本文只对船舶航速规律进行统计，如何将此与进出港各阶段航向分布相结合，以提高船舶进出港的效率和安全，将是下一步的主要研究方向.
参考文献：
[1]吴兆麟， 朱军. 海上交通工程[M]. 2版. 大连： 大连海事大学出版社， 2004： 1011.
[2]孙方力， 孙文强. 船载自动识别系统[M]. 大连： 大连海事大学出版社， 2004： 110.
[3]国际海事组织. 通用船载自动识别系统国际标准汇编[G]. 袁安存， 张淑芳， 编译. 大连： 大连海事大学出版社， 2005： 110.
