CN112693578B - 基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法 - Google Patents

基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法 Download PDF

Info

Publication number
CN112693578B
CN112693578B CN202011619590.8A CN202011619590A CN112693578B CN 112693578 B CN112693578 B CN 112693578B CN 202011619590 A CN202011619590 A CN 202011619590A CN 112693578 B CN112693578 B CN 112693578B
Authority
CN
China
Prior art keywords
heave
acceleration
semi
submersible
term
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
Application number
CN202011619590.8A
Other languages
English (en)
Other versions
CN112693578A (zh
Inventor
刘福顺
高树健
纪翔
田哲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ocean University of China
Original Assignee
Ocean University of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Ocean University of China filed Critical Ocean University of China
Priority to CN202011619590.8A priority Critical patent/CN112693578B/zh
Publication of CN112693578A publication Critical patent/CN112693578A/zh
Priority to PCT/CN2021/129514 priority patent/WO2022142746A1/zh
Priority to US17/775,321 priority patent/US11760449B2/en
Application granted granted Critical
Publication of CN112693578B publication Critical patent/CN112693578B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B63—SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63B—SHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B1/00—Hydrodynamic or hydrostatic features of hulls or of hydrofoils
    • B63B1/02—Hydrodynamic or hydrostatic features of hulls or of hydrofoils deriving lift mainly from water displacement
    • B63B1/10—Hydrodynamic or hydrostatic features of hulls or of hydrofoils deriving lift mainly from water displacement with multiple hulls
    • B63B1/107—Semi-submersibles; Small waterline area multiple hull vessels and the like, e.g. SWATH
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B63—SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63B—SHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B35/00—Vessels or similar floating structures specially adapted for specific purposes and not otherwise provided for
    • B63B35/44—Floating buildings, stores, drilling platforms, or workshops, e.g. carrying water-oil separating devices
    • B63B35/4413—Floating drilling platforms, e.g. carrying water-oil separating devices
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B63—SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63B—SHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B79/00—Monitoring properties or operating parameters of vessels in operation
    • B63B79/10—Monitoring properties or operating parameters of vessels in operation using sensors, e.g. pressure sensors, strain gauges or accelerometers
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B63—SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63B—SHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B79/00—Monitoring properties or operating parameters of vessels in operation
    • B63B79/20—Monitoring properties or operating parameters of vessels in operation using models or simulation, e.g. statistical models or stochastic models
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B63—SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63B—SHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B79/00—Monitoring properties or operating parameters of vessels in operation
    • B63B79/30—Monitoring properties or operating parameters of vessels in operation for diagnosing, testing or predicting the integrity or performance of vessels

Landscapes

  • Engineering & Computer Science (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • Ocean & Marine Engineering (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Fluid Mechanics (AREA)
  • Architecture (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Other Liquid Machine Or Engine Such As Wave Power Use (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明涉及一种基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法，在垂荡运动中，基于线性势流理论，忽略附加质量和辐射阻尼的耦合影响，对垂荡加速度进行表征；考虑实测海洋环境的噪声影响、环境缓慢变化引起的低频影响以及加速度传感器自身基线漂移误差引起的影响，引入噪声项、低频变化项、基线漂移误差项，并对噪声项、低频变化项、基线漂移误差项进行统一的普罗尼序列归一化表征；对归一化表征后的垂荡加速度进行漂移项去除，通过去除漂移项后剩余的普罗尼序列建立垂荡加速度与垂荡运动参量之间的关系式，对垂荡运动参量进行估计。相比于传统基于滤波器的方法，本发明通过普罗尼序列进行垂荡运动参量的估计，具有更高的计算精度。

Description

基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法
技术领域
本发明属于半潜式平台的垂荡运动技术领域，尤其涉及一种基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法。
背景技术
在半潜式平台的运行过程中，许多海洋技术的应用都是基于结构的垂荡运动开展的，因此对半潜式平台的垂荡运动进行监测是至关重要的。例如，半潜式平台的垂荡运动是立管系统的重要外力之一，直接关系到立管系统的整体动力响应分析，并严重影响垂直运输系统的稳定性。另一个应用的例子是在海洋结构物的安装过程中，高精度的垂荡运动监测可以有效地提高海上结构物的安装效率。此外，对于半潜式平台来说，恶劣海况下平台的运动和水面的变化会引起平台气隙值的变化，而负气隙可能会对平台造成破坏，甚至是人员伤亡。
随着全球定位系统的发展，目前大多数半潜式平台的垂荡运动的定位和监测都是基于全球定位系统开展的。然而，全球定位系统的采样率较低，一般不超过20Hz，且精度较差。此外，在某些极端恶劣的条件下，可能会导致必要的运动信息的丢失。
虽然在理论上可以通过对垂荡加速度进行积分得到结构的垂荡速度信息，然后再次积分得到结构的垂荡运动参量信息，但是在实际测试时，结构的初始速度和初始位移通常是未知的，这两个未知项会使积分结果产生漂移的问题。同时，现场测试时的加速度传感器由于不可避免地基线误差，也会给结果带来较大的误差。为了从结构的垂荡加速度中估计结构的垂荡运动信息，Richter等人基于惯性测量单元，通过使用自适应的升沉滤波器，提出了三种减小积分误差的相位修正法。通过对每个滤波器进行误差分析，从而导出误差函数，然后将误差函数最小化从而获得每个滤波器的最佳参数。Kch l er等人将惯性测量单元作为独立的运动传感器推导了一种用于垂荡运动参量估计的观测器。在该方法中，将垂荡运动近似为叠加的正弦波，然后通过快速傅里叶变换将精确近似的正弦波与相应的频率一起识别出来，再利用识别的参数建立估计垂荡运动的观测器模型。但是上述方法通过滤波对漂移项进行滤除，不可避免的会造成垂荡运动参量中的信息丢失，从而导致估计结果精度较差。
发明内容
本发明在上述现有方法不足的基础上提供了一种基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法，基于线性势流理论推导了结构垂荡方向上的运动方程，并通过普罗尼序列建立了半潜式海洋平台垂荡加速度与垂荡运动参量之间的关系，避免了传统基于滤波器方法带来的误差问题，具有较高的计算精度和实用性。
为了实现上述目的，本发明提供了一种基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法，包括：
在半潜式海洋平台垂荡运动中，基于线性势流理论，忽略附加质量和辐射阻尼的耦合影响，对半潜式海洋平台的垂荡加速度进行表征，确定垂荡加速度理论值；
考虑半潜式海洋平台垂荡运动实测海洋环境的噪声影响、环境缓慢变化引起的低频影响以及加速度传感器自身基线漂移误差引起的影响，引入噪声项、低频变化项、基线漂移误差项，确定垂荡加速度实测值；
对垂荡加速度实测值中的垂荡加速度理论值项、噪声项、低频变化项、基线漂移误差项进行统一的普罗尼序列归一化表征；
对归一化表征后的垂荡加速度进行漂移项去除，通过去除漂移项后剩余的普罗尼序列建立半潜式海洋平台垂荡加速度与垂荡运动参量之间的关系式，对半潜式海洋平台的垂荡运动参量进行估计。
优选的，在半潜式海洋平台垂荡运动中，基于线性势流理论，忽略附加质量和辐射阻尼的耦合影响，同时考虑在流体中受到的波浪力、回复力和辐射力的影响，将垂荡运动表示为：
Figure GDA0003307425910000032
式中，m代表半潜式海洋平台的质量，
Figure GDA0003307425910000031
代表半潜式海洋平台的垂荡加速度，fw(t)为半潜式海洋平台受到的波浪荷载，fm(t)代表半潜式海洋平台受到的系泊力，fs(t)代表半潜式海洋平台受到的回复力，fr(t)代表半潜式海洋平台受到的辐射力；
其中，回复力fs(t)表示为：
fs(t)＝-czzo(t)＝-ρgAwzo(t) (2)
式中，zo(t)代表半潜式海洋平台的垂向位移；cz是半潜式海洋平台在垂荡方向的回复刚度，与水线面面积Aw，流体密度ρ以及重力加速度g有关；
辐射力fr(t)表示为：
Figure GDA0003307425910000041
式中，
Figure GDA0003307425910000042
代表半潜式海洋平台垂荡方向的速度，m∞和kz是无限频率处的垂荡方向的附加质量和脉冲响应函数；
由式(1)-(3)，将半潜式海洋平台垂荡运动表示为：
Figure GDA0003307425910000043
式中，f0(t)＝fw(t)+fm(t)；
则确定半潜式海洋平台垂荡加速度理论值表示为：
Figure GDA0003307425910000044
理论上，半潜式海洋平台的垂向加速度可以建模为一组谐波的叠加的形式，则由式(5)，将垂荡加速度理论值表征为：
Figure GDA0003307425910000045
式中，Ai、fi和θi分别代表垂向加速度中第i个分量的幅值、频率和相位，
Figure GDA0003307425910000046
和
Figure GDA0003307425910000047
为使用普罗尼序列对半潜式海洋平台垂荡加速度理论值进行拟合时使用的参数。
优选的，考虑半潜式海洋平台垂荡运动实测海洋环境的噪声影响、环境缓慢变化引起的低频影响以及加速度传感器自身基线漂移误差引起的影响，引入噪声项、低频变化项、基线漂移误差项，确定垂荡加速度实测值为
Figure GDA0003307425910000048
式中，n(t)代表噪声项，v(t)代表低频变化项，b代表基线漂移误差项。
优选的，引入普罗尼序列对垂荡加速度实测值中的噪声项、低频变化项、基线漂移误差项分别进行表征，即：
Figure GDA0003307425910000051
式中，
Figure GDA0003307425910000052
其中，An、fn、ξn和θn分别代表噪声项中各个成分的幅值、频率、阻尼和相位；
Figure GDA0003307425910000053
式中，
Figure GDA0003307425910000054
Dv＝-ξv+j2πfv，其中，Av、fv、ξv和θv分别代表低频变化项中各个成分的幅值、频率、阻尼和相位；
b＝EeFt (10)
式中，E和F为基线漂移误差项进行拟合时使用的参数；
由式(6)-(10)，将垂荡加速度实测值中的垂荡加速度理论值项、噪声项、低频变化项、基线漂移误差项进行统一的普罗尼序列表征，得到：
Figure GDA0003307425910000055
进一步对垂荡加速度实测值进行归一化表征：
Figure GDA0003307425910000056
式中，Np＝Ni+Nn+Nv+1,
Figure GDA0003307425910000057
和
Figure GDA0003307425910000058
为使用普罗尼序列对半潜式海洋平台的垂荡加速度进行归一化表征时使用的普罗尼序列参数。
优选的，根据计算出的普罗尼参数
Figure GDA0003307425910000059
确定归一化表征后的垂荡加速度各成分的频率，即：
Figure GDA00033074259100000510
对求解出的频率进行排序，并去除其中的极小频率成分，即漂移项，则得到去除漂移项的归一化表征的垂荡加速度实测值为：
Figure GDA0003307425910000061
式中，
Figure GDA0003307425910000062
和
Figure GDA0003307425910000063
为使用普罗尼序列对去除漂移项的归一化表征的垂荡加速度实测值表征的普罗尼序列参数。
优选的，根据去除漂移项的归一化表征的垂荡加速度实测值，确定垂荡运动响应：
Figure GDA0003307425910000064
即垂荡加速度与垂荡运动参量之间的关系式为：
Figure GDA0003307425910000065
则将半潜式海洋平台的真实垂荡运动参量表征为：
Figure GDA0003307425910000066
与现有技术相比，本发明的优点和积极效果在于：
本发明提供的基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法，基于线性势流理论，忽略了垂荡运动中附加质量和辐射阻尼的耦合影响，推导了半潜式平台垂荡运动参量的运动方程，建立了结构在波浪作用下的垂荡加速度的数学模型。同时，从多方面考虑了环境、设备等对于测试的半潜式海洋平台垂荡加速度的影响，包括复杂海洋环境的噪声影响、潮流缓慢变化引起的低频影响以及加速度传感器自身基线漂移误差引起的影响，使得计算结果将更符合于实际海况测试，具有更高的实际应用价值。同时，通过对半潜式海洋平台垂荡加速度、环境噪声、潮流变化及传感器基线漂移进行统一的普罗尼序列表征，并通过对各成分的频率进行筛选和剔除，从而通过剩余的普罗尼序列建立了垂荡加速度和垂荡运动参量之间的转换关系，实现了半潜式平台垂荡运动参量的预报，具有较高的计算精度和实用性，避免了传统基于滤波器方法带来的误差问题。
附图说明
图1为本发明的基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法整体流程图；
图2为试验布置简图；
图3为使用加速度传感器和光学六自由度仪测试的半潜式平台垂荡加速度和位移时域图，其中，(a)为垂荡加速度时域图，(b)为垂荡运动参量时域图；
图4为使用普罗尼参数对实测垂荡加速度的拟合结果图，其中，(a)为使用普罗尼信号对垂荡加速度的拟合结果，(b)为100到110秒的局部加速度信号的拟合结果；
图5为使用本发明方法重构的半潜式平台垂荡运动参量结果图，其中，(a)为本发明方法重构得到的结构位移与测试位移的对比图，(b)为100到110秒的局部垂荡信号的估计结果。
具体实施方式
以下结合附图对本发明的具体实施方式进行进一步的描述。
本发明提供了一种基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法，如图1所示，具体包括：
(1)在半潜式海洋平台垂荡运动中，基于线性势流理论，忽略附加质量和辐射阻尼的耦合影响，对半潜式海洋平台的垂荡加速度进行表征，确定垂荡加速度理论值。具体为：
在半潜式海洋平台垂荡运动中，基于线性势流理论，忽略附加质量和辐射阻尼的耦合影响，同时考虑在流体中受到的波浪力、回复力和辐射力的影响，将垂荡运动表示为：
Figure GDA0003307425910000081
式中，m代表半潜式海洋平台的质量，
Figure GDA0003307425910000082
代表半潜式海洋平台的垂荡加速度，fw(t)为半潜式海洋平台受到的波浪荷载，fm(t)代表半潜式海洋平台受到的系泊力，fs(t)代表半潜式海洋平台受到的回复力，fr(t)代表半潜式海洋平台受到的辐射力。
其中，回复力fs(t)表示为：
fs(t)＝-czzo(t)＝-ρgAwzo(t) (2)
式中，zo(t)代表半潜式海洋平台的垂向位移；cz是半潜式海洋平台在垂荡方向的回复刚度，与水线面面积Aw，流体密度ρ以及重力加速度g有关。
辐射力fr(t)表示为：
Figure GDA0003307425910000091
式中，
Figure GDA0003307425910000092
代表半潜式海洋平台垂荡方向的速度，m∞和kz是无限频率处的垂荡方向的附加质量和脉冲响应函数。
由式(1)-(3)，将半潜式海洋平台垂荡运动表示为：
Figure GDA0003307425910000093
式中，f0(t)＝fw(t)+fm(t)；
则确定半潜式海洋平台垂荡加速度理论值表示为：
Figure GDA0003307425910000094
理论上，半潜式海洋平台的垂向加速度可以建模为一组谐波的叠加的形式，则由式(5)，将垂荡加速度理论值表征为：
Figure GDA0003307425910000095
式中，Ai、fi和θi分别代表垂向加速度中第i个分量的幅值、频率和相位，
Figure GDA0003307425910000096
和
Figure GDA0003307425910000097
为使用普罗尼序列对半潜式海洋平台垂荡加速度理论值进行拟合时使用的参数。
因此，对于波浪作用下的半潜式海洋平台，本实施例中同时考虑了在流体中受到的波浪力、回复力和辐射力的影响，并基于线性势流理论，忽略垂荡运动中附加质量和辐射阻尼的耦合影响，建立了半潜式平台垂荡加速度的理论模型。
(2)在实际的半潜式海洋平台运行环境中，除了结构的运动外，由于复杂的海洋环境及机械运转会产生大量的噪声干扰，实测的垂荡加速度中还包含有大量的环境噪声和流体缓慢变化引起的效应。此外，受加速度传感器自身的限制，在测试中不可避免地会带来加速度传感器基线漂移带来的误差。因此，进一步考虑半潜式海洋平台垂荡运动实测海洋环境的噪声影响、环境缓慢变化引起的低频影响以及加速度传感器自身基线漂移误差引起的影响，引入噪声项、低频变化项、基线漂移误差项，确定垂荡加速度实测值为：
Figure GDA0003307425910000101
式中，n(t)代表噪声项，v(t)代表低频变化项，b代表基线漂移误差项。
因此，本实施例中对于波浪作用下的半潜式海洋平台的垂荡加速度响应，除了在波浪运动下的结构运动外，还考虑了由于海洋环境及机械运转带来的噪声项，以及由于潮差等效性引起的缓慢变化效应，同时也考虑了使用加速度传感器进行测试时不可避免地引入的基线漂移问题，相比于仅使用谐波的叠加对半潜式平台垂荡运动参量的表示，对垂荡加速度的表示方式更符合实际海洋环境中结构的运行状态。
(3)对垂荡加速度实测值中的垂荡加速度理论值项、噪声项、低频变化项、基线漂移误差项进行统一的普罗尼序列归一化表征，具体为：
引入普罗尼序列对式(7)垂荡加速度实测值中的噪声项、低频变化项、基线漂移误差项分别进行表征，即：
Figure GDA0003307425910000102
式中，
Figure GDA0003307425910000103
其中，An、fn、ξn和θn分别代表噪声项中各个成分的幅值、频率、阻尼和相位；
Figure GDA0003307425910000104
式中，
Figure GDA0003307425910000111
Dv＝-ξv+j2πfv，其中，Av、fv、ξv和θv分别代表低频变化项中各个成分的幅值、频率、阻尼和相位；
b＝EeFt (10)
式中，E和F为基线漂移误差项进行拟合时使用的参数；
由式(6)-(10)，将垂荡加速度实测值中的垂荡加速度理论值项、噪声项、低频变化项、基线漂移误差项进行统一的普罗尼序列表征，得到：
Figure GDA0003307425910000112
进一步对垂荡加速度实测值进行归一化表征：
Figure GDA0003307425910000113
式中，Np＝Ni+Nn+Nv+1,
Figure GDA0003307425910000114
和
Figure GDA0003307425910000115
为使用普罗尼序列对半潜式海洋平台的垂荡加速度进行归一化表征时使用的普罗尼序列参数。
因此，本实施例中基于普罗尼序列可以对直流信号、谐波信号、增幅(衰减)振荡信号进行拟合的优势，分别对垂荡加速度中的噪声成分，潮流变化等引起的缓慢变化以及加速度传感器引起的基线漂移误差项分别进行表征，从而将半潜式海洋平台的垂荡加速度归一化的表征为了普罗尼序列的形式。
(4)对归一化表征后的垂荡加速度进行漂移项去除，通过去除漂移项后剩余的普罗尼序列建立半潜式海洋平台垂荡加速度与垂荡运动参量之间的关系式，对半潜式海洋平台的垂荡运动参量进行估计，具体为：
根据计算出的普罗尼参数
Figure GDA0003307425910000116
确定归一化表征后的垂荡加速度各成分的频率，即：
Figure GDA0003307425910000121
对求解出的频率进行排序，并去除其中的极小频率成分(低频噪声以及加速度传感器带来的基线漂移)，即漂移项，则得到去除漂移项的归一化表征的垂荡加速度实测值为：
Figure GDA0003307425910000122
式中，
Figure GDA0003307425910000123
和
Figure GDA0003307425910000124
为使用普罗尼序列对去除漂移项的归一化表征的垂荡加速度实测值表征的普罗尼序列参数。
根据去除漂移项的归一化表征的垂荡加速度实测值，确定垂荡运动响应：
Figure GDA0003307425910000125
即垂荡加速度与垂荡运动参量之间的关系式为：
Figure GDA0003307425910000126
则将半潜式海洋平台的真实垂荡运动参量表征为：
Figure GDA0003307425910000127
因此，本实施例中通过对归一化表征后的垂荡加速度进行漂移项去除，去除了其中会引起垂荡运动参量漂移的低频项成分，然后使用剩余的普罗尼序列建立了半潜式海洋平台垂荡加速度与垂荡运动参量之间的关系，克服了传统方法基于积分和滤波器的病态问题。
综上可知，本发明提供的基于垂荡加速度的半潜式平台垂荡运动参量预报方法，主要基于线性势流理论推导了结构垂荡方向上的运动方程，并通过普罗尼序列建立了半潜式海洋平台垂荡加速度与垂荡运动参量之间的数学关系。该方法首先对垂荡加速度中的噪声信号、潮差等引起的缓慢信号以及加速度传感器引起的基线漂移成分进行归一化的表征，然后通过频率筛选去除了其中引起漂移的低频项成分，从而去除了其中加速度传感器带来的基线漂移成分以及低频噪声，最后利用剩余的普罗尼序列通过推导半潜式平台垂荡加速度的普罗尼序列与垂荡运动响应之间的数学关系，建立了结构垂荡加速度与垂荡运动参量之间的关系。本发明方法不同于传统基于滤波器的方法，通过普罗尼序列建立了半潜式平台垂荡加速度与位移之间的转换关系，而未使用传统积分和滤波器对垂荡运动参量数据进行校正，具有更高的预报精度。同时，在本发明方法中充分考虑了多方面的因素影响，包括海洋环境影响及传感器硬件设施影响等，使得本发明方法更具有实际应用价值。
下面以半潜式海洋平台一具体试验实例对上述方法进行验证：
本例采用了放置于波浪水槽中的半潜式海洋平台的运动响应数据进行计算分析，测试中使用造波机造波，并使用加速度传感器对半潜式海洋平台的垂荡加速度响应进行记录。同时，为了验证转换结果的正确性，使用光学六自由度仪对结构的垂荡运动参量进行了记录。试验平台搭建如图2所示。测试中激光位移传感器和加速度传感器的采样频率均设置为50Hz。
本例对安装的加速度传感器记录的半潜式海洋平台的垂荡加速度进行分析，实验中测试得到的波浪作用下的平台垂荡加速度响应如图3(a)所示。同时，为了验证本发明方法通过对加速度分析后得到的垂荡运动参量的正确性，实验中使用了光学六自由度仪进行半潜式平台的垂荡运动参量记录，测试的垂荡运动参量如图3(b)所示。从图中可以看到，半潜式平台由静止开始运动，为了考虑初始速度和位移的影响，在之后的分析中，选择图3中30秒到180秒的信号进行分析。
在分析中，首先由公式(12)对垂荡加速度实测值中的垂荡加速度理论值项、噪声项、低频变化项、基线漂移误差项进行统一的普罗尼序列表征，表征结果与测试加速度如图4(a)所示，图4(b)对100秒到110秒的表征结果进行局部放大，可以看到使用普罗尼序列可以较好的表征测试加速度信号。然后使用公式(13)对普罗尼序列进行筛选，去除其中的低频成分，最终得到剩余的普罗尼序列，如公式(14)。最后，将剩余的普罗尼序列代入公式(15)求解得到对应于结构垂荡加速度的垂荡运动参量。
然后，通过使用经过筛选后滤除漂移项的普罗尼序列，可以重构半潜式平台的真实垂荡运动参量，转换的结果如图5(a)所示，图5(b)对100秒到110秒的重构结果进行局部放大，与上述测试结果对比可知，可以看到使用剩余的普罗尼序列与光学六自由度仪测试的半潜式平台的垂荡运动参量具有较好的一致性，这也证明了本发明方法的正确性。
以上所述，仅是本发明的较佳实施例而已，并非是对本发明作其它形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例应用于其它领域，但是凡是未脱离本发明技术方案内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型，仍属于本发明技术方案的保护范围。

Claims (6)

1.一种基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法，其特征在于，包括：
在半潜式海洋平台垂荡运动中，基于线性势流理论，忽略附加质量和辐射阻尼的耦合影响，对半潜式海洋平台的垂荡加速度进行表征，确定垂荡加速度理论值；
考虑半潜式海洋平台垂荡运动实测海洋环境的噪声影响、环境缓慢变化引起的低频影响以及加速度传感器自身基线漂移误差引起的影响，引入噪声项、低频变化项、基线漂移误差项，确定垂荡加速度实测值；
对垂荡加速度实测值中的垂荡加速度理论值项、噪声项、低频变化项、基线漂移误差项进行统一的普罗尼序列归一化表征；
对归一化表征后的垂荡加速度进行漂移项去除，通过去除漂移项后剩余的普罗尼序列建立半潜式海洋平台垂荡加速度与垂荡运动参量之间的关系式，对半潜式海洋平台的垂荡运动参量进行估计。
2.根据权利要求1所述的基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法，其特征在于：
在半潜式海洋平台垂荡运动中，基于线性势流理论，忽略附加质量和辐射阻尼的耦合影响，同时考虑在流体中受到的波浪力、回复力和辐射力的影响，将垂荡运动表示为：
Figure FDA0003307425900000011
式中，m代表半潜式海洋平台的质量，
Figure FDA0003307425900000021
代表半潜式海洋平台的垂荡加速度，fw(t)为半潜式海洋平台受到的波浪荷载，fm(t)代表半潜式海洋平台受到的系泊力，fs(t)代表半潜式海洋平台受到的回复力，fr(t)代表半潜式海洋平台受到的辐射力；
其中，回复力fs(t)表示为：
fs(t)＝-czzo(t)＝-ρgAwzo(t) (2)
式中，zo(t)代表半潜式海洋平台的垂向位移；cz是半潜式海洋平台在垂荡方向的回复刚度，与水线面面积Aw，流体密度ρ以及重力加速度g有关；
辐射力fr(t)表示为：
Figure FDA0003307425900000022
式中，
Figure FDA0003307425900000023
代表半潜式海洋平台垂荡方向的速度，m∞和kz是无限频率处的垂荡方向的附加质量和脉冲响应函数；
由式(1)-(3)，将半潜式海洋平台垂荡运动表示为：
Figure FDA0003307425900000024
式中，f0(t)＝fw(t)+fm(t)；
则确定半潜式海洋平台垂荡加速度理论值表示为：
Figure FDA0003307425900000025
理论上，半潜式海洋平台的垂向加速度可以建模为一组谐波的叠加的形式，则由式(5)，将垂荡加速度理论值表征为：
Figure FDA0003307425900000031
式中，Ai、fi和θi分别代表垂向加速度中第i个分量的幅值、频率和相位，
Figure FDA0003307425900000032
和vi为使用普罗尼序列对半潜式海洋平台垂荡加速度理论值进行拟合时使用的参数。
3.根据权利要求2所述的基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法，其特征在于，考虑半潜式海洋平台垂荡运动实测海洋环境的噪声影响、环境缓慢变化引起的低频影响以及加速度传感器自身基线漂移误差引起的影响，引入噪声项、低频变化项、基线漂移误差项，确定垂荡加速度实测值为
Figure FDA0003307425900000033
式中，n(t)代表噪声项，v(t)代表低频变化项，b代表基线漂移误差项。
4.根据权利要求3所述的基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法，其特征在于：
引入普罗尼序列对垂荡加速度实测值中的噪声项、低频变化项、基线漂移误差项分别进行表征，即：
Figure FDA0003307425900000034
式中，
Figure FDA0003307425900000035
其中，An、fn、ξn和θn分别代表噪声项中各个成分的幅值、频率、阻尼和相位；
Figure FDA0003307425900000036
式中，
Figure FDA0003307425900000037
Dv＝-ξv+j2πfv，其中，Av、fv、ξv和θv分别代表低频变化项中各个成分的幅值、频率、阻尼和相位；
b＝EeFt (10)
式中，E和F为基线漂移误差项进行拟合时使用的参数；
由式(6)-(10)，将垂荡加速度实测值中的垂荡加速度理论值项、噪声项、低频变化项、基线漂移误差项进行统一的普罗尼序列表征，得到：
Figure FDA0003307425900000041
进一步对垂荡加速度实测值进行归一化表征：
Figure FDA0003307425900000042
式中，Np＝Ni+Nn+Nv+1,
Figure FDA0003307425900000043
和
Figure FDA0003307425900000044
为使用普罗尼序列对半潜式海洋平台的垂荡加速度进行归一化表征时使用的普罗尼序列参数。
5.根据权利要求4所述的基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法，其特征在于：
根据计算出的普罗尼参数
Figure FDA0003307425900000045
确定归一化表征后的垂荡加速度各成分的频率，即：
Figure FDA0003307425900000046
对求解出的频率进行排序，并去除其中的极小频率成分，即漂移项，则得到去除漂移项的归一化表征的垂荡加速度实测值为：
Figure FDA0003307425900000047
式中，
Figure FDA0003307425900000048
和
Figure FDA0003307425900000049
为使用普罗尼序列对去除漂移项的归一化表征的垂荡加速度实测值表征的普罗尼序列参数。
6.根据权利要求5所述的基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法，其特征在于：
根据去除漂移项的归一化表征的垂荡加速度实测值，确定垂荡运动响应：
Figure FDA0003307425900000051
即垂荡加速度与垂荡运动参量之间的关系式为：
Figure FDA0003307425900000052
则将半潜式海洋平台的真实垂荡运动参量表征为：
Figure FDA0003307425900000053
CN202011619590.8A 2020-12-30 2020-12-30 基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法 Active CN112693578B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202011619590.8A CN112693578B (zh) 2020-12-30 2020-12-30 基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法
PCT/CN2021/129514 WO2022142746A1 (zh) 2020-12-30 2021-11-09 基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法
US17/775,321 US11760449B2 (en) 2020-12-30 2021-11-09 Method for predicting heaving motion parameters of semi-submersible offshore platform based on heaving acceleration

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011619590.8A CN112693578B (zh) 2020-12-30 2020-12-30 基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法

Publications (2)

Publication Number Publication Date
CN112693578A CN112693578A (zh) 2021-04-23
CN112693578B true CN112693578B (zh) 2021-11-30

Family

ID=75512867

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011619590.8A Active CN112693578B (zh) 2020-12-30 2020-12-30 基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法

Country Status (3)

Country Link
US (1) US11760449B2 (zh)
CN (1) CN112693578B (zh)
WO (1) WO2022142746A1 (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112693578B (zh) 2020-12-30 2021-11-30 中国海洋大学 基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法
CN114194347B (zh) * 2022-01-05 2022-12-27 广东海洋大学 半潜式海洋平台的动力定位方法、装置、设备及介质
CN115659131B (zh) * 2022-09-27 2023-07-04 中国海洋大学 一种海工结构近场波浪传播方向识别方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2234060C1 (ru) * 2003-02-10 2004-08-10 Тульский государственный университет Устройство для измерения вертикальной качки судна
GB2440088B (en) * 2005-04-15 2009-12-30 Mitsui Shipbuilding Eng Automatic vessel position holding control method and controller
CN203349841U (zh) * 2013-08-07 2013-12-18 广州航海高等专科学校 一种舰船波浪运动检测装置
WO2017086482A1 (ja) * 2015-11-20 2017-05-26 流体テクノ株式会社 海象推定装置及び海象推定方法
CN107600352B (zh) * 2017-08-22 2019-03-05 哈尔滨工程大学 一种基于船舶水动力在线预报的船舶减摇控制系统
CN108549616B (zh) * 2018-02-28 2021-11-23 哈尔滨工程大学 基于g-n波浪模型的船舶在恶劣海况中大幅运动的预报方法
CN109446634B (zh) 2018-10-23 2022-07-15 哈尔滨工程大学 基于泰勒展开边界元方法的船舶运动预报方法
CN111336996A (zh) * 2020-03-28 2020-06-26 哈尔滨工程大学 一种海洋环境与船舶运动监测预报系统及监测预报方法
CN112693578B (zh) 2020-12-30 2021-11-30 中国海洋大学 基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法

Also Published As

Publication number Publication date
US11760449B2 (en) 2023-09-19
US20230129913A1 (en) 2023-04-27
WO2022142746A1 (zh) 2022-07-07
WO2022142746A9 (zh) 2022-11-10
CN112693578A (zh) 2021-04-23

Similar Documents

Publication Publication Date Title
CN112693578B (zh) 基于垂荡加速度的半潜式海洋平台垂荡运动参量预报方法
He et al. Optimized minimum generalized Lp/Lq deconvolution for recovering repetitive impacts from a vibration mixture
Gao et al. Improvement study of modal analysis for offshore structures based on reconstructed displacements
CN101221066A (zh) 工程非线性振动检测方法
CN109556700B (zh) 一种海底悬跨管道振动检测方法
CN113820003B (zh) 一种适用于桥梁振动监测的加速度实时重构动态位移方法
US8893550B2 (en) Non-invasive vibro-acoustic analysis
CN105654062A (zh) 海洋结构弱模态识别与时域重构方法
CN111323227A (zh) 一种航空发动机转子故障特征的提取方法
CN116975608A (zh) 一种基于窄带滤波的主成分分析的结构振型识别方法
JP5101327B2 (ja) 高架橋の揺れ方によるリアルタイム地震被害推定方法及びその装置
Cheng et al. AR model-based crosstalk cancellation method for operational transfer path analysis
Gelman The new second and higher order spectral technique for damage monitoring of structures and machinery
Randall et al. Applied digital signal processing
Zhang et al. Data Processing Based on Low-Precision IMU Equipment to Predict Wave Height and Wave Period
Kim Real-time Monitoring by Using Kalman Filter in the Ocean
CN112198911B (zh) 一种随机线谱自适应跟踪消除方法和系统
Sheen et al. A quantified index for bearing vibration analysis based on the resonance modes of mechanical system
Mosher Results from a new separation algorithm for planetary gear system vibration measurements
Sharkh et al. A novel Kalman filter based technique for calculating the time history of vertical displacement of a boat from measured acceleration
Schmittner et al. Application of focusing wave groups in model testing practice
Chai et al. Practical Application of ARMA Model Method for Blade Tip Timing
Zhang Hybrid wave models and their applications for steep ocean waves
Lu Removal of transfer function effects for gear condition monitoring using cepstrum-based operational modal analysis
Matthaiou On robust statistical outlier analysis for damage identification

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