CN116050720B - 模拟仿真作业过程策略的航道疏浚船舶配比系统及方法 - Google Patents
模拟仿真作业过程策略的航道疏浚船舶配比系统及方法 Download PDFInfo
- Publication number
- CN116050720B CN116050720B CN202210530592.2A CN202210530592A CN116050720B CN 116050720 B CN116050720 B CN 116050720B CN 202210530592 A CN202210530592 A CN 202210530592A CN 116050720 B CN116050720 B CN 116050720B
- Authority
- CN
- China
- Prior art keywords
- dredger
- ship
- engineering
- dredgers
- channel
- 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 83
- 230000008569 process Effects 0.000 title claims abstract description 60
- 238000004088 simulation Methods 0.000 claims abstract description 75
- 238000007664 blowing Methods 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 17
- 230000008859 change Effects 0.000 claims description 9
- 125000004122 cyclic group Chemical group 0.000 claims description 8
- 238000011166 aliquoting Methods 0.000 claims description 4
- 239000000969 carrier Substances 0.000 claims description 4
- 239000010802 sludge Substances 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 15
- 238000004422 calculation algorithm Methods 0.000 description 10
- 238000010276 construction Methods 0.000 description 4
- 238000012876 topography Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000013049 sediment Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000008439 repair process Effects 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/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
- G06Q10/06313—Resource planning in a project environment
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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/40—Business processes related to the transportation industry
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- General Business, Economics & Management (AREA)
- Computer Hardware Design (AREA)
- Marketing (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Biodiversity & Conservation Biology (AREA)
- Automation & Control Theory (AREA)
- Game Theory and Decision Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种模拟仿真作业过程策略的航道疏浚船舶配比系统及方法,包括:模拟仿真模块:依据输入信息模拟仿真作业过程,计算每个时间单位下船舶分布情况、等待情况,以及任务推进情况,估算总工程完成时间,输出该船舶配比下的总工程成本;最优配比计算模块:将船舶配比输入模拟仿真模块,得到在各个配比下的总工程成本,并输出最小工程成本对应的船舶配比作为当前工程环境下的解决方案,计算每种最佳船舶配比方案。本发明创新性地设计了一种航道疏浚工程中的船舶作业运输模拟模型;使得模型对各类工程环境均有良好的仿真能力与准确性;相比于经验公式的方法能够寻找到更优的配比方案与更低的成本。
Description
技术领域
本发明涉及内河航道疏浚工程规划领域,具体地,涉及一种模拟仿真作业过程策略的航道疏浚船舶配比系统及方法。
背景技术
内河的航道疏浚在河流较多的城市中一直是一项重要的工程作业任务。航道疏浚作业需要挖泥船、运泥船和吹泥船协同配合完成,挖泥船在疏浚处挖泥,并将淤泥装至运泥船上,由运泥船运至指定排放地点,再由吹泥船将运泥船上的淤泥吹填至指定位置,最后运泥船返回疏浚处开始下一次循环过程。
为了减少疏浚工程的作业成本,在施工前会依据各类船舶的施工效率、装载容量和租用成本设计三类船舶的配比,以在满足工期要求的情况下尽可能实现成本的最低化。当前船舶配比一般采取依据施工效率和运输效率之间的比例形成配比,并按照经验添加作业时常、闸门开放时间以及各类突发因素对疏浚作业速度的影响参数。但是在实际工程中,依据经验对设计时确定的船舶配比进行修正一般并不准确,且并不能依据各类工程环境的变化而灵活有效地调整影响参数。
专利文献CN114117609A(申请号:CN202111435257.6)公开了一种航道整治方法及装置,装置包括:采用RTK配合多波束测深系统测量河口水下地形,根据地形确定拦门沙坎地理位置坐标,分析水下地形特性;利用水沙二相数值模拟系统,根据河口水下地形搭建二维水动力学模型;利用二维水动力学模型计算本底流场,确定水深和主流分布特征;根据水深和主流分布特征,结合航道设计指标确定拦门沙坎泥沙的疏浚范围和疏浚规模;根据径流动力条件和机械扰动参数、工程布局,利用二维水动力学模型模拟疏浚效果,优化疏浚范围和疏浚规模;根据优化后的疏浚范围和疏浚规模选定泥沙机械扰动采用的疏浚船型以及船只数量,以通过船只进行疏浚整治。但该发明没有通过设计航道疏浚工程中的船舶作业运输模拟模型的方法进行航道疏浚船舶配比计算。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种模拟仿真作业过程策略的航道疏浚船舶配比系统及方法。
根据本发明提供的一种模拟仿真作业过程策略的航道疏浚船舶配比系统,包括:
模拟仿真模块:依据输入信息模拟仿真作业过程,计算每个时间单位下船舶分布情况、等待情况,以及任务推进情况,估算总工程完成时间,输出该船舶配比下的总工程成本;
最优配比计算模块:将船舶配比输入模拟仿真模块,得到在各个配比下的总工程成本,并输出最小工程成本对应的船舶配比作为当前工程环境下的解决方案,计算每种最佳船舶配比方案。
优选地,在所模拟仿真模块中:
将总共疏浚工程量、运泥船、挖泥船和吹泥船的数量和租用成本信息、航道长度和满载、空载航行速度信息以及航道上闸门开放时间信息输入模块中,依据输入信息模拟仿真挖泥、运泥和吹泥的作业过程,并计算每个时间单位下运泥船在航道上的分布情况、在挖泥处和吹泥处的等待情况以及随着时间变化工程任务的推进情况,估算总工程完成时间,并依据船舶租金输出该船舶配比下的总工程成本。
优选地,运泥船在航道上分布的模拟过程:
以预设时间段为一个时间单位,设航道长度为L,满载航行速度为v1,空载航行速度为v2,满载航行时间t1为:
t1=L/v1
空载航行时间t2为:
t2=L/v2
航行时间t1和t2的计量单位为预设时间段,且向上取整;
通过将航道按照航行时间分段的方式来表示整条航道上不同位置运泥船的分布情况,已知满载航行时间为t1,空载航行时间为t2,将去程航道t1等分,将回程航道t2等分,每一等分即表示预设时间段内出发的运泥船的数量,各个等分能够表示运泥船当前时刻在航道上对应位置的分布情况。
优选地,运泥船在航道上航行的模拟过程:
模拟随着时间推进的航道疏浚作业情况,通过循环的方式不断更新时间,每一次循环代表度过了一个预设时间段,在上一次循环结束后,若此时闸门是开启状态,则所有运泥船向下一格推进一个位置,新的位置表示下一个预设时间段运泥船所处的位置;
去程第一等分新添加挖泥船在过去预设时间段内装满运泥船的数量,回程第一等分新添加吹泥船在过去预设时间段内卸空运泥船的数量,具体数量为
n1=m1p1/v
n2=m2p2/v
其中,n1为去程新加入的运泥船数量,n2为回程新加入的运泥船数量,v为运泥船容积,m1和m2分别为挖泥船和吹泥船的数量,p1和p2分别为挖泥船的挖泥速度和吹泥船的吹泥速度,由于运泥船数量必须为整数,因此n1和n2的计算结果自动作向下取整,由于向下取整导致缺少统计的工程量会被叠加,当累计叠加量达到一艘运泥船的容量时,则该时刻的n1和n2在计算结果的基础上加一,并将累计叠加量清零;最后一等分在下一预设时间段自动被添加进入等待区队列中。
优选地,运泥船在闸门处的模拟过程:
当闸门关闭时,行驶到此处的运泥船需要在闸门处等待其下一次开启,设闸门在往返航道上所处位置分别为x1和x2,则除了闸门所处位置外,其余均向下一分位推进一个位置,x1和x2位置的运泥船数量在下一时刻由于闸门关闭等待所以不推进,而x1-1和x2-1分位由于继续向前推进,因此在下一时刻,x1和x2位置的运泥船数量分别为上一时刻x1和x1-1两个位置运泥船数量之和,而x1+1和x2+1位置的运泥船数量为0。
优选地,挖泥处和吹泥处的模拟过程:
挖泥处和吹泥处用来完成对各自等待区船舶进行装泥或吹泥处理,并更新剩余工程作业量,设挖泥处等待船舶数量为N1,吹泥处等待船舶数量为N2,且预设时间段内挖泥量和吹泥量分别为m1p1和m2p2,若
m1p1>N1v
或
m2p2>N2v
则在单位时间内,挖泥船或吹泥船的作业速度大于运泥船在挖泥或吹泥等待区的等待工程量,在下一时刻等待区运泥船会全部清零,当前等待区的运泥船在下一时刻会全部被卸泥或装泥完毕,故在此条件下
n1=N1
或
n2=N2
随后更新当前时刻的N1和N2分别为回程和去程上一时刻最后等分位中运泥船的数量;
若等待队列中的船舶数量过多,那么单位时间内挖泥船和运泥船的作业速度并不足以将等待中的船舶全部清空,故仍然有
n1=m1p1/v
n2=m2p2/v
同时N1和N2增加量为上一时刻最后等分位的运泥船数量n1和n2;
对剩余工程量进行更新,更新按照挖泥处剩余需要疏浚的方量为准,当每一时刻挖泥处的挖泥船由等待区装泥完成进入去程后,更新当前时刻疏浚完成的方量,剩余方量为M,则更新后的剩余方量为M-n1v,当更新后剩余方量小于0时,疏浚作业结束,跳出循环并输出作业完成总时间和总成本。
优选地,在所最优配比计算模块中:
将挖泥船、运泥船和吹泥船的配比通过循环的方式依次输入模拟仿真模块,得到在各个配比下的总工程成本,并输出最小工程成本对应的船舶配比作为当前工程环境下的解决方案,通过改变闸门开放时间,计算每种最佳船舶配比方案。
优选地,循环寻找最优配比:
依据实际工程中挖泥船、吹泥船和运泥船的可租用数量,为满足基本的作业需求,三者的可行数量下限均为1艘,循环寻找最低成本下最优的船舶配比。
优选地,闸门开放时间与船舶最优配比的关系:
基于模拟仿真模块,通过模拟的方式得到在改变闸门开放时间的情况下的三种船舶的最佳配比情况,最佳船舶配比不能通过闸门开放时间精确计算得到,总工程租金y与闸门开放时间x成负相关关系。
根据本发明提供的一种模拟仿真作业过程策略的航道疏浚船舶配比方法,采用所述的模拟仿真作业过程策略的航道疏浚船舶配比系统,执行包括:
步骤S1:设计可调节变量,包括船舶数量、租金成本、航行速度、作业效率、航道长度、闸门开关时间;
步骤S2:设计模拟船舶作业运输过程的模型,模型能够模拟某一时刻运泥船在航道上的分布情况以及在挖泥处、吹泥处的等待情况;
步骤S3:基于船舶作业运输过程的模拟结果估计给定工程量下的总成本,并按照成本最小化的原则通过循环迭代的方式寻找并给出最佳船舶配比。
与现有技术相比,本发明具有如下的有益效果:
1、本发明创新性地设计了一种航道疏浚工程中的船舶作业运输模拟模型;
2、本发明将航道长度、航速变化、闸门开关时间等作为输入参数,使得模型对各类工程环境均有良好的仿真能力与准确性;
3、本发明通过模拟仿真的方式估算某一船舶配比下的工程成本,并通过循环迭代寻找最低工程成本下的船舶配比,相比于经验公式的方法能够寻找到更优的配比方案与更低的成本;
4、本发明详细考虑了闸门开放时间对最优船舶配比的影响,并给出二者和总工程租金的具体关系;
5、本发明计算的配比成本相比于当前基准算法的配比成本降低。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明的算法框架图;
图2为本发明模拟仿真模块的实现示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
实施例1:
本发明涉及内河航道疏浚工程规划领域,具体地,涉及以节约工程成本为目标的用于内河航道疏浚的挖泥船、吹泥船和运泥船配比设计方案,更为具体地,涉及基于模拟仿真船舶作业运输过程的船舶配比设计算法。
本发明公开了一种基于模拟仿真作业过程策略的航道疏浚船舶配比算法。为了尽可能节约航道疏浚的工程成本,需要充分利用挖泥船、运泥船和吹泥船。由于航道上闸门的存在,运泥船在特定时间内无法顺利通航,导致基于作业速度配比船舶的策略并不是最优解。本发明设计模型模拟运泥船、挖泥船和吹泥船的作业过程,并得到总工程时间和总成本;而后本发明通过循环迭代的思想利用模拟模型得到成本最低的船舶配比,经过验证,本发明计算的配比成本相比于当前基准算法的配比成本降低大约15%,实际工程中考虑挖泥船数据不足可大约降低10%。最后,本发明着重讨论了航道上闸门开放时间对船舶配比的影响。
根据本发明提供的一种模拟仿真作业过程策略的航道疏浚船舶配比系统,如图1-图2所示,包括:
模拟仿真模块:依据输入信息模拟仿真作业过程,计算每个时间单位下船舶分布情况、等待情况,以及任务推进情况,估算总工程完成时间,输出该船舶配比下的总工程成本;
具体地,在所模拟仿真模块中:
将总共疏浚工程量、运泥船、挖泥船和吹泥船的数量和租用成本信息、航道长度和满载、空载航行速度信息以及航道上闸门开放时间信息输入模块中,依据输入信息模拟仿真挖泥、运泥和吹泥的作业过程,并计算每个时间单位下运泥船在航道上的分布情况、在挖泥处和吹泥处的等待情况以及随着时间变化工程任务的推进情况,估算总工程完成时间,并依据船舶租金输出该船舶配比下的总工程成本。
具体地,运泥船在航道上分布的模拟过程:
以预设时间段为一个时间单位,设航道长度为L,满载航行速度为v1,空载航行速度为v2,满载航行时间t1为:
t1=L/v1
空载航行时间t2为:
t2=L/v2
航行时间t1和t2的计量单位为预设时间段,且向上取整;
通过将航道按照航行时间分段的方式来表示整条航道上不同位置运泥船的分布情况,已知满载航行时间为t1,空载航行时间为t2,将去程航道t1等分,将回程航道t2等分,每一等分即表示预设时间段内出发的运泥船的数量,各个等分能够表示运泥船当前时刻在航道上对应位置的分布情况。
具体地,运泥船在航道上航行的模拟过程:
模拟随着时间推进的航道疏浚作业情况,通过循环的方式不断更新时间,每一次循环代表度过了一个预设时间段,在上一次循环结束后,若此时闸门是开启状态,则所有运泥船向下一格推进一个位置,新的位置表示下一个预设时间段运泥船所处的位置;
去程第一等分新添加挖泥船在过去预设时间段内装满运泥船的数量,回程第一等分新添加吹泥船在过去预设时间段内卸空运泥船的数量,具体数量为
n1=m1p1/v
n2=m2p2/v
其中,n1为去程新加入的运泥船数量,n2为回程新加入的运泥船数量,v为运泥船容积,m1和m2分别为挖泥船和吹泥船的数量,p1和p2分别为挖泥船的挖泥速度和吹泥船的吹泥速度,由于运泥船数量必须为整数,因此n1和n2的计算结果自动作向下取整,由于向下取整导致缺少统计的工程量会被叠加,当累计叠加量达到一艘运泥船的容量时,则该时刻的n1和n2在计算结果的基础上加一,并将累计叠加量清零;最后一等分在下一预设时间段自动被添加进入等待区队列中。
具体地,运泥船在闸门处的模拟过程:
当闸门关闭时,行驶到此处的运泥船需要在闸门处等待其下一次开启,设闸门在往返航道上所处位置分别为x1和x2,则除了闸门所处位置外,其余均向下一分位推进一个位置,x1和x2位置的运泥船数量在下一时刻由于闸门关闭等待所以不推进,而x1-1和x2-1分位由于继续向前推进,因此在下一时刻,x1和x2位置的运泥船数量分别为上一时刻x1和x1-1两个位置运泥船数量之和,而x1+1和x2+1位置的运泥船数量为0。
具体地,挖泥处和吹泥处的模拟过程:
挖泥处和吹泥处用来完成对各自等待区船舶进行装泥或吹泥处理,并更新剩余工程作业量,设挖泥处等待船舶数量为N1,吹泥处等待船舶数量为N2,且预设时间段内挖泥量和吹泥量分别为m1p1和m2p2,若
m1p1>N1v
或
m2p2>N2v
则在单位时间内,挖泥船或吹泥船的作业速度大于运泥船在挖泥或吹泥等待区的等待工程量,在下一时刻等待区运泥船会全部清零,当前等待区的运泥船在下一时刻会全部被卸泥或装泥完毕,故在此条件下
n1=N1
或
n2=N2
随后更新当前时刻的N1和N2分别为回程和去程上一时刻最后等分位中运泥船的数量;
若等待队列中的船舶数量过多,那么单位时间内挖泥船和运泥船的作业速度并不足以将等待中的船舶全部清空,故仍然有
n1=m1p1/v
n2=m2p2/v
同时N1和N2增加量为上一时刻最后等分位的运泥船数量n1和n2;
对剩余工程量进行更新,更新按照挖泥处剩余需要疏浚的方量为准,当每一时刻挖泥处的挖泥船由等待区装泥完成进入去程后,更新当前时刻疏浚完成的方量,剩余方量为M,则更新后的剩余方量为M-n1v,当更新后剩余方量小于0时,疏浚作业结束,跳出循环并输出作业完成总时间和总成本。
最优配比计算模块:将船舶配比输入模拟仿真模块,得到在各个配比下的总工程成本,并输出最小工程成本对应的船舶配比作为当前工程环境下的解决方案,计算每种最佳船舶配比方案。
具体地,在所最优配比计算模块中:
将挖泥船、运泥船和吹泥船的配比通过循环的方式依次输入模拟仿真模块,得到在各个配比下的总工程成本,并输出最小工程成本对应的船舶配比作为当前工程环境下的解决方案,通过改变闸门开放时间,计算每种最佳船舶配比方案。
具体地,循环寻找最优配比:
依据实际工程中挖泥船、吹泥船和运泥船的可租用数量,为满足基本的作业需求,三者的可行数量下限均为1艘,循环寻找最低成本下最优的船舶配比。
具体地,闸门开放时间与船舶最优配比的关系:
基于模拟仿真模块,通过模拟的方式得到在改变闸门开放时间的情况下的三种船舶的最佳配比情况,最佳船舶配比不能通过闸门开放时间精确计算得到,总工程租金y与闸门开放时间x成负相关关系。
根据本发明提供的一种模拟仿真作业过程策略的航道疏浚船舶配比方法,采用所述的模拟仿真作业过程策略的航道疏浚船舶配比系统,执行包括:
步骤S1:设计可调节变量,包括船舶数量、租金成本、航行速度、作业效率、航道长度、闸门开关时间;
步骤S2:设计模拟船舶作业运输过程的模型,模型能够模拟某一时刻运泥船在航道上的分布情况以及在挖泥处、吹泥处的等待情况;
步骤S3:基于船舶作业运输过程的模拟结果估计给定工程量下的总成本,并按照成本最小化的原则通过循环迭代的方式寻找并给出最佳船舶配比。
实施例2:
实施例2为实施例1的优选例,以更为具体地对本发明进行说明。
为解决工程环境变化多样导致船舶配比估算困难的问题,发明了一种基于模拟船舶作业运输过程的船舶配比设计算法,有效满足了对工程成本的估算与优化。本发明提出一种基于模拟船舶作业运输过程的船舶配比设计算法,着重提高船舶配比设计结果的可靠性与准确性。首先本发明设计了一个模型,该模型能够直接仿真船舶在航道上的分布变化与作业进度情况,并输出疏浚作业全部完成时的总工程成本;其次,本发明通过循环迭代过程获取到在给定作业量和船舶使用成本的条件下的最优船舶配比以及该配比下的最低工程成本。最后,本发明通过运行模拟仿真模型,着重讨论闸门开放时间与船舶配比的关系。
本发明提出了一种航道疏浚作业中的船舶配比设计算法,目的是通过模拟船舶作业运输过程,计算在给定船舶配比的情况下的工程成本,有效解决由于航道和工程环境多变导致的船舶配比设计估算困难以及总工程成本难以进一步降低的问题。首先,本发明设计了尽可能多的可调节变量,包括船舶数量、租金成本、航行速度、作业效率、航道长度、闸门开关时间等作为输入,以尽可能满足各种工程环境;其次,本发明设计了一种模拟船舶作业运输过程的模型,该模型可以模拟某一时刻运泥船在航道上的分布情况以及在挖泥处、吹泥处的等待情况;最后,本发明基于船舶作业运输过程的模拟结果估计给定工程量下的总成本,并按照成本最小化的原则通过循环迭代的方式寻找并给出最佳船舶配比。
本发明的核心算法框架如图1所示,主要由两个模块构成。
第一部分是模拟仿真模块。将总共疏浚工程量、运泥船、挖泥船和吹泥船的数量和租用成本信息、航道长度和满载、空载航行速度信息以及航道上闸门开放时间信息输入模块中,其将会依据输入信息模拟仿真挖泥、运泥和吹泥的作业过程,并计算每个时间单位下运泥船在航道上的分布情况、在挖泥处和吹泥处的等待情况,以及随着时间变化工程任务的推进情况,最后估算总工程完成时间,并依据船舶租金输出该船舶配比下的总工程成本。
第二部分为最优配比计算模块。将可能存在的挖泥船、运泥船和吹泥船的配比通过循环的方式依次输入模拟仿真模块,得到在各个配比下的总工程成本,并输出最小工程成本对应的船舶配比作为当前工程环境下的解决方案。与此同时,通过改变闸门开放时间,计算每种最佳船舶配比方案,讨论二者的对应关系。
1.模拟仿真模块
(1)运泥船在航道上分布的模拟过程
本发明假定挖泥船和吹泥船作业速度恒定、运泥船的船速恒定,但由于运泥船在从挖泥处到吹泥处是满载状态,从吹泥处返回挖泥处是空载状态,因此往返船速有细微差别。本系统以1小时为一个时间单位,设航道长度为L,满载航行速度为v1,空载航行速度为v2,因此满载航行时间
t1=L/v1
空载航行时间
t2=L/v2
由于本模块的时间单位为1小时,因此航行时间t1和t2的计量单位均为小时,且向上取整,以此即能保证随后用数组记录航道上运泥船的分布情况。
本发明通过将航道按照航行时间分段的方式来表示整条航道上不同位置运泥船的分布情况,具体实现示意图如图2所示。已知满载航行时间为t1,空载航行时间为t2,故将去程航道t1等分,将回程航道t2等分,每一等分即表示一小时内出发的运泥船的数量,即各个等分能够表示运泥船当前时刻在航道上对应位置的分布情况。
(2)运泥船在航道上航行的模拟过程
本发明模拟了随着时间推进的航道疏浚作业情况,其通过循环的方式不断更新时间,每一次循环代表度过了一个小时。在上一次循环结束后,若此时闸门是开启状态,则所有运泥船向下一格推进一个位置,新的位置表示下一小时运泥船所处的位置,
第一等分和最后一等分的处理方式较为特殊。去程和回程第一等分新添加的分别是挖泥船和吹泥船在过去一小时内装满或卸空运泥船的数量,具体数量为
n1=m1p1/v
n2=m2p2/v
其中n1和n2即分别为去程和回程新加入的运泥船数量,v为运泥船容积,m1和m2分别为挖泥船和吹泥船的数量,p1和p2分别为挖泥船的挖泥速度和吹泥船的吹泥速度,由于运泥船数量必须为整数,因此n1和n2的计算结果会自动作向下取整,由于向下取整导致缺少统计的工程量会被叠加,当累计叠加量达到一艘运泥船的容量时,则该时刻的n1和n2会在计算结果的基础上+1,并将累计叠加量清零。最后一等分在下一小时会自动被添加进入等待区队列中。
(3)运泥船在闸门处的模拟过程
前文所述为在闸门正常开启时的运泥船分布随时间的位置变化过程,当闸门关闭时,行驶到此处的运泥船需要在闸门处等待其下一次开启。设闸门在往返航道上所处位置分别为x1和x2,则除了闸门所处位置外,其余均向下一分位推进一个位置。x1和x2位置的运泥船数量在下一时刻由于闸门关闭等待所以不推进,而x1-1和x2-1分位由于继续向前推进,因此在下一时刻,x1和x2位置的运泥船数量分别为上一时刻x1和x1-1两个位置运泥船数量之和,而x1+1和x2+1位置的运泥船数量为0。
(4)挖泥处和吹泥处的模拟过程
挖泥处和吹泥处用来完成对各自等待区船舶进行装泥或吹泥处理,并更新剩余工程作业量。设挖泥处等待船舶数量为N1,吹泥处等待船舶数量为N2,且一小时内挖泥量和吹泥量分别为m1p1和m2p2。若
m1p1>N1v
或
m2p2>N2v
则在单位时间内,挖泥船或吹泥船的作业速度大于其在等待区的等待工程量,即在下一时刻等待区运泥船会全部清零,当前等待区的运泥船在下一时刻会全部被卸泥或装泥完毕,故在此条件下
n1=N1
或
n2=N2
随后更新当前时刻的N1和N2分别为回程和去程上一时刻最后等分位中运泥船的数量。
反之,若等待队列中的船舶数量过多,那么单位时间内挖泥船和运泥船的作业速度并不足以将等待中的船舶全部清空,故仍然有
n1=m1p1/v
n2=m2p2/v
同时N1和N2增加量为上一时刻最后等分位的运泥船数量n1和n2。
为了衡量当前工程完成度情况,最后需要对剩余工程量进行更新,更新按照挖泥处剩余需要疏浚的方量为准,即当每一时刻挖泥处的挖泥船由等待区装泥完成进入去程后,即更新当前时刻疏浚完成的方量,设剩余方量为M,则更新后的剩余方量为M-n1v,当更新后剩余方量小于0时,则认为疏浚作业结束,跳出循环并输出作业完成总时间和总成本。
2.最优配比计算模块
本发明的模拟仿真模块的功能是依据输入的三类船舶数量配比,输出其作业完成总时间和总成本,最优配比计算模块即是通过循环迭代运行模拟仿真模块的方式,寻找到总成本最低的船舶配比方案。
(1)循环寻找最优配比
首先,需要明确一点,总工程成本只与船舶配比有关,与各类船舶具体数量无关。若两类船舶配比方案完全相同,但数量相差k倍,则显然运泥船在航道上的数量分布、等待区的数量比和挖泥船、吹泥船的作业速度也相差k倍,总工程完工时间缩短为1/k,故总租金成本由于数量相差k倍但时间缩短至1/k因此保持不变,只需要考虑配比即可。
具体流程为:依据实际工程中挖泥船、吹泥船和运泥船的可租用数量,例如,设定挖泥船和吹泥船的可行数量上限为10艘,运泥船的可行数量上限为50艘;为满足基本的作业需求,三者的可行数量下限均为1艘,循环寻找最低成本下最优的船舶配比。
(2)闸门开放时间与船舶最优配比的关系
在传统方案中,不考虑闸门开关拦截对运泥船的影响,只需要考虑运泥船运量、挖泥船和吹泥船作业速度以及航道长度之间的比例关系即可大致推算出配比。然而,由于闸门开放时间的限制,导致运泥船必须要在特定的时间在闸门处等待开放通行,闸门开放时间的变化会极大且持续地影响运泥船在航道上的分布情况,即传统方法并不能准确计算在给定闸门开放时间的条件下的工程进展情况。本发明基于其所设计的模拟仿真模块,通过模拟的方式得到在改变闸门开放时间的情况下的三种船舶的最佳配比情况。例如,设定挖泥船和吹泥船的可行数量上限为10艘,运泥船的可行数量上限为50艘,迭代计算过程如下表1所示。
通过对上表的分析可知,最佳船舶配比不能通过闸门开放时间精确计算得到,然而显然总工程租金y与闸门开放时间x成明显的负相关关系,其拟合结果为
y=-12.639x+529.79
其中拟合R2=0.84,相关关系明显。
本发明采用Python3.10编写。设定挖泥船、吹泥船和运泥船每艘每天的租金成本为5174元、3736元和2700元,挖泥船和吹泥船的作业速度分别为154m3/h和500m3/h,运泥船的运载量为250m3。设定总工程量为225000m3,闸门开放时间为每天12h,空载航行时间为9h,满载航行时间为12h。本发明通过运算得到上述状态下挖泥船、吹泥船和运泥船的最佳配比为7:2:43,总成本为326.90万元,通过传统算法计算得到的配比为2:1:21,总成本为384.59万元,成本节约比例为15.0%。在实际工程中,由于挖泥船数量可能不足,因此如果配置比例同除以2向下取整变为3:1:21时,工程成本约为345.29万元,也能节约大致10.2%的租金成本。
本领域技术人员知道,除了以纯计算机可读程序代码方式实现本发明提供的系统、装置及其各个模块以外,完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统、装置及其各个模块以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同程序。所以,本发明提供的系统、装置及其各个模块可以被认为是一种硬件部件,而对其内包括的用于实现各种程序的模块也可以视为硬件部件内的结构;也可以将用于实现各种功能的模块视为既可以是实现方法的软件程序又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (5)
1.一种模拟仿真作业过程策略的航道疏浚船舶配比系统,其特征在于,包括:
模拟仿真模块:依据输入信息模拟仿真作业过程,计算每个时间单位下船舶分布情况、等待情况,以及任务推进情况,估算总工程完成时间,输出该船舶配比下的总工程成本;
最优配比计算模块:将船舶配比输入模拟仿真模块,得到在各个配比下的总工程成本,并输出最小总工程成本对应的船舶配比作为当前工程环境下的解决方案,计算每种最佳船舶配比方案;
在所模拟仿真模块中:
将总共疏浚工程量、运泥船、挖泥船和吹泥船的数量和租用成本信息、航道长度和满载、空载航行速度信息以及航道上闸门开放时间信息输入模块中,依据输入信息模拟仿真挖泥、运泥和吹泥的作业过程,并计算每个时间单位下运泥船在航道上的分布情况、在挖泥处和吹泥处的等待情况以及随着时间变化工程任务的推进情况,估算总工程完成时间,并依据船舶租金输出该船舶配比下的总工程成本;
运泥船在航道上分布的模拟过程:
以预设时间段为一个时间单位,设航道长度为L,满载航行速度为v1,空载航行速度为v2,满载航行时间t1为:
t1=L/v1
空载航行时间t2为:
t2=L/v2
航行时间t1和t2的计量单位为预设时间段,且向上取整;
通过将航道按照航行时间分段的方式来表示整条航道上不同位置运泥船的分布情况,已知满载航行时间为t1,空载航行时间为t2,将去程航道t1等分,将回程航道t2等分,每一等分即表示预设时间段内出发的运泥船的数量,各个等分能够表示运泥船当前时刻在航道上对应位置的分布情况;
运泥船在航道上航行的模拟过程:
模拟随着时间推进的航道疏浚作业情况,通过循环的方式不断更新时间,每一次循环代表度过了一个预设时间段,在上一次循环结束后,若此时闸门是开启状态,则所有运泥船向下一格推进一个位置,新的位置表示下一个预设时间段运泥船所处的位置;
去程第一等分新添加挖泥船在过去预设时间段内装满运泥船的数量,回程第一等分新添加吹泥船在过去预设时间段内卸空运泥船的数量,具体数量为
n1=m1p1/v
n2=m2p2/v
其中,n1为去程新加入的运泥船数量,n2为回程新加入的运泥船数量,v为运泥船容积,m1和m2分别为挖泥船和吹泥船的数量,p1和p2分别为挖泥船的挖泥速度和吹泥船的吹泥速度,由于运泥船数量必须为整数,因此n1和n2的计算结果自动作向下取整,由于向下取整导致缺少统计的工程量会被叠加,当累计叠加量达到一艘运泥船的容量时,则该时刻的n1和n2在计算结果的基础上加一,并将累计叠加量清零;最后一等分在下一预设时间段自动被添加进入等待区队列中;
运泥船在闸门处的模拟过程:
当闸门关闭时,行驶到此处的运泥船需要在闸门处等待其下一次开启,设闸门在往返航道上所处位置分别为x1和x2,则除了闸门所处位置外,其余均向下一分位推进一个位置,x1和x2位置的运泥船数量在下一时刻由于闸门关闭等待所以不推进,而x1-1和x2-1分位由于继续向前推进,因此在下一时刻,x1和x2位置的运泥船数量分别为上一时刻x1和x1-1两个位置运泥船数量之和,而x1+1和x2+1位置的运泥船数量为0;
挖泥处和吹泥处的模拟过程:
挖泥处和吹泥处用来完成对各自等待区船舶进行装泥或吹泥处理,并更新剩余工程作业量,设挖泥处等待船舶数量为N1,吹泥处等待船舶数量为N2,且预设时间段内挖泥量和吹泥量分别为m1p1和m2p2,若
m1p1>N1v
或
m2p2>N2v
则在单位时间内,挖泥船或吹泥船的作业速度大于运泥船在挖泥或吹泥等待区的等待工程量,在下一时刻等待区运泥船会全部清零,当前等待区的运泥船在下一时刻会全部被卸泥或装泥完毕,故在此条件下
n1=N1
或
n2=N2
随后更新当前时刻的N1和N2分别为回程和去程上一时刻最后等分位中运泥船的数量;
若等待队列中的船舶数量过多,那么单位时间内挖泥船和运泥船的作业速度并不足以将等待中的船舶全部清空,故仍然有
n1=m1p1/v
n2=m2p2/v
同时N1和N2增加量为上一时刻最后等分位的运泥船数量n1和n2;
对剩余工程量进行更新,更新按照挖泥处剩余需要疏浚的方量为准,当每一时刻挖泥处的挖泥船由等待区装泥完成进入去程后,更新当前时刻疏浚完成的方量,剩余方量为M,则更新后的剩余方量为M-n1v,当更新后剩余方量小于0时,疏浚作业结束,跳出循环并输出作业完成总时间和总成本。
2.根据权利要求1所述的模拟仿真作业过程策略的航道疏浚船舶配比系统,其特征在于,在所最优配比计算模块中:
将挖泥船、运泥船和吹泥船的配比通过循环的方式依次输入模拟仿真模块,得到在各个配比下的总工程成本,并输出最小工程成本对应的船舶配比作为当前工程环境下的解决方案,通过改变闸门开放时间,计算每种最佳船舶配比方案。
3.根据权利要求2所述的模拟仿真作业过程策略的航道疏浚船舶配比系统,其特征在于:
循环寻找最优配比:
依据实际工程中挖泥船、吹泥船和运泥船的可租用数量,为满足基本的作业需求,三者的可行数量下限均为1艘,循环寻找最低成本下最优的船舶配比。
4.根据权利要求2所述的模拟仿真作业过程策略的航道疏浚船舶配比系统,其特征在于:
闸门开放时间与船舶最优配比的关系:
基于模拟仿真模块,通过模拟的方式得到在改变闸门开放时间的情况下的三种船舶的最佳配比情况,最佳船舶配比不能通过闸门开放时间精确计算得到,总工程租金y与闸门开放时间x成负相关关系。
5.一种模拟仿真作业过程策略的航道疏浚船舶配比方法,其特征在于,采用权利要求1-4任一项所述的模拟仿真作业过程策略的航道疏浚船舶配比系统,执行包括:
步骤S1:设计可调节变量,包括船舶数量、租金成本、航行速度、作业效率、航道长度、闸门开关时间;
步骤S2:设计模拟船舶作业运输过程的模型,模型能够模拟某一时刻运泥船在航道上的分布情况以及在挖泥处、吹泥处的等待情况;
步骤S3:基于船舶作业运输过程的模拟结果估计给定工程量下的总成本,并按照成本最小化的原则通过循环迭代的方式寻找并给出最佳船舶配比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210530592.2A CN116050720B (zh) | 2022-05-16 | 2022-05-16 | 模拟仿真作业过程策略的航道疏浚船舶配比系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210530592.2A CN116050720B (zh) | 2022-05-16 | 2022-05-16 | 模拟仿真作业过程策略的航道疏浚船舶配比系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116050720A CN116050720A (zh) | 2023-05-02 |
CN116050720B true CN116050720B (zh) | 2023-09-05 |
Family
ID=86127790
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210530592.2A Active CN116050720B (zh) | 2022-05-16 | 2022-05-16 | 模拟仿真作业过程策略的航道疏浚船舶配比系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116050720B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104657538A (zh) * | 2015-01-04 | 2015-05-27 | 河海大学 | 一种多节点航道调度的仿真方法 |
CN106156428A (zh) * | 2016-07-06 | 2016-11-23 | 大连海事大学 | 一种基于样本分析的船舶交通流模拟生成方法 |
KR20220057164A (ko) * | 2020-10-29 | 2022-05-09 | 주식회사 에이스이앤티 | 설계데이터를 이용한 작업순서 및 작업배분 시뮬레이션 방법 |
-
2022
- 2022-05-16 CN CN202210530592.2A patent/CN116050720B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104657538A (zh) * | 2015-01-04 | 2015-05-27 | 河海大学 | 一种多节点航道调度的仿真方法 |
CN106156428A (zh) * | 2016-07-06 | 2016-11-23 | 大连海事大学 | 一种基于样本分析的船舶交通流模拟生成方法 |
KR20220057164A (ko) * | 2020-10-29 | 2022-05-09 | 주식회사 에이스이앤티 | 설계데이터를 이용한 작업순서 및 작업배분 시뮬레이션 방법 |
Non-Patent Citations (1)
Title |
---|
考虑船舶航速的船队规划非线性模型;杨秋平;谢新连;裴光石;;华南理工大学学报(自然科学版)(10);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116050720A (zh) | 2023-05-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110083983B (zh) | 一种船舶分段航速优化方法和智能管理系统 | |
CN111651838B (zh) | 船舶舱容测量方法、系统以及终端 | |
CN106202265A (zh) | 流域大尺度复杂河网海量精细规则网格的绘编方法 | |
CN113902172A (zh) | 一种污水处理方法、系统、装置及介质 | |
CN104153332B (zh) | 一种基于水面漂浮物迁移轨迹的清污工程布置方法 | |
Xia et al. | Quadtree-based domain decomposition for parallel map-matching on GPS data | |
CN116912459B (zh) | 一种变网格多尺度混合式资料同化方法 | |
CN116050720B (zh) | 模拟仿真作业过程策略的航道疏浚船舶配比系统及方法 | |
CN101573729A (zh) | 用于多级光线追踪的方法和装置 | |
Tassi et al. | GAIA-a unified framework for sediment transport and bed evolution in rivers, coastal seas and transitional waters in the TELEMAC-MASCARET modelling system | |
CN110096828A (zh) | 基于bim技术的楼层框架梁上部通长筋或架立筋的建模及施工方法 | |
Wijnberg et al. | Beach-dune modelling in support of Building with Nature for an integrated spatial design of urbanized sandy shores | |
CN102912760B (zh) | 水库下游河道最低通航水位生成方法 | |
CN110659342A (zh) | 基于信息提取特征的规则库自动制图综合处理方法及系统 | |
CN113887010A (zh) | 基于实景建模与bim的水利水电场内交通规划组织方法 | |
CN116776464B (zh) | 一种特定航线江海直达集装箱船型线谱系生成方法及系统 | |
CN102799712A (zh) | 一种混凝土结构三维模型的快速配筋方法 | |
Bai et al. | Global time optimization method for dredging construction cycles of trailing suction hopper dredger based on Grey System Model | |
CN116432552A (zh) | 一种水库减淤调度方法、装置、设备及介质 | |
CN112288200A (zh) | 内陆河流溢油应急处置方案确定方法及装置 | |
Gladkov et al. | Hydromorphological substantiation of channel stability of navigable rivers in engineering water transport regulation of river runoff | |
CN113469405B (zh) | 基于子图分割的跨城市尾气迁移预测方法、系统及介质 | |
CN114482174A (zh) | 一种疏浚工程水下精准抛填施工方法 | |
KR100929927B1 (ko) | 선박 화물창 벽면주행 이동로봇의 경로 계획방법 | |
CN116341283B (zh) | 一种基于深水淤泥采样的水下地形仿真方法和系统 |
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 |