Sea state estimation based on vessel responses

Sea state estimation from vessel motion responses: Smoothness and robustness using Bezier surface and L1 optimization

The complexity of the ocean environment is challenging for safe voyages and marine operations. An important consideration is thereby given to the sea states at the planning and operating stages for an offshore project. The wave spectra can be measured and detected by several measurement instruments. This paper presents nonparametric sea state estimation based on vessel motion responses (wave buoy analogy).

Because a vessel is asymmetric, the cross spectra of 6DOF vessel motions $S_{ij}^{wave}$ are related to wave spectra S_wave. The relationship between the vessel motions and wave spectra are characterized by response amplitude operators (RAO), which are influenced by the hull geometry and can be calculated by well-known hydrodynamic formulas. After discretizing the directional spectrum concerning incoming wave direction ${\beta }_k$ and wave frequency ${\omega }_m$, respectively, a linear equation $b=Af$ is formulated according to the hydrodynamics with a tall matrix $A$; shown in Figure 2. The vector of directional wave spectrum f can be solved as an inverse problem. However, the best-fit solution in this application does not always give a practical meaning. To solve the problem, typical least square (LS) methods are adopted with additional implicit smoothness conditions existing in a spectrum. The research innovatively augments this classic setup in two aspects.

First, two-dimensional smoothness is achieved using the Bézier surface. Instead of considering the constant slopes between a node and its two neighbor nodes in the wave direction and frequency, respectively, the size of the smoothness scope in the discretized network can be arbitrarily selected (see Figure 3). The new smoothness condition is a more general form, where a Bézier surface is the Cartesian product of two orthogonal Bezier curves. The classic constant-slope condition is a particular example of a Bézier curve with three nodes.

Second, the potential of sparseness in the ocean wave spectrum is realized through L1 optimization. The sparsity pattern of a wave spectrum with a threshold of 0.1 reveals the powerful components in the wave spectrum; see Figure 4. Converting from time-domain sampling, disturbances are noticed in the discretized motion cross-spectra, resulting in extra risks to the estimate results. Since L1 optimization is more robust to wildpoints, sparse regression is adopted in the inverse problem. The simulation results (e.g., in Figure 5) verify that the L1 norm $(1,1,1,1)$ has much higher robustness against the disturbances and errors in the vessel response spectrum than the LS methods $(2,2,2,2)$. The estimate given by the LS cost function $(p_1,r_1,p_2,r_2)=(2,2,2,2)$ is risky. The estimates of L1 optimization $(p_1,r_1,p_2,r_2)=(1,1,1,1)$ is “cleaner” showing a higher sparsity, and the energy in the estimated shows a better concentration.

Ref: Zhengru Ren, Xu Han, Amrit Shankar Verma, and Johann Alexander Dirdal, Roger Skjetne. Sea state estimation based on vessel motion responses: improved smoothness and robustness using Bézier surface and L1 optimization. Marine Structures, 2021. In press. Available online. [PDF]

Report: Norwegian SciTech News, EurekAlert, TU.no, Press-News, One American Northwest Arkansas, WorldNewsEra PopularMechanics