用一维水力学方程求取泻湖水位及潮汐汊道断面流速——以海南陵水新村港为例

本文通过一维水力学运动方程及连续性方程,得出一个反映泻湖水位变化的二阶非线性有阻尼的强迫振动方程,采用龙格-库特四步格式进行求解,并以海南陵水新村港为例,计算了泻湖水位及潮汐汊道的断面流速.研究结果表明:(1)一维水力学方程可容易地求得泻湖-潮汐汊道体系中泻湖内的水位及汊道的断面流速,计算简单快捷.(2)一维水力学方程可分析潮汐汊道体系中的一些基本现象.如潮汐汊道对潮波的"过滤器"特征;由于泻湖-潮汐汊道体系的阻尼系数较大所导致的潮波自外向泻湖传播时衰减剧烈;径流对泻湖水位及汊道流速产生有利于落潮优势的影响;风可使泻湖水位抬升或下降等.(3)无论是否考虑泻湖与潮汐汊道面积随水位的变化,一维水力学方程都具有相当大的缺陷.由于其不能反映泻湖内地形的空间变化,也不能解决漫滩与归槽等过程,在外海潮波无次生波的情况下,它本身不能产生次生潮波,因而不能正确地反映潮汐汊道体系的涨落潮历时与流速不对称,从而不能分析泥沙的净输运及潮汐汊道的稳定性问题.(4)新村港泻湖-潮汐汊道体系中潮汐汊道的断面面积大小对整个体系的稳定性及演变趋势起着控制作用.     

环肋椭圆柱壳-流场耦合系统自由振动的解析求解方法

提出了分析环肋椭圆柱壳-流场耦合系统自由振动的一种解析求解方法。环肋通过"刚度均摊"到均质椭圆柱壳体上,考虑流场对壳体振动特性的影响,通过Helmholtz方程描述声介质的波动,建立环肋椭圆柱壳-流场耦合系统的自由振动方程,进而求得耦合系统的自由振动固有频率。为了验证方法的准确性,将水下环肋椭圆柱壳退化成水下环肋圆柱壳和不加肋的水下椭圆柱光壳,两个退化模型的计算结果与已有文献结果吻合良好。并详细讨论了椭圆度、环肋间距等主要参数对耦合系统自由振动特性的影响。     

海底浅层介质切变波的初步研究

初步分析了中国沿海海底浅层介质切变波的传播特性,讨论了多种海底浅层介质切变波与压缩波、物理力学性质之间的关系,并与Hamilton和Chen等人的结果进行了比较。认为海底浅层介质存在切变模量,切变波可以被测量出来,切变波速为50—600m·s-1,测量频率在50—200kHz之间,用波速比可以评价海底浅层介质的应力 应变性质。研究工作有助于描述和评价中国沿海海底浅层介质声学物理性状,并为最终建立中国沿海海底地声模型提供依据。     

非均匀各向同性介质中弹性波传播的盒式积分法模拟

为方便积分差值法的精确处理,利用积分差值法,将运动方程、自由边界条件及人工边界上的吸收边界条件进行积分,根据高斯定理得出了二维非均匀介质中弹性波传播的二阶显式差分格式。并针对几个典型的地质模型,计算了其弹性波传播的数值模拟结果,结果表明,该方法有较好的精度和稳定性。     

半球形多向波聚合波道振荡水柱气室优化设计

为了解决振动水柱式波浪能转换装置收集多向波浪问题,本文设计了半球形多向聚合波道振荡水柱气室结构,以适合远海单点波浪能采集和发电。在规则波正向入射条件下,基于流体仿真分析软件(FLUENT)、流体动力学连续性假设和粘性不可压缩流体动量守恒的运动方程(Navier-Stokes方程)建立半球形振荡气室和三维数值波浪水槽模型。仿真结果表明:增设气室后壁,合理设计波道开口角度实现多向迎波捕获波浪能,优化前壁形状可降低波浪触底反射带来的能量耗散,同时提高了气室内空气压强和出气口速度,有效提升波浪能俘获效率,为后续发电的二次能量转换提供高效的空气动力。     

海底管线柔性导流板的变形方程研究

海底输油管线是海洋油气田开发的生命线工程,在复杂的海洋环境中易发生破坏.为使管线实现自埋防护,提出用柔性材料代替刚性导流板安装在管线顶部,不仅可以增加阻水面积,加大冲刷深度和范围,而且可以减缓刚性导流板对管线上部水流的扰动强度和尾流涡旋的扰动强度,减小管线的振动.利用力学关系由变形方程推导出海底管线与海底间距离为零以及不为零情况下柔性导流板变形的二次曲面方程和自由端最大挠度.在单向流条件下将不同材料不同长度的柔性导流板安装在管线中轴上方开展试验研究,测量柔性导流板在水流作用下的变形曲面.对比分析利用公式计算柔性导流板的自由端最大挠度和曲面变形数据与试验结果,发现两者吻合较好.     

局部微分求积法的深水包络孤立波数值模拟

利用局部微分求积法(LDQ)对非线性薛定谔(Schr dinger)方程进行数值求解,分别模拟了单深水孤立波运动,同向双深水孤立波追赶碰撞耦合运动,高阶孤立波振动和孤立波的反射与透射现象,得到各情况下的数值结果。从数值模拟及图像中揭示非线性薛定谔方程的性质和特点,阐述深水孤立波形成的物理意义、运动方式和运动规律,分析在不同初值条件下波形的变化特点,验证了LDQ法对该类问题的有效性。     

分层流水槽箱内垂荡板式内波造波模拟研究

提出了一种新型内波造波方法——箱内垂荡板式造波.应用FLUENT软件的有限体积法,采用k-ε两方程湍流模型封闭RANS方程作为求解控制方程,建立二维数值内波水槽.通过模仿物理造波机运动的数值模拟结果比较表明:箱内垂荡板式造波机生成的内波与理论解十分吻合,与双板联动式造波和平板拍击式造波相比,控制简单,对自由水面的扰动小...     

声波在水-多孔介质海底界面上的反射与透射

声波在海底界面的反射和透射是海底散射、海底混响、海底目标探测的重要问题。利用Biot多孔介质声传播理论对声波在水-多孔海底界面上的反射和透射进行了分析,具体给出反射声波的反射系数,3种透射声波的透射系数以及声能透射系数随入射波入射角和频率(10~40 kHz)的变化关系,分析了各种透射波对透射声能的贡献。多孔海底介质参数分别采用Stoll和Chotiros给出的2组参数进行理论计算。     

浅地层剖面的自由界面多次波预测与衰减

本文基于仅包含一次波和自由界面多次波的原始浅剖资料地震波场,推导出适应于单道数据的自由界面多次波预测方程,并提出基于同相轴追踪的多次波剔除与多道维纳滤波相结合的组合匹配衰减方法,实现了高精度的浅剖数据自由界面多次波预测与衰减,理论模型与实际浅地层剖面的处理实验均证明了该多次波预测与匹配衰减方法的有效性。     

Wave dispersion equation in a porous seabed

The wave dispersion equation has played a very important role in the development of ocean surface wave theories. The evaluation of the length of a water wave is an essential example of solving the dispersion relation. Conventional ocean wave theories have been based on an assumption of a rigid impermeable seabed. Thus, the conventional wave dispersion equation can only be used in the case of a wave propagating over a rigid impermeable seabed. For waves propagating over a porous seabed (such as a sandy bed), the conventional dispersion relation is no longer valid because of the absence of the characteristics of the porous seabed. The objective of this study is to establish a new wave dispersion equation for waves propagating over a porous seabed. Based on the new relation, the effects of a porous seabed on wave characteristics (such as the wavelength and wave profile) are discussed in detail.     

On the response of a Coulomb‐damped poroelastic bed to water waves

Abstract

The problem of forced vibration of a slightly inelastic porous bed by water waves is treated analytically on the basis of a linearized expression of the nonlinear damping term for the grain‐to‐grain friction in bed soils and the linear theory by Biot (1962a [Jour. Appl. Physics, 33:1482–1498]) on the elastic wave propagation in porous media. A dispersion relation of water waves is obtained as a function of wave frequency, water depth, permeability, Poisson's ratio, rigidity, and specific loss of bed soil. Three types of elastic waves are induced in a bed by water waves: a shear wave and a compressional wave in the skeletal frame of soil, and a compressional wave in the pore fluid. The compressional wave, due to the motion of the pore fluid relative to the skeletal frame of soil, is highly damped by the viscosity of pore fluid and only a short range effect near the boundaries of discontinuity, such as a sea‐seabed interface. The seabed response to water waves is characterized by the two Mach numbers, i.e., the ratio of water‐wave speed to shear‐wave speed in soil and the ratio of water‐wave speed to compressional‐wave speed in soil. Most of the water‐wave propagation problems fall into the subsonic flow condition, where elastic waves in the bed travel faster than water waves.

For sandy beds, generally the speeds of compressional and shear waves are much higher than the phase velocity of the water wave. For this case, the solution of the Coulomb‐damped poroelastic bed response presented in this paper approaches the solution of the massless poroelastic bed response in Yamamoto et al. (1978 [Jour. Fluid Mech., 87(1): 193–206]). The damping of water waves due to internal grain‐to‐grain friction is equally or more significant than the damping due to percolation in sand beds.

For clay beds, the speed of the shear wave in soil becomes low and comparable to the phase speed of the water wave. The bed motion for this case is considerably amplified due to the near‐resonance vibration of shear mode of bed vibration. The water wavelength on a clay bed is significantly shortened compared to the water wavelength over a rigid bed. The water wave damping due to internal grain‐to‐grain friction in soil becomes much larger compared to the water wave damping due to percolation in clay beds. Long water waves over a soft clayey bed attenuate within several wavelengths of travel distance.     

Research on Fatigue invalidation of A Submarine Pipeline Span Subjected to Waves and Currents

1 .IntroductionRecently, more and more importance is attached to the development of offshore oil and gas alloverthe world and many submarine pipelines have been laid.Submarine pipelines are the lifeline ofthe transportation of offshore oil and gas .Submar…     

Vortex-Induced Vibrations of A Free-Spanning Pipe Based on A Nonlinear Hysteretic Soil Model at the Shoulders

The pipe-soil interactions at shoulders can significantly affect the vortex-induced vibrations (VIV) of free-spanning pipes in the subsea. In this paper, the seabed soil reacting force on the pipe is directly calculated with a nonlinear hysteretic soil model. For the VIV in the middle span, a classic van der Pol wake oscillator is adopted. Based on the Euler-Bernoulli beam theory, the vibration equations of the pipe are obtained which are different in the middle span and at the two end shoulders. The static configuration of the pipe is firstly calculated and then the VIV is simulated.The present model is validated with the comparisons of VIV experiment, pipe-soil interaction experiment and the simulation results of VIV of free-spanning pipes in which the seabed soil is modelled with spring-dashpots. With the present model, the influence of seabed soil on the VIV of a free-spanning pipe is analyzed. The parametric studies show that when the seabed soil has a larger suction area, the pipe vibrates with smaller bending stresses and is safer.While with the increase of the shear strength of the seabed soil, the bending stresses increase and the pipe faces more danger.     

不规则波下船舶运动响应特征数值分析

针对不规则波浪作用下Wigley型船的运动响应问题进行了系统的研究,采用统计学方法深入探讨了船舶不规则运动幅值和响应周期的分布规律,并通过傅里叶变换对船舶运动响应进行了频谱特征分析。结果表明,船舶横摇方向与升沉和纵摇方向随机运动的响应特征有显著差异。在升沉与纵摇方向,波浪谱峰频率远离自振频率,前十分之一大振幅运动对应周期离散性较小,基本稳定在波浪谱峰周期附近,但小振幅运动周期分布离散性较大,频谱分析指出船舶升沉与纵摇运动响应频谱在波浪谱峰频率附近出现明显峰值。而在横摇方向,波浪谱峰频率与自振频率相耦合,不同振幅的横摇运动响应周期均稳定在自振周期附近,且周期离散性较小,频谱分析也表明横摇运动响应频谱主要集中于船舶运动自振频率附近。     

Characteristics of Vibrational Wave Propagation and Attenuation in Submarine Fluid-Filled Pipelines

As an important part of lifeline engineering in the development and utilization of marine resources, the submarine fluid-filled pipeline is a complex coupling system which is subjected to both internal and external flow fields. By utilizing Kennard’s shell equations and combining with Helmholtz equations of flow field, the coupling equations of submarine fluid-filled pipeline for n=0 axisymmetrical wave motion are set up. Analytical expressions of wave speed are obtained for both s=1 and s=2 waves, which correspond to a fluid-dominated wave and an axial shell wave, respectively. The numerical results for wave speed and wave attenuation are obtained and discussed subsequently. It shows that the frequency depends on phase velocity, and the attenuation of this mode depends strongly on material parameters of the pipe and the internal and the external fluid fields. The characteristics of PVC pipe are studied for a comparison. The effects of shell thickness/radius ratio and density of the contained fluid on the model are also discussed. The study provides a theoretical basis and helps to accurately predict the situation of submarine pipelines, which also has practical application prospect in the field of pipeline leakage detection.     

Cnoidal water wave induced seepage in a permeable seabed with a defined thickness

Extreme waves can induce seepage in a seabed and cause problems to marine structures in coastal regions. In this study, the seepage under cnoidal waves was studied using the transient seepage equation. An analytical solution is presented for the pore pressure in a seabed of defined thickness. Parametric studies were carried out to examine the influence of air content in the pore water, and of the soil hydraulic conductivity on the seepage. It has been shown that the air content and the soil hydraulic conductivity can affect the pore pressure response significantly. An increase in the air content or a decrease in the soil hydraulic conductivity will increase the magnitude of the pore pressure gradient and results in the pore pressure varying sharply. The liquefaction potential of a seabed under cnoidal waves is discussed. Consequently, comparative studies are carried out to show that the soil shear modulus and Poisson constant can influence the difference between the transient seepage equation and Biot's equation, and the transient seepage equation is a limit of Biot's equation.     

An integral description of the propagation of steady-state perturbations of the heavy liquid surface

The system of equations of motion describing the gravity wave propagation in a perfect heavy liquid layer is transformed into a new integral equation for the free surface elevations. In the limit cases, this integral equation describes the linear and nonlinear periodic waves as well as the known types of solitary waves. In this case a dispersion equation arises because perturbations of the second and higher orders of smallness are neglected. The integral equation allows for the propagation of invariable surface perturbations of arbitrary forms if their spatial spectrum is concentrated near small wave numbers (compared to the inverse wave amplitude). Several examples of solutions are presented.     

Experimental Study on Vortex-Induced Vibrations of Submarine Pipeline near Seabed Boundary in Ocean Currents

1 .IntroductionThe submarine pipeline is a commonfacility widely usedfor offshore oil and gastransport . Whena pipeline is installed on a seabed and not buried,unsupportedspans may exist insomelocations ,es-peciallyinthe uneven zones of the seabed.The spa…     

渗透海床上矩形Bragg防波堤对波浪反射的研究

与海床不可渗透的情况相比,波浪在可渗透海床上传播时会发生波能衰减。本文将基于可渗透海床上一维修正型缓坡方程,建立方程求解的有限差分模型。将通过与不可渗透海床上矩形Bragg防波堤对波浪反射系数解析解的对比,验证有限差分模型的正确性和适用性。将进一步研究海床可渗透情况下,海床的渗透性参数、坝体的相对宽度、数量、浸没度对波浪反射系数的影响及其与海床不可渗透情况下的差异。本文研究发现,Bragg共振发生时的反射系数随坝体数量的增多而增大,随海床渗透性参数和坝体浸没度的增大而减小,并且存在一个坝体相对宽度值会使Bragg共振反射达到最大。相较于海床不可渗透的情况,发生Bragg共振反射的波浪频率几乎相同,但反射系数减小,而且零反射（或全透射）现象不再存在。