CN109146179A - 沿海港口船舶作业条件监测预警方法 - Google Patents
沿海港口船舶作业条件监测预警方法 Download PDFInfo
- Publication number
- CN109146179A CN109146179A CN201810968239.6A CN201810968239A CN109146179A CN 109146179 A CN109146179 A CN 109146179A CN 201810968239 A CN201810968239 A CN 201810968239A CN 109146179 A CN109146179 A CN 109146179A
- Authority
- CN
- China
- Prior art keywords
- wind
- data
- ship
- wave
- force
- 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.)
- Granted
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 23
- 230000004044 response Effects 0.000 claims abstract description 19
- 238000005259 measurement Methods 0.000 claims abstract description 14
- 238000013178 mathematical model Methods 0.000 claims abstract description 12
- 238000004422 calculation algorithm Methods 0.000 claims description 21
- 230000002068 genetic effect Effects 0.000 claims description 13
- 238000007667 floating Methods 0.000 claims description 12
- 238000011068 loading method Methods 0.000 claims description 11
- 230000033001 locomotion Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 7
- 238000004873 anchoring Methods 0.000 claims description 5
- 230000007613 environmental effect Effects 0.000 claims description 5
- 230000009471 action Effects 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 239000003651 drinking water Substances 0.000 claims description 3
- 235000020188 drinking water Nutrition 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 231100001261 hazardous Toxicity 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 238000005096 rolling process Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 206010008190 Cerebrovascular accident Diseases 0.000 claims description 2
- 208000006011 Stroke Diseases 0.000 claims description 2
- 210000000349 chromosome Anatomy 0.000 description 5
- 108090000623 proteins and genes Proteins 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 210000004027 cell Anatomy 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000010353 genetic engineering Methods 0.000 description 2
- 239000003921 oil Substances 0.000 description 2
- 239000004677 Nylon Substances 0.000 description 1
- 239000004952 Polyamide Substances 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009395 breeding Methods 0.000 description 1
- 230000001488 breeding effect Effects 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013277 forecasting method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003278 mimic effect Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 229920001778 nylon Polymers 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 229920002647 polyamide Polymers 0.000 description 1
- 238000010186 staining Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 238000012360 testing method Methods 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"
-
- G06Q50/40—
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
本发明提供了沿海港口船舶作业条件监测预警方法，具体包括如下步骤：通过历史资料搜集以及现场实测数据，建立风资料数据库，得到现场监测的风数据以及预报的风数据；建立波浪资料数据库，得到现场波浪监测数据、以及预测的外海波浪数据、预测的泊位处波浪数据；建立潮流资料数据库，首先建立制定港区的潮流数学模型，给出码头泊位的潮流预测数据；建立船舶动力响应资料数据库，针对港区内的特定船型建立船舶运动量响应数学模型，基于上述步骤中预报的风浪流条件，预测对应的船舶运动量、系缆力及护舷撞击力数值；根据PIANC规范和相关标准，评估系泊船舶的运动量、系缆力和护舷撞击力，给出对应的预警等级。
Description
技术领域
港口航道及近海工程技术领域，尤其是涉及沿海港口船舶作业综合条件预报预警技术。
背景技术
船舶作业条件作为海港码头运营中的重要参考因素，它会影响船舶靠泊操作、码头装卸调度等，进而影响码头运营的经济效益，因此准确预测码头泊位处的风浪、潮流以及船舶动力响应对保障码头系泊安全，合理进行生产调度等有着非常重要的经济效益和安全意义。
经调研，我国沿海港口年吞吐量已超过120亿吨，原油吞吐量达4.2亿吨，沿海港口建设投资已多年维持在850亿元以上，港口大国地位凸显。调研发现，由于建设条件的限制，港口运营期间，不同位置处的泊位作业条件有较大不同，生产效率不一，运营安全在某些泊位还存在较大隐患。例如，一些中小港口的离岸泊位和口门附近的大型泊位，因掩护条件相对较差，受外海涌浪影响较大。
目前，海洋系统的海浪预报主要针对某一海区，尚未直接服务于具体的港口和码头，且波浪预报只有波高，没有波周期，对港口生产作业缺乏针对性和有效性。迫切需要建立一套集波浪预报、潮位潮流预报、船舶系靠泊于一体的预报预警系统，为港口运营和生产提供技术支撑。
发明内容
有鉴于此，本发明旨在提出沿海港口船舶作业条件监测预警方法，对港内泊位处未来3天的作业条件进行预测，并针对系泊船舶在外界环境下的动力响应进行预测，根据对应的运动量、缆绳拉力与规范比较，给出船舶系泊条件下的作业安全等级，码头运营人员可根据以上的预测信息，更为合理地安排船舶靠泊操作以及码头装卸作业。
为达到上述目的，本发明的技术方案是这样实现的：
沿海港口船舶作业条件监测预警方法，具体包括如下步骤：
(1)通过历史资料搜集以及现场实测数据，建立风资料数据库，得到现场监测的风数据以及预报的风数据；
(2)建立波浪资料数据库，得到现场波浪监测数据、以及预测的外海波浪数据、预测的泊位处波浪数据；与现场实测数据进行对比，在原有算法的基础上进行自学习；
(3)建立潮流资料数据库，首先建立制定港区的潮流数学模型，给出码头泊位的潮流预测数据；与现场实测数据进行对比，在原有算法的基础上进行自学习；
(4)建立船舶动力响应资料数据库，针对港区内的特定船型建立船舶运动量响应数学模型，基于上述步骤(1)(2)(3)中预报的风浪流条件，预测对应的船舶运动量、系缆力及护舷撞击力数值；与现场实测数据进行对比，在原有算法的基础上进行自学习；
(5)根据PIANC规范和相关标准，评估系泊船舶的运动量、系缆力和护舷撞击力，给出对应的预警等级，其中包括安全作业、基本满足作业、不可作业但可系泊以及不可系泊需逃离码头等信息。
进一步的，所述步骤(1)中，风资料数据库包括常规风以及台风等灾害条件下的数据。
进一步的，所述步骤(2)中，根据风资料数据库中的数据，运用遗传算法得到外海波浪数据库，然后通过抛物型缓坡方程推算工程区的设计波浪条件及小范围计算波浪边界条件；进而采用丹麦DHI的Mike21计算软件中的BW波浪计算模块得出码头泊位处波浪预测数据。
进一步的，所述步骤(3)中，潮流数学模型采用大小模型嵌套的方法进行计算，选用连续方程、X方向动量方程、Y向动量方程等控制方程计算模型中各个网格单元的流速流向，给出码头泊位的潮流预测数据。
进一步的，所述步骤(4)中，所述的船舶运动量响应数学模型首先基于势流理论对船舶的水动力参数进行计算，包括附加质量、附加阻尼、波浪荷载、单位幅值响应算子(RAO)等参数，然后考虑系泊系统和碰撞护舷的恢复刚度，建立系泊船舶在风、浪、流等环境载荷和系泊载荷的共同作用下的计算模型；在时域内，考虑系泊系统之后的平衡方程如公式(1)所示：
其中，Mij为浮体质量矩阵；Aij为浮体附加质量矩阵；Xj为浮体的位移向量(包括纵移、横移、升沉、横摇、纵摇和艏摇)；Cij为浮体恢复力系数矩阵；为系统的时延函数；Fm为系泊力，n为锚泊线的数量；Fw为波浪力，Fc为流载荷；Fwind为风载荷；
Fm缆绳拉力根据输入的缆绳拉伸力与变形曲线计算，与缆绳的相对伸长(ΔL/L)和缆绳的最小破断力BL有关；
如果带缆桩的坐标为(X1,Y1,Z1)，船舶在初始位置对应的导缆孔坐标(X2,Y2,Z2)，缆绳原长可由下式确定：
如果船舶发生运动后，导缆孔的坐标变为(X'2,Y'2,Z'2)，则缆绳长度应为：
计算出缆绳的相对伸长根据曲线则可得出缆绳对应的拉力；
FW波浪力的计算根据水动力软件HYDROSTAR计算一阶波浪载荷和二阶波浪载荷，并考虑浅水效应，计算二阶差频波浪载荷传递函数(QTF)，得到作用在港内系泊船舶上的波浪作用力，而以往在数值模拟港内系泊船舶时，多未考虑波浪对船舶的激励作用，只考虑了系泊缆的静张力；
Fc根据OCIMF规范的有关规定，船舶受到的纵荡方向流力Fc1、横荡方向流力Fc2和艏摇流力矩Mc分别为：
Fc1＝1/2CxcρwVc 2LBP T；
Fc2＝1/2CycρwVc 2LBP T；
Mc＝1/2Cxycρw Vc 2LBP 2T
其中，系数Cxc、Cyc和Cxyc为流体系数，可根据不同装载、不同船型进行选取；ρw为空气密度；Vc为风速；LBP为船舶的垂线间长；T为吃水；
Fwind根据OCIMF规范的有关规定计算，船舶受到的纵荡方向风力Fwind1、横荡方向风力Fwind2和艏摇风力矩Mwind分别为：
Fwind1＝1/2CxwρwVw 2AT；
Fwind2＝1/2CywρwVw 2AL；
Mwind＝1/2Cxywρw Vw 2AL LBP
其中，系数Cxw、Cyw和Cxyw为风力系数，可根据不同装载、不同船型进行选取；ρw为空气密度；Vw为风速；为船舶的横向受风面积AT，AL为船舶的纵向受风面积；LBP为船舶的垂线间长。
相对于现有技术，本发明所述的沿海港口船舶作业条件监测预警方法具有以下优势：
应用本发明方法，在目前仅有外海波浪条件预报的前提下，能够较为科学地针对工程港区的外界环境条件进行分析，为预测船舶的动力响应提供了可靠的依据；同时结合具体的船型，将外界环境条件反馈到船舶的动态响应中，建立港口船舶作业安全方法，降低了船舶系泊作业期间船舶断缆以及碰撞的风险，避免了由于天气和海况信息的掌握不全面、调整指令时效性以及仅凭个人经验进行系泊作业的弊端，使结果更为合理和准确，对实际操作更具有指导意义。
本发明的另一目的在于提出一种沿海港口船舶作业条件监测预警系统，对码头作业条件进行预测预警，保障码头系泊安全，提高码头经济效益和重要的安全意义。
为了达到上述目的，本发明的技术方案是这样实现的：
一种沿海港口船舶作业条件监测预警系统，包括风资料数据库模块、波浪资料数据模块、潮流资料数据库模块以及船舶动力响应资料数据库模块，其中风资料数据库模块包括风数据现场监测模块和风数据预报模块，波浪资料数据库模块包括波浪数据现场监测模块、外海波浪数据预测模块和泊位处波浪数据预测模块，潮流资料数据库模块包括潮流数据现场监测模块、潮流数据预报模块，船舶动力响应资料数据库模块包括运动量、缆力及护舷撞击力现场监测模块以及运动量、缆力及护舷撞击力数据预报模块。
本发明所述的一种沿海港口船舶作业条件监测预警系统与上述的一种沿海港口船舶作业条件监测预警方法的有益效果相同，在此不再赘述。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为本发明实施例所述的沿海港口船舶作业条件预测预警方法的分析示意图。
具体实施方式
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”等仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”等的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中，除非另有说明，“多个”的含义是两个或两个以上。
在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以通过具体情况理解上述术语在本发明中的具体含义。
下面将参考附图并结合实施例来详细说明本发明。
本发明的目的在于克服现有技术的不足，对港内泊位处未来3天的作业条件进行预测，并针对系泊船舶在外界环境下的动力响应进行预测，根据对应的运动量、缆绳拉力与规范比较，给出船舶系泊条件下的作业安全等级。码头运营人员可根据以上的预测信息，更为合理地安排船舶靠泊操作以及码头装卸作业。
如图1所示，本发明的沿海港口船舶作业条件监测预警方法，它包括以下步骤：
(1)通过历史资料搜集以及现场实测数据，建立风资料数据库，其中包括常规风以及台风等灾害条件下的数据；
(2)建立波浪资料数据库，首先根据上述的风数据，运用遗传算法得到外海波浪数据库，然后通过抛物型缓坡方程推算工程区的设计波浪条件及小范围计算波浪边界条件；进而采用丹麦DHI的Mike21计算软件中的BW波浪计算模块得出码头泊位处波浪预测数据。与现场实测数据进行对比，在原有算法的基础上进行自学习，提高预测精度；
其中，遗传算法如下：遗传算法是一种全局优化算法，它借用了生物遗传学的观点，通过自然选择、遗传、变异等作用机制，实现各个个体的适应性提高。遗传算法类似于自然进化，通过作用于染色体上的基因寻找好的染色体来求解问题。与自然界相似，遗传算法对求解问题的本身一无所知，它所需要的仅是对算法所产生的每个染色体进行评价，并基于适应值来选择染色体，使适应性好的染色体有更多的繁殖机会。在遗传算法中，通过随机方式产生若干个要求解问题的数字编码，即染色体，形成初始群体：通过适应度函数给每个个体一数值评价，淘汰低适应度的个体，选择高适应度的个体参加遗传操作，经过遗传操作后的个体集合形成下一代新的种群。对这个新种群进行下一轮进化。这就是遗传算法的基本原理。
选择即从当前群体中选择适应值高的个体以生成交配池的过程，目的是为了从当前群体中选出优良的个体，使它们有机会作为父代产生后代个体，遗传算法通过选择操作模仿生物中的优胜劣汰操作。适应度越高的个体被遗传到下一代的概率较高。选择的目的是为了避免基因损失，提高全局收敛性。有很多常用的选择算子。
由样本集合Φ＝{(ui,vi)},i＝1,2,…,m,m＜n，其中m为样本个数，u，v为遗传因子，ui，vi为期望值。定义适应度函数为：
具体实现步骤如下：
根据预报区域过去30年风浪数据，生成一个初始群体。
风浪响应时间分析，例如在渤海区域取较长风距500km，对于风可取风速作为是否产生影响的因子，如风速10m/s，则有效影响时间为50000s，为13小时；对于波浪取波长周期比为影响因子，如3s的波浪，速度约为5m/s，则影响时间为26小时。风浪的成长需要一定的过程，综合考虑取一个较长的时间3天，即72小时。
相似基因获取技术，对于渤海的风，有4个主要因子：风速的大小及过程、风向及过程、风速的空间分布、风向的空间分布。在渤海中取n个点，将72小时分隔为m个均匀时间段，预测期内的数据为Ui,j、Vi,j，i＝1,…n；j＝1,…m。在历史数据中相同的位置，相同的时间间隔，任意的时间起点t，得到ui,j、vi,j，i＝1,…n；j＝1,…m。
构建基因符合程度参数kt i,j＝((ui,j-Ui,j)^2+(vi,j-Vi,j)^2)/(Ui,j^2+Vi,j^2)，kt为时间和空间上的总和，在接近30年的数据中取kt的最小值，就是基因差异最小的遗传样本。
风浪预报技术，风与浪之间存在复杂的相关关系，与地形潮位边界等因素有关。应用实测资料进行相关分析是有效且高效的方法。在风浪预报中首先将历史资料作为实测历史数据进行处理，由此可以建立起风与浪的关系。对72小时的风速大小建立不同权重，离当前越近权重越大，得到代表风速；同样对风的方向进行处理，得到代表风向。权重系数可取三角形分布。
(3)建立潮流资料数据库，首先建立工程港区的潮流数学模型，潮流数学模型采用大小模型嵌套的方法进行计算，选用连续方程、X方向动量方程、Y向动量方程等控制方程计算模型中各个网格单元的流速流向，给出码头泊位的潮流预测数据；与现场实测数据进行对比，在原有算法的基础上进行自学习，提高预测精度；
(4)建立船舶动力响应资料数据库，针对港区内的特定船型建模分析，基于上述预报的风浪流条件，预测对应的船舶运动量、系缆力及护舷撞击力数值；与现场实测数据进行对比，在原有算法的基础上进行自学习，提高预测精度；
船舶运动量响应数学模型：首先基于势流理论对船舶的水动力参数进行计算，包括附加质量、附加阻尼、波浪荷载、单位幅值响应算子(RAO)等参数，然后考虑系泊系统和碰撞护舷的恢复刚度，建立系泊船舶在风、浪、流等环境载荷和系泊载荷的共同作用下的计算模型。在时域内，考虑系泊系统之后的平衡方程如公式(1)所示：
其中，Mij为浮体质量矩阵；Aij为浮体附加质量矩阵；Xj为浮体的位移向量(包括纵移、横移、升沉、横摇、纵摇和艏摇)；Cij为浮体恢复力系数矩阵；为系统的时延函数；Fm为缆绳拉力，n为锚泊线的数量；Fw为波浪力，Fc为流载荷；Fwind为风载荷。其中系泊力Fm、波浪力、流载荷Fc和风载荷Fwind。
Fm缆绳拉力根据输入的缆绳拉伸力与变形曲线计算，与缆绳的相对伸长(ΔL/L)和缆绳的最小破断力BL(不同缆绳会有对应的BL值)有关。
如果带缆桩的坐标为(X1,Y1,Z1)，船舶在初始位置对应的导缆孔坐标(X2,Y2,Z2)，缆绳原长可由下式确定：
如果船舶发生运动后，导缆孔的坐标变为(X'2,Y2',Z'2)，则缆绳长度应为：
计算出缆绳的相对伸长根据曲线则可得出缆绳对应的拉力。
FW波浪力的计算根据水动力软件HYDROSTAR计算一阶波浪载荷和二阶波浪载荷，并考虑浅水效应，计算二阶差频波浪载荷传递函数(QTF)，得到作用在港内系泊船舶上的波浪作用力，而以往在数值模拟港内系泊船舶时，多未考虑波浪对船舶的激励作用，只考虑了系泊缆的静张力。
Fc根据OCIMF规范的有关规定，船舶受到的纵荡方向流力Fc1、横荡方向流力Fc2和艏摇流力矩Mc分别为：
Fc1＝1/2CxcρwVc 2LBP T；
Fc2＝1/2CycρwVc 2LBP T；
Mc＝1/2Cxycρw Vc 2LBP 2T
其中，系数Cxc、Cyc和Cxyc为流体系数，可根据不同装载、不同船型进行选取；ρw为空气密度；Vc为风速；LBP为船舶的垂线间长；T为吃水。
Fwind根据OCIMF规范(Oil Companies International Marine Forum)的有关规定计算，船舶受到的纵荡方向风力Fwind1、横荡方向风力Fwind2和艏摇风力矩Mwind分别为：
Fwind1＝1/2CxwρwVw 2AT；
Fwind2＝1/2CywρwVw 2AL；
Mwind＝1/2Cxywρw Vw 2AL LBP
其中，系数Cxw、Cyw和Cxyw为风力系数，可根据不同装载、不同船型进行选取；ρw为空气密度；Vw为风速；为船舶的横向受风面积AT，AL为船舶的纵向受风面积；LBP为船舶的垂线间长。
(5)根据PIANC规范和相关标准(见下)，评估系泊船舶的运动量、系缆力和护舷撞击力，给出对应的预警等级，其中包括安全作业、基本满足作业、不可作业但可系泊以及不可系泊需逃离码头等信息。
(6)进行系统架构搭建以及系统单元测试，通过系统集成实现系统的各项功能，达到预测预警的目的。
其中，相关标准如下表所示：
(1)运动量标准：
表1国际航运协会不同船舶安全作业推荐允许运动范围标准(PIANC，1995)
(2)系缆力标准：
根据我国《港口工程荷载规范》(JTS 144-1-2010)10.2.4条规定：“系缆力的标准值不应大于缆绳的破断力。缆绳破断力应按缆绳材质和规格确定。当缺乏资料时，可按附录G确定”。
根据石油公司国际海事论坛(OCIMF)“Mooring Equipment Guidelines(2008)”的规定：“对于钢缆(steel wire)，其缆绳所受拉力不应大于其最小破断力(MBL)的55％；对于合成纤维缆(synthetic rope)，其缆绳所受拉力不应大于其最小破断力的50％；对于尼龙缆(polyamide rope)，其缆绳所受拉力不应大于其最小破断力的45％”。
(3)撞击力标准：
对于护舷的撞击力和撞击能量，当计算所得到的撞击力和撞击能量超过护舷的设计撞击力和撞击能量时，则认为护舷型号不满足要求。
其中运动量数值为表中60％则为作业条件良好，运动量数值接近表中数值为基本可满足作业，运动量数值大于表中数值为不可作业但可系泊。运动量超过表中数值，且系缆力大于55％的最小破断力，则船舶需要逃离码头。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (6)
1.沿海港口船舶作业条件监测预警方法，其特征在于：具体包括如下步骤：
(1)建立风资料数据库，得到现场监测的风数据以及预报的风数据，其中预报的风数据来源于GFS数值预报模式；
(2)基于ECMWF和搜集的工程实测数据，建立波浪资料数据库，联合预报的风数据，通过遗传算法得到预测的外海波浪数据、预测的泊位处波浪数据；与现场实测数据进行对比，在原有算法的基础上进行自学习；
(3)建立潮流资料数据库，首先建立制定港区的潮流数学模型，给出码头泊位的潮流预测数据；与现场实测数据进行对比，在原有算法的基础上进行自学习；
(4)建立船舶动力响应资料数据库，针对港区内的特定船型建立船舶运动量响应数学模型，基于上述步骤(1)(2)(3)中预报的风浪流条件，预测对应的船舶运动量、系缆力及护舷撞击力数值；与现场实测数据进行对比，在原有算法的基础上进行自学习；
(5)根据PIANC规范和相关标准，评估系泊船舶的运动量、系缆力和护舷撞击力，给出对应的预警等级，其中包括安全作业、基本满足作业、不可作业但可系泊以及不可系泊需逃离码头等信息。
2.根据权利要求1所述的沿海港口船舶作业条件监测预警方法，其特征在于：所述步骤(1)中，风资料数据库包括常规风以及台风等灾害条件下的数据。
3.根据权利要求1所述的沿海港口船舶作业条件监测预警方法，其特征在于：所述步骤(2)中，根据风资料数据库中的数据，运用遗传算法得到外海波浪数据库，然后通过抛物型缓坡方程推算工程区的设计波浪条件及小范围计算波浪边界条件；进而采用丹麦DHI的Mike21计算软件中的BW波浪计算模块得出码头泊位处波浪预测数据。
4.根据权利要求1所述的沿海港口船舶作业条件监测预警方法，其特征在于：所述步骤(3)中，潮流数学模型采用大小模型嵌套的方法进行计算，选用连续方程、X方向动量方程、Y向动量方程等控制方程计算模型中各个网格单元的流速流向，给出码头泊位的潮流预测数据。
5.根据权利要求1所述的沿海港口船舶作业条件监测预警方法，其特征在于：所述步骤(4)中，所述的船舶运动量响应数学模型首先基于势流理论对船舶的水动力参数进行计算，包括附加质量、附加阻尼、波浪荷载、单位幅值响应算子(RAO)等参数，然后考虑系泊系统和碰撞护舷的恢复刚度，建立系泊船舶在风、浪、流等环境载荷和系泊载荷的共同作用下的计算模型，该模型的运动方程如公式(1)所示：
其中，Mij为浮体质量矩阵；Aij为浮体附加质量矩阵；Xj为浮体的位移向量(包括纵移、横移、升沉、横摇、纵摇和艏摇)；Cij为浮体恢复力系数矩阵；为系统的时延函数；Fm为缆绳拉力，n为锚泊线的数量；Fw为波浪力，Fc为流载荷；Fwind为风载荷；其中系泊力Fm、波浪力、流载荷Fc和风载荷Fwind；
Fm缆绳拉力根据输入的缆绳拉伸力与变形曲线计算，与缆绳的相对伸长(ΔL/L)和缆绳的最小破断力BL有关；
如果带缆桩的坐标为(X1,Y1,Z1)，船舶在初始位置对应的导缆孔坐标(X2,Y2,Z2)，缆绳原长可由下式确定：
如果船舶发生运动后，导缆孔的坐标变为(X'2,Y2',Z'2)，则缆绳长度应为：
计算出缆绳的相对伸长根据曲线则可得出缆绳对应的拉力；
FW波浪力的计算根据水动力软件HYDROSTAR计算一阶波浪载荷和二阶波浪载荷，并考虑浅水效应，计算二阶差频波浪载荷传递函数(QTF)，得到作用在港内系泊船舶上的波浪作用力，而以往在数值模拟港内系泊船舶时，多未考虑波浪对船舶的激励作用，只考虑了系泊缆的静张力；
Fc根据OCIMF规范的有关规定，船舶受到的纵荡方向流力Fc1、横荡方向流力Fc2和艏摇流力矩Mc分别为：
Fc1＝1/2CxcρwVc 2LBPT；
Fc2＝1/2CycρwVc 2LBPT；
Mc＝1/2CxycρwVc 2LBP 2T
其中，系数Cxc、Cyc和Cxyc为流体系数，可根据不同装载、不同船型进行选取；ρw为空气密度；Vc为风速；LBP为船舶的垂线间长；T为吃水；
Fwind根据OCIMF规范的有关规定计算，船舶受到的纵荡方向风力Fwind1、横荡方向风力Fwind2和艏摇风力矩Mwind分别为：
Fwind1＝1/2CxwρwVw 2AT；
Fwind2＝1/2CywρwVw 2AL；
Mwind＝1/2CxywρwVw 2ALLBP
其中，系数Cxw、Cyw和Cxyw为风力系数，可根据不同装载、不同船型进行选取；ρw为空气密度；Vw为风速；为船舶的横向受风面积AT，AL为船舶的纵向受风面积；LBP为船舶的垂线间长。
6.一种沿海港口船舶作业条件监测预警系统，其特征在于：包括风资料数据库模块、波浪资料数据模块、潮流资料数据库模块以及船舶动力响应资料数据库模块，其中风资料数据库模块包括风数据现场监测模块和风数据预报模块，波浪资料数据库模块包括波浪数据现场监测模块、外海波浪数据预测模块和泊位处波浪数据预测模块，潮流资料数据库模块包括潮流数据现场监测模块、潮流数据预报模块，船舶动力响应资料数据库模块包括运动量、缆力及护舷撞击力现场监测模块以及运动量、缆力及护舷撞击力数据预报模块。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810968239.6A CN109146179B (zh) | 2018-08-23 | 2018-08-23 | 沿海港口船舶作业条件监测预警方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810968239.6A CN109146179B (zh) | 2018-08-23 | 2018-08-23 | 沿海港口船舶作业条件监测预警方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109146179A true CN109146179A (zh) | 2019-01-04 |
| CN109146179B CN109146179B (zh) | 2021-07-23 |
Family
ID=64791239
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810968239.6A Active CN109146179B (zh) | 2018-08-23 | 2018-08-23 | 沿海港口船舶作业条件监测预警方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109146179B (zh) |
Cited By (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108399279A (zh) * | 2018-01-25 | 2018-08-14 | 交通运输部天津水运工程科学研究所 | 波浪影响下码头系泊船舶损失作业天数分析与控制方法 |
| CN110175696A (zh) * | 2019-04-22 | 2019-08-27 | 中国水产科学研究院渔业工程研究所 | 一种基于多元回归的渔港船舶进出港动态预测方法及系统 |
| CN110443426A (zh) * | 2019-08-12 | 2019-11-12 | 上海海洋中心气象台 | 基于浅滩效应和陡度联合预测船舶风险的方法及装置 |
| CN110532685A (zh) * | 2019-08-29 | 2019-12-03 | 山东交通学院 | 浮式结构物摇荡运动响应预报方法 |
| CN110852619A (zh) * | 2019-11-12 | 2020-02-28 | 山东交通学院 | 一种在船舶性能的评估中对海浪阻力的修正方法 |
| CN111125881A (zh) * | 2019-11-21 | 2020-05-08 | 天津海运职业学院 | 一种船舶靠离泊过程中系缆作业仿真的方法 |
| CN111311964A (zh) * | 2020-03-04 | 2020-06-19 | 上海楼语智能科技有限公司 | 一种智能安全监控平台 |
| CN111507586A (zh) * | 2020-02-19 | 2020-08-07 | 中交疏浚技术装备国家工程研究中心有限公司 | 通航水域水上施工安全风险预警分级方法及系统 |
| CN111532398A (zh) * | 2020-04-26 | 2020-08-14 | 杭州智海人工智能有限公司 | 一种关于无人船的模拟与观测算法 |
| CN111797476A (zh) * | 2020-07-06 | 2020-10-20 | 华电重工股份有限公司 | 浮吊起吊天气窗口判定方法 |
| CN112050740A (zh) * | 2020-09-03 | 2020-12-08 | 交通运输部天津水运工程科学研究所 | 系泊安全监测及预警系统 |
| CN112484957A (zh) * | 2020-08-06 | 2021-03-12 | 重庆交通大学 | 滑坡涌浪对斜坡式码头-趸船-船舶模型系统作用力的试验装置及试验方法 |
| CN112700125A (zh) * | 2020-12-29 | 2021-04-23 | 宁波市气象服务中心 | 一种港口气象灾害风险及气象服务效益评估的方法 |
| CN113032714A (zh) * | 2021-03-05 | 2021-06-25 | 大连中远海运重工有限公司 | 液化天然气运输船船对船靠泊用护舷的选型计算方法 |
| CN113553785A (zh) * | 2021-07-14 | 2021-10-26 | 海博泰科技(青岛)有限公司 | 一种开敞式码头及港池波浪预报方法 |
| CN113553775A (zh) * | 2021-09-01 | 2021-10-26 | 交通运输部天津水运工程科学研究所 | 一种基于机器学习的港内码头系泊船舶运动量预报方法 |
| CN113704965A (zh) * | 2021-07-22 | 2021-11-26 | 中国海洋大学 | 浮托安装荷载转移过程动态模拟分析方法 |
| CN113761668A (zh) * | 2021-09-28 | 2021-12-07 | 江南造船(集团)有限责任公司 | 船舶系泊控制方法、系统、介质及电子设备 |
| CN113779698A (zh) * | 2021-09-14 | 2021-12-10 | 浙江数智交院科技股份有限公司 | 一种水流作用下的船舶系泊系统简化设计方法 |
| CN114162266A (zh) * | 2021-12-08 | 2022-03-11 | 交通运输部天津水运工程科学研究所 | 一种基于现场实时反馈的岸基智能系泊系统及方法 |
| CN115410419A (zh) * | 2022-08-23 | 2022-11-29 | 交通运输部天津水运工程科学研究所 | 一种船舶系泊预警方法、系统、电子设备及存储介质 |
| CN116341424A (zh) * | 2023-05-30 | 2023-06-27 | 交通运输部天津水运工程科学研究所 | 一种作用于船舶上的水流力的综合计算方法 |
| CN116910481A (zh) * | 2023-07-27 | 2023-10-20 | 中国舰船研究设计中心 | 一种基于遗传算法的船舶任务系统装载弹量优化方法 |
| CN117786604A (zh) * | 2024-02-23 | 2024-03-29 | 交通运输部水运科学研究所 | 一种基于港区灾害数据的融合安全韧性度量分析系统 |
| CN117786604B (zh) * | 2024-02-23 | 2024-04-26 | 交通运输部水运科学研究所 | 一种基于港区灾害数据的融合安全韧性度量分析系统 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101329696A (zh) * | 2007-11-16 | 2008-12-24 | 武汉理工大学 | 进港航道通过能力计算模型及预测仿真方法 |
| CN107942860A (zh) * | 2017-11-27 | 2018-04-20 | 大连交通大学 | 开敞式码头系泊缆力预警方法及系统 |
-
2018
- 2018-08-23 CN CN201810968239.6A patent/CN109146179B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101329696A (zh) * | 2007-11-16 | 2008-12-24 | 武汉理工大学 | 进港航道通过能力计算模型及预测仿真方法 |
| CN107942860A (zh) * | 2017-11-27 | 2018-04-20 | 大连交通大学 | 开敞式码头系泊缆力预警方法及系统 |
Non-Patent Citations (3)
| Title |
|---|
| 岳俊涛: "开敞式码头安全保障与预警系统网络查询平台设计", 《计算机与现代化》 * |
| 张勇军等: "《热轧电气自动化与计算机控制技术》", 31 October 2016, 冶金工业出版社 * |
| 邹志利等: "风浪流作用下系泊船系缆力和碰撞力的数值模拟", 《中国海洋平台》 * |
Cited By (35)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108399279B (zh) * | 2018-01-25 | 2021-05-25 | 交通运输部天津水运工程科学研究所 | 波浪影响下码头系泊船舶损失作业天数分析与控制方法 |
| CN108399279A (zh) * | 2018-01-25 | 2018-08-14 | 交通运输部天津水运工程科学研究所 | 波浪影响下码头系泊船舶损失作业天数分析与控制方法 |
| CN110175696A (zh) * | 2019-04-22 | 2019-08-27 | 中国水产科学研究院渔业工程研究所 | 一种基于多元回归的渔港船舶进出港动态预测方法及系统 |
| CN110175696B (zh) * | 2019-04-22 | 2020-03-27 | 中国水产科学研究院渔业工程研究所 | 一种基于多元回归的渔港船舶进出港动态预测方法及系统 |
| CN110443426A (zh) * | 2019-08-12 | 2019-11-12 | 上海海洋中心气象台 | 基于浅滩效应和陡度联合预测船舶风险的方法及装置 |
| CN110532685A (zh) * | 2019-08-29 | 2019-12-03 | 山东交通学院 | 浮式结构物摇荡运动响应预报方法 |
| CN110852619A (zh) * | 2019-11-12 | 2020-02-28 | 山东交通学院 | 一种在船舶性能的评估中对海浪阻力的修正方法 |
| CN111125881A (zh) * | 2019-11-21 | 2020-05-08 | 天津海运职业学院 | 一种船舶靠离泊过程中系缆作业仿真的方法 |
| CN111125881B (zh) * | 2019-11-21 | 2023-04-21 | 天津海运职业学院 | 一种船舶靠离泊过程中系缆作业仿真的方法 |
| CN111507586A (zh) * | 2020-02-19 | 2020-08-07 | 中交疏浚技术装备国家工程研究中心有限公司 | 通航水域水上施工安全风险预警分级方法及系统 |
| CN111507586B (zh) * | 2020-02-19 | 2024-03-22 | 中交疏浚技术装备国家工程研究中心有限公司 | 通航水域水上施工安全风险预警分级方法及系统 |
| CN111311964A (zh) * | 2020-03-04 | 2020-06-19 | 上海楼语智能科技有限公司 | 一种智能安全监控平台 |
| CN111532398A (zh) * | 2020-04-26 | 2020-08-14 | 杭州智海人工智能有限公司 | 一种关于无人船的模拟与观测算法 |
| CN111797476A (zh) * | 2020-07-06 | 2020-10-20 | 华电重工股份有限公司 | 浮吊起吊天气窗口判定方法 |
| CN112484957A (zh) * | 2020-08-06 | 2021-03-12 | 重庆交通大学 | 滑坡涌浪对斜坡式码头-趸船-船舶模型系统作用力的试验装置及试验方法 |
| CN112050740A (zh) * | 2020-09-03 | 2020-12-08 | 交通运输部天津水运工程科学研究所 | 系泊安全监测及预警系统 |
| CN112700125A (zh) * | 2020-12-29 | 2021-04-23 | 宁波市气象服务中心 | 一种港口气象灾害风险及气象服务效益评估的方法 |
| CN113032714A (zh) * | 2021-03-05 | 2021-06-25 | 大连中远海运重工有限公司 | 液化天然气运输船船对船靠泊用护舷的选型计算方法 |
| CN113553785A (zh) * | 2021-07-14 | 2021-10-26 | 海博泰科技(青岛)有限公司 | 一种开敞式码头及港池波浪预报方法 |
| CN113553785B (zh) * | 2021-07-14 | 2023-12-26 | 海博泰科技(青岛)有限公司 | 一种开敞式码头及港池波浪预报方法 |
| CN113704965A (zh) * | 2021-07-22 | 2021-11-26 | 中国海洋大学 | 浮托安装荷载转移过程动态模拟分析方法 |
| CN113553775B (zh) * | 2021-09-01 | 2022-03-08 | 交通运输部天津水运工程科学研究所 | 一种基于机器学习的港内码头系泊船舶运动量预报方法 |
| CN113553775A (zh) * | 2021-09-01 | 2021-10-26 | 交通运输部天津水运工程科学研究所 | 一种基于机器学习的港内码头系泊船舶运动量预报方法 |
| CN113779698A (zh) * | 2021-09-14 | 2021-12-10 | 浙江数智交院科技股份有限公司 | 一种水流作用下的船舶系泊系统简化设计方法 |
| CN113761668A (zh) * | 2021-09-28 | 2021-12-07 | 江南造船(集团)有限责任公司 | 船舶系泊控制方法、系统、介质及电子设备 |
| CN114162266A (zh) * | 2021-12-08 | 2022-03-11 | 交通运输部天津水运工程科学研究所 | 一种基于现场实时反馈的岸基智能系泊系统及方法 |
| CN114162266B (zh) * | 2021-12-08 | 2022-11-15 | 交通运输部天津水运工程科学研究所 | 一种基于现场实时反馈的岸基智能系泊系统及方法 |
| CN115410419B (zh) * | 2022-08-23 | 2024-02-02 | 交通运输部天津水运工程科学研究所 | 一种船舶系泊预警方法、系统、电子设备及存储介质 |
| CN115410419A (zh) * | 2022-08-23 | 2022-11-29 | 交通运输部天津水运工程科学研究所 | 一种船舶系泊预警方法、系统、电子设备及存储介质 |
| CN116341424B (zh) * | 2023-05-30 | 2023-08-15 | 交通运输部天津水运工程科学研究所 | 一种作用于船舶上的水流力的综合计算方法 |
| CN116341424A (zh) * | 2023-05-30 | 2023-06-27 | 交通运输部天津水运工程科学研究所 | 一种作用于船舶上的水流力的综合计算方法 |
| CN116910481A (zh) * | 2023-07-27 | 2023-10-20 | 中国舰船研究设计中心 | 一种基于遗传算法的船舶任务系统装载弹量优化方法 |
| CN116910481B (zh) * | 2023-07-27 | 2024-02-02 | 中国舰船研究设计中心 | 一种基于遗传算法的船舶任务系统装载弹量优化方法 |
| CN117786604A (zh) * | 2024-02-23 | 2024-03-29 | 交通运输部水运科学研究所 | 一种基于港区灾害数据的融合安全韧性度量分析系统 |
| CN117786604B (zh) * | 2024-02-23 | 2024-04-26 | 交通运输部水运科学研究所 | 一种基于港区灾害数据的融合安全韧性度量分析系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109146179B (zh) | 2021-07-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109146179A (zh) | 沿海港口船舶作业条件监测预警方法 | |
| AU2020102354A4 (en) | Morning and early warning method for coastal port ship operation conditions | |
| Aryai et al. | Reliability of multi-purpose offshore-facilities: Present status and future direction in Australia | |
| CN102221389A (zh) | 结合统计模型与动力模型的乘潮水位预报方法 | |
| CN105836066A (zh) | 开敞式码头系泊船姿态预测方法及系统 | |
| Clemente et al. | Proof-of-concept study on a wave energy converter based on the roll oscillations of multipurpose offshore floating platforms | |
| Sulaiman et al. | Mooring analysis for very large offshore aquaculture ocean plantation floating structure | |
| CN203745863U (zh) | 沉管隧道管节海上浮运和沉放施工作业监控系统 | |
| CN103698823A (zh) | 沉管隧道管节海上浮运沉放施工气象窗口预报系统与方法 | |
| JP2002149767A (ja) | 港湾における船舶の運用管理システム及びその運用管理方法 | |
| CN113283138A (zh) | 一种基于深度学习的深海养殖平台动力响应分析方法 | |
| Shi | Offshore Wind Turbine-Ice Interactions | |
| Oladokun | Hydrodynamic Model Test and Simulation for the Design of Offshore Aquaculture System for Macro Algae Farming Journal of Aquaculture & Livestock Production | |
| Feng | Autonomous Underwater Vehicle (AUV) | |
| Aps et al. | Bayesian inference for oil spill related Net Environmental Benefit Analysis | |
| Özgül et al. | Two experimental fish aggregating systems (fads) in the Aegean sea: their design and application | |
| Peng et al. | Experimental Investigation on Element Immersing Process of Immersed Tube Tunnel of Hong Kong-Zhuhai-Macao Bridge | |
| Stemland | Assessment of Service Vessel Operability In Exposed Aquaculture-An exploratory approach combining vessel response and discrete-event simulation | |
| Tee et al. | Design and analysis of wave energy buoy integrated with seaweed farming | |
| Gan | Alternative Fuel for Ship Propulsion | |
| Zhang | Application of Image Processing in Ice–Structure Interaction | |
| Jiang et al. | AUV/ROV/HOV Resistance | |
| Jiang et al. | Analysis of Renewable Energy Devices | |
| Chen et al. | AUV/ROV/HOV Propulsion System | |
| Schreier et al. | Committee V. 6: Ocean Space Utilization |
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 |