浅海环境下溢油海面的仿真

在考虑水深因素的条件下,提出了一种适合于浅海环境下溢油海面的仿真方法,利用TMA谱模型和Marangoni溢油理论模型,计算有限水深下溢油海面的海浪谱;然后根据海浪的色散关系和Longuet-Higgins海浪模型,计算有限水深下溢油海面的铅直位移;并分析了Kitaigorodskii深度函数的特性。仿真结果表明,有限水深下海面的铅直位移比深水的小,而有限水深下溢油海面的粗糙度要比有限水深下和深水清洁海面的都小,结果与Marangoni溢油理论模型相吻合。     

GPS浮标数据反演海浪谱的理论仿真与试验验证

GPS浮标作为一种新型的海洋测量设备,近年来在海面高度现场测量和星载高度计定标方面取得了重要应用。通过仿真试验对反演海浪谱的方法和流程进行研究,旨在探索从GPS浮标测量的海面高度序列中提取海浪谱的方法。首先,使用Longuest模型生成了海浪波面位移时间序列,并通过Pierson-Moscowitz风浪谱对波面位移的统计特性进行约束,其随机性由相位引入。结合典型潮汐和GPS浮标仪器噪声的仿真时间序列,合成了仿真时间长度1h的1Hz(每秒1次)随机海面高度序列。然后,利用自相关函数法,进行高通滤波和数据压缩,得到了仿真的海浪谱。该仿真结果和理论海浪谱非常接近,可满足海浪谱反演的需求。最后,通过山东石岛外海的GPS浮标现场试验,验证了本文提出的反演方法的适用性。本文的研究解决了GPS浮标反演海浪谱的关键问题,丰富了海浪谱反演的手段,拓展了GPS浮标的应用领域,有望为未来我国的星载波谱仪定标服务。     

波面位移非线性特征数值研究

基于Longuest-Higgins(1963)非线性海浪模型,在有限水深且存在均匀背景流的条件下,根据Song(2006)给出的波面位移二阶表达式,采用Combi海浪频谱计算了海表面定点波面位移时间序列和波面位移概率统计分布。分析了波面位移统计分布随风速、水深、反波龄和均匀背景流的变化特征和规律以及不同海况条件下二阶非线性项对波面位移统计分布的影响。结果表明:二阶非线性项使波面位移分布偏离正态分布,二阶非线性作用受风速、水深、反波龄和均匀背景流的影响。风速增大、水深降低、反波龄减小或者均匀背景流和风速传播方向相反均使波面位移二阶非线性项的作用加强,无因次波面位移概率密度分布的偏度和峰度随之增大,反之则二阶非线性项作用减弱。当均匀背景流和风速相同时,虽然使非线性项的作用减弱,但平均波面位移反而比静止水平面降低。当均匀背景流和风速相反时,虽然使非线性作用增强,但平均波面位移反而趋于静止水平面。得到如下结论:二阶非线性项对于波面位移有显著影响,数值模拟波面位移需要增加二阶非线性项。通过以上研究,提高了数值模拟波面位移的准确性,而波面位移是海浪最基本的特征量,从而增强了海浪模拟和预报的准确性,对海洋工程、海–气相互作用、上层海洋动力学等具有重要意义。     

海浪波面极大值分布的非线性影响

依非线性海浪模型,在二阶近似下,利用特征函数展开技术和直接求矩的方法,导出了定点波面位移及其对时间的一阶和二阶导数的联合分布。由此联合分布,导出了二阶非线性近似下的波面极大值分布。它由线性意义下的海浪频谱及所考虑的二阶级。波相互作用所确定。当忽略非线性影响时,文中给出的波面极大值分布退化为Cartwright和Longuet-Higgins所导出的分布。     

基于GPS高频数据的海浪测量方法

研究了基于GPS高频数据进行海浪测量的方法。利用TRACK解算得海上载体高精度的垂向位移,经浪潮分离提取海浪信号,分别采用相关函数法和周期图法估计海浪信号的功率谱,并计算海浪要素。利用实测数据进行试验,结果表明,周期图法推算的平均波高和平均周期精度较高,与测波仪结果差异分别小于2cm和0.25s,基于GPS高频数据的海浪测量方法可有效反演海浪要素。     

基于谱分析的GPS浮标测量参数提取方法

介绍采用谱分析原理提取海浪信息的方法。根据GPS浮标测量得到的海浪高度数据,先通过滑动平均的方法分离海浪信息和潮流潮位信息,然后对海浪信息采用AR模型法进行功率谱估计,并对海浪信息进行连续小波变换,对AR模型法功率谱估计结果和小波变换结果进行综合研究,提取海浪特征参数———周期。最后通过一个GPS浮标试验对上述方法进行了验证。     

基于三阶斯托克斯波的谱估计及海浪数值模拟研究

为了拓宽建立深水非线性海浪频谱模型的研究途径,基于三阶斯托克斯波理论,提出了以改进的三阶斯托克斯波为组成波,以修正的自回归模型(AR)谱估计方法得到的估计谱(简称新谱)为靶谱的非线性海浪模型。通过基于模拟频谱和相位谱反演海浪波面高度时间序列(以下简称海浪时历反演方法)的谱比较方法验证了新谱的优越性,为从根本上解决精确评估海浪谱提供了技术支撑,对于增强未来掌控海战场具有重要意义。     

一种实体嵌入的三维海浪可视化方法

针对实体嵌入对海面的随机干扰问题,提出一种海浪三维可视化方法。该方法采用扰动的思想改进了海浪的数学模型,并利用信息更加全面的方向谱描述海浪,进而将方向谱和改进海浪模型共同应用于海浪的三维可视化,形象展示实体嵌入时对海浪的随机性影响,并给出了三维实体的建模方法。通过多组实体嵌入的海浪可视化仿真试验的比较与分析,验证了本文提出方法的有效性。     

数值模式与统计模型相耦合的近岸海浪预报方法

针对数值模式和统计模型预报近岸海浪存在的局限性,构建了数值模式和统计模型相耦合的近岸海浪预报框架,在模式计算格点和近岸预报目标点之间定义一个海浪能量密度谱传递系数,通过经验正交函数分解和卡尔曼滤波方法建立传递系数的统计预报模型并与数值模式进行耦合。经过对近岸波浪观测站1a的预报试验表明:该方法能够提高近岸海浪有效波高预报精度,有效波高的均方根误差降低了约0.16m,平均相对误差降低约9%。进一步试验和分析发现,该方法的预报有效时间小于24h,将海浪能量密度谱经过分解后得到的基本模态反映了近岸波侯的主要特征,海浪能量密度谱传递系数的变化体现了波侯的季节变化特点。     

用MC算法优化基于SPH的海浪粒子模型

结合风浪场复杂的动力学、时空特性以及对目前海浪仿真方法的研究总结,针对海浪仿真模拟中存在的真实感不足,难以模拟海浪破碎等问题,根据流体的物理特性,采用光滑粒子流体动力学(SPH)的方法实现了海浪的动态演变仿真。然后针对基于SPH方法的海浪粒子模型存在的离散性和真实感较差的问题,采用一种基于移动立方体法(MC)的海浪运动的自由表面抽取算法,完成了海浪场的表面建模,优化了基于粒子系统的海浪模拟效果。实验表明,用MC算法优化的方法既满足了海浪动态演变的仿真要求,又使水体模型更加连贯,更接近真实水体效果,优化效果较为明显,是综合利用体绘制与面绘制法实现三维可视化模拟的一次成功的探索。     

波谱仪调制功率谱估计算法性能的仿真研究

基于星载波谱仪海浪方向谱探测原理,仿真了不同海况、风速下的海浪波谱仪接收信号,并采用周期图法、Welch法、AR模型法以及最小方差法共4种不同的调制谱估计方法反演出海浪谱,比较各种调制谱估计方法的海浪方向谱反演性能。仿真结果表明:对于一定方位向下的一维海浪谱反演,不同调制谱估计方法反演海浪谱性能优劣没有绝对的顺序。对于二维海浪谱反演,在成长中海浪条件下周期图法反演性能最差,其他3种方法的反演性能没有绝对优劣顺序;对于成熟风浪,最小方差法在积分能量误差、有效波高误差两个指标上的反演性能最好,在主波波向、主波波长误差上,周期图法反演性能最差,其他3种方法没有绝对优劣顺序。在涌浪条件下,AR模型法反演性能优于其他3种方法。在不同海况下,随着波谱仪中心入射角的下降,反演性能会下降。基于这些仿真结果,本文推荐最小方差法作为充分成长海浪的海浪方向谱反演的调制谱估计方法,AR模型法作为涌浪海浪方向谱反演的调制谱估计方法。     

球坐标系下MASNUM海浪数值模式的建立及其应用

为开展海浪对海洋上混合层的搅拌混合作用及其对海气界面通量的影响等研究,在LAGFD WAM区域海浪数值模式基础上建立了球坐标系下的全球海浪数值模式.重点导出了球坐标系下的海浪能量谱平衡方程及其复杂特征线方程,该组方程包含了背景流场对波动传播的调整、波动沿大圆传播的折射等.数值积分则采用复杂特征线嵌入计算格式.初步数值模拟结果表明,该海浪全球数值模式能够较为精确地刻画海浪的动力过程.     

Nonlinear Dynamic Analysis of Deepwater Drilling Risers Subjected to Random Loads

Excited by ocean currents, random wave and vessel motion, deepwater drilling risers exhibit significant dynamic response. In time domain, a method is proposed to calculate the nonlinear dynamic response of deepwater drilling risers subjected to random wave and dynamic large displacement vessel motion boundary condition. Structural and functional loads, external and internal pressure, free surface effect of irregular wave, hydrodynamic forces induced by current and wave, as well as wave and low frequency (drift) motion of the drilling vessel are all accounted for. An example is presented which illustrates the application of the proposed method. The study shows that long term drift motion of the vessel has profound effect on the envelopes of bending stress and lateral displacement, as well as the range of lower flex joint angle of the deepwater riser. It can also be concluded that vessel motion is the principal dynamic loading of nonlinear dynamic response for the deepwater risers rather than wave force.     

表面波浪对海洋上层混合的作用

A new three-dimensional numerical model is derived through a wave average on the primitive N-S equations, in which both the＂Coriolis-Stokes forcing＂ and the＂Stokes-Vortex force＂ are considered. Three ideal experiments are run using the new model applied to the Princeton ocean model （POM）. Numerical results show that surface waves play an important role on the mixing of the upper ocean. The mixed layer is enhanced when wave effect is considered in conjunction with small Langmuir numbers. Both surface wave breaking and Stokes production can strengthen the turbulent mixing near the surface. However, the influence of wave breaking is limited to a thin layer, but Stokes drift can affect the whole mixed layer. Furthermore, the vertical mixing coefficients clearly rise in the mixed layer, and the upper ocean mixed layer is deepened especially in the Antarctic Circumpolar Current when the model is applied to global simulations. It indicates that the surface gravity waves are indispensable in enhancing the mixing in the upper ocean, and should be accounted for in ocean general circulation models.     

Polarimetric entropy of the ocean surface with a two-scale scattering model

The relationships among an ocean wave spectrum,a fully polarimetric coherence matrix,and radar parameters are deduced with an electromagnetic wave theory.Furthermore,the relationship between the polarimetric entropy and ocean wave spectrum is established based on the definition of entropy and a twoscale scattering model of the ocean surface.It is the first time that the polarimetric entropy of the ocean surface is presented in theory.Meanwhile,the relationships among the fully polarimetric entropy and the parameters related to radar and ocean are discussed.The study is the basis of further monitoring targets on the ocean surface and deriving oceanic information with the entropy from the ocean surface.The contrast enhancement between human-made targets and the ocean surface with the entropy is presented with quad-pol airborne synthetic aperture radar（AIRSAR） data.     

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

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.     

一种改进的随机起伏海面的仿真方法

根据线性滤波器模型的基本原理,提出了一种改进的基于高阶滤波器的随机起伏海面仿真方法.利用Pierson的半经验海浪谱模型,得到了各种不同风速条件下海面随机起伏的仿真结果,并且对仿真结果进行了比较和讨论,结果表明,利用高阶滤波器方法能更加准确地仿真不同风速下的实际海面、更好地再现随机起伏海面的海浪谱特征.这一方法的提出为海浪的物理模拟实验提供了更为准确的数值计算方法,也为研究目标与海面共存条件下高频地波雷达信号的回波特性提供了更好的数值仿真平台.     

A spectral barotropic model of the wind-driven world ocean

A global spectral barotropic ocean model is introduced to describe the depth-averaged flow. The equations are based on vorticity and divergence (instead of horizontal momentum); continents exert a nearly infinite drag on the fluid. The coding follows that of spectral atmospheric general circulation models using triangular truncation and implicit time integration to provide a first step for seamless coupling to spectral atmospheric global circulation models and an efficient method for filtering of ocean wave dynamics. Five experiments demonstrate the model performance: (i) Bounded by an idealized basin geometry and driven by a zonally uniform wind stress, the ocean circulation shows close similarity with Munk’s analytical solution. (ii) With a real land–sea mask the model is capable of reproducing the spin-up, location and magnitudes of depth-averaged barotropic ocean currents. (iii) The ocean wave-dynamics of equatorial waves, excited by a height perturbation at the equator, shows wave dispersion and reflection at eastern and western coastal boundaries. (iv) The model reproduces propagation times of observed surface gravity waves in the Pacific with real bathymetry. (v) Advection of tracers can be simulated reasonably by the spectral method or a semi-Langrangian transport scheme. This spectral barotropic model may serve as a first step towards an intermediate complexity spectral atmosphere–ocean model for studying atmosphere–ocean interactions in idealized setups and long term climate variability beyond millennia.