双相各向异性介质中弹性波方程的有限元解法及波场模拟

基于双相各向异性介质模型,首先推导了双相各向异性介质中弹性波传播的动力学方程及其Galerkin变分方程和有限元运动方程,然后给出了孔隙弹性波方程的有限元数值解法以及二维双相PTL介质中波场模拟的人为吸收边界条件. 最后,利用本文给出的有限元方法对双相PTL介质和双相各向同性介质中的弹性波传播进行了数值模拟. 结果表明:有限元方法和吸收边界条件有效、可行,在理想相界条件下,不论是从固体位移,还是从流体位移的波场快照都能看到明显的慢速拟P波;在黏滞相界情况下,能否观察到慢速拟P波,与含流体地层介质的耗散性质有关.对实际含流体介质,从流体位移分量的波场快照比从固体位移波场快照更容易观察到慢速拟P波.     

快速有效的转换波共转换点叠加技术

转换波共转换点(Common Conversion Point简称CCP)叠加的关键在于CCP抽道集和非双曲线正常时差校正NMO(Normal MoveOut). 目前方法的精度限制了其在中-浅层或大炮检距情况下的应用. 我们对CCP叠加技术进行了系统研究,导出了新的CCP位置计算公式和非双曲线时距关系式,并给出了具体的CCP抽道集方法. 理论模型试验和实验资料处理表明本文方法精度高、简便易行,特别是对于中-浅层和大炮检距情况也能得到良好的叠加效果.     

Diagnosis of earth-fill dams by synthesised approach of sounding and surface wave method

The spatial distribution of the strength inside the earth-fill is identified by the sounding tests. In this research, the Swedish weight sounding (SWS) is employed, and the spatial high-density test is performed to identify the spatial correlation structure. Furthermore, the synthesised approach of the SWS and surface wave method, which is one of the geophysical method, is proposed to compensate the shortage of each approach. Consequently, the correlation structure of an earth-fill could be identified accurately, and the high resolution of the spatial distribution could be visualised based on the survey results.     

浅层高分辨率反射地震试验研究：以雄安新区第四系为例

雄安新区城市规划建设对地下空间的高精度勘探提出了高要求,作为具有高分辨率的地面物探方法,反射波地震勘探技术具有良好应用前景。本文探究了厚覆盖区条件下的高精度有效反射地震探测技术,并利用反射纵波技术在雄安新区进行提高浅层分辨率的地震勘探方法试验。根据实际探测目标,通过对比分析不同检波器、震源类型、道距大小及接收排列长度等不同采集因素试验效果,探究已有条件下适宜的震源类型和最佳采集参数。结合实际资料处理和分层解释,试验研究表明,采用高频检波器、选择可控震源激发效果最佳,针对试验目标,结合小道距、适当排列接收等措施,可有效提高厚覆盖区50~300 m深度内反射纵波地震勘探的纵向分辨率,试验获得了丰富的有效波组信息,为内部地层精细划分和解释提供了依据。该方法可为类似环境下城市地下空间的精细探测提供有益参考借鉴。     

浅水波长计算问题

浅水波周期保持常值,而波长缩减时其计算需进行迭代。提出了两个精度较高的经验公式供计算应用。水深波长比差值ΔLD0的变化规律反映出浅水波长缩减为一复杂过程,并与波能密度加大相关。LD0≤0.14时,特别是LD0=0.056前后,其特征会更加显著,由此可探讨波长变化的内在机理。     

斜坡平台波浪破碎数值模拟

波浪在斜坡地形上破碎,破波后稳定波高多采用物理模型试验方法进行研究,利用近岸波浪传播变形的抛物型缓坡方程和波能流平衡方程,导出了适用于斜坡上波浪破碎的数值模拟方法。首先根据波能流平衡方程和缓坡方程基本型式分析波浪在破波带内的波能变化和衰减率,推导了波浪传播模型中波能衰减因子和破波能量流衰减因子之间的关系;其次,利用陡坡地形上的高阶抛物型缓坡方程建立了波浪传播和波浪破碎数学模型;最后,根据物理模型试验实测数据对数值模拟的效果进行验证。数值计算与试验资料比较表明,该模型可以较好地模拟斜坡地形的波浪传播波高变化。     

Temporal variation of modulated-wave-train geometries and their influence on vertical bending moments of a container ship

A towing experiment was conducted using a modulated wave train to investigate the vertical bending responses of a hydro-structural container ship model. In the experiment, a spatially periodic modulated wave train, as a model of a freak wave in successive high waves mimicking the so-called three sisters, was generated by the recently established higher-order spectral method wave generation (HOSM-WG) method. HOSM-WG enables us to control the location and timing of the maximum crest height in a wave tank. With precise control of the towing carriage, an experiment was conducted in which the timing of the encounters between the ship model and the modulated wave train was accurately determined. The maximum sagging moment (SM) was found to increase in proportion with the encounter wave height. However, because of differences in the relative depth of the fore and aft troughs, the maximum SM is highly variable for a given wave height. The temporal wave-geometry evolution caused the relative trough-depth to vary significantly within a wave period in the vicinity of the maximum crest height. As a result, depending on the encounter timing, the SM varied considerably for a given wave height. The temporal variation of the wave geometry is a robust feature of a modulated wave train and is common between the spatially periodic and temporally periodic modulated wave trains.     

Obtaining the distribution of quiescent periods directly from the power spectral densities of Sea waves

There is a growing practical interest in the ability to increase the sea states at which marine operations can be safely undertaken by exploiting the quiescent periods that are well known to exist under a wide range of sea conditions. While the actual prediction of quiescent periods at sea for the control of operations is a deterministic process, the long term planning of future maritime tasks that rely on these quiescent periods is a statistical process involving the anticipated quiescence properties of the forecasted sea conditions in the geographical region of interest. It is in principle possible to obtain such data in tabular form either large scale simulation or from field data. However, such simulations are computationally intensive and libraries of appropriate field data are not common. Thus, it is clearly attractive to develop techniques that exploit standard wave spectral models for describing the quiescence statistics directly from such spectra. The present study focuses upon such techniques and is a first step towards the production of a computationally low-cost quiescence prediction tool and compares its efficacy against simulations. Two significant properties emerge for a large class of wave spectral models that encompasses the ubiquitous Neumann and Pierson Moskowitz or Bretschneider forms. Firstly, the auto-correlation function of the wave profile that are required to produce the quiescence property can be obtained analytically in terms of standard special functions. This considerably reduces the computational cost making desktop computer-based planning tools a reality. Secondly, for each class of these parametric spectra, the probability of a given number of consecutive wave heights (normalised to the significant wave heights) less than some critical value is in fact independent of absolute wave height. Thus, for a broad class of practically interesting wave spectra all that is required to obtain the statistical distribution of the quiescent periods is simple rescaling.     

Semi-analytical solution for the water wave diffraction by arrays of truncated circular cylinders in front of a vertical wall

The present study investigates the combined wave field that is induced by the continuous interaction of plane waves with an array of truncated circular cylinders in front of a rigid wall. The long-term goal of the study is the investigation of possible increase in the efficiency of cylindrical Wave Energy Converters (WECs) by putting in the vicinity of the array a barrier to propagation, a wall, that could assist the reflection of the incoming waves. The main task is to develop a generic solution method that is free of conceptual simplifications employed, e.g. by the method of images and the assumption of “pure” wave reflection. To cope with the set task, the proposed method relies on the semi-analytical formulation of the velocity potentials, while the solution is sought by combined expressions that involve polar and elliptical harmonics. The wall is represented as an elliptical cylinder with zero semi-minor axis. This assumption has eventually a beneficial effect to the underlying formulation given that it simplifies significantly the expansions of the involved diffraction potentials.     

涌浪对风浪能量影响的实验研究

在实验室风浪水槽中进行纯风浪和混合浪波面位移观测,研究波长较长的规则波对风浪能量的影响.本文用混合浪和纯风浪中的风浪显著波的零阶谱矩之比代表混合浪中的风浪与纯风浪能量之比,并以此表征涌浪对风浪能量的影响.研究了该能量比随涌浪波陡S、风区x、波龄倒数u/C、涌浪频率与纯风浪谱峰频率之比fs/fwp的变化规律.结果表明,涌浪对风浪能量的抑制作用随涌浪波陡的增加、波龄倒数的增大及涌浪频率与纯风浪谱峰频率之比的增大而增强.发现该能量比依赖于无因次量R=(1+80(πS)2)1.9(fs/fwp)0.9(u/C)0.27,并拟合得到2者的经验关系.此外,本文实验还发现,在某些情况下,涌浪的存在使风浪能量增加.