基于SVR和LSTM算法反演散射计下有效波高研究

有效波高反演对于海洋工程及海洋环境安全具有重要意义.我国海洋二号(HY-2A)卫星载有散射计和高度计等获取海洋要素的仪器.散射计可获取海洋风场数据但无法直接获取有效波高数据,高度计可获取海洋有效波高数据但覆盖区域狭小.本文将散射计与高度计各自优势结合,利用支持向量回归(SVR)和长短期记忆(LSTM)智能算法反演散射计...     

A Spectral Approach for Determining Altimeter Wind Speed Model Functions

We propose a new analytical algorithm for the estimation of wind speeds from altimeter data using the mean square slope of the ocean surface, which is obtained by integration of a widely accepted wind-wave spectrum including the gravity-capillary wave range. It indicates that the normalized radar cross section depends not only on the wind speed but also on the wave age. The wave state effect on the altimeter radar return becomes remarkable with increasing wind speed and cannot be neglected at high wind speeds. A relationship between wave age and nondimensional wave height based on buoy observational data is applied to compute the wave age using the significant wave height of ocean waves, which could be simultaneously obtained from altimeter data. Comparison with actual data shows that this new algorithm produces more reliable wind speeds than do empirical algorithms. This revised version was published online in July 2006 with corrections to the Cover Date.     

An inverse method for reconstructing the density and sound speed profiles of a layered ocean bottom

A standard inverse problem in underwater acoustics is the reconstruction of the ocean subbottom structure (e.g., the density and sound speed profiles) from an aperture- and bandlimited knowledge of the reflection coefficient. In this paper we describe an inverse solution method due to Candel et al. [12] which is based on the scattering of acoustic plane waves by a one-dimensional inhomogeneous medium. As a consequence of applying the forward scattering approximation to a local wave representation of the acoustic field, they obtain an expression for the reflection coefficient in the form of a nonlinear Fourier transform of the logarithmic derivative of the local admittance. Inversion of this integral transform enables the recovery of the admittance profile via the numerical integration of two first-order differential equations which require as reflection data a single impulse response of the medium. Separate recovery of both the density and sound speed profiles requires two impulse responses for two different grazing angles. In this case, four differential equations need to be integrated instead of two. To illustrate the capability of the method, we present numerical reconstructions which are based on synthetic reflection data for a geoacoustic model that represents the acoustic properties of the surficial sediments for a site in the Hatteras Abyssal Plain.     

A semiempirical algorithm for SAR wave height retrieval and its validation using Envisat ASAR wave mode data

Many synthetic aperture radar(SAR) wave height retrieval algorithms have been developed.However,the wave height retrievals from most existing methods either depend on other input as the first guess or are restricted to the long wave regime.A semiempirical algorithm is presented,which has the objective to estimate the wave height from SAR imagery without any prior knowledge.The proposed novel algorithm was developed based on the theoretical SAR ocean wave imaging mechanism and the empirical relation between two types of wave period.The dependency of the proposed model on radar incident and wave direction was analyzed.For Envisat advanced synthetic aperture radar(ASAR) wave mode data,the model can be reduced to the simple form with two input parameters,i.e.,the cutoff wavelength and peak wavelength of ocean wave,which can be retrieved from SAR imagery without any prior knowledge of wind or wave.Using Envisat ASAR wave mode data and the collocated buoy measurements from NDBC,the semiempirical algorithm is validated and compared with the Envisat ASAR level 2 products.The root-mean-square-error(RMSE) and scatter index(SI) in respect to the in situ measurements are 0.52 m and 19% respectively.Validation results indicate that,for Envisat ASAR wave mode data,the proposed method works well.     

FY-3A/MERSI, ocean color algorithm, products and demonstrative applications

A medium resolution spectral imager (MERSI) on-board the first spacecraft of the second generation of China’s polar-orbit meteorological satellites FY-3A, is a MODIS-like sensor with 20 bands covering visible to thermal infrared spectral region. FY-3A/MERSI is capable of making continuous global observations, and ocean color application is one of its main targets. The objective is to provide information about the ocean color products of FY-3A/MERSI, including sensor calibration, ocean color algorithms, ocean color products validation and applications. Although there is a visible on-board calibration device, it cannot realize the on-board absolute radiometric calibration in the reflective solar bands. A multisite vicarious calibration method is developed, and used for monitoring the in-flight response change and providing post-launch calibration coefficients updating. FY-3A/MERSI ocean color products consist of the water-leaving reflectance retrieved from an atmospheric correction algorithm, a chlorophyll a concentration (CHL1) and a pigment concentration (PIG1) from global empirical models, the chlorophyll a concentration (CHL2), a total suspended mater concentration (TSM) and the absorption coefficient of CDOM and NAP (YS443) from China’s regional empirical models. The atmospheric correction algorithm based on lookup tables and ocean color components concentration estimation models are described. By comparison with in situ data, the FY-3A/MERSI ocean color products have been validated and preliminary results are presented. Some successful ocean color applications such as algae bloom monitoring and coastal suspended sediment variation have demonstrated the usefulness of FY-3A/MERSI ocean color products.     

"Principal component inverse" algorithm for detection in thepresence of reverberation

Detection in the presence of reverberation is often difficult in active sonar, due to the reflection/diffusion/diffraction of the transmitted signal by the ocean surface, ground, and volume. A modelization of reverberation is often used to improve detection because classical algorithms are inefficient. A commonly used reverberation model is colored and nonstationary noise. This model leads to elaborate detection algorithms which normalize and whiten reverberation. In this paper, we focus on a more deterministic model which considers reverberation as a sum of echoes issued from the transmitted signal. The Principal Component Inverse (PCI) algorithm is used with this model to estimate and delete the reverberation echoes. A rank analysis of the observation matrix shows that PCI is efficient in this configuration under some conditions, such as when the transmitted signal is Frequency Modulated. Both methods are validated with real sonar surface reverberation noise. We show that whitening has poor performance when reverberation and target echo have the same properties, while PCI maintains the same performance whatever the reverberation characteristics. Further, we extend the algorithms to spatio-temporal data. We propose a new algorithm for PCI which allows better echo separation. This new method is shown to be more efficient on real spatio-temporal data     

基于理论的气旋条件下双极化C波段合成孔径雷达波参数检测

Theoretical-based ocean wave retrieval algorithms are applied by inverting a synthetic aperture radar(SAR)intensity spectrum into a wave spectrum, that has been developed based on a SAR wave mapping mechanism. In our previous studies, it was shown that the wave retrieval algorithm, named the parameterized first-guess spectrum method(PFSM), works for C-band and X-band SAR at low to moderate sea states. In this work, we investigate the performance of the PFSM algorithm when it is applied for dual-polarization c-band sentinel-1(S-1) SAR acquired in extra wide-swath(EW) and interferometric wide-swath(IW) mode under cyclonic conditions.Strong winds are retrieved from six vertical-horizontal(VH) polarization S-1 SAR images using the c-band crosspolarization coupled-parameters ocean(C-3 PO) model and then wave parameters are obtained from the image at the vertical-vertical(VV) polarization channel. significant wave height(SWH) and mean wave period(MWP) are compared with simulations from the WAVEWATCH-III(WW3) model. The validation shows a 0.69 m root mean square error(RMSE) of SWH with a –0.01 m bias and a 0.62 s RMSE of MWP with a –0.17 s bias. Although the PFSM algorithm relies on a good quality SAR spectrum, this study confirms the applicability for wave retrieval from an S-1 SAR image. Moreover, it is found that the retrieved results have less accuracy on the right sector of cyclone eyes where swell directly affects strong wind-sea, while the PFSM algorithm works well on the left and rear sectors of cyclone eyes where the interaction of wind-sea and swell is relatively poor.     

Comparison of inversion algorithms for HF radar wave measurements

All ocean wave components contribute to the second-order scattering of a high-frequency (HF) radio wave by the sea surface. It is therefore theoretically possible to estimate the ocean wave spectrum from the radar backscatter. To extract the wave information, it is necessary to solve the nonlinear integral equation that describes the relationship between the backscatter spectrum and the ocean wave directional spectrum. Different inversion techniques have been developed for this problem by different researchers, but there is at present no accepted “best” method. This paper gives an assessment of the current status of two methods for deriving sea-state information from HF radar observations of the sea surface. The methods are applied to simulated data and to an experimental data set with sea-truth being provided by a directional wave buoy     

Measurement of ocean wave spectra using narrow-beam HE radar

A data interpretation algorithm is developed to extract ocean wave information from HF radar backscatter observed by a narrow-beam antenna system. The basis of this measurement is the inversion of the integral equation representing the second-order radar cross section of the ocean surface. This equation is numerically inverted by approximating it as a matrix equation and pseudoinverting the kernel matrix using a singular value decomposition. As a test of this algorithm, comparisons are made between wave spectrum estimates obtained from a WAVEC buoy and a pair of 25.4-MHz ground wave radars, using data collected during the 1986 Canadian Atlantic Storms Program (CASP). Overall, the results of this experiment have been positive and have demonstrated both the basic feasibility of the inversion algorithm and the wave sensing capability of HF radar. For example, significant wave height estimates deduced by two radars differed from the buoy, in an absolute value sense, by only 0.12 m on average. When using only one radar, the mean difference of this important parameter from the buoy was a reasonable 0.33 m     

A simplified pseudo-fluid model derived from biot theory through low grazing angle approximation

We propose an approximation technique with which the porous Biot model can be converted into a pseudo-fluid medium model, i.e., a medium represented by its sound speed and density. This technique begins from an analytic derivation of the reflection coefficient on a planar interface of fluid and porous ocean sediment. Invoking the low grazing angle approximation, useful for solving long range propagation problems, the pseudo-fluid medium is obtained. Alternate pseudo-fluid medium can be obtained through the weak frame approximation (Williams, J. Acoust. Soc. Am. 110, 2276-2281, 2001). In this paper, we discuss the accuracy and limitations of the low grazing angle approximation by numerically comparing the reflection coefficient to the full Biot model cases as functions of various Biot parameters, frequency, and water sound speed. The usefulness of the present low grazing angle approximated pseudo-fluid medium in ocean wave propagation modeling is demonstrated by comparing the transmission loss results with that of the full Biot model.     

天然气水合物赋存地层的地震波场模拟及对比

通过分析天然气水合物在海洋中的6种主要赋存状态类型,总结了每种赋存状态之间的相互转化关系及其物性参数计算方法,并应用到地震波场的正演模拟中.对比研究了声波模型、弹性波模型和双相介质模型对各种水合物赋存地层的响应特征,结果表明:1)当地层中存在孔隙充填型水合物且下伏地层不含游离气时,双相介质模拟的含水合物层底界表现负极性...     

Model-based identification: an adaptive approach to ocean-acousticprocessing

A model-based approach is developed to solve an adaptive ocean-acoustic signal-processing problem. Model-based signal processing is a well-defined methodology enabling the inclusion of propagation models, measurement models, and noise models into sophisticated processing algorithms. Here, we investigate the design of a so-called model-based identifier (MBID) for a general nonlinear state-space structure and apply it to a shallow water ocean-acoustic problem characterized by the normal-mode model. In this problem, we assume that the structure of the model is known and we show how this parameter-adaptive processor can be configured to jointly estimate both the modal functions and the horizontal wave numbers directly from the measured pressure-field and sound speed. We first design the model-based identifier using a model developed from a shallow-water ocean experiment and then apply it to a corresponding set of experimental data demonstrating the feasibility of this approach. It is also shown that one of the benefits of this adaptive approach is a solution to the so-called “mismatch” problem in matched-field processing (MFP)     

Static and dynamic modeling of a SAR imaged ocean scene

A number of models exist that attempt to explain wave imagery obtained with a synthetic aperture radar (SAR). These models are of two types; static models that depend on instantaneous surface features and dynamic models that employ surface velocities. Radar backscatter values (sigma_{0}) were calculated from 1.3- and 9.4-GHz SAR data collected off Marineland, FL. Thesigma_{0}data (averaged over many wave trains) collected at Marineland can best be modeled by the Bragg-Rice-Phillips model which is based on roughness variation and the complex dielectric constant of oceans. This result suggests that capillaries on the surface of oceanic waves are the primary cause for the surface return observed by a SAR. Salinity and temperature of the sea at small and medium incidence angles produce little effect upon sea-surface reflection coefficients atX-band, for either of the linear polarizations. The authors' observation of moving ocean, imaged by the SAR and studied in the SAR optical correlator, support a theory that the ocean surface appears relatively stationary in the absence of currents. The reflecting surface is most likely moving slowly (i.e., capillaries) relative to the phase velocity of the large gravity waves.     

A method for estimating the physical and acoustic properties of the sea bed using chirp sonar data

This paper proposes a method, based on the Biot model, for estimating the physical and acoustic properties of surficial ocean sediments from normal incidence reflection data acquired by a chirp sonar. The inversion method estimates sediment porosity from reflection coefficient measurements and, using the estimated porosity and the measured change in fast wave attenuation with frequency, estimates the permeability of the top sediment layer. The spectral ratio of echoes from the interface at the base of the upper sediment layer and from the sediment-water interface provides a measure of the change in attenuation with frequency. Given the porosity and permeability estimates, the Kozeny-Carman equation provides the mean grain size and the inversion method yields the acoustic properties of top sediment layer. The inversion technique is tested using chirp sonar data collected at the 1999 Sediment Acoustics Experiment (SAX-99) site. Remote estimates of porosity, grain size, and permeability agree with direct measurements of those properties.     

卫星海洋遥感数据的无损压缩及解压算法

在分析卫星海洋遥感单轨三级数据产品特征的基础上,针对游程编码无损压缩算法特征,提出了优化游程编码压缩算法,运用该算法实现了海洋遥感三级数据产品的无损压缩,并与流行的WinRAR压缩软件进行了比较测试,结果表明上述算法对海洋遥感数据均有较高的压缩率,优化游程编码压缩算法在空间复杂度及时间复杂度上均有明显的比较优势,可极大节约海洋遥感数据的存储与共享发布空间。     

Results from the Geoacoustic Inversion Techniques Workshop using modal inverse methods

This paper presents the results of a perturbative inverse approach applied to the range-dependent acoustic data provided as part of the Geoacoustic Inversion Techniques Workshop. The method is based on the Hankel transform relationship between complex-pressure field data measured on a horizontal array and the depth-dependent Green's function for a horizontally stratified medium. The input data to the inversion algorithm are discrete values of horizontal wavenumbers estimated for locally range-independent subapertures of the acoustic data. Inversion results are presented with emphasis placed on inverting for compressional wave speed as a function of both range and depth in the bottom.     

基于光流法的X波段雷达浪向反演研究

浪向代表着海浪传播方向,它是海上交通安全以及海岸资源管理的重要海洋环境参数之一。本文针对X波段测波雷达对海杂波的连续成像,提出了一种基于光流法的海浪传播方向反演新方法。该方法利用X波段测波雷达接收的海面回波图像序列直接进行光流运动估计,将得到的每个像素点的运动矢量进行加权平均,最后得到实际的海浪传播方向。与传统基于频域的X波段测波雷达浪向反演方法相比,本文提出的方法直接从时域来进行回波信号处理,无需提前得到调制传递函数以及精确的海流,减少了由于海流估算的不准确以及调制传递函数选取的误差而给雷达浪向反演带来的影响。同时,该方法简单高效,占用资源少,将其应用至仿真雷达回波以及现场实测数据来进行浪向反演,反演得到的浪向与仿真设定浪向值以及浮标实时观测浪向结果都有很好的吻合度,变化趋势也完全一致,进而验证了该方法的有效性以及准确性。     

Analysis and simulation of wave records through fast Hartley transform

The Hartley transform, a real-valued alternative to the complex Fourier transform, is presented as an efficient tool for the analysis and synthesis of ocean surface wave records. Basic theoretical properties of this real-valued transform are briefly reviewed. Similarities and differences between Fourier and Hartley integral transforms, as well as computational benefits and disadvantages between numerical algorithms used to evaluate their discrete versions, are presented. The fast Hartley transform algorithm is used to simulate stationary Gaussian time series of the sea surface elevation and to estimate the spectral density function, the Hilbert transform and the envelope function of wave records.     

Wave reflection from vertical breakwater with porous structure

This paper presents numerical solutions for the wave reflection from submerged porous structures in front of the impermeable vertical breakwater. A new time-dependent mild-slope equation involves the parameters of the porous medium including the porosity, the friction factor and the inertia coefficient, etc. is derived for solving the boundary value problem. A comprehensive comparison between the present model and the existing analytical solution for the case of simple rectangular geometries of the submerged structure is performed first. Then, more complicated cases such as the inclined and trapezoidal submerged porous structures in front of the vertical breakwater with sloping bottom are considered. This study also examines the effects of the permeable properties and the geometric configurations of the porous structure to the wave reflection. It is found that the submerged porous structure with trapezoidal shape has more efficiency to reduce the wave reflection than that of triangular shape. The numerical results show that the minimum wave reflection is occurred when the breakwater is located at the intermediate water depth.     

Detection of free gas and gas hydrate based on 3D seismic data and cone penetration testing: An example from the Nigerian Continental Slope

We present a new method to characterize free gas, gas hydrates and carbonate concretions occurrence which are considered as high-risk factors for sub-sea developments in the Niger delta. This method is based on the combination of 3D seismic data to the geotechnical site characterizations using piezocone CPTU tests (Cone Penetration Test with additional measurement of the pore water pressure). A special processing of the 3D seismic data has enabled the determination of the interval compressional velocity. Using the effective-medium theory, velocity anomalies (negative and positive) within the first 15 m were translated in gas hydrate and free gas distribution. The calibration of the P wave velocity anomalies was done thanks to in-situ geotechnical testing carried out during two oceanographic surveys (2003 and 2004). Comparison between in-situ testing, recovered cores and the prediction of the gas and the gas hydrate distribution based on the compressional wave velocity have shown that 3D seismic data is a valuable tool to identify heterogeneous areas but the use of the piezocone was essential to discriminate between gas hydrate occurrences and carbonate concretions' presence. Furthermore, in-situ compressional wave velocity (V_p) measurements have clearly demonstrated what it was suspected from the 3D seismic data, the co-existence in the study area between gas hydrate and free gas.