基于LIC的海流动态模拟与可视化

线积分卷积(LIC)方法是一种针对矢量场数据的可视化技术,本文主要研究了LIC方法应用于海洋流场的可视化质量改进问题。主要包括以下几点改进:根据海流运动方向改变稀疏背景噪声的指向以使流线变得更加细长明亮;针对稀疏噪声中海流变化剧烈的区域流线会出现分叉的问题,采用忽略稀疏背景噪声流线追踪时非主流线部分的纹理点的积分计算方法来解决此问题,同时可以省去不必要纹理点的积分计算的开支,加快显示速度;采用赋予海流随机的初始相位的方法来解决稀疏噪声循环动画突兀的问题;通过流线的运动速度快慢或者流线的颜色来表示海流速度大小;通过对海流数据进行插值并改善边界流线追踪条件的方法可形成完整细致的全球海流图像。实验结果表明,该方法可以生成纹理对比强烈、细节清晰、动画连贯的高质量海流运动图像,具有良好的可视化效果。     

层化海洋中的二阶波散射作用问题

引入Stokes有限振幅水波的多色波作用机制 ,采用二层层化海洋模型 ,分析了层化海洋中非线性波对结构的散射作用问题。提出了二层海洋中二阶水波散射条件的数学表达式 ,进而推导了二层海洋中二阶波对一般型圆柱散射作用的格林积分解式 ,并给出了对解式中自由面及分层界面无限积分的递推算式。     

Effects of Second Order Scattered Wave in Stratified Ocean

By using the theories on Stokes multicolored water waves and taking the two- layer ocean as a basic model of stratified ocean, the paper analyzes the problems related to the effects of the nonlinear water wave on offshore structures. A mathematical expression is presented to describe second order wave radiation conditions. Using integral principle, the analytical integral solutions are given to evaluate second order scattered wave loads on general vertical circular cylinders in the two-layer ocean, and the special recurrence formulas for infinite integrals over free and stratified surfaces are derived.     

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     

Analysis of the response of a non-linear data assimilation algorithm in a one-dimensional model of the upper ocean

A one-dimensional (1D) integral dynamico-stochastic model of the upper ocean with a non-linear assimilation algorithm is considered. The accuracy of computing the characteristics of the upper layer depends essentially on the values of the empirical coefficients. A numerical experiment was carried out which verified the efficiency of the model's adaptive mechanism operation when different values of the empirical coefficients and their variances were preset. Recommendations on assigment of the model's initial parameters are derived.Translated by Mikhail M. Trufanov.     

西北太平洋浪流相互作用对有效波高的影响研究

西北太平洋强流区会对海浪的特征和分布产生显著的影响,尤其是研究台风过程中海流与海浪的相互作用具有重要的研究意义。本文以ROMS海洋模式和SWAN海浪模式为基础,构建了浪流耦合模式系统,对2013年10月6-17日间的台风“丹娜丝”、“百合”、“韦帕”过程中西北太平洋浪流相互作用中海流对有效波高的影响进行了研究。通过对比模式模拟有效波高与浮标观测资料,发现耦合后的有效波高比非耦合结果更接近观测值,耦合模式中海流的存在对有效波高的分布有明显的影响。研究表明,特别是在有效波高峰值处,海流引起的有效波高增大最大可达1 m。海浪浪向及流向的空间分布以及中国近海浮标处浪向与流向的时间序列表明,流向与浪向反向时,海流的影响造成有效波高增大;二者同向时,有效波高减小。海流对有效波高的调整会沿着海浪传播的方向传播相当一段距离。在西北太平洋的海浪场计算中,引入海流的耦合模式计算结果对改善强流区海浪预报具有重要意义,并且海流的模拟精度对于高精度的海浪预报非常重要。     

底部呈二维缓变情况下超大型浮体水弹性响应的两种计算方法比较

超大型浮体在海洋资源开发和海洋空间利用方面有重要应用前景.非均匀海洋环境中的水弹性响应是其应用中的一个重要问题.在近海中最典型的非均匀海洋环境当属由于底部变化引起的非均匀现象.本文分别采用多重尺度法（零阶近似）和常规的有限水深势流格林函数边界积分法,对底部呈二维缓变情况下超大型浮体的水弹性响应问题进行了研究和对比,并与实验工况进行了对照.两种方法与试验结果吻合较好,证明非均匀海洋环境确实对超大型浮体的水弹性响应具有一定的影响.     

Extraction of ocean wave parameters from HF backscatter received bya four-element array: analysis and application

An algorithm that would extend the capabilities of a four-element square array known as the Coastal Oceans Dynamics Applications Radar (CODAR) to include the yielding of directional wave-height spectra from backscattered radiation is addressed. General expressions for the first- and second-order broadbeam radar cross-sections of the ocean surface are applied to the array. A Fourier-basis-function approach allows the broadbeam cross-sections to be written as a system of integral equations. The second-order radar return involves a double integral whose integrand contains nonlinear combinations of the unknowns, namely, the Fourier coefficients of the ocean wave directional spectrum. The first-order portion of the radar spectrum is used to linearize this integral. The matrix system then formulated is solved using a singular value decomposition (SVD) approach, and the resulting ocean spectral coefficients are used to give the directional spectrum. Test results for the algorithm are reported and discussed     

基于海面多普勒谱模型的海浪反演精度分析

Microwave remote sensing is one of the most useful methods for observing the ocean parameters. The Doppler frequency or interferometric phase of the radar echoes can be used for an ocean surface current speed retrieval,which is widely used in spaceborne and airborne radars. While the effect of the ocean currents and waves is interactional. It is impossible to retrieve the ocean surface current speed from Doppler frequency shift directly. In order to study the relationship between the ocean surface current speed and the Doppler frequency shift, a numerical ocean surface Doppler spectrum model is established and validated with a reference. The input parameters of ocean Doppler spectrum include an ocean wave elevation model, a directional distribution function, and wind speed and direction. The suitable ocean wave elevation spectrum and the directional distribution function are selected by comparing the ocean Doppler spectrum in C band with an empirical geophysical model function(CDOP). What is more, the error sensitivities of ocean surface current speed to the wind speed and direction are analyzed. All these simulations are in Ku band. The simulation results show that the ocean surface current speed error is sensitive to the wind speed and direction errors. With VV polarization, the ocean surface current speed error is about 0.15 m/s when the wind speed error is 2 m/s, and the ocean surface current speed error is smaller than 0.3 m/s when the wind direction error is within 20° in the cross wind direction.     

Radar cross section analysis of marine targets using a combining method of physical optics/geometric optics and a Monte-Carlo simulation

In this paper, the radar cross section of flat plates on ocean surfaces is statistically investigated. A combining method of physical optics and geometric optics is applied to establish an effective backscattering analysis procedure. This method is a high-frequency analysis method originally derived from a simplified Stratton-Chu integral equation, assuming that the radar is far away from the target so that Kirchhoff approximation is valid. A Monte-Carlo simulation method is adopted to statistically analyze the effects of undulated ocean surfaces. The ocean surfaces are randomly generated by Pierson-Moskowitz ocean wave spectrum and a directional distribution function. Numerical investigations are carried out for flat plates, with the same height and width but with different inclined angles, on ocean surfaces of various significant wave heights.     

岛屿岛礁海域海浪能谱模型研究进展

波浪能谱模型在岛屿岛礁海域的波浪预报研究和海洋工程中应用广泛,但存在模式计算格点无法充分体现岛屿岛礁的复杂地形特征和很难刻画波浪受到岛屿岛礁影响发生变形物理过程等两个关键问题。多重网格嵌套方案、岛屿次网格地形效应计算方案以及非结构网格、无网格、动态自适应四叉树网格等技术在体现岛屿岛礁复杂地形方面取得了较好的效果;将相位解析模型与波浪能谱模型优势互补是提高能谱模型对岛屿近岸波浪变形物理过程计算能力的一个有效方法。开展球坐标系下波作用密度谱方程的自适应四叉树网格求解方法研究,借鉴相位解析模型最新成果完善能谱模式的绕射、反射、底摩擦等物理过程,是提高岛屿岛礁海域海浪精细预报技术水平的前沿性、探索性研究方向。     

On the problem of heat balance in integral models of the upper ocean

The heat balance equations for the mixed layer, the jump layer modelled by the temperature discontinuity surface, the heat flux, current velocity, and the thermocline, where the temperature is parameterized vertically, are derived in the spatial unsteady hydrodynamic integral model of the upper ocean.Translated by Mikhail M. Trufanov.     

Radial velocity of ocean surface current estimated from SAR Doppler frequency measurements—a case study of Kuroshio in the East China Sea

Ocean currents are a key element in ocean processes and in meteorology, affecting material transport and modulating climate change patterns. The Doppler frequency shift information of the synthetic aperture radar (SAR) echo signal can reflect the dynamic characteristics of the sea surface, and has become an essential sea surface dynamic remote sensing parameter. Studies have verified that the instantaneous Doppler frequency shift can realize the SAR detection of the sea surface current. However, the validation of SAR-derived ocean current data and a thorough analysis of the errors associated with them remain lacking. In this study, we derive high spatial resolution flow measurements for the Kuroshio in the East China Sea from SAR data using a theoretical model of shifts in Doppler frequency driven by ocean surface current. Global ocean multi observation (MOB) products and global surface Lagrangian drifter (GLD) data are used to validate the Kuroshio flow retrieved from the SAR data. Results show that the central flow velocity for the Kuroshio derived from the SAR is 0.4–1.5 m/s. The error distribution between SAR ocean currents and MOB products is an approximate standard normal distribution, with the 90％ confidence interval concentrated between –0.1 m/s and 0.1 m/s. Comparative analysis of SAR ocean current and GLD products, the correlation coefficient is 0.803, which shows to be significant at a confidence level of 99％. The cross-validation of different ocean current dataset illustrate that the SAR radial current captures the positions and dynamics of the Kuroshio central flow and the Kuroshio Counter Current, and has the capability to monitor current velocity over a wide range of values.     

Numerical modeling of currents in the Drake Passage with assimilation of the experimental data of 2003

Computation of ocean currents in the Drake Passage is carried out with variational assimilation of the data of the hydrographical section across the Drake Passage carried out on December 11–15, 2003, and other data. A stream-eddy structure of the easterly Antarctic Circumpolar Current and a westerly current on the Antarctic Slope are obtained. Water transports by the different current field components and the integral transport across the Drake Passage are estimated. The necessity of direct current measurements in the Drake Passage for a correct estimate of the transport is confirmed.     

HF radar ocean current algorithm based on MUSIC and the validation experiments

Wuhan University's ocean state measuring and analyzing radar (OSMAR2000), working at around 7.5 MHz in the low region of the HF band with a 120-m-long linear receiving antenna array, can measure ocean surface current at ranges of up to 200 km. An ocean surface current algorithm based on direction finding (DF) using the multiple signal classification (MUSIC) method is developed for the OSMAR2000 radar. This paper describes the OSMAR2000 ocean surface current algorithm based on MUSIC and the validation experiments in the East China Sea. The results of the ocean surface current measurements demonstrate that the OSMAR2000 ocean surface current algorithm based on MUSIC is feasible for the long range of ocean surface current mapping with a sufficient bearing resolution.     

Oceanic sound speed profile inversion

A modal (full-wave) method has been developed to predict ocean sound speed profiles from propagated acoustic field data. The method assumes a point source of sound in the ocean and uses as data the values of the transmitted acoustic field at an array. The formalism for depth-dependent sound speeds consists of the standard Hankel integral transform of the depth solution. In the travel length coordinate, the latter is written exactly, using the Green's function, in terms of an integral equation whose kernel includes the sound speed profile correction. A Born approximation to this equation is used. This is just the WKB solution, and permits the use of a nontrivial input (or guess) profile, here chosen as bilinear. The use of asymptotic methods enables us to write the data as an integral transform over the profile correction. The transform can be inverted. An example is presented for full-bandwidth inversion.     

层化海洋中二阶多色波对可渗透结构的绕射问题

可渗透结构具有使波浪作用减弱的效应,而海水的层化及水波的非线性使结构的波绕射产生多层复杂机制。将可渗透结构应用于复杂海况条件中,海水的层化性、波浪的非线性及结构的透空性构成了波绕射的一个十分复杂的数学问题。该问题存在理论研究的必要性,而文章则着重探讨其数学分析的可能性。通过引入二层海的层化海模式及Stokes二阶波的非线性波模式,给出了二阶多色波对透空结构的波绕射的定解问题提法,提出了复合形式的二阶多色波辐射条件式及可渗透结构的二阶物面条件式,应用特征函数解法与积分法推导了多色波对结构绕射的一阶势解与二阶作用的耦合积分解式,并讨论了解式所涉及无穷积分的算法。     

Extraction of Ocean Wave Spectra From Simulated Noisy Bistatic High-Frequency Radar Data

An algorithm is developed for the inversion of bistatic high-frequency (HF) radar sea echo to give the nondirectional wave spectrum. The bistatic HF radar second-order cross section of patch scattering, consisting of a combination of four Fredholm-type integral equations, contains a nonlinear product of ocean wave directional spectrum factors. The energy inside the first-order cross section is used to normalize this integrand. The unknown ocean wave spectrum is represented by a truncated Fourier series. The integral equation is then converted to a matrix equation and a singular value decomposition (SVD) method is invoked to pseudoinvert the kernel matrix. The new algorithm is verified with simulated radar Doppler spectrum for varying water depths, wind velocities, and radar operating frequencies. To make the simulation more realistic, zero-mean Gaussian noise from external sources is also taken into account     

A Geometrically Exact Formulation for Three-Dimensional Numerical Simulation of the Umbilical Cable in A Deep-Sea ROV System

This paper proposes a geometrically exact formulation for three-dimensional static and dynamic analyses of the umbilical cable in a deep-sea remotely operated vehicle (ROV) system. The presented formulation takes account of the geometric nonlinearities of large displacement, effects of axial load and bending stiffness for modeling of slack cables. The resulting nonlinear second-order governing equations are discretized spatially by the finite element method and solved temporally by the generalized- implicit time integration algorithm, which is adapted to the case of varying coefficient matrices. The ability to consider three-dimensional union action of ocean current and ship heave motion upon the umbilical cable is the key feature of this analysis. The presented formulation is firstly validated, and then three numerical examples for the umbilical cable in a deep-sea ROV system are demonstrated and discussed, including the steady configurations only under the action of depth-dependent ocean current, the dynamic responses in the case of the only ship heave motion, and in the case of the combined action of the ship heave motion and ocean current.     

A Geometrically Exact Formulation for Three-Dimensional Numerical Simulation of the Umbilical Cable in A Deep-Sea ROV System

This paper proposes a geometrically exact formulation for three-dimensional static and dynamic analyses of the umbilical cable in a deep-sea remotely operated vehicle(ROV) system. The presented formulation takes account of the geometric nonlinearities of large displacement, effects of axial load and bending stiffness for modeling of slack cables. The resulting nonlinear second-order governing equations are discretized spatially by the finite element method and solved temporally by the generalized-a implicit time integration algorithm, which is adapted to the case of varying coefficient matrices. The ability to consider three-dimensional union action of ocean current and ship heave motion upon the umbilical cable is the key feature of this analysis. The presented formulation is firstly validated, and then three numerical examples for the umbilical cable in a deep-sea ROV system are demonstrated and discussed, including the steady configurations only under the action of depth-dependent ocean current, the dynamic responses in the case of the only ship heave motion, and in the case of the combined action of the ship heave motion and ocean current.