本发明涉及海事空中巡航管理技术领域，具体公开了一种面向综合监管的海事巡航路径动态开辟方法及系统，其中，该方法包括：根据监管辖区的范围和实时动态特征、监管路径、巡航可视范围综合确定目标函数；根据单个监管设施的单次监管的最短距离、固定成本分别确定约束条件；基于目标函数和约束条件，根据模拟退火算法动态开辟飞机监管路径，确定最佳路径。本发明根据辖区水域的实时动态监管需求，能够实时计算辖区监管力量的需求满足能力量化数值。并在给定投资和安全航行距离等约束条件下，开辟最佳的巡航推荐路径。
1.一种面向综合监管的海事巡航路径动态开辟方法，其特征在于，包括：步骤S101：根据监管辖区的范围和实时动态特征、监管路径、巡航可视范围确定目标函数；
步骤S102：根据巡航飞机或船舶的单次监管的最短距离、固定成本分别确定约束条件；
步骤S103：基于所述目标函数和所述约束条件，根据模拟退火算法开辟所述巡航飞机或船舶的监管路径，确定最佳路径；包括以下步骤S1031‑S1035：步骤S1031：设定预期路径个数，并初始化监管路径并设定初始参数；
步骤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)
再使用snew替换sj，生成新的航线pnewi，替换Pcur中的路径pi，生成新的监管路线方案Pnew；
所述步骤S1034还包括：根据下式依次计算每一个目标对s2的接受概率ac：利用计算机程序随机生成一个0‑1之间的随机数rand，如果rand
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)
再使用snew替换sj，生成新的航线pnewi，替换Pcur中的路径pi，生成新的监管路线方案Pnew；
所述路线选择单元还用于：根据下式依次计算每一个目标对s2的接受概率ac：利用计算机程序随机生成一个0‑1之间的随机数rand，如果rand
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为海巡船舶每公里代价。
一种面向综合监管的海事巡航路径动态开辟方法及系统
技术领域
[0001] 本发明属于海事空中巡航管理技术领域，具体涉及一种面向综合监管的海事巡航路径动态开辟方法及系统。
背景技术
[0002] 如何动态生成最佳路径，满足在有限投入条件下，同时兼顾监管水域覆盖更广阔、到达更深远等多个管理目标，是有效巡航的关键。
[0003] 如公开号为CN111428916A的发明专利《一种海上救援船舶巡航路径规划方法》，针对应急场景，通过求解线性规划问题，提出了一种基于海洋环境数据可以确定救援船最优巡航路径的方法。公开号为CN113361765A的发明专利《一种船载无人机协同监测港口船舶大气污染路径规划的方法》，针对海上船舶对港口环境的污染，提供了一种船载无人机协同监测港口船舶大气污染路径规划的方法。曲小同等（2020年）以台州为实验水域，基于K‑means算法选择了28个海上任务点，然后运用蚁群算法开展了水域日常和应急立体协同巡航监管路径规划研究。王力等（2019年）以烟台为实验水域，运用聚类获取重点巡航范围后，基于混合整数线性规划方法开展了无人机协同巡航模型研究。
[0004] 已有算法均先基于需求或者聚类等算法确定巡航路径节点，将优化路径设计简化为传统的旅行商问题，即在可能连接节点的所有巡航路径中，选择代价最低且满足限制性条件的路径，亦即以最低的代价遍历所有接节点。这种方法忽略了巡航的过程性，点之间只有直线连接一种方式，是在有限解集对最优方案的选择。而面向综合监管的巡航路径规划中，为了最大化覆盖监管需求，往往需要根据全域需求覆盖情况，开辟新的路径，需要在优化路径设计时同步生成路径节点。目前的已有方法均不能实现节点的自动选择与新路径的开辟。
[0005] 因此，如何在一定的约束条件下（如航程/航时约束、传感器性能约束、飞行安全距离约束等），根据水域实时监管需求，集合各监管力量（如飞机、船艇等）的巡航能力，将监管业务分配给不同的监管力量，是亟需解决问题。
发明内容
[0006] 为实现本发明目的提供的一种面向综合监管的海事巡航路径动态开辟方法，包括：
[0007] 步骤S101：根据监管辖区的范围和实时动态特征、监管路径、巡航可视范围综合确定目标函数；
[0008] 步骤S102：根据巡航飞机或船舶的单次监管的最短距离、固定成本分别确定约束条件；
[0009] 步骤S103：基于目标函数和约束条件，在模拟退火算法的基础上开辟最佳巡航飞机或船舶监管路径，确定最大化满足监管需求的路径；
[0010] 包括以下步骤S1031‑S1035：
[0011] 步骤S1031：设定预期路径个数，并初始化监管路径并设定初始参数；
[0012] 步骤S1032：计算目标函数的初始值；
[0013] 步骤S1033：根据空间扰动产生新的监管路线；
[0014] 步骤S1034：根据Metropolis准则计算对所述新的监管路线的接受概率，确定下一次选择的路线；
[0015] 步骤S1035：判断当前路径是否是最佳路径，若是，则终止计算，并输出最佳路径；若不是，则重复步骤S1033‑步骤S1035。
[0016] 在其中一些具体实施例中，步骤S101还包括：对监管辖区的范围内的监管点进行动态网格离散处理，得到n个监管点，综合确定目标函数为：
[0017] ；
[0018] ；
[0019] 式中，pathes=（p1，p2，…，pk）为航线集合，表示监管辖区范围内海巡飞机或船舶分别的监管线路，其中，pi为第i条航线；Dn为监管辖区的范围内的监管点需求总量；dn为在确定路径pathes的情况下n个监管点中未被满足的需求量；N为航线上航迹点的总数。
[0020] 在其中一些具体实施例中，步骤S102中基于单个监管设施的单次监管的最短距离确定的约束条件l(p)为：
[0021] ；
[0022] 式中，A为单个飞机监管路径长度的公里数阈值；xj、yj分别为该点的经纬度坐标，xj‑1、yj‑1分别为当前点的前一个点的经纬度坐标，K为航线条数。
[0023] 在其中一些具体实施例中，步骤S102中基于固定成本确定的约束条件cost为：
[0024] ；
[0025] 式中，B为固定成本；L1为飞机航线个数；L2为海巡船舶航线个数；pa为第a个飞机航线；pc为第c个船舶航线；C1为巡航飞机每公里代价，C2为海巡船舶每公里代价。
[0026] 在其中一些具体实施例中，步骤S1033包括：
[0027] 选择当前路径Pcur中的一条路径pi；其次，在pi的k‑2个中间节点中，随机选择pi中任意一个监管点sj，其空间坐标为（xj，yj），以随机的距离r和角度v扰动sj的坐标位置产生新的空间坐标snew，所述新的空间坐标计算如下式所示：
[0028] ；
[0029] ；
[0030] 再使用snew替换sj，生成新的航线pnewi，替换Pcur中的路径pi，生成新的监管路线方案Pnew。
[0031] 为实现本发明目的，本申请还提供了一种面向综合监管的海事巡航路径动态开辟系统，包括：
[0032] 目标函数定义模块：用于根据监管辖区的范围和实时动态特征、监管路径、巡航可视范围综合确定目标函数；
[0033] 约束函数生成模块：用于根据单个监管设施的单次监管的最短距离、固定成本分别确定约束条件；
[0034] 最佳路径生成模块：用于基于所述目标函数和约束条件，根据模拟退火算法开辟巡航飞机或船舶监管路径，生成最佳路径；
[0035] 所述最佳路径生成模块包括：
[0036] 初始化单元：用于设定预期路径个数，并初始化监管路径并设定初始参数；
[0037] 初始值计算单元：用于计算目标函数的初始值；
[0038] 路线生成单元：用于根据空间扰动产生新的监管路线；
[0039] 路线选择单元：用于根据Metropolis准则计算对所述新的监管路线的接受概率，确定下一次选择的路线；
[0040] 最佳路径确定单元：用于判断当前路径是否是最佳路径，若是，则终止计算，并输出最佳路径；若不是，则跳转至所述路径生成单元重新生成新的监管路线。
[0041] 在其中一些具体实施例中，目标函数定义模块还用于：对监管辖区的范围内的监管点进行动态网格离散处理，得到n个监管点，综合确定目标函数为：
[0042] ；
[0043] ；
[0044] 式中，pathes=（p1，p2，…，pk）为航线集合，表示监管辖区范围内海巡飞机或船舶分别的监管线路，其中，pi为第i条航线；Dn为监管辖区的范围内的监管点需求总量；dn为在确定路径pathes的情况下n个监管点中未被满足的需求量；N为航线上航迹点的总数。
[0045] 在其中一些具体实施例中，约束函数生成模块用于基于单个监管设施的单次监管的最短距离确定的约束条件l(p)为：
[0046] ；
[0047] 式中，A为单个飞机监管路径长度的公里数阈值；xj、yj分别为该点的经纬度坐标，xj‑1、yj‑1分别为当前点的前一个点的经纬度坐标，K为航线条数。
[0048] 在其中一些具体实施例中，约束函数生成模块用于基于固定成本确定的约束条件cost为：
[0049] ；
[0050] 式中，B为固定成本；L1为飞机航线个数；L2为海巡船舶航线个数；pa为第a个飞机航线；pc为第c个船舶航线；C1为巡航飞机每公里代价，C2为海巡船舶每公里代价。
[0051] 在其中一些具体实施例中，路线生成单元用于按照以下方式产生新的监管路线：
[0052] 选择当前路径Pcur中的一条路径pi；其次，在pi的k‑2个中间节点中，随机选择pi中任意一个监管点sj，其空间坐标为（xj，yj），以随机的距离r和角度v扰动sj的坐标位置产生新的空间坐标snew，所述新的空间坐标计算如下式所示：
[0053] ；
[0054] ；
[0055] 再使用snew替换sj，生成新的航线pnewi，替换Pcur中的路径pi，生成新的监管路线方案Pnew。
[0056] 本发明的有益效果：
[0057] 本发明的一种面向综合监管的海事巡航路径动态开辟方法及系统，根据辖区水域的实时动态监管需求，能够实时计算辖区监管力量的需求满足能力量化数值。并在给定投资和安全航行距离等约束条件下，开辟最佳的巡航推荐路径。
附图说明
[0058] 为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0059] 图1是本发明一种面向综合监管的海事巡航路径动态开辟方法一些具体实施例的流程示意图；
[0060] 图2是本发明一种面向综合监管的海事巡航路径动态开辟方法的生成监管路径的流程示意图；
[0061] 图3是本发明一种面向综合监管的海事巡航路径动态开辟系统一些具体实施例的结构示意图；
[0062] 图4是本发明一种面向综合监管的海事巡航路径动态开辟系统的生成监管路径的结构示意图。
具体实施方式
[0063] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。
[0064] 所述实施例的示例在附图中示出，其中自始至终相同或类似的符号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
[0065] 参照图1所示，一种面向综合监管的海事巡航路径动态开辟方法，包括：
[0066] 步骤S101：根据监管辖区的范围和实时动态特征、监管路径、巡航可视范围综合确定目标函数。
[0067] 在本发明一些具体实施例中，步骤S101还包括：对监管辖区的范围内的监管点进行动态网格离散处理，得到n个监管点，综合确定目标函数为：
[0068] ；
[0069] ；
[0070] 式中，pathes=（p1，p2，…，pk）为航线集合，表示监管辖区范围内海巡飞机或船舶分别的监管线路，其中，pi为第i条航线；Dn为监管辖区的范围内的监管点需求总量；dn为在确定路径pathes的情况下n个监管点中未被满足的需求量；N为航线上航迹点的总数。
[0071] 具体的，巡航力量集合，L=（1,2,…,l），表示辖区拥有的海巡船舶或飞机总数。
[0072] 航线集合，pathes=（p1,p2,…,pk），表示辖区各海巡船舶或飞机总数分别的监管线路，是k条航线的集合。其中，pi为第i条航线。
[0073] 航迹点集合， ，s1为其出发基地和返回的目的地。
[0074] 航迹点sj=（xj，yj）。航迹点sj为巡航区域中所开辟的路径节点，xj，yj分别为该点的经纬度坐标。
[0075] N为区域内离散格网的数目，n为遍历中的某一个离散格网。
[0076] 步骤S102：根据巡航飞机或船舶的单次监管的最短距离、固定成本分别确定约束条件。
[0077] 在本发明一些具体实施例中，令起始点与终止点均为监管点，即s1=sk=s0。根据飞机性能和执法人员疲劳等因素，尽量减少远航程的航线条数和航次，因此单个飞机监管路径长度要求小于A公里，基于单个监管设施的单次监管的最短距离确定的约束条件l(p)为：
[0078] ；
[0079] 式中，A为单个飞机监管路径长度的公里数阈值；xj、yj分别为该点的经纬度坐标，xj‑1、yj‑1分别为当前点的前一个点的经纬度坐标，K为航线条数，根据历史统计数据设计总监管费用cost约束条件为：
[0080] ；
[0081] 其中，L1为飞机航线个数，L2为海巡船舶航线个数，pa为第a个飞机航线，pc为第c个船舶航线，C1为巡航飞机每公里代价，C2为海巡船舶巡航每公里代价。
[0082] 步骤S103：基于目标函数和约束条件，根据模拟退火算法开辟巡航飞机或船舶监管路径，确定最佳路径。
[0083] 在本发明一些具体实施例中，参照图2所示，
[0084] 步骤S1031：设定预期路径个数，并初始化监管路径并设定初始参数。
[0085] 具体的，设定预期路径个数，生成P1条初始的飞机路径和P2条初始的船舶路径，所有路径的方案为pathes；其次，以40公里为间隔，加密所有路径中的节点，设置当前路径Pcur=pathes；最后优化初始参数，包括初始温度T0=1，当前温度t=T0退火速率c=0.99，最低温度‑7Tmin=1×10 最大迭代次数为itermax=10000，连续拒绝次数为conRejmax=100，当前迭代次数为iter=0，当前连续拒绝次数conRej=0。
[0086] 步骤S1032：计算目标函数的初始值。
[0087] 具体的，Pcur中的每一条路径，按监管手段d对于监管任务m的满足程度 建立视域缓冲区（为简化处理，飞机设为800公里，船舶为600公里），生成单条航线轨迹的可视域。其次，将缓冲区与监管需求栅格进行空间叠置，将缓冲区内部监管需求置位0，获得更新后的需求值dn(pathes)；
[0088] ；
[0089] 计算海巡船和飞机覆盖后监管需求的剩余值Ocur，作为初始目标函数值，如下式所示：
[0090] 。
[0091] 步骤S1033：根据空间扰动产生新的监管路线。
[0092] 具体的，随机选择Pcur中的一条路径pi；其次，在pi的k‑2个中间节点中，随机选择pi中任意一个监管点sj，其空间坐标为（xj，yj），以随机的距离r和角度v扰动sj的坐标位置，产生snew，其坐标计算如下式所示：
[0093] ；
[0094] ；
[0095] 然后，使用snew替换sj，生成新的航线pnewi，替换Pcur中的路径pi，生成新的监管路线方案Pnew；进一步的，对Pnew中的每一条路径，按照视域做缓冲区，其中飞机为800公里，船舶为600公路，将缓冲区与监管需求栅格进行空间叠置，将缓冲区内部监管需求置位0，获得更新后的需求值Onew，如下所示：
[0096] 。
[0097] 步骤S1034：根据Metropolis准则计算对新的监管路线的接受概率，确定下一次选择的路线。
[0098] 具体的，根据下式依次计算每一个目标对s2的接受概率ac：
[0099] ；
[0100] 利用计算机程序随机生成一个0‑1之间的随机数rand，如果rand
[0101] 步骤S1035：判断当前路径是否是最佳路径，若是，则终止计算；若不是，则重复步骤S1033‑步骤S1035。
[0102] 具体的，如果conRe≥conRejmax，或者iter≥itermax，或者t
[0103] 为实现本发明目的，本申请还提供了一种面向综合监管的海事巡航路径动态开辟系统，参照图3所示，包括：
[0104] 目标函数定义模块10：用于根据监管辖区的范围和实时动态特征、监管路径、巡航可视范围综合确定目标函数。
[0105] 在本发明一些具体实施例中，对监管辖区的范围内的监管点进行动态网格离散处理，得到n个监管点，综合确定目标函数为：
[0106] ；
[0107] ；
[0108] 式中，pathes=（p1，p2，…，pk）为航线集合，表示监管辖区范围内海巡飞机或船舶分别的监管线路，其中，pi为第i条航线；Dn为监管辖区的范围内的监管点需求总量；dn为在确定路径pathes的情况下n个监管点中未被满足的需求量；N为航线上航迹点的总数。
[0109] 具体的，巡航力量集合，L=（1,2,…,l），表示辖区拥有的海空力量个数。
[0110] 航线集合，pathes=（p1,p2,…,pk），表示辖区各个海空力量（船舶和飞机）分别的监管线路，是k条航线的集合。其中，pi为第i条航线。
[0111] 航迹点集合， ，s1为其出发基地和返回的目的地。
[0112] 航迹点sj=（xj，yj）。航迹点sj为巡航区域中所开辟的路径节点，xj，yj分别为该点的经纬度坐标。
[0113] N为区域内离散格网的数目，n为遍历中的某一个离散格网。
[0114] 约束函数生成模块20：用于根据巡航飞机或船舶的单次监管的最短距离、固定成本分别确定约束条件。
[0115] 在本发明一些具体实施例中，令起始点与终止点均为监管点，即s1=sk=s0。根据飞机性能和执法人员疲劳等因素，尽量减少远航程的航线条数和航次，因此单个飞机监管路径长度要求小于A公里，基于单个监管设施的单次监管的最短距离确定的约束条件l(p)为：
[0116] ；
[0117] 式中，A为单个飞机监管路径长度的公里数阈值；xj、yj分别为该点的经纬度坐标，xj‑1、yj‑1分别为当前点的前一个点的经纬度坐标，K为航线条数。
[0118] 根据历史统计数据设计总监管费用cost约束条件：
[0119] ；
[0120] 其中，L1为飞机航线个数，L2为海巡船舶航线个数，pa为第a个飞机航线，pc为第c个船舶航线，C1为巡航飞机每公里代价，C2为海巡船舶每公里代价。
[0121] 最佳路径生成模块30：用于基于目标函数和约束条件，根据模拟退火算法开辟巡航飞机或船舶监管路径，确定最佳路径。
[0122] 在本发明一些具体实施例中，参照图4所示，
[0123] 初始化单元31：用于设定预期路径个数，并初始化监管路径并设定初始参数。
[0124] 具体的，设定预期路径个数，生成P1条初始的飞机路径和P2条初始的船舶路径，所有路径的方案为pathes；其次，以40公里为间隔，加密所有路径中的节点，设置当前路径Pcur=pathes；最后优化初始参数，包括初始温度T0=1，当前温度t=T0退火速率c=0.99，最低温度‑7Tmin=1×10 最大迭代次数为itermax=10000，连续拒绝次数为conRejmax=100，当前迭代次数为iter=0，当前连续拒绝次数conRej=0。
[0125] 初始值计算单元32：用于计算目标函数的初始值。
[0126] 具体的，Pcur中的每一条路径，按监管手段d对于监管任务m的满足程度 建立视域缓冲区（为简化处理，飞机设为800公里，船舶为600公里），生成单条航线轨迹的可视域。其次，将缓冲区与监管需求栅格进行空间叠置，将缓冲区内部监管需求置位0，获得更新后的需求值dn(pathes)；
[0127] ；
[0128] 计算海巡船和飞机覆盖后监管需求的剩余值Ocur，作为初始目标函数值，如下式所示：
[0129] 。
[0130] 路线生成单元33：用于根据空间扰动产生新的监管路线。
[0131] 具体的，首先，随机选择Pcur中的一条路径pi；其次，在pi的k‑2个中间节点中，随机选择pi中任意一个监管点sj，其空间坐标为（xj，yj），以随机的距离r和角度v扰动sj的坐标位置，产生snew，其坐标计算如下式所示：
[0132] ；
[0133] ；
[0134] 然后，使用snew替换sj，生成新的航线pnewi，替换Pcur中的路径pi，生成新的监管路线方案Pnew；进一步的，对Pnew中的每一条路径，按照视域做缓冲区，其中飞机为800公里，船舶为600公路，将缓冲区与监管需求栅格进行空间叠置，将缓冲区内部监管需求置位0，获得更新后的需求值Onew，如下所示：
[0135] 。
[0136] 路线选择单元34：用于根据Metropolis准则计算对新的监管路线的接受概率，确定下一次选择的路线。
[0137] 具体的，根据下式依次计算每一个目标对s2的接受概率ac：
[0138] ；
[0139] 利用计算机程序随机生成一个0‑1之间的随机数rand，如果rand
[0140] 最佳路径确定单元35：用于判断当前路径是否是最佳路径，若是，则终止计算，并输出最佳路径；若不是，则跳转至路径生成单元重新生成新的监管路线。
[0141] 具体的，如果conRe≥conRejmax，或者iter≥itermax，或者t
[0142] 本发明的一种面向综合监管的海事巡航路径动态开辟方法及系统根据辖区水域的实时动态监管需求或风险值，能够实时计算辖区监管力量的需求满足能力量化数值，并在给定投资和安全航行距离等约束条件下，给出最优的巡航推荐路径。
[0143] 在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、“一个具体实施例”或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
[0144] 最后应说明的是：以上实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不驱使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。
