CN116777096B - 一种面向综合监管的海事巡航路径动态开辟方法及系统 - Google Patents
一种面向综合监管的海事巡航路径动态开辟方法及系统 Download PDFInfo
- Publication number
- CN116777096B CN116777096B CN202311001139.3A CN202311001139A CN116777096B CN 116777096 B CN116777096 B CN 116777096B CN 202311001139 A CN202311001139 A CN 202311001139A CN 116777096 B CN116777096 B CN 116777096B
- Authority
- CN
- China
- Prior art keywords
- supervision
- path
- new
- route
- aircraft
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 11
- 238000002922 simulated annealing Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 8
- 230000001105 regulatory effect Effects 0.000 claims description 7
- 230000005641 tunneling Effects 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 4
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 230000009191 jumping Effects 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 8
- 238000013139 quantization Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 3
- 235000006679 Mentha X verticillata Nutrition 0.000 description 2
- 235000002899 Mentha suaveolens Nutrition 0.000 description 2
- 235000001636 Mentha x rotundifolia Nutrition 0.000 description 2
- 238000000137 annealing Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000000779 smoke Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
- G06Q10/047—Optimisation of routes or paths, e.g. travelling salesman problem
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Feedback Control In General (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明涉及海事空中巡航管理技术领域,具体公开了一种面向综合监管的海事巡航路径动态开辟方法及系统,其中,该方法包括:根据监管辖区的范围和实时动态特征、监管路径、巡航可视范围综合确定目标函数;根据单个监管设施的单次监管的最短距离、固定成本分别确定约束条件;基于目标函数和约束条件,根据模拟退火算法动态开辟飞机监管路径,确定最佳路径。本发明根据辖区水域的实时动态监管需求,能够实时计算辖区监管力量的需求满足能力量化数值。并在给定投资和安全航行距离等约束条件下,开辟最佳的巡航推荐路径。
Description
技术领域
本发明属于海事空中巡航管理技术领域,具体涉及一种面向综合监管的海事巡航路径动态开辟方法及系统。
背景技术
如何动态生成最佳路径,满足在有限投入条件下,同时兼顾监管水域覆盖更广阔、到达更深远等多个管理目标,是有效巡航的关键。
如公开号为CN111428916A的发明专利《一种海上救援船舶巡航路径规划方法》,针对应急场景,通过求解线性规划问题,提出了一种基于海洋环境数据可以确定救援船最优巡航路径的方法。公开号为CN113361765A的发明专利《一种船载无人机协同监测港口船舶大气污染路径规划的方法》,针对海上船舶对港口环境的污染,提供了一种船载无人机协同监测港口船舶大气污染路径规划的方法。曲小同等(2020年)以台州为实验水域,基于K-means算法选择了28个海上任务点,然后运用蚁群算法开展了水域日常和应急立体协同巡航监管路径规划研究。王力等(2019年)以烟台为实验水域,运用聚类获取重点巡航范围后,基于混合整数线性规划方法开展了无人机协同巡航模型研究。
已有算法均先基于需求或者聚类等算法确定巡航路径节点,将优化路径设计简化为传统的旅行商问题,即在可能连接节点的所有巡航路径中,选择代价最低且满足限制性条件的路径,亦即以最低的代价遍历所有接节点。这种方法忽略了巡航的过程性,点之间只有直线连接一种方式,是在有限解集对最优方案的选择。而面向综合监管的巡航路径规划中,为了最大化覆盖监管需求,往往需要根据全域需求覆盖情况,开辟新的路径,需要在优化路径设计时同步生成路径节点。目前的已有方法均不能实现节点的自动选择与新路径的开辟。
因此,如何在一定的约束条件下(如航程/航时约束、传感器性能约束、飞行安全距离约束等),根据水域实时监管需求,集合各监管力量(如飞机、船艇等)的巡航能力,将监管业务分配给不同的监管力量,是亟需解决问题。
发明内容
为实现本发明目的提供的一种面向综合监管的海事巡航路径动态开辟方法,包括:
步骤S101:根据监管辖区的范围和实时动态特征、监管路径、巡航可视范围综合确定目标函数;
步骤S102:根据巡航飞机或船舶的单次监管的最短距离、固定成本分别确定约束条件;
步骤S103:基于目标函数和约束条件,在模拟退火算法的基础上开辟最佳巡航飞机或船舶监管路径,确定最大化满足监管需求的路径;
包括以下步骤S1031-S1035:
步骤S1031:设定预期路径个数,并初始化监管路径并设定初始参数;
步骤S1032:计算目标函数的初始值;
步骤S1033:根据空间扰动产生新的监管路线;
步骤S1034:根据Metropolis准则计算对所述新的监管路线的接受概率,确定下一次选择的路线;
步骤S1035:判断当前路径是否是最佳路径,若是,则终止计算,并输出最佳路径;若不是,则重复步骤S1033-步骤S1035。
在其中一些具体实施例中,步骤S101还包括:对监管辖区的范围内的监管点进行动态网格离散处理,得到n个监管点,综合确定目标函数为:
;
;
式中,pathes=(p1,p2,…,pk)为航线集合,表示监管辖区范围内海巡飞机或船舶分别的监管线路,其中,p i 为第i条航线;D n 为监管辖区的范围内的监管点需求总量;d n 为在确定路径pathes的情况下n个监管点中未被满足的需求量;N为航线上航迹点的总数。
在其中一些具体实施例中,步骤S102中基于单个监管设施的单次监管的最短距离确定的约束条件l(p)为:
;
式中,A为单个飞机监管路径长度的公里数阈值;xj、yj分别为该点的经纬度坐标,xj-1、yj-1分别为当前点的前一个点的经纬度坐标,K为航线条数。
在其中一些具体实施例中,步骤S102中基于固定成本确定的约束条件cost为:
;
式中,B为固定成本;L 1为飞机航线个数;L 2为海巡船舶航线个数;pa为第a个飞机航线;pc为第c个船舶航线;C1为巡航飞机每公里代价,C2为海巡船舶每公里代价。
在其中一些具体实施例中,步骤S1033包括:
选择当前路径P cur 中的一条路径p i ;其次,在p i 的k-2个中间节点中,随机选择p i 中任意一个监管点s j,其空间坐标为(x j,y j),以随机的距离r和角度v扰动s j的坐标位置产生新的空间坐标snew,所述新的空间坐标计算如下式所示:
;
;
再使用snew替换s j,生成新的航线p newi ,替换P cur 中的路径p i ,生成新的监管路线方案P new。
为实现本发明目的,本申请还提供了一种面向综合监管的海事巡航路径动态开辟系统,包括:
目标函数定义模块:用于根据监管辖区的范围和实时动态特征、监管路径、巡航可视范围综合确定目标函数;
约束函数生成模块:用于根据单个监管设施的单次监管的最短距离、固定成本分别确定约束条件;
最佳路径生成模块:用于基于所述目标函数和约束条件,根据模拟退火算法开辟巡航飞机或船舶监管路径,生成最佳路径;
所述最佳路径生成模块包括:
初始化单元:用于设定预期路径个数,并初始化监管路径并设定初始参数;
初始值计算单元:用于计算目标函数的初始值;
路线生成单元:用于根据空间扰动产生新的监管路线;
路线选择单元:用于根据Metropolis准则计算对所述新的监管路线的接受概率,确定下一次选择的路线;
最佳路径确定单元:用于判断当前路径是否是最佳路径,若是,则终止计算,并输出最佳路径;若不是,则跳转至所述路径生成单元重新生成新的监管路线。
在其中一些具体实施例中,目标函数定义模块还用于:对监管辖区的范围内的监管点进行动态网格离散处理,得到n个监管点,综合确定目标函数为:
;
;
式中,pathes=(p1,p2,…,pk)为航线集合,表示监管辖区范围内海巡飞机或船舶分别的监管线路,其中,p i 为第i条航线;D n 为监管辖区的范围内的监管点需求总量;d n 为在确定路径pathes的情况下n个监管点中未被满足的需求量;N为航线上航迹点的总数。
在其中一些具体实施例中,约束函数生成模块用于基于单个监管设施的单次监管的最短距离确定的约束条件l(p)为:
;
式中,A为单个飞机监管路径长度的公里数阈值;xj、yj分别为该点的经纬度坐标,xj-1、yj-1分别为当前点的前一个点的经纬度坐标,K为航线条数。
在其中一些具体实施例中,约束函数生成模块用于基于固定成本确定的约束条件cost为:
;
式中,B为固定成本;L 1为飞机航线个数;L 2为海巡船舶航线个数;p a 为第a个飞机航线;p c 为第c个船舶航线;C1为巡航飞机每公里代价,C2为海巡船舶每公里代价。
在其中一些具体实施例中,路线生成单元用于按照以下方式产生新的监管路线:
选择当前路径P cur 中的一条路径p i ;其次,在p i 的k-2个中间节点中,随机选择p i 中任意一个监管点s j,其空间坐标为(x j,y j),以随机的距离r和角度v扰动s j的坐标位置产生新的空间坐标snew,所述新的空间坐标计算如下式所示:
;
;
再使用snew替换s j,生成新的航线p newi ,替换P cur 中的路径p i ,生成新的监管路线方案P new。
本发明的有益效果:
本发明的一种面向综合监管的海事巡航路径动态开辟方法及系统,根据辖区水域的实时动态监管需求,能够实时计算辖区监管力量的需求满足能力量化数值。并在给定投资和安全航行距离等约束条件下,开辟最佳的巡航推荐路径。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明一种面向综合监管的海事巡航路径动态开辟方法一些具体实施例的流程示意图;
图2是本发明一种面向综合监管的海事巡航路径动态开辟方法的生成监管路径的流程示意图;
图3是本发明一种面向综合监管的海事巡航路径动态开辟系统一些具体实施例的结构示意图;
图4是本发明一种面向综合监管的海事巡航路径动态开辟系统的生成监管路径的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。
所述实施例的示例在附图中示出,其中自始至终相同或类似的符号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
参照图1所示,一种面向综合监管的海事巡航路径动态开辟方法,包括:
步骤S101:根据监管辖区的范围和实时动态特征、监管路径、巡航可视范围综合确定目标函数。
在本发明一些具体实施例中,步骤S101还包括:对监管辖区的范围内的监管点进行动态网格离散处理,得到n个监管点,综合确定目标函数为:
;
;
式中,pathes=(p1,p2,…,pk)为航线集合,表示监管辖区范围内海巡飞机或船舶分别的监管线路,其中,p i 为第i条航线;D n 为监管辖区的范围内的监管点需求总量;d n 为在确定路径pathes的情况下n个监管点中未被满足的需求量;N为航线上航迹点的总数。
具体的,巡航力量集合,L=(1,2,…,l),表示辖区拥有的海巡船舶或飞机总数。
航线集合,pathes=(p1,p2,…,pk),表示辖区各海巡船舶或飞机总数分别的监管线路,是k条航线的集合。其中,p i 为第i条航线。
航迹点集合,,s1为其出发基地和返回的目的地。
航迹点sj=(xj,yj)。航迹点sj为巡航区域中所开辟的路径节点,xj,yj分别为该点的经纬度坐标。
N为区域内离散格网的数目,n为遍历中的某一个离散格网。
步骤S102:根据巡航飞机或船舶的单次监管的最短距离、固定成本分别确定约束条件。
在本发明一些具体实施例中,令起始点与终止点均为监管点,即s1=sk=s0。根据飞机性能和执法人员疲劳等因素,尽量减少远航程的航线条数和航次,因此单个飞机监管路径长度要求小于A公里,基于单个监管设施的单次监管的最短距离确定的约束条件l(p)为:
;
式中,A为单个飞机监管路径长度的公里数阈值;xj、yj分别为该点的经纬度坐标,xj-1、yj-1分别为当前点的前一个点的经纬度坐标,K为航线条数,根据历史统计数据设计总监管费用cost约束条件为:
;
其中,L1为飞机航线个数,L2为海巡船舶航线个数,pa为第a个飞机航线,pc为第c个船舶航线,C1为巡航飞机每公里代价,C2为海巡船舶巡航每公里代价。
步骤S103:基于目标函数和约束条件,根据模拟退火算法开辟巡航飞机或船舶监管路径,确定最佳路径。
在本发明一些具体实施例中,参照图2所示,
步骤S1031:设定预期路径个数,并初始化监管路径并设定初始参数。
具体的,设定预期路径个数,生成P 1条初始的飞机路径和P 2条初始的船舶路径,所有路径的方案为pathes;其次,以40公里为间隔,加密所有路径中的节点,设置当前路径P cur =pathes;最后优化初始参数,包括初始温度T 0=1,当前温度t=T 0退火速率c=0.99,最低温度Tmin=1×10-7最大迭代次数为itermax=10000,连续拒绝次数为conRejmax=100,当前迭代次数为iter=0,当前连续拒绝次数conRej=0。
步骤S1032:计算目标函数的初始值。
具体的,P cur 中的每一条路径,按监管手段d对于监管任务m的满足程度建立视域缓冲区(为简化处理,飞机设为800公里,船舶为600公里),生成单条航线轨迹的可视域。其次,将缓冲区与监管需求栅格进行空间叠置,将缓冲区内部监管需求置位0,获得更新后的需求值d n (pathes);
;
计算海巡船和飞机覆盖后监管需求的剩余值O cur ,作为初始目标函数值,如下式所示:
。
步骤S1033:根据空间扰动产生新的监管路线。
具体的,随机选择P cur 中的一条路径p i ;其次,在p i 的k-2个中间节点中,随机选择p i 中任意一个监管点s j,其空间坐标为(x j,y j),以随机的距离r和角度v扰动s j的坐标位置,产生snew,其坐标计算如下式所示:
;
;
然后,使用snew替换s j,生成新的航线p newi ,替换P cur 中的路径p i ,生成新的监管路线方案P new;进一步的,对P new中的每一条路径,按照视域做缓冲区,其中飞机为800公里,船舶为600公路,将缓冲区与监管需求栅格进行空间叠置,将缓冲区内部监管需求置位0,获得更新后的需求值O new,如下所示:
。
步骤S1034:根据Metropolis准则计算对新的监管路线的接受概率,确定下一次选择的路线。
具体的,根据下式依次计算每一个目标对s2的接受概率ac:
;
利用计算机程序随机生成一个0-1之间的随机数rand,如果rand<ac,则接受新的扰动路线,设置P cur =P new,O cur =O new,conRej=0;否则,拒绝新的路径P new,设置conRej=conRej+1;无论接受新路线与否,更新迭代次数iter=iter+1,温度t=t×c。
步骤S1035:判断当前路径是否是最佳路径,若是,则终止计算;若不是,则重复步骤S1033-步骤S1035。
具体的,如果conRe≥conRejmax,或者iter≥itermax,或者t<minT,则当前路径为最佳路径,终止计算;否则转步骤S1033继续迭代,直至输出当前路径为最佳路径。若当前路径为最佳路径,输出当前路径P cur 、当前目标函数取值O cur 以及迭代终止参数conRej、iter和t。
为实现本发明目的,本申请还提供了一种面向综合监管的海事巡航路径动态开辟系统,参照图3所示,包括:
目标函数定义模块10:用于根据监管辖区的范围和实时动态特征、监管路径、巡航可视范围综合确定目标函数。
在本发明一些具体实施例中,对监管辖区的范围内的监管点进行动态网格离散处理,得到n个监管点,综合确定目标函数为:
;
;
式中,pathes=(p1,p2,…,pk)为航线集合,表示监管辖区范围内海巡飞机或船舶分别的监管线路,其中,p i 为第i条航线;D n 为监管辖区的范围内的监管点需求总量;d n 为在确定路径pathes的情况下n个监管点中未被满足的需求量;N为航线上航迹点的总数。
具体的,巡航力量集合,L=(1,2,…,l),表示辖区拥有的海空力量个数。
航线集合,pathes=(p1,p2,…,pk),表示辖区各个海空力量(船舶和飞机)分别的监管线路,是k条航线的集合。其中,p i 为第i条航线。
航迹点集合,,s1为其出发基地和返回的目的地。
航迹点sj=(xj,yj)。航迹点sj为巡航区域中所开辟的路径节点,xj,yj分别为该点的经纬度坐标。
N为区域内离散格网的数目,n为遍历中的某一个离散格网。
约束函数生成模块20:用于根据巡航飞机或船舶的单次监管的最短距离、固定成本分别确定约束条件。
在本发明一些具体实施例中,令起始点与终止点均为监管点,即s1=sk=s0。根据飞机性能和执法人员疲劳等因素,尽量减少远航程的航线条数和航次,因此单个飞机监管路径长度要求小于A公里,基于单个监管设施的单次监管的最短距离确定的约束条件l(p)为:
;
式中,A为单个飞机监管路径长度的公里数阈值;xj、yj分别为该点的经纬度坐标,xj-1、yj-1分别为当前点的前一个点的经纬度坐标,K为航线条数。
根据历史统计数据设计总监管费用cost约束条件:
;
其中,L1为飞机航线个数,L2为海巡船舶航线个数,pa为第a个飞机航线,pc为第c个船舶航线,C1为巡航飞机每公里代价,C2为海巡船舶每公里代价。
最佳路径生成模块30:用于基于目标函数和约束条件,根据模拟退火算法开辟巡航飞机或船舶监管路径,确定最佳路径。
在本发明一些具体实施例中,参照图4所示,
初始化单元31:用于设定预期路径个数,并初始化监管路径并设定初始参数。
具体的,设定预期路径个数,生成P 1条初始的飞机路径和P 2条初始的船舶路径,所有路径的方案为pathes;其次,以40公里为间隔,加密所有路径中的节点,设置当前路径P cur =pathes;最后优化初始参数,包括初始温度T 0=1,当前温度t=T 0退火速率c=0.99,最低温度Tmin=1×10-7最大迭代次数为itermax=10000,连续拒绝次数为conRejmax=100,当前迭代次数为iter=0,当前连续拒绝次数conRej=0。
初始值计算单元32:用于计算目标函数的初始值。
具体的,P cur 中的每一条路径,按监管手段d对于监管任务m的满足程度建立视域缓冲区(为简化处理,飞机设为800公里,船舶为600公里),生成单条航线轨迹的可视域。其次,将缓冲区与监管需求栅格进行空间叠置,将缓冲区内部监管需求置位0,获得更新后的需求值d n (pathes);
;
计算海巡船和飞机覆盖后监管需求的剩余值O cur ,作为初始目标函数值,如下式所示:
。
路线生成单元33:用于根据空间扰动产生新的监管路线。
具体的,首先,随机选择P cur 中的一条路径p i ;其次,在p i 的k-2个中间节点中,随机选择p i 中任意一个监管点s j,其空间坐标为(xj,yj),以随机的距离r和角度v扰动s j的坐标位置,产生snew,其坐标计算如下式所示:
;
;
然后,使用snew替换s j,生成新的航线p newi ,替换P cur 中的路径p i ,生成新的监管路线方案P new;进一步的,对P new中的每一条路径,按照视域做缓冲区,其中飞机为800公里,船舶为600公路,将缓冲区与监管需求栅格进行空间叠置,将缓冲区内部监管需求置位0,获得更新后的需求值O new,如下所示:
。
路线选择单元34:用于根据Metropolis准则计算对新的监管路线的接受概率,确定下一次选择的路线。
具体的,根据下式依次计算每一个目标对s2的接受概率ac:
;
利用计算机程序随机生成一个0-1之间的随机数rand,如果rand<ac,则接受新的扰动路线,设置P cur =P new,O cur =O new,conRej=0;否则,拒绝新的路径P new,设置conRej=conRej+1;无论接受新路线与否,更新迭代次数iter=iter+1,温度t=t×c。
最佳路径确定单元35:用于判断当前路径是否是最佳路径,若是,则终止计算,并输出最佳路径;若不是,则跳转至路径生成单元重新生成新的监管路线。
具体的,如果conRe≥conRejmax,或者iter≥itermax,或者t<minT,则当前路径为最佳路径,终止计算;否则转路线生成单元继续迭代,直至输出当前路径为最佳路径。若当前路径为最佳路径,输出当前路径P cur 、当前目标函数取值O cur 以及迭代终止参数conRej、iter和t。
本发明的一种面向综合监管的海事巡航路径动态开辟方法及系统根据辖区水域的实时动态监管需求或风险值,能够实时计算辖区监管力量的需求满足能力量化数值,并在给定投资和安全航行距离等约束条件下,给出最优的巡航推荐路径。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、“一个具体实施例”或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
最后应说明的是:以上实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不驱使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。
Claims (6)
1.一种面向综合监管的海事巡航路径动态开辟方法,其特征在于,包括:
步骤S101:根据监管辖区的范围和实时动态特征、监管路径、巡航可视范围确定目标函数;
步骤S102:根据巡航飞机或船舶的单次监管的最短距离、固定成本分别确定约束条件;
步骤S103:基于所述目标函数和所述约束条件,根据模拟退火算法开辟所述巡航飞机或船舶的监管路径,确定最佳路径;包括以下步骤S1031-S1035:
步骤S1031:设定预期路径个数,并初始化监管路径并设定初始参数;
步骤S1032:计算目标函数的初始值;
步骤S1033:根据空间扰动产生新的监管路线;
步骤S1034:根据Metropolis准则计算对所述新的监管路线的接受概率,确定下一次选择的路线;
步骤S1035:判断当前路径是否是最佳路径,若是,则终止计算,并输出最佳路径;若不是,则重复步骤S1033-步骤S1035;
所述步骤S101还包括:对所述监管辖区的范围内的监管点进行动态网格离散处理,得到n个监管点,综合确定目标函数为:
式中,pathes=(p1,p2,…,pk)为航线集合,表示监管辖区范围内海巡飞机或船舶分别的监管线路,其中,pi为第i条航线;Dn为监管辖区的范围内的监管点需求总量;dn为在确定路径pathes的情况下n个监管点中未被满足的需求量;N为航线上航迹点的总数;
所述步骤S1033还包括:选择当前路径Pcur中的一条路径pi;其次,在pi的k-2个中间节点中,随机选择pi中任意一个监管点sj,其空间坐标为(xj,yj),以随机的距离r和角度v扰动sj的坐标位置产生新的空间坐标snew,所述新的空间坐标计算如下式所示:
xnew=xj+r*cos(v)
ynew=yj+r*sin(v)
再使用snew替换sj,生成新的航线pnewi,替换Pcur中的路径pi,生成新的监管路线方案Pnew;
所述步骤S1034还包括:根据下式依次计算每一个目标对s2的接受概率ac:
利用计算机程序随机生成一个0-1之间的随机数rand,如果rand<ac,则接受新的扰动路线,设置Pcur=Pnew,Ocur=Onew,conRej=0;否则,拒绝新的路径Pnew,设置conRej=conRej+1;无论接受新路线与否,更新迭代次数iter=iter+1,温度t=t*c。
2.根据权利要求1所述的面向综合监管的海事巡航路径动态开辟方法,其特征在于,所述步骤S102中包括:基于单个监管设施的单次监管的最短距离确定如下约束条件l(p):
式中,A为单个飞机监管路径长度的公里数阈值;xj、yj分别为该点的经纬度坐标,xj-1、yj-1分别为当前点的前一个点的经纬度坐标,K为航线条数。
3.根据权利要求1所述的面向综合监管的海事巡航路径动态开辟方法,其特征在于,所述步骤S102中还包括:基于固定成本确定如下约束条件cost:
式中,B为固定成本;L1为飞机航线个数;L2为海巡船舶航线个数;pa为第a个飞机航线;pc为第c个船舶航线;C1为巡航飞机每公里代价,C2为海巡船舶每公里代价。
4.一种面向综合监管的海事巡航路径动态开辟系统,其特征在于,包括:
目标函数定义模块:用于根据监管辖区的范围和实时动态特征、监管路径、巡航可视范围确定目标函数;
约束函数生成模块:用于根据巡航飞机或船舶的单次监管的最短距离、固定成本分别确定约束条件;
最佳路径生成模块:用于基于所述目标函数和约束条件,根据模拟退火算法开辟巡航飞机或船舶监管路径,生成最佳路径;
所述最佳路径生成模块包括:
初始化单元:用于设定预期路径个数,并初始化监管路径并设定初始参数;
初始值计算单元:用于计算目标函数的初始值;
路线生成单元:用于根据空间扰动产生新的监管路线;
路线选择单元:用于根据Metropolis准则计算对所述新的监管路线的接受概率,确定下一次选择的路线;
最佳路径确定单元:用于判断当前路径是否是最佳路径,若是,则终止计算,并输出最佳路径;若不是,则跳转至所述路径生成单元重新生成新的监管路线;
所述目标函数定义模块还用于:对所述监管辖区的范围内的监管点进行动态网格离散处理,得到n个监管点,综合确定目标函数为:
式中,pathes=(p1,p2,…,pk)为航线集合,表示监管辖区范围内海巡飞机或船舶分别的监管线路,其中,pi为第i条航线;Dn为监管辖区的范围内的监管点需求总量;dn为在确定路径pathes的情况下n个监管点中未被满足的需求量;N为航线上航迹点的总数;
所述路线生成单元还用于:选择当前路径Pcur中的一条路径pi;其次,在pi的k-2个中间节点中,随机选择pi中任意一个监管点sj,其空间坐标为(xj,yj),以随机的距离r和角度v扰动sj的坐标位置产生新的空间坐标snew,所述新的空间坐标计算如下式所示:
xnew=xj+r*cos(v)
ynew=yj+r*sin(v)
再使用snew替换sj,生成新的航线pnewi,替换Pcur中的路径pi,生成新的监管路线方案Pnew;
所述路线选择单元还用于:根据下式依次计算每一个目标对s2的接受概率ac:
利用计算机程序随机生成一个0-1之间的随机数rand,如果rand<ac,则接受新的扰动路线,设置Pcur=Pnew,Ocur=Onew,conRej=0;否则,拒绝新的路径Pnew,设置conRej=conRej+1;无论接受新路线与否,更新迭代次数iter=iter+1,温度t=t*c。
5.根据权利要求4所述的面向综合监管的海事巡航路径动态开辟系统,其特征在于,所述约束函数生成模块用于基于海巡船舶或飞机的单次监管的最短距离确定如下约束条件l(p):
式中,A为单个飞机监管路径长度的公里数阈值;xj、yj分别为该点的经纬度坐标,xj-1、yj-1分别为当前点的前一个点的经纬度坐标,K为航线条数。
6.根据权利要求4所述的面向综合监管的海事巡航路径动态开辟系统,其特征在于,所述约束函数生成模块用于基于固定成本确定的如下约束条件cost:
式中,B为固定成本;L1为飞机航线个数;L2为海巡船舶航线个数;pa为第a个飞机航线;pc为第c个船舶航线;C1为巡航飞机每公里代价,C2为海巡船舶每公里代价。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311001139.3A CN116777096B (zh) | 2023-08-10 | 2023-08-10 | 一种面向综合监管的海事巡航路径动态开辟方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311001139.3A CN116777096B (zh) | 2023-08-10 | 2023-08-10 | 一种面向综合监管的海事巡航路径动态开辟方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116777096A CN116777096A (zh) | 2023-09-19 |
CN116777096B true CN116777096B (zh) | 2023-10-31 |
Family
ID=87994695
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311001139.3A Active CN116777096B (zh) | 2023-08-10 | 2023-08-10 | 一种面向综合监管的海事巡航路径动态开辟方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116777096B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106679666A (zh) * | 2016-12-26 | 2017-05-17 | 北京联合大学 | 无人机大范围巡航路径规划模型及方法 |
CN115565093A (zh) * | 2022-12-07 | 2023-01-03 | 交通运输部水运科学研究所 | 一种海事飞机空中巡航可达域的量化评估方法与系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2668609A4 (en) * | 2011-01-25 | 2014-10-08 | Nextgen Aerosciences Llc | METHOD AND APPARATUS FOR DYNAMIC MANAGEMENT OF AIRCRAFT TRAJECTORIES |
-
2023
- 2023-08-10 CN CN202311001139.3A patent/CN116777096B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106679666A (zh) * | 2016-12-26 | 2017-05-17 | 北京联合大学 | 无人机大范围巡航路径规划模型及方法 |
CN115565093A (zh) * | 2022-12-07 | 2023-01-03 | 交通运输部水运科学研究所 | 一种海事飞机空中巡航可达域的量化评估方法与系统 |
Non-Patent Citations (3)
Title |
---|
基于A-Star和改进模拟退火算法的航迹规划;岳秀;张超峰;张伟;王亚刚;;控制工程;27(08);第1365-1371页 * |
基于模拟退火算法的航迹规划方法研究;高晓静, 李俊山, 赵宗涛, 李新社;微电子学与计算机(05);第10-14页 * |
无人机海事巡航路径优化研究;魏蕾;李云斌;;天津航海(04);第61-64页 * |
Also Published As
Publication number | Publication date |
---|---|
CN116777096A (zh) | 2023-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Huang et al. | Non-dominated immune-endocrine short feedback algorithm for multi-robot maritime patrolling | |
CN111428916B (zh) | 一种海上救援船舶巡航路径规划方法 | |
Zhao et al. | Multi-objective weather routing algorithm for ships based on hybrid particle swarm optimization | |
CN111538349B (zh) | 面向多任务的长航程auv自主决策方法 | |
Ma et al. | Optimized BP neural network algorithm for predicting ship trajectory | |
CN112819255A (zh) | 基于粒子群-遗传算法的多准则船舶航线确定方法、装置、计算机设备及可存储介质 | |
CN116777096B (zh) | 一种面向综合监管的海事巡航路径动态开辟方法及系统 | |
Jin et al. | An overview: Maritime spatial-temporal trajectory mining | |
CN109726910B (zh) | 一种海运航线选择方法 | |
CN117636691B (zh) | 船舶航行的安全治理方法、装置、电子设备及存储介质 | |
Pan et al. | Safe and efficient UAV navigation near an airport | |
CN102298855A (zh) | 基于仿真的危险化学品船舶交通组织优化方法 | |
CN115981377A (zh) | 一种无人机动态避障方法及系统 | |
Yoo et al. | Ship route optimization considering on-time arrival probability under environmental uncertainty | |
CN115713872A (zh) | 一种基于soa的环境自适应航路规划方法 | |
Cai et al. | Optimization of the crossing waypoints in air route network | |
CN114882736A (zh) | 可容忍风险船舶航路与海上区域建筑安全间距计算方法 | |
Zhang et al. | UAVs vs. Pirates: an anticipatory swarm monitoring method using an adaptive pheromone map | |
Xue | Assessment of Ship‐Overtaking Situation Based on Swarm Intelligence Improved KDE | |
Liu et al. | Optimization design of ship collision avoidance strategy based on genetic algorithm | |
Zhou et al. | Theory and application of vessel speed dynamic control considering safety and environmental factors | |
Bye et al. | A RECEDING HORIZON GENETIC ALGORITHM FOR DYNAMIC MULTI-TARGET ASSIGNMENT AND TRACKING-A Case Study on the Optimal Positioning of Tug Vessels along the Northern Norwegian Coast | |
Moon et al. | Develop and evaluate of intelligent autonomous-ship framework | |
Liu et al. | Aircraft rerouting decision-making model under severe weather | |
Yang et al. | An Offshore Maritime Edge Computing Network: Design and Experiment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |