# A Novel Method to Estimate the Sea State for Recycling Work on the Sea Surface

Authors: Haixu Zhang; Xiaosu Xu; Guangcai Wang; Yongyun Zhu
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

The recycling of marine exploration equipment after it has surfaced is greatly affected by sea state. In order to estimate the sea state in real time, this paper proposes a method for measuring wave elevation, which modifies the integrated results of GNSS/SINS in the up direction by virtual horizontal lines to extract wave fluctuation information. From these wave information, the significant wave heights (SWH) can be calculated as the only input parameter of P-M spectrum, and a series of wave height data can be further simulated. When the GNSS is interrupted due to severe sea state, the simulated data can be integrated with the SINS to deal with the data distortion problem. The simulation results show that the application of wave spectrum in the GNSS intermittent situation has obvious improvement effect and important significance.


---

## Body

## 1. Introduction

From the experience of engineering practice, the recycling of marine detection equipment is often carried out in the good sea state, as shown in Figure1, because the safety of the crane operation will be seriously threatened by the shaking of the ship in the bad sea state. Although the recycling work can be arranged in advance according to weather forecast, the complex sea state still cannot be estimated very well, especially the impact of swell [1, 2]. Missing the best recycling time may cause the device to run out of power, lose contact, and disappear into the sea. Therefore, the acquisition of real-time sea state information is of great significance in the process of equipment recycling.Figure 1
Equipment to be recycled in the good sea state.In order to obtain wave elevation data under different application backgrounds, many scholars have carried out research studies on it. Using buoys to collect wave data is the most common and effective method. Crosby et al. [3] provide a new method for estimating offshore boundary conditions according to the buoy observations at 16 near-shore buoy sites. Buoys are also used to validate Haiyang-2A altimeter SWH data [4 and 5]. In addition to buoys, high-frequency surface wave radar, X-band marine radar, and satellite can also work for wave data acquisition [6–9].However, for detection equipment such as AUVs and underwater gliders, their position after surfacing is uncertain beforehand, and buoys cannot be preset to measure the sea conditions around them. The use of satellite or radar to measure wave heights requires additional hardware support. As the only subsystem that can sense wave elevation that they are equipped with, the strap-down inertial navigation system (SINS) cannot be simply used to obtain the sea state due to its error accumulation.The SINS is irreplaceable in the aspect of navigation with the restrictions of other sensors because of its high precision in short time [10–13]. Especially, the combination of the SINS and GNSS has been widely used [14–19], which not only limits the divergence of the SINS but also can provide high-positioning accuracy. If applied to the detection of wave elevation on the sea surface in the up direction, the following two problems need to be solved: (1)
The error of the integrated method is too large, which covers the real wave elevation information(2)
Limited by the output frequency of the GNSS, the real-time performance of the results needs to be improvedIn order to solve the above problems and provide the real-time sea state of the target equipment, a novel fusion method is proposed to estimate the sea surface wave elevation. In this way, the crew can get real-time information about the sea conditions around the ocean detection equipment to avoid working under bad sea conditions and arrange the recycling work reasonably. The method makes full use of the high precision of the SINS in short time to sensitive wave characteristic and establishes the wave model of current time to fuse information in the next moment.The remainder of the paper is organized as follows. The mathematical model is described in Section2. Section 3 introduces the fusion about SINS, GNSS, and simulated waves. The simulations and semiphysics experiment are discussed in Section 4. Finally, Section 5 concludes this paper.


## 2. System Model

### 2.1. SINS State Equation

The physical parameters measured by gyroscopes and accelerometers of the SINS are used to calculate the trajectory of the carrier in the up direction when floating on the sea surface. The state equation is given as follows according to Wang et al [20]:(1)X˙=FX+GW,where the state vector can be expressed as follows:(2)X=ϕxϕyϕzδVEnδVNnδVUnδλδLδhεxεyεz∇x∇y∇zT,where ϕx, ϕy, and ϕz denote misalignment angles, δVEn, δVNn, and δVUn denote velocity errors, δλ, δL, and δh denote the longitude, latitude, and height errors, respectively, εx, εy, and εz denote gyroscope biases in three directions of the frame b, and ∇x, ∇y, and ∇z denote accelerometer biases.The state transition matrixF is expressed as follows:(3)F=MaaMavMap−Cbn03×3MvaMvvMvp03×3Cbn03×3MpvMpp03×303×306×306×306×306×306×3,where direction cosine matrix Cbn denotes the rotational transform from the b-frame to the n-frame. And, we have(4)Maa=0ωiesinL+VERN+htanL−ωiecosL−VERN+h−ωiesinL−VERN+htanL0−VNRM+hωiecosL+VERN+hVNRM+h0,(5)Mav=0−1RM+h01RN+h00tanLRN+h00,(6)Map=00VNRM+h2−ωiesinL0−VERN+h2ωiecosL+VEsec2LRN+h0−VEtanLRN+h2,(7)Mva=0−fUnfNnfUn0−fEn−fNnfEn0,(8)Mvv=VNtanL−VURN+h2ωiesinL+VEtanLRN+h−2ωiecosL−VERN+h−2ωiesinL−2VEtanLRN+h−VURM+h−VNRM+h2ωiecosL+2VERN+h2VNRM+h0,(9)Mvp=Vn×2M1+M3+M4,(10)M1=000−ωiesinL00ωiecosL00,(11)M3=00VNRM+h200−VERN+h2VEsec2LRN+h0−VEtanLRN+h2,(12)M4=000000−gesin2Lβ−4β1cos2L0β2,(13)Mpv=01RM+h0secLRN+h00001,(14)Mpp=00−VNRM+h2VEsecLtanLRN+h0−VEsecLRN+h2000,where ωie is the rotational angular velocity of the Earth, L and h are the position components in latitude and altitude, VE, VN, and VU are the velocity components in the east, north, and up directions, fNn, fEn, and fUn are the measurements of specific force in the east, north, and up directions, ge is the acceleration of gravity, RM and RN are the radii of curvature in meridian and prime vertical [21], β is gravity flattening, and f is flattening of ellipsoid, and we have(15)β1=2βf+f28,(16)β2=3.08×10−6s−2.The system noise drive matrixG can be expressed as follows:(17)G=−Cbn03×303×3Cbn09×309×3.And, the noise vectorW is written as follows:(18)W=wgbwab,where wgb and wab denote angular velocity measurement noise of gyroscopes and specific force measurement noise of accelerometers, respectively.


### 2.2. Measurement Equation

Satellite positioning is very effective in dealing with the divergence of the SINS, so the GNSS in the detector can be used to correct the accumulated errors in the up direction. The measurement equation of the GNSS/SINS integrated method is defined as follows:(19)Z=HX+V,where the measurement vector is expressed as follows:(20)Z=hSINS−hGNSS,where hSINS and hGNSS denote the measurements in the up direction of the SINS and GNSS, respectively.The observation matrixH and observation error matrix V are expressed as follows:(21)H=01×8I1×101×6,(22)V=vz,where the variable vz is the measurement error of the GNSS in the up direction. Although corrected by the GNSS, the error is still unacceptable relative to the wave heights.


## 2.1. SINS State Equation

The physical parameters measured by gyroscopes and accelerometers of the SINS are used to calculate the trajectory of the carrier in the up direction when floating on the sea surface. The state equation is given as follows according to Wang et al [20]:(1)X˙=FX+GW,where the state vector can be expressed as follows:(2)X=ϕxϕyϕzδVEnδVNnδVUnδλδLδhεxεyεz∇x∇y∇zT,where ϕx, ϕy, and ϕz denote misalignment angles, δVEn, δVNn, and δVUn denote velocity errors, δλ, δL, and δh denote the longitude, latitude, and height errors, respectively, εx, εy, and εz denote gyroscope biases in three directions of the frame b, and ∇x, ∇y, and ∇z denote accelerometer biases.The state transition matrixF is expressed as follows:(3)F=MaaMavMap−Cbn03×3MvaMvvMvp03×3Cbn03×3MpvMpp03×303×306×306×306×306×306×3,where direction cosine matrix Cbn denotes the rotational transform from the b-frame to the n-frame. And, we have(4)Maa=0ωiesinL+VERN+htanL−ωiecosL−VERN+h−ωiesinL−VERN+htanL0−VNRM+hωiecosL+VERN+hVNRM+h0,(5)Mav=0−1RM+h01RN+h00tanLRN+h00,(6)Map=00VNRM+h2−ωiesinL0−VERN+h2ωiecosL+VEsec2LRN+h0−VEtanLRN+h2,(7)Mva=0−fUnfNnfUn0−fEn−fNnfEn0,(8)Mvv=VNtanL−VURN+h2ωiesinL+VEtanLRN+h−2ωiecosL−VERN+h−2ωiesinL−2VEtanLRN+h−VURM+h−VNRM+h2ωiecosL+2VERN+h2VNRM+h0,(9)Mvp=Vn×2M1+M3+M4,(10)M1=000−ωiesinL00ωiecosL00,(11)M3=00VNRM+h200−VERN+h2VEsec2LRN+h0−VEtanLRN+h2,(12)M4=000000−gesin2Lβ−4β1cos2L0β2,(13)Mpv=01RM+h0secLRN+h00001,(14)Mpp=00−VNRM+h2VEsecLtanLRN+h0−VEsecLRN+h2000,where ωie is the rotational angular velocity of the Earth, L and h are the position components in latitude and altitude, VE, VN, and VU are the velocity components in the east, north, and up directions, fNn, fEn, and fUn are the measurements of specific force in the east, north, and up directions, ge is the acceleration of gravity, RM and RN are the radii of curvature in meridian and prime vertical [21], β is gravity flattening, and f is flattening of ellipsoid, and we have(15)β1=2βf+f28,(16)β2=3.08×10−6s−2.The system noise drive matrixG can be expressed as follows:(17)G=−Cbn03×303×3Cbn09×309×3.And, the noise vectorW is written as follows:(18)W=wgbwab,where wgb and wab denote angular velocity measurement noise of gyroscopes and specific force measurement noise of accelerometers, respectively.


## 2.2. Measurement Equation

Satellite positioning is very effective in dealing with the divergence of the SINS, so the GNSS in the detector can be used to correct the accumulated errors in the up direction. The measurement equation of the GNSS/SINS integrated method is defined as follows:(19)Z=HX+V,where the measurement vector is expressed as follows:(20)Z=hSINS−hGNSS,where hSINS and hGNSS denote the measurements in the up direction of the SINS and GNSS, respectively.The observation matrixH and observation error matrix V are expressed as follows:(21)H=01×8I1×101×6,(22)V=vz,where the variable vz is the measurement error of the GNSS in the up direction. Although corrected by the GNSS, the error is still unacceptable relative to the wave heights.


## 3. Improved Method for Estimating Sea State

The solution of the SINS in the up direction represents the characteristic of the target moving along with waves on the sea surface; thus, it can provide a reference for staffs to know the sea state around the remote target. However, the SINS will diverge quickly in the up direction without any other sensors to perform correction. The GNSS plays a certain role in the correction of the SINS, but considering its own error and real-time performance, the results obtained by means of the integrated method are still not enough to represent the characteristics of wave heights. Therefore, we need to find another way to extract the waveform from the existing conditions.


### 3.1. Corrected Waves from SINS/GNSS

The SINS can observe the wave heights with good precision and high frequency. However, the wave heights observed by the SINS in the up direction always deviate from the horizontal line due to error accumulation, as shown in Figure2.Figure 2
The time history curve of wave elevation observed by the SINS.Through the comparison in Figure2, we can also get that the wave fluctuation is retained in the initial period of time, while gradually it disappears in the following time. Therefore, in order to obtain wave information, the divergence problem must be solved first.The combination of the GNSS and SINS is a classic way to weigh divergence and accuracy in navigation. It is also effective in solving the positioning problem in the up direction, as shown in Figure3. However, although the divergence problem is solved, the magnitude of the integrated method errors is still unacceptable relative to the wave heights.Figure 3
Comparison of the integrated results and true values in the up direction.Waves always float about the horizontal line up and down, but the error accumulation of the SINS seems to make the horizontal line no longer horizontal. To obtain the heights more accurately, a simulation of the virtual horizontal line for the waves in the data fragment is needed first, as shown in Figure4.Figure 4
The schematic diagram of the virtual horizontal line.In the data segment with the virtual horizontal line as the reference, different from the observed valueh, wave height vh is the difference value in the up direction in a short time, and the problem of error accumulation can be less concerned. In addition, t1 and t2 are the time corresponding to the measured wave heights, while vt1 and vt2 are to the virtual horizontal line. Therefore, the trajectory of the waves can be corrected after finding the virtual horizontal line. The least squares method is an effective tool to solve this problem, whose minimum mean square error function can be expressed as follows:(23)Jθ=B−TθTB−Tθ,where(24)B=h1h2⋯hnT,(25)T=t11⋮⋮tn1,(26)θ=ab,where n is the number of measured values in a data segment. The virtual horizontal line can be approximated as follows:(27)h=at+b.Waves can be considered as the result of the superposition of countless random simple cosine waves, so the average height of the resultant wave surface is zero. According to this principle, the sum of all wave heights to the virtual horizontal line should also be zero. Considering the errors, the sum of wave heights can be denoted as follows:(28)SUMd1⋯di⋯dn=ε,where di is the distance from the ith observed wave height to the virtual horizontal line, whose sign depends on the position relative to the virtual horizontal line. In addition, ε comes from system errors. Then, the sample values referring to the true horizontal line can be corrected as follows:(29)chi=di−εn.The corrected sample valuechi can be used to get the wave information referring to the real horizontal line.


### 3.2. Estimating Waves by P-M Spectrum

The above method to calculate wave heights is based on the SINS and GNSS integrated method in good condition. However, GNSS signal interruption is not rare for the detector floating on the sea surface in the harsh marine environment. In order to extract the wave data more accurately under finite conditions, the model of wave mechanics for ocean engineering is introduced in this paper.Waves are unpredictable due to the influence of many factors. Although wave elevation is a random value that changes with time, it can still be characterized by mathematical simulation due to its stationarity and ergodic property of various states [22]. As a single input parameter spectrum, the P-M spectrum, according to Shouhua et al. [23], can be used for wave elevation simulation based on the acquired SWH by the SINS, which is defined as follows:(30)Sω=Aω−5exp−Bω−4,where(31)A=αg2,(32)B=0.74gV19.44.Phillips constantα in the formula above is defined as follows:(33)α=8.1×10−3.The wave spectrum in equation (30) contains only one parameter, namely, V19.4, which represents the wind speed at 19.4 meters above the sea surface. However, IMU cannot provide wind speed information, so further conversion is needed.According to the following definition of the average wave height of partial big waves,(34)H1/p=8m0lnp1/2+p2πm01/21−erflnp1/2.We obtain the average wave elevation of the first third of the waves as follows:(35)H1/3=Hs≈4.0m01/2,where Hs represents the SWH and m0 is the zeroth moment of the spectrum expressed as follows:(36)m0=A4B.The valueHs can be approximated as below according to equations (35) and (36):(37)Hs≈0.209V19.42g.It can be seen that the wave elevation is proportional to the square of the wind speed, so the wave spectrum represented by the single parameter of the wave elevation can be written as follows:(38)Sω=0.78ω5exp−3.11ω4Hs2.The constituent waves of high frequency and small frequency provide little energy, and the bulk of the energy comes from the constituent waves with the narrow band. If the part of the total energy is allowed to be omitted on the high and low sides of the spectrum, the upper and lower bounds of frequency of the constituent waves can be written as follows:(39)ωL=−3.11Hs2lnμ1/4,(40)ωH=−3.11Hs2ln1−μ1/4.In engineering applications, the time history of a series of irregular waves can be simulated by the target spectrum after getting the characteristic wave parameters. After taking P-M spectrum as the target spectrum and dividing the total energy intoN parts according to the equipartition energy method, the waves’ time history can be expressed as follows:(41)ηt=2m0N1/2∑i=1Ncosω^it+φi,where(42)ω^i=ωi+ωi−12.In addition, the valueφi is the initial phase uniformly distributed in the interval of [0, 2π], while ωi is the boundary frequency defined as follows:(43)ωi=BlnN/i1/4.According to equations (38) and (41), a part of the waves can be simulated, as shown in Figure 5, with SWH of 3 meters as the only parameter of P-M spectrum.Figure 5
Wave simulation by P-M spectrum.The waves simulated by the spectrum satisfy the stationarity and ergodicity of the general waves statistically, which fluctuates around the horizontal line, and the mean value is close to zero. As long as the SWH of the previous moment is collected, the approximate value of the wave heights in the next short moment can be estimated by this method.Compared with other spectra, P-M spectrum is more suitable for the engineering application proposed in this paper because its only input parameter is more easy to be obtained by the SINS. However, the integral results of the SINS will diverge without the correction of the GNSS, which indirectly determines that the wave heights simulated by P-M spectrum cannot be applied independently.


### 3.3. Get Wave Heights with Intermittent GNSS

It is not uncommon for the surfaced detection equipment to lose GNSS signals in the bad sea state. Without effective correction, the divergence of the SINS in the up direction will weaken its ability to extract wave features, as shown in Figure6.Figure 6
The wave height results with intermittent GNSS signals.
(a)(b)As can be seen from the above figures, (a) is the measurement results at the beginning of the intermittent state, while (b) represents that the state lasted for a period. Although the measurement effect is good at the beginning, as time goes by, the divergence of the SINS cannot be effectively corrected, and the wave height data calculated according to equation (29) gradually distorts until the characteristics of wave heights cannot be recognized at last.In order to improve this problem, we tried to use wave spectrum to simulate wave height data to assist the SINS to maintain a considerable height measurement in the intermittent GNSS state. The only input parameter of P-M spectrum is SWH according to equation (38). Therefore, it is necessary to analyze the latest measured data fragment to calculate SWH before the GNSS signal is interrupted. After obtaining the corrected wave height information by equation (29), the zero-up crossing method, as shown in Figure 7, is used to calculate wave heights, where time t1 and t2 are the zero-up crossing points and h1 and h2 are the wave heights, respectively.Figure 7
The wave height results with intermittent GNSS signals.After taking the average wave elevation of the first third of the waves as SWH, the wave elevation can be immediately simulated by P-M spectrum as follows according to equation (41):(44)ηt=η1η2η3…ηn,where the time interval can be set small enough to obtain a suitable sampling frequency to match the SINS. In addition, the measurement equation of the SINS and simulated wave elevation integrated method is defined as follows:(45)Zs=HX+Vs,where the observation error matrix Vs represents the error of the simulated wave elevation, and the measurement vector is expressed as follows:(46)Zs=hSINS−ηt.The integrated method can fill the gap of GNSS interruption in bad sea conditions. In addition, under the correction of the connected GNSS, the wave elevation information in the environment of the surface detection equipment can be output to provide the sea state information continuously and stably.


## 3.1. Corrected Waves from SINS/GNSS

The SINS can observe the wave heights with good precision and high frequency. However, the wave heights observed by the SINS in the up direction always deviate from the horizontal line due to error accumulation, as shown in Figure2.Figure 2
The time history curve of wave elevation observed by the SINS.Through the comparison in Figure2, we can also get that the wave fluctuation is retained in the initial period of time, while gradually it disappears in the following time. Therefore, in order to obtain wave information, the divergence problem must be solved first.The combination of the GNSS and SINS is a classic way to weigh divergence and accuracy in navigation. It is also effective in solving the positioning problem in the up direction, as shown in Figure3. However, although the divergence problem is solved, the magnitude of the integrated method errors is still unacceptable relative to the wave heights.Figure 3
Comparison of the integrated results and true values in the up direction.Waves always float about the horizontal line up and down, but the error accumulation of the SINS seems to make the horizontal line no longer horizontal. To obtain the heights more accurately, a simulation of the virtual horizontal line for the waves in the data fragment is needed first, as shown in Figure4.Figure 4
The schematic diagram of the virtual horizontal line.In the data segment with the virtual horizontal line as the reference, different from the observed valueh, wave height vh is the difference value in the up direction in a short time, and the problem of error accumulation can be less concerned. In addition, t1 and t2 are the time corresponding to the measured wave heights, while vt1 and vt2 are to the virtual horizontal line. Therefore, the trajectory of the waves can be corrected after finding the virtual horizontal line. The least squares method is an effective tool to solve this problem, whose minimum mean square error function can be expressed as follows:(23)Jθ=B−TθTB−Tθ,where(24)B=h1h2⋯hnT,(25)T=t11⋮⋮tn1,(26)θ=ab,where n is the number of measured values in a data segment. The virtual horizontal line can be approximated as follows:(27)h=at+b.Waves can be considered as the result of the superposition of countless random simple cosine waves, so the average height of the resultant wave surface is zero. According to this principle, the sum of all wave heights to the virtual horizontal line should also be zero. Considering the errors, the sum of wave heights can be denoted as follows:(28)SUMd1⋯di⋯dn=ε,where di is the distance from the ith observed wave height to the virtual horizontal line, whose sign depends on the position relative to the virtual horizontal line. In addition, ε comes from system errors. Then, the sample values referring to the true horizontal line can be corrected as follows:(29)chi=di−εn.The corrected sample valuechi can be used to get the wave information referring to the real horizontal line.


## 3.2. Estimating Waves by P-M Spectrum

The above method to calculate wave heights is based on the SINS and GNSS integrated method in good condition. However, GNSS signal interruption is not rare for the detector floating on the sea surface in the harsh marine environment. In order to extract the wave data more accurately under finite conditions, the model of wave mechanics for ocean engineering is introduced in this paper.Waves are unpredictable due to the influence of many factors. Although wave elevation is a random value that changes with time, it can still be characterized by mathematical simulation due to its stationarity and ergodic property of various states [22]. As a single input parameter spectrum, the P-M spectrum, according to Shouhua et al. [23], can be used for wave elevation simulation based on the acquired SWH by the SINS, which is defined as follows:(30)Sω=Aω−5exp−Bω−4,where(31)A=αg2,(32)B=0.74gV19.44.Phillips constantα in the formula above is defined as follows:(33)α=8.1×10−3.The wave spectrum in equation (30) contains only one parameter, namely, V19.4, which represents the wind speed at 19.4 meters above the sea surface. However, IMU cannot provide wind speed information, so further conversion is needed.According to the following definition of the average wave height of partial big waves,(34)H1/p=8m0lnp1/2+p2πm01/21−erflnp1/2.We obtain the average wave elevation of the first third of the waves as follows:(35)H1/3=Hs≈4.0m01/2,where Hs represents the SWH and m0 is the zeroth moment of the spectrum expressed as follows:(36)m0=A4B.The valueHs can be approximated as below according to equations (35) and (36):(37)Hs≈0.209V19.42g.It can be seen that the wave elevation is proportional to the square of the wind speed, so the wave spectrum represented by the single parameter of the wave elevation can be written as follows:(38)Sω=0.78ω5exp−3.11ω4Hs2.The constituent waves of high frequency and small frequency provide little energy, and the bulk of the energy comes from the constituent waves with the narrow band. If the part of the total energy is allowed to be omitted on the high and low sides of the spectrum, the upper and lower bounds of frequency of the constituent waves can be written as follows:(39)ωL=−3.11Hs2lnμ1/4,(40)ωH=−3.11Hs2ln1−μ1/4.In engineering applications, the time history of a series of irregular waves can be simulated by the target spectrum after getting the characteristic wave parameters. After taking P-M spectrum as the target spectrum and dividing the total energy intoN parts according to the equipartition energy method, the waves’ time history can be expressed as follows:(41)ηt=2m0N1/2∑i=1Ncosω^it+φi,where(42)ω^i=ωi+ωi−12.In addition, the valueφi is the initial phase uniformly distributed in the interval of [0, 2π], while ωi is the boundary frequency defined as follows:(43)ωi=BlnN/i1/4.According to equations (38) and (41), a part of the waves can be simulated, as shown in Figure 5, with SWH of 3 meters as the only parameter of P-M spectrum.Figure 5
Wave simulation by P-M spectrum.The waves simulated by the spectrum satisfy the stationarity and ergodicity of the general waves statistically, which fluctuates around the horizontal line, and the mean value is close to zero. As long as the SWH of the previous moment is collected, the approximate value of the wave heights in the next short moment can be estimated by this method.Compared with other spectra, P-M spectrum is more suitable for the engineering application proposed in this paper because its only input parameter is more easy to be obtained by the SINS. However, the integral results of the SINS will diverge without the correction of the GNSS, which indirectly determines that the wave heights simulated by P-M spectrum cannot be applied independently.


## 3.3. Get Wave Heights with Intermittent GNSS

It is not uncommon for the surfaced detection equipment to lose GNSS signals in the bad sea state. Without effective correction, the divergence of the SINS in the up direction will weaken its ability to extract wave features, as shown in Figure6.Figure 6
The wave height results with intermittent GNSS signals.
(a)(b)As can be seen from the above figures, (a) is the measurement results at the beginning of the intermittent state, while (b) represents that the state lasted for a period. Although the measurement effect is good at the beginning, as time goes by, the divergence of the SINS cannot be effectively corrected, and the wave height data calculated according to equation (29) gradually distorts until the characteristics of wave heights cannot be recognized at last.In order to improve this problem, we tried to use wave spectrum to simulate wave height data to assist the SINS to maintain a considerable height measurement in the intermittent GNSS state. The only input parameter of P-M spectrum is SWH according to equation (38). Therefore, it is necessary to analyze the latest measured data fragment to calculate SWH before the GNSS signal is interrupted. After obtaining the corrected wave height information by equation (29), the zero-up crossing method, as shown in Figure 7, is used to calculate wave heights, where time t1 and t2 are the zero-up crossing points and h1 and h2 are the wave heights, respectively.Figure 7
The wave height results with intermittent GNSS signals.After taking the average wave elevation of the first third of the waves as SWH, the wave elevation can be immediately simulated by P-M spectrum as follows according to equation (41):(44)ηt=η1η2η3…ηn,where the time interval can be set small enough to obtain a suitable sampling frequency to match the SINS. In addition, the measurement equation of the SINS and simulated wave elevation integrated method is defined as follows:(45)Zs=HX+Vs,where the observation error matrix Vs represents the error of the simulated wave elevation, and the measurement vector is expressed as follows:(46)Zs=hSINS−ηt.The integrated method can fill the gap of GNSS interruption in bad sea conditions. In addition, under the correction of the connected GNSS, the wave elevation information in the environment of the surface detection equipment can be output to provide the sea state information continuously and stably.


## 4. Experiment and Simulation

In order to demonstrate the effectiveness of the integrated method by wave spectrum for the surfaced detection equipment to estimate the sea state, experiment and simulation will be carried out to verify the following aspects:(1)
Advantages of P-M spectrum in wave elevation simulation(2)
The performance of the virtual horizontal line in wave elevation extraction(3)
The response effect of the novel integrated method to GNSS interruptionTherefore, each experiment or simulation is carried out in the following order.


### 4.1. P-M Spectrum for Simulating Wave Elevation

Before using wave spectrum to integrate the SINS in the up direction, it is necessary to analyze the statistical characteristic of simulated waves to verify its feasibility. A set of wave elevation data collected in the offshore area of the Yellow Sea in China is selected as the reference sample to verify the advantages of wave spectrum in wave elevation simulation. The comparisons of true and simulated waves (SW) are shown in Figure8.Figure 8
Comparison of true and simulated values of wave elevation.Since the waves can be regarded as a random process, the simulation based on the wave spectrum cannot reproduce the motion trajectory, but it can be used to approximate the true values statistically. Furthermore, 100 simulations were carried out according to the same series of true waves, and the statistical parameters are shown in Figure9.Figure 9
Comparison of statistical parameters.The waves could be regarded as the result of superposition of an infinite number of random simple cosine waves, so the average height of the waves should be zero in theory. In this respect, wave simulation works well as shown above. The variance of wave elevation can be regarded as the wave energy. When the waves are in a stable stage, the variance changes slowly. Limited by the sampling frequency of the wave buoy, the variance of the true values in Figure9 is slightly higher than the simulated values. At the same time, the sensitivity of wave simulation to SWH is within the acceptable range. In general, the simulated waves have good statistical characteristics, and more importantly, their sampling frequency is not limited, which will make their positioning assistance in the up direction more meaningful.


### 4.2. Advantage of the Virtual Horizontal Line in the Up Direction

After verifying the feasibility of wave spectrum simulation, the motion state simulation of the target in the direction of up floating on the sea surface is shown in Figure10. The target moves up and down according to the trajectory of wave elevation, with uniform distribution of velocity and acceleration, which conforms to the characteristics of the stable wave motion.Figure 10
The simulated motion state.The positioning results of the GNSS/SINS in the up direction are shown in Figure11. It can be seen that, with the help of the GNSS, the divergence of the SINS in the up direction is solved, but it is still not enough to represent the characteristic of wave elevation.Figure 11
GNSS/SINS positioning results in the up direction.In view of the problem in Figure11, the virtual horizontal line of the original data was calculated according to the idea in Figure 4, and the wave information is extracted, as shown in Figure 12, according to equation (29).Figure 12
The corrected results from the GNSS/SINS integrated method.As can be seen from the figure, the revised measurements are closer to the longitudinal trajectories of the waves without divergence, which is sufficient to estimate the sea state. Thus, this approach, which uses the virtual horizontal line to correct GNSS/SINS integrated data, performs well at extracting information in the Up direction.


### 4.3. Integrated Method When GNSS Breaks

The GNSS is an important means to correct the divergence of the SINS, and its combination with the SINS can figure wave height information well. However, when the severe sea conditions cause the GNSS to be intermittent for a long time, the results of the SINS will diverge because they cannot be corrected in time. The wave height data simulated by wave spectrum can predict the waves in a short period statistically. Applying it to the measurement correction during GNSS interruption can provide a significant improvement effect, as shown in Figure13.Figure 13
Advantages of SW applied in the integrated method during GNSS interruption.
(a)(b)As can be seen from Figure13(a), compared with the reference value waves, the results of the GNSS/SINS integrated method did not perform well and could not obtain the wave height information. That is, because the interruption of the GNSS caused the SINS unable to be effectively corrected, the divergence of the SINS in the up direction gradually weakened the wave height information. The integrated method of GNSS/SINS/SW can fetch the wave height information very well, and the effect is better compared with the reference value. This is attributed to the good performance of SW in the Up direction during GNSS interruption, which made up for the deficiency of hardware caused by the GNSS failure with effective correction. The error comparison in Figure 13(b) further illustrates the effectiveness of the GNSS/SINS/SW method, whose error is far less than that of the GNSS/SINS method, which can be controlled at about 0.3 meters or less.In order to show the excellent performance of the method in a more comprehensive way, the average wave heights (AWH) and SWH are calculated every five minutes for ten hours, and the results are shown in Figure14.Figure 14
The performance of GNSS/SINS/SW in AWH and SWH.
(a)(b)AWH and SWH are two important parameters to estimate the sea state. Especially, in the application background of this paper, the staff can estimate the sea state based on this information and then make a judgment on whether it is suitable for the recycling work. According to the statistical results of 120 samples in Figure14(a), AWH and SWH calculated by the GNSS/SINS/SW method has a certain difference between the calculated results and reference value at some sampling points, but the errors are still within the acceptable range. In general, AWH and SWH calculated by this method perform well overall, as can be seen intuitively in Figure 14(b). Their errors can be controlled at about 0.1 meters, and in some better performances, the errors are so small that can almost be ignored in the engineering application of this paper. The results can be used to estimate the sea state well and provide an effective reference for the recycling work.


## 4.1. P-M Spectrum for Simulating Wave Elevation

Before using wave spectrum to integrate the SINS in the up direction, it is necessary to analyze the statistical characteristic of simulated waves to verify its feasibility. A set of wave elevation data collected in the offshore area of the Yellow Sea in China is selected as the reference sample to verify the advantages of wave spectrum in wave elevation simulation. The comparisons of true and simulated waves (SW) are shown in Figure8.Figure 8
Comparison of true and simulated values of wave elevation.Since the waves can be regarded as a random process, the simulation based on the wave spectrum cannot reproduce the motion trajectory, but it can be used to approximate the true values statistically. Furthermore, 100 simulations were carried out according to the same series of true waves, and the statistical parameters are shown in Figure9.Figure 9
Comparison of statistical parameters.The waves could be regarded as the result of superposition of an infinite number of random simple cosine waves, so the average height of the waves should be zero in theory. In this respect, wave simulation works well as shown above. The variance of wave elevation can be regarded as the wave energy. When the waves are in a stable stage, the variance changes slowly. Limited by the sampling frequency of the wave buoy, the variance of the true values in Figure9 is slightly higher than the simulated values. At the same time, the sensitivity of wave simulation to SWH is within the acceptable range. In general, the simulated waves have good statistical characteristics, and more importantly, their sampling frequency is not limited, which will make their positioning assistance in the up direction more meaningful.


## 4.2. Advantage of the Virtual Horizontal Line in the Up Direction

After verifying the feasibility of wave spectrum simulation, the motion state simulation of the target in the direction of up floating on the sea surface is shown in Figure10. The target moves up and down according to the trajectory of wave elevation, with uniform distribution of velocity and acceleration, which conforms to the characteristics of the stable wave motion.Figure 10
The simulated motion state.The positioning results of the GNSS/SINS in the up direction are shown in Figure11. It can be seen that, with the help of the GNSS, the divergence of the SINS in the up direction is solved, but it is still not enough to represent the characteristic of wave elevation.Figure 11
GNSS/SINS positioning results in the up direction.In view of the problem in Figure11, the virtual horizontal line of the original data was calculated according to the idea in Figure 4, and the wave information is extracted, as shown in Figure 12, according to equation (29).Figure 12
The corrected results from the GNSS/SINS integrated method.As can be seen from the figure, the revised measurements are closer to the longitudinal trajectories of the waves without divergence, which is sufficient to estimate the sea state. Thus, this approach, which uses the virtual horizontal line to correct GNSS/SINS integrated data, performs well at extracting information in the Up direction.


## 4.3. Integrated Method When GNSS Breaks

The GNSS is an important means to correct the divergence of the SINS, and its combination with the SINS can figure wave height information well. However, when the severe sea conditions cause the GNSS to be intermittent for a long time, the results of the SINS will diverge because they cannot be corrected in time. The wave height data simulated by wave spectrum can predict the waves in a short period statistically. Applying it to the measurement correction during GNSS interruption can provide a significant improvement effect, as shown in Figure13.Figure 13
Advantages of SW applied in the integrated method during GNSS interruption.
(a)(b)As can be seen from Figure13(a), compared with the reference value waves, the results of the GNSS/SINS integrated method did not perform well and could not obtain the wave height information. That is, because the interruption of the GNSS caused the SINS unable to be effectively corrected, the divergence of the SINS in the up direction gradually weakened the wave height information. The integrated method of GNSS/SINS/SW can fetch the wave height information very well, and the effect is better compared with the reference value. This is attributed to the good performance of SW in the Up direction during GNSS interruption, which made up for the deficiency of hardware caused by the GNSS failure with effective correction. The error comparison in Figure 13(b) further illustrates the effectiveness of the GNSS/SINS/SW method, whose error is far less than that of the GNSS/SINS method, which can be controlled at about 0.3 meters or less.In order to show the excellent performance of the method in a more comprehensive way, the average wave heights (AWH) and SWH are calculated every five minutes for ten hours, and the results are shown in Figure14.Figure 14
The performance of GNSS/SINS/SW in AWH and SWH.
(a)(b)AWH and SWH are two important parameters to estimate the sea state. Especially, in the application background of this paper, the staff can estimate the sea state based on this information and then make a judgment on whether it is suitable for the recycling work. According to the statistical results of 120 samples in Figure14(a), AWH and SWH calculated by the GNSS/SINS/SW method has a certain difference between the calculated results and reference value at some sampling points, but the errors are still within the acceptable range. In general, AWH and SWH calculated by this method perform well overall, as can be seen intuitively in Figure 14(b). Their errors can be controlled at about 0.1 meters, and in some better performances, the errors are so small that can almost be ignored in the engineering application of this paper. The results can be used to estimate the sea state well and provide an effective reference for the recycling work.


## 5. Conclusions

Before recycling the detection equipment on the sea surface, the staff needs to know the sea state of its surrounding environment to determine whether it is suitable for operation. Without additional sensors, the sea state can be estimated intuitively by wave height measurement with the SINS and GNSS of the equipment such as AUVs and underwater gliders.Under this application background, this paper proposes a wave correction method, which amends the results of the GNSS/SINS integrated method in the up direction by the virtual horizontal line. In addition, an integrated method of SINS and P-M spectrum is presented, which solves the problem of data divergence caused by GNSS interruption under bad sea conditions. The simulation results show that the problem is obviously improved, and the wave spectrum is of great significance in locating the sea-surface target.



---

*Source: 6647765-2021-01-18.xml*