一种模拟非均匀介质中弹性波传播新的k-space方法          下载免费PDF全文

黄继伟  刘洪 《地球物理学报》2020,63(8):3091-3104

传统的伪谱(PS)方法,采用傅里叶变换(FT)计算空间导数具有很高的精度,每个波长仅需要两个采样点,而时间导数采用有限差分(FD)近似因而精度较低.当采用大时间步长时,由于时空精度不平衡,PS法存在不稳定性问题.原始的k-space方法可以有效地克服这些问题但是却无法适用于非均匀介质.为了提高原始k-space方法模拟非均匀介质波动方程的精度,我们提出了一种新的k-space算子族.它是用非均匀介质的变速度代替原k-space算子中的常数补偿速度构造得到,引入低秩近似可以高效求解.我们将构造的新的k-space算子应用于耦合的二阶位移波动方程,而不是交错网格一阶速度应力波动方程,使模拟弹性波的计算存储量减少.我们从数学上证明了基于二阶波动方程的k-space方法与基于一阶波动方程的k-space方法是等价的.数值模拟实验表明,与传统的PS、交错网格PS和原始的k-space方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.  相似文献   

基于分区间泛化Kriging近似模型的船舶局部结构形状优化     

张干锋  王德禹 《海洋工程》2018,36(6):46-57

为提高Kriging近似模型在船舶结构性能多维度响应预测方面的适用性,对常规Kriging近似模型进行分区间泛化改进:一是对设计样本点的各个维度(分量)进行划区,并在每一个分区间内采用最优拉丁超立方(OLhd)取样;二是引入比例系数w1组合高斯与指数型相关函数提高模型对数据的泛化能力,从而在每个划分的区间内建立泛化的Kriging近似模型。通过SCH测试函数,验证了构造的分区间泛化Kriging模型有效性。结合参数化建模和改进Kriging近似模型对某过渡肘板和舱口角隅边界进行形状优化,优化第一阶段由改进的近似模型通过多岛遗传算法得到全局初步的最优解,第二阶段在初步解的基础上缩小优化变量范围,由少量的FEM计算即可在小范围内搜寻到精确的最优设计变量。结果表明:分区间泛化Kriging近似模型在预测多维度响应时较常规Kriging模型预测精度更高;分阶段的形状优化流程在保证极小误差和缩小计算成本的情况下能够得到理想的应力分布和重量优化结果,有助于船舶结构的轻量化研究。  相似文献   

On effective characteristic of Rayleigh surface wave propagation in porous fluid-saturated media at low frequencies     

《Soil Dynamics and Earthquake Engineering》2014

  相似文献   

基于人工鱼群算法和SCA等效理论的井中横波速度预测方法     

魏忠宇  刘财  郭智奇  张冰  刘喜武 《世界地质》2018,37(1):282-288

在建立页岩岩石物理模型的基础上,根据等效自相容近似(SCA)岩石物理模型,构建出岩石的纵波速度、横波速度与岩石密度、组分和孔隙度等的定量关系,得出使理论纵波速度和实际纵波速度最接近的孔隙纵横比,进而将该孔隙纵横比作为约束条件来实现横波速度预测。反演算法利用人工鱼群算法来计算最佳孔隙纵横比,并将预测的横波速度与实际测得的横波速度对比,证明了人工鱼群算法的有效性。  相似文献   

Generation of oceanic internal gravity waves by a cyclonic surface stress disturbance     

《Dynamics of Atmospheres and Oceans》2019

Atmospheric cyclones with strong winds significantly impact ocean circulation, regional sea surface temperature, and deep water formation across the global oceans. Thus they are expected to play a key role in a variety of energy transport mechanisms. Even though wind-generated internal gravity waves are thought to contribute significantly to the energy balance of the deep ocean, their excitation mechanisms are only partly understood.The present study investigates the generation of internal gravity waves during a geostrophic adjustment process in a Boussinesq model with axisymmetric geometry. The atmospheric disturbance is set by an idealized pulse of cyclonic wind stress with a Rankine vortex structure. Strength, radius and duration of the forcing are varied. The effect upon wave generation of stratification with variable mixed-layer depth is also examined.Results indicate that internal gravity waves are generated after approximately one inertial period. The outward radial energy flux is dominated by waves having structure close to vertical mode-1 and with frequency close to the inertial frequency. Less energetic higher mode waves are observed to be generated close to the sea floor underneath the storm. The total radiated energy corresponds to approximately 0.02% of the wind input. Deeper mixed-layer conditions as well as weaker stratification reduce this fraction.The low energy transfer rates suggest that other processes that drive vertical motion like surface heat fluxes, turbulent motion, mixed region collapse and storm translation are essential for significant energy extraction by internal gravity waves to occur.  相似文献   

Bragg reflection of water waves by multiple floating horizontal flexible membranes with submerged rectangular bars on the seabed     

《Applied Ocean Research》2019

Bragg reflection of water waves by multiple floating horizontal flexible membranes is investigated based on the linear wave theory and the assumption of small membrane response. Under the floating horizontal membranes, periodical submerged rectangular bars are arranged on the flat seabed. The total reflection and transmission coefficients are obtained by using the eigenfunction expansion method and the wide spacing approximation. The calculated coefficients are validated with the results available in the literature, which shows that the present method is applicable. The characteristics of Bragg reflection are systematically investigated by changing various parameters including the height of the rectangular bars, the number, the tension, the spacing, and the length of the flexible membranes. The results can help designing multiple floating horizontal flexible membranes as effective floating breakwaters by taking advantage of Bragg reflection.  相似文献   

Characteristics of eigen-system for flexural gravity wave problems     

S. Mandal  A. Chakrabarti 《地球物理与天体物理流体动力学》2017,111(4):249-281

For a class of boundary value problems associated with flexural gravity waves, characteristics of the eigen-system are analysed which are used to study the convergence of the expansion formulae in both cases of single and two layer fluids in water of finite and infinite depths. Using Green’s function technique, the spectral representations of the vertical eigenfunctions are obtained which are used along with the orthogonal mode-coupling relations to prove the convergence of the expansion formulae. Using the expansion formulae flexural gravity wave scattering due to multiple articulations in the presence of compression are investigated. The problem is studied in both the cases of single-layer and two-layer fluids in finite water depth. Effects of compressive force, stiffness of the connectors, length of the elastic plates and water depth, position of interface on wave scattering by articulated plates are studied by analysing the reflection and transmission coefficients. This procedure of proving the convergence of expansion formulae is independent of water depth. The concept and methodology for dealing with wave-structure interaction problems discussed here can be generalised to deal with problems of similar nature arising in the broad area of mathematical physics and engineering.  相似文献   

Estimation of poromechanical and thermal conductivity properties of unsaturated isotropically microcracked cement pastes     

Benoît Bary 《国际地质力学数值与分析法杂志》2011,35(14):1560-1586

  相似文献   

Incorporating velocity shear into the magneto-Boussinesq approximation     

Jordan A. Bowker  Evy Kersalé 《地球物理与天体物理流体动力学》2014,108(5):553-567

Motivated by consideration of the solar tachocline, we derive, via an asymptotic procedure, a new set of equations incorporating velocity shear and magnetic buoyancy into the Boussinesq approximation. We demonstrate, by increasing the magnetic field scale height, how these equations are linked to the magneto-Boussinesq equations of Spiegel and Weiss (Magnetic buoyancy and the Boussinesq approximation. Geophys. Astrophys. Fluid Dyn. 1982, 22, 219–234).  相似文献   

Common‐shot Fresnel beam migration based on wave‐field approximation in effective vicinity under complex topographic conditions          下载免费PDF全文

Jianping Huang  Wenyuan Liao  Xin Wang  Zhenchun Li 《Geophysical Prospecting》2016,64(3):554-570

Gaussian beam migration is a versatile imaging method for geologically complex land areas, which overcomes the limitation of Kirchhoff migration in imaging multiple arrivals and has no steep‐dip limits of one‐way wave‐equation migration. However, its imaging accuracy depends on the geometry of Gaussian beam that is determined by the initial parameter of dynamic ray tracing. As a result, its applications in exploration areas with strong variations in topography and near‐surface velocity are limited. Combined with the concept of Fresnel zone and the theory of wave‐field approximation in effective vicinity, we present a more robust common‐shot Fresnel beam imaging method for complex topographic land areas in this paper. Compared with the conventional Gaussian beam migration for irregular topography, our method improves the beam geometry by limiting its effective half‐width with Fresnel zone radius. Moreover, through a quadratic travel‐time correction and an amplitude correction that is based on the wave‐field approximation in effective vicinity, it gives an accurate method for plane‐wave decomposition at complex topography, which produces good imaging results in both shallow and deep zones. Trials of two typical models and its application in field data demonstrated the validity and robustness of our method.  相似文献