一种模拟非均匀介质中弹性波传播新的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方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.  相似文献   

A nearly analytic symplectic partitioned Runge–Kutta method based on a locally one-dimensional technique for solving two-dimensional acoustic wave equations     

Chol Sim  Chunyou Sun  Nam Yun 《Geophysical Prospecting》2020,68(4):1253-1269

In this paper, we develop a new nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique for numerically solving two-dimensional acoustic wave equations. We first split two-dimensional acoustic wave equation into the local one-dimensional equations and transform each of the split equations into a Hamiltonian system. Then, we use both a nearly analytic discrete operator and a central difference operator to approximate the high-order spatial differential operators, which implies the symmetry of the discretized spatial differential operators, and we employ the partitioned second-order symplectic Runge–Kutta method to numerically solve the resulted semi-discrete Hamiltonian ordinary differential equations, which results in fully discretized scheme is symplectic unlike conventional nearly analytic symplectic partitioned Runge–Kutta methods. Theoretical analyses show that the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique exhibits great higher stability limits and less numerical dispersion than the nearly analytic symplectic partitioned Runge–Kutta method. Numerical experiments are conducted to verify advantages of the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique, such as their computational efficiency, stability, numerical dispersion and long-term calculation capability.  相似文献   

Response of a scale-model pile group for a jacket foundation of an offshore wind turbine in liquefiable ground during shaking table tests     

Yung-Yen Ko  Yi-Ting Li 《地震工程与结构动力学》2020,49(15):1682-1701

To investigate the seismic response of a pile group during liquefaction, shaking table tests on a 1/25 scale model of a 2 × 2 pile group were conducted, which were pilot tests of a test project of a scale-model offshore wind turbine with jacket foundation. A large laminar shear box was utilized as the soil container to prepare a liquefiable sandy ground specimen. The pile group model comprising four slender aluminum piles with their pile heads connected by a rigid frame was designed with similitude considerations focusing on soil–pile interaction. The input motions were 2-Hz sinusoids with various acceleration amplitudes. The excess pore water pressure generation indicated that the upper half of the ground specimen reached initial liquefaction under the 50-gal-amplitude excitation, whereas in the 75-gal-amplitude test, almost entire ground was liquefied. Accelerations in soil, on the movable frames composing the laminar boundary of the shear box, and along the pile showed limited difference at the same elevation before liquefaction. After liquefaction, the soil and the movable-frame accelerations that represented the ground response considerably reduced, whereas both the movable frames and the piles exhibited high-frequency jitters other than 2-Hz sinusoid, and meantime, remarkable phase difference between the responses of the pile group and the ground was observed, all probably due to the substantial degradation of liquefied soil. Axial strains along the pile implied its double-curvature bending behavior, and the accordingly calculated moment declined significantly after liquefaction. These observations demonstrated the interaction between soil and piles during liquefaction.  相似文献   

黄土山区煤层群分层开采地表位移规律分析     

沈宇恒 《北京测绘》2020,(4):490-494

为了研究黄土山区地形对煤层群开采地表位移的影响,以榆林某矿为地质原型,利用数值软件模拟地表为近水平、正坡、负坡、凸面和凹面五种地形条件下煤层群开采。结果表明:黄土山区地形的变化,对煤层群开采顶板垂直应力影响不明显,对地表位移影响较大;单向坡开采,斜坡滑移引起的下沉量在坡顶与开采沉陷量形成“叠加”,在坡底形成“抵消”;在组合坡中,煤层群开采地表位移量相差很大,说明在变坡点附近地表位移较为复杂。  相似文献   

Experimental investigation on a cabin-suspended catamaran in terms of motion reduction and wave energy harvesting by means of a semi-active motion control system     

《Applied Ocean Research》2019

A novel concept catamaran equipped with a suspended cabin, named Wave Harmonizer Type 4 (WHzer-4), is proposed and evaluated. The mass-spring-mass system is constructed by mounting four sets of suspensions in-between the cabin and the twin-hull. Two sets of dual motor/generators (M/Gs) are attached on the center beam of the cabin's deck fore and aft. Each shaft-end of the dual M/Gs is connected to the twin-hull through a rack-pinion gear unit. In this way the vertical relative motion between the cabin and the twin-hull can be transferred into the rotational motion of the M/Gs, and vice versa. A semi-active motion control system, which contains a proportional-integral (PI) controller, is designed and applied to each of the dual M/Gs for the aim of absorbing wave energy under the condition of suppressing the local vertical velocity of the cabin as much as possible. A 1/5 scale model ship with a length of 1.6 m is built, and a forced-oscillation bench test is implemented to validate the performance of the control system. Then, a series of towing tank tests is carried out in regular head waves. The heave and pitch responses of the cabin, those of the twin-hull and the corresponding wave energy capture width ratio (CWR) at five control scenarios and two reference scenarios are investigated. Discussion on the results of the tank test shows that the motion reduction of the cabin and the wave energy harvesting can be achieved simultaneously at a few wave conditions. However, at other conditions, although noticeable amount of wave energy is harvested, motion reduction of the heave and pitch of the cabin could not be obtained at the same time. It is suggested that varying the gain settings of the PI controllers according to the location of the controllers may improve the effectiveness of the proposed control system.  相似文献   

Hydrodynamic performance of a T-shaped floating breakwater     

《Applied Ocean Research》2019

The comprehensive utilization of floating breakwaters, specially acting as a supporting structure for offshore marine renewable energy explorations, has received more and more attention recently. Based on linear water-wave theory, the hydrodynamic performance of a T-shaped floating breakwater is semi-analytically investigated through the matched eigenfunction expansion method (MEEM). Auxiliary functions, to speed up the convergence and improve the accuracy in the numerical computations, are introduced to represent the singular behavior of fluid field near the lower salient corners of the structure. The effects of the height and installation position of the vertical screen on the reflection and transmission coefficients, dynamic response and wave forces are examined. It is found that the presence of the screen shifts the resonance frequency of RAO for both surge and pitch modes to the low-frequency area, while has no effect on heave mode. The identical added masses, damping and transmission coefficients can be obtained in the cases where the screen holds the same distance away from the longitudinal central axis of the upper box-type structure. Moreover, a relatively small pitch response can be achieved in a wide wave–frequency range, when the breakwater is Γ-shaped.  相似文献   

GPS/PWV时间序列特征提取方法的研究     

胡广保  叶世榕  张彦祥  夏朋飞  夏凤雨 《大地测量与地球动力学》2019,39(1):36-40

针对GAMIT/GLOBK软件解算得到的4 a GPS/PWV时间序列的特征提取问题，提出基于滤波辅助的部分集成经验模态分解（PEEMD）与Hilbert谱分析相结合的特征提取方法。首先，在PEEMD方法的基础上,结合滤波辅助的PEEMD方法与Hilbert谱分析，建立GPS/PWV时间序列特征提取模型；然后，将所提出的方法应用于TNML测站4 a的GPS/PWV长时间序列和7 d的GPS/PWV短时间序列分析中，并将滤波辅助的PEEMD结果与传统的小波分解结果进行对比。结果表明，该特征提取方法能准确有效地提取出 GPS/PWV时间序列中的周年周期和日周期特征分量，滤波辅助的PEEMD分解结果与小波分解结果一致，且提取的特征分量与原始信号更加吻合。  相似文献   

Modelling of viscoacoustic wave propagation in transversely isotropic media using decoupled fractional Laplacians     

Zhihao Qiao  Chengyu Sun  Jie Tang 《Geophysical Prospecting》2020,68(8):2400-2418

A new wave equation is derived for modelling viscoacoustic wave propagation in transversely isotropic media under acoustic transverse isotropy approximation. The formulas expressed by fractional Laplacian operators can well model the constant-Q (i.e. frequency-independent quality factor) attenuation, anisotropic attenuation, decoupled amplitude loss and velocity dispersion behaviours. The proposed viscoacoustic anisotropic equation can keep consistent velocity and attenuation anisotropy effects with that of qP-wave in the constant-Q viscoelastic anisotropic theory. For numerical simulations, the staggered-grid pseudo-spectral method is implemented to solve the velocity–stress formulation of wave equation in the time domain. The constant fractional-order Laplacian approximation method is used to cope with spatial variable-order fractional Laplacians for efficient modelling in heterogeneous velocity and Q media. Simulation results for a homogeneous model show the decoupling of velocity dispersion and amplitude loss effects of the constant-Q equation, and illustrate the influence of anisotropic attenuation on seismic wavefields. The modelling example of a layered model illustrates the accuracy of the constant fractional-order Laplacian approximation method. Finally, the Hess vertical transversely isotropic model is used to validate the applicability of the formulation and algorithm for heterogeneous media.  相似文献   

结合PS特征点的SBAS地表形变监测     

聂运菊  熊佳诚  程朋根  罗跃 《测绘通报》2020,(4):91-95

针对SARscape SBAS技术存在人工选取地面控制点影响监测结果的问题，本文提出了一种结合PS特征点的SBAS地表形变监测方法。通过设置PS技术中的振幅离差指数与相干性阈值，提取PS监测结果中的特征点作为地面控制点进行轨道精炼与形变反演。以2017年3月至2018年12月覆盖南昌地区的Sentinel-1影像为数据源，进行了地表形变监测与验证。结果表明：①该方法与PS方法监测结果均显示南昌市城区在监测时段内整体呈小幅度上升趋势，年平均沉降速率相关性达到了0.959，具有高度一致性。②将该方法、PS方法、人工选取地面控制点的SBAS方法的监测结果与水准数据进行对比，该方法与水准数据的误差不超过3.5 mm，说明了该方法的有效性。  相似文献   

基于深度残差网络的机载LiDAR点云分类     

赵传  郭海涛  卢俊  余东行  张保明 《测绘学报》2020,49(2):202-213

机载LiDAR点云的分类是利用其进行城市场景三维重建的关键步骤之一。为充分利用现有的图像领域性能较好的深度学习网络模型,提高点云分类精度,并降低训练时间和对训练样本数量的要求,本文提出一种基于深度残差网络的机载LiDAR点云分类方法。首先提取归一化高程、表面变化率、强度和归一化植被指数4种具有较高区分度的点云低层次特征;然后通过设置不同的邻域大小和视角,利用所提出的点云特征图生成策略,得到多尺度和多视角点云特征图;再将点云特征图输入到预训练的深度残差网络,提取多尺度和多视角深层次特征;最后构建并训练神经网络分类器,利用训练的模型对待分类点云进行预测,经后处理得到分类结果。利用ISPRS三维语义标记竞赛的公开标准数据集进行试验,结果表明,本文方法可有效区分建筑物、地面、车辆等8类地物,分类结果的总体精度为87.1%,可为城市场景三维重建提供可靠的信息。  相似文献