一种船舶冲突解脱的规划方法技术

本发明专利技术涉及一种船舶冲突解脱的规划方法，包括如下几个步骤：首先通过海上交通控制中心获得其在每一采样时刻推测的各船舶在未来时段内的船舶轨迹；然后在每一采样时刻，基于船舶当前的运行状态和历史位置观察序列，获取海域风场变量的数值；然后在每一采样时刻，基于各船舶的运行状态和设定的船舶在海域内运行时需满足的安全规则集，当船舶间有可能出现违反安全规则的状况时，对其动态行为实施监控并为海上交通控制中心提供及时的告警信息；最后当告警信息出现时，在满足船舶物理性能和海域交通规则的前提下，通过设定优化指标函数以及融入风场变量数值，采用自适应控制理论方法对船舶避撞轨迹进行滚动规划，并将规划结果传输给各船舶执行。

A planning method for ship collision avoidance

The invention relates to a ship collision free planning method includes the following steps: firstly, through the maritime traffic control center obtained at each sampling time that each ship in the future period of ship trajectory; and then at each sampling time, the running condition of the ship and the historical position of the observation sequence based on variable gain the numerical sea wind field; and then at each sampling time, safety rules for the operation state of each ship and set the ship running in the area need to be based on the set, when the ship is likely to occur in violation of safety rules situation, to monitor the implementation of its dynamic behavior and provide timely warning information for maritime traffic control center when the alarm information; finally appeared, in order to meet the demand of ship physical properties and sea traffic rules, by setting the optimization index function and thawing The adaptive control theory is used to plan the Collision Avoidance Trajectory of the ship, and the planning results are transmitted to each ship.

【技术实现步骤摘要】
本申请是申请号为：201410844662.7，专利技术创造名称为《船舶冲突解脱的规划方法》，申请日为：2014年12月30日的专利技术专利申请的分案申请。
本专利技术涉及一种海域交通管制方法，尤其涉及一种基于滚动规划策略的船舶冲突解脱的规划方法。
技术介绍
随着全球航运业的快速发展，部分繁忙海域内的交通愈加拥挤。在船舶交通流密集复杂海域，针对船舶间的冲突情形仍然采用航行计划结合人工间隔调配的管制方式已不能适应航运业的快速发展。为保证船舶间的安全间隔，实施有效的冲突调配就成为海域交通管制工作的重点。船舶冲突解脱是航海领域中的一项关键技术，安全高效的解脱方案对于增加海域船舶流量以及确保海运安全具有重大意义。为了提高船舶的航行效率，船用雷达自动标绘仪目前已经被广泛应用到船舶监控和避碰中，该设备通过提取船舶相关信息为船舶间冲突情形的判定提供参考依据。尽管此类设备极大降低了人工监控的负荷，但它并不具备船舶自动冲突解脱功能。针对船舶冲突解脱问题，目前的处理方式主要包括几何式确定性算法和启发式智能算法两大类方案，相关文献研究主要集中在无约束条件下两船舶间的冲突避让规划算法并且多以“离线形式”为存在冲突的船舶规划解脱轨迹，由此造成各个船舶解脱轨迹的动态适应性和鲁棒性较差。
技术实现思路
本专利技术要解决的技术问题是提供一种鲁棒性较好的船舶冲突解脱的规划方法，该方法可有效防止船舶运行冲突。实现本专利技术目的的技术方案是提供一种船舶冲突解脱的规划方法，包括如下几个步骤：①通过海上交通控制中心获得其在每一采样时刻推测的各船舶在未来时段内的船舶轨迹；②在每一采样时刻，基于船舶当前的运行状态和历史位置观察序列，获取海域风场变量的数值；③在每一采样时刻，基于各船舶的运行状态和设定的船舶在海域内运行时需满足的安全规则集，当船舶间有可能出现违反安全规则的状况时，对其动态行为实施监控并为海上交通控制中心提供及时的告警信息；④当告警信息出现时，在满足船舶物理性能和海域交通规则的前提下，通过设定优化指标函数以及融入风场变量数值，采用模型预测控制理论方法对船舶避撞轨迹进行滚动规划，并将规划结果传输给各船舶执行，其具体过程如下：4.1)设定船舶避撞轨迹规划的终止参考点位置P、避撞策略控制时域Θ、轨迹预测时域W；4.2)设定在给定优化指标函数的前提下，基于合作式避撞轨迹规划思想，通过给各个船舶赋予不同的权重以及融入实时风场变量滤波数值，得到各个船舶的避撞轨迹和避撞控制策略并将规划结果传输给各船舶执行，且各船舶在滚动规划间隔内仅实施其第一个优化控制策略；4.3)在下一采样时刻，重复步骤4.2直至各船舶均到达其解脱终点。进一步的，所述步骤②获取海域风场变量的数值的具体过程如下：2.1)设定船舶的停靠位置为轨迹参考坐标原点并在水平面上建立横坐标轴和纵坐标轴；2.2)在船舶处于直线运行状态和匀速转弯运行状态时，构建海域风场线性滤波模型x1(t+Δt)＝F(t)x1(t)+w(t)和z(t)＝H(t)x1(t)+v(t)获取风场变量数值，其中Δt表示采样间隔，x1(t)表示t时刻的状态向量，z(t)表示t时刻的观测向量，且x1(t)＝[x(t),y(t),vx(t),vy(t),wx(t),wy(t)]T，其中x(t)和y(t)分别表示t时刻船舶位置在横坐标轴和纵坐标轴上的分量，vx(t)和vy(t)分别表示t时刻船舶速度在横坐标轴和纵坐标轴上的分量，wx(t)和wy(t)分别表示t时刻风场数值在横坐标轴和纵坐标轴上的分量，F(t)和H(t)分别表示状态转移矩阵和输出测量矩阵，w(t)和v(t)分别表示系统噪声向量和测量噪声向量：在船舶处于变速转弯运行状态时，构建海域风场非线性滤波模型x1(t+Δt)＝Ψ(t,x1(t),u(t))+w(t)、z(t)＝Ω(t,x1(t))+v(t)和u(t)＝[ωa(t),γa(t)]T，其中Ψ(·)和Ω(·)分别表示状态转移矩阵和输出测量矩阵，ωa(t)和γa(t)分别表示转弯率和加速率：其中：Δt表示采样时间间隔，2.3)根据所构建的滤波模型获取风场变量的数值。进一步的，所述步骤③中对各船舶的动态行为实施监控并为海上交通控制中心提供及时的告警信息的具体过程如下：3.1)构造船舶在海域内运行时需满足的安全规则集Dmr(t)≥Dmin，其中Dmr(t)表示任意两个船舶m和船舶r在t时刻的距离，Dmin表示船舶间的最小安全距离；3.2)依据采样时间，建立由船舶连续运行状态至离散采样状态的观测器Λ:Γ→Ξ，其中Γ表示船舶的连续运行状态，Ξ表示船舶的离散采样状态；3.3)当船舶m和r的观测器Λm和Λr的离散观测数值Ξm和Ξr在t时刻表明该向量不在安全规则集中时，即关系式Dmr(t)≥Dmin不成立时，立刻向海上交通控制中心发出告警信息。进一步的，步骤④中，步骤4.2)的具体过程是：令其中表示t时刻船舶R当前所在位置和下一航道点间的距离的平方，PR(t)＝(xRt,yRt)，那么t时刻船舶R的优先级指数可设定为：其中zt表示t时刻海域内存在冲突的船舶数目，由优先级指数的含义可知，船舶距离其下一航道点越近，其优先级越高；设定优化指标，其中R∈I(t)表示船舶代码且I(t)＝{1,2,...,Zt本文档来自技高网...

【技术保护点】
一种船舶冲突解脱的规划方法，其特征在于包括如下几个步骤：①通过海上交通控制中心获得其在每一采样时刻推测的各船舶在未来时段内的船舶轨迹；②在每一采样时刻，基于船舶当前的运行状态和历史位置观察序列，获取海域风场变量的数值；③在每一采样时刻，基于各船舶的运行状态和设定的船舶在海域内运行时需满足的安全规则集，当船舶间有可能出现违反安全规则的状况时，对其动态行为实施监控并为海上交通控制中心提供及时的告警信息；④当告警信息出现时，在满足船舶物理性能和海域交通规则的前提下，通过设定优化指标函数以及融入风场变量数值，采用模型预测控制理论方法对船舶避撞轨迹进行滚动规划，并将规划结果传输给各船舶执行，其具体过程如下：4.1)设定船舶避撞轨迹规划的终止参考点位置P、避撞策略控制时域Θ、轨迹预测时域W；4.2)设定在给定优化指标函数的前提下，基于合作式避撞轨迹规划思想，通过给各个船舶赋予不同的权重以及融入实时风场变量滤波数值，得到各个船舶的避撞轨迹和避撞控制策略并将规划结果传输给各船舶执行，且各船舶在滚动规划间隔内仅实施其第一个优化控制策略；4.3)在下一采样时刻，重复步骤4.2直至各船舶均到达其解脱终点；所述步骤②获取海域风场变量的数值的具体过程如下：2.1)设定船舶的停靠位置为轨迹参考坐标原点并在水平面上建立横坐标轴和纵坐标轴；2.2)在船舶处于直线运行状态和匀速转弯运行状态时，构建海域风场线性滤波模型x1(t+△t)＝F(t)x1(t)+w(t)和z(t)＝H(t)x1(t)+v(t)获取风场变量数值，其中△t表示采样间隔，x1(t)表示t时刻的状态向量，z(t)表示t时刻的观测向量，且x1(t)＝[x(t),y(t),vx(t),vy(t),wx(t),wy(t)]T，其中x(t)和y(t)分别表示t时刻船舶位置在横坐标轴和纵坐标轴上的分量，vx(t)和vy(t)分别表示t时刻船舶速度在横坐标轴和纵坐标轴上的分量，wx(t)和wy(t)分别表示t时刻风场数值在横坐标轴和纵坐标轴上的分量，F(t)和H(t)分别表示状态转移矩阵和输出测量矩阵，w(t)和v(t)分别表示系统噪声向量和测量噪声向量：F(t)=10sin(ωa(t)Δt)ωa(t)1-cos(ωa(t)Δt)ωa(t)Δt001cos(ωa(t)Δt)-1ωa(t)sin(ωa(t)Δt)ωa(t)0Δt00cos(ωa(t)Δt)sin(ωa(t)Δt)0000-sin(ωa(t)Δt)cos(ωa(t)Δt)00000010000001]]>H(k)=100000010000001000000100;]]>在船舶处于变速转弯运行状态时，构建海域风场非线性滤波模型x1(t+△t)＝Ψ(t,x1(t),u(t))+w(t)、z(t)＝Ω(t,x1(t))+v(t)和u(t)＝[ωa(t),γa(t)]T，其中Ψ(·)和Ω(·)分别表示状态转移矩阵和输出测量矩阵，ωa(t)和γa(t)分别表示转弯率和加速率：Ψ=x(t)+vx(t)(sin(ωa(t)Δt)ωa(t)+γa(t)Δtvx2(t)+vy2(t)C5)+vy(t)(1-cos(ωa(t)Δt)ωa(t)+γa(t)Δtvx2(t)+vy2(t)C6)+wx(t)y(t)-vx(t)(1-cos(ωa(t)Δt)ωa(t)+γa(t)Δtvx2(t)+vy2(t)C6)+vy(t)(sin(ωa(t)Δt)ωa(t)+γa(t)Δtvx2(t)+vy2(t)C5)+wy(t)((1+γa(t)Δtvx2(t)+vy2(t))(vx(t)cos(ωa(t)Δt)+vy(t)sin(ωa(t)Δt)))((1+γa(t)Δtvx2(t)+vy2(t))(vy(t)cos(ωa(t)Δt)-vx(t)sin(ωa(t)Δt)))wx(t)wy(t),]]>其中：△t表示采样时间间隔，C5=(sin(ωa(t)Δt)ωa(t)-1-cos(ω...

【技术特征摘要】
1.一种船舶冲突解脱的规划方法，其特征在于包括如下几个步骤：①通过海上交通控制中心获得其在每一采样时刻推测的各船舶在未来时段内的船舶轨迹；②在每一采样时刻，基于船舶当前的运行状态和历史位置观察序列，获取海域风场变量的数值；③在每一采样时刻，基于各船舶的运行状态和设定的船舶在海域内运行时需满足的安全规则集，当船舶间有可能出现违反安全规则的状况时，对其动态行为实施监控并为海上交通控制中心提供及时的告警信息；④当告警信息出现时，在满足船舶物理性能和海域交通规则的前提下，通过设定优化指标函数以及融入风场变量数值，采用模型预测控制理论方法对船舶避撞轨迹进行滚动规划，并将规划结果传输给各船舶执行，其具体过程如下：4.1)设定船舶避撞轨迹规划的终止参考点位置P、避撞策略控制时域Θ、轨迹预测时域W；4.2)设定在给定优化指标函数的前提下，基于合作式避撞轨迹规划思想，通过给各个船舶赋予不同的权重以及融入实时风场变量滤波数值，得到各个船舶的避撞轨迹和避撞控制策略并将规划结果传输给各船舶执行，且各船舶在滚动规划间隔内仅实施其第一个优化控制策略；4.3)在下一采样时刻，重复步骤4.2直至各船舶均到达其解脱终点；所述步骤②获取海域风场变量的数值的具体过程如下：2.1)设定船舶的停靠位置为轨迹参考坐标原点并在水平面上建立横坐标轴和纵坐标轴；2.2)在船舶处于直线运行状态和匀速转弯运行状态时，构建海域风场线性滤波模型x1(t+△t)＝F(t)x1(t)+w(t)和z(t)＝H(t)x1(t)+v(t)获取风场变量数值，其中△t表示采样间隔，x1(t)表示t时刻的状态向量，z(t)表示t时刻的观测向量，且x1(t)＝[x(t),y(t),vx(t),vy(t),wx(t),wy(t)]T，其中x(t)和y(t)分别表示t时刻船舶位置在横坐标轴和纵坐标轴上的分量，vx(t)和vy(t)分别表示t时刻船舶速度在横坐标轴和纵坐标轴上的分量，wx(t)和wy(t)分别表示t时刻风场数值在横坐标轴和纵坐标轴上的分量，F(t)和H(t)分别表示状态转移矩阵和输出测量矩阵，w(t)和v(t)分别表示系统噪声向量和测量噪声向量：F(t)=10sin(ωa(t)Δt)ωa(t)1-cos(ωa(t)Δt)ωa(t)Δt001cos(ωa(t)Δt)-1ωa(t)sin(ωa(t)Δt)ωa(t)0Δ...

【专利技术属性】
技术研发人员：韩云祥，赵景波，李广军，
申请(专利权)人：江苏理工学院，
类型：发明
国别省市：江苏;32