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

南海东北部内潮的模态结构和传播     

刘倩  谢晓辉  尚晓东  陈桂英  汪宏 《海洋学报(英文版)》2019,38(9):12-23

The evolution of energy, energy flux and modal structure of the internal tides(ITs) in the northeastern South China Sea is examined using the measurements at two moorings along a cross-slope section from the deep continental slope to the shallow continental shelf. The energy of both diurnal and semidiurnal ITs clearly shows a～14-day spring-neap cycle, but their phases lag that of barotropic tides, indicating that ITs are not generated on the continental slope. Observations of internal tidal energy flux suggest that they may be generated at the Luzon Strait and propagate west-northwest to the continental slope in the northwestern SCS. Because the continental slope is critical-supercritical with respect to diurnal ITs, about 4.6 kJ/m~2 of the incident energy and 8.7 kW/m of energy flux of diurnal ITs are reduced from the continental slope to the continental shelf. In contrast, the semidiurnal internal tides enter the shelf because of the sub-critical topography with respect to semidiurnal ITs.From the continental slope to the shelf, the vertical structure of diurnal ITs shows significant variation, with dominant Mode 1 on the deep slope and dominant higher modes on the shelf. On the contrary, the vertical structure of the semidiurnal ITs is stable, with dominant Mode 1.  相似文献   

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.  相似文献   

煤体静载破坏微震、电磁辐射及裂纹扩展特征研究     

艾迪昊  李成武  赵越超  李光耀 《岩土力学》2020,41(6):2043-2051

为研究型煤在单轴压缩破裂过程中产生的微震、电磁辐射信号与裂纹演化特征的对应关系，利用自主设计的低噪声静态加载试验系统，对0～0.25、0.25～0.5、0.5～1.0、1～2 mm共4种不同粒度的型煤进行了单轴压缩破坏试验，同步采集了煤样破坏过程中的微震、电磁辐射信号及破坏视频图像，提出了一种煤体裂纹快速提取方法并计算了型煤裂纹面积的变化规律。研究结果表明：型煤在单轴压缩过程中产生的微震、电磁辐射信号及裂纹面积在时域上具有良好的同步性。型煤破坏过程中裂纹面积随时间变化曲线可分为4个阶段。第1阶段为压实阶段，煤样所受应力值较小，其表面的裂纹面积以极为缓慢的速率增加。第2阶段为初始破裂阶段，随着应力的增加和内部弹性势能的积聚，型煤表面裂纹面积的增加速率较压实阶段有明显提高，伴随产生许多细小裂纹。第3阶段为加速破裂阶段，随着应力继续增加以及材料内部积聚弹性能的释放，试样变形过程加速，导致裂纹面积增速进一步增大。第4阶段为卸压阶段，试样的裂纹面积达到最大值，承载能力急剧降低，发生失稳破坏。  相似文献   

三峡库区龙门寨危岩体崩塌产生涌浪研究   总被引：1，自引：1，他引：0  

郑嘉豪  黄波林  张全  赵海林  冯万里  王健  陈小婷 《地质力学学报》2020,26(4):533-543

长江两岸高耸的危岩体对航道、沿岸居民带来巨大安全隐患。大宁河属于长江一级支流，龙门寨危岩体位于大宁河上，距离巫山县城仅1 km。利用FLOW-3D软件，模拟了145 m、175 m两种水位工况下龙门寨危岩体崩塌产生涌浪过程和涌浪传播过程。模拟结果表明，涌浪在145 m水位工况下最大浪高约为17.9 m，175 m水位工况下最大浪高约为11.6 m；在巫山县的五个码头处，两种水位工况最大涌浪爬高分别约为10.9 m、3.8 m；根据涌浪高度，对大宁河进行危险分区，145 m水位工况下极高危险区长度约4.4 km，很高危险区长度约1.9 km；175 m水位工况下极高危险区长度约3.0 km，很高危险区长度约1.0 km。研究结果有助于防控龙门寨危岩体潜在涌浪灾害危害，保障大宁河航道和巫山县码头安全，同时也为三峡库区滑坡涌浪灾害提供了预警依据。   相似文献   

Detection thresholds in structural health monitoring     

《Soil Dynamics and Earthquake Engineering》2014

  相似文献   

Sampling efficiency in Monte Carlo based uncertainty propagation strategies: Application in seawater intrusion simulations     

《Advances in water resources》2014

The implementation of Monte Carlo simulations (MCSs) for the propagation of uncertainty in real-world seawater intrusion (SWI) numerical models often becomes computationally prohibitive due to the large number of deterministic solves needed to achieve an acceptable level of accuracy. Previous studies have mostly relied on parallelization and grid computing to decrease the computational time of MCSs. However, another approach which has received less attention in the literature is to decrease the number of deterministic simulations by using more efficient sampling strategies. Sampling efficiency is a measure of the optimality of a sampling strategy. A more efficient sampling strategy requires fewer simulations and less computational time to reach a certain level of accuracy. The efficiency of a sampling strategy is highly related to its space-filling characteristics.This paper illustrates that the use of optimized Latin hypercube sampling (OLHS) strategies instead of the widely employed simple random sampling (SRS) and Latin hypercube sampling (LHS) strategies, can significantly improve sampling efficiency and hence decrease the simulation time of MCSs. Nine OLHS strategies are evaluated including: improved Latin hypercube sampling (IHS); optimum Latin hypercube (OLH) sampling; genetic optimum Latin hypercube (GOLH) sampling; three sampling strategies based on the enhanced stochastic evolutionary (ESE) algorithm namely φ_p-ESE which employs the φ_p space-filling criterion, CLD-ESE which utilizes the centered L₂-discrepancy (CLD) space-filling criterion, and SLD-ESE which uses the star L₂-discrepancy (SLD) space-filling criterion; and three sampling strategies based on the simulated annealing (SA) algorithm namely φ_p-SA which employs the φ_p criterion, CLD-SA which uses the CLD criterion, and SLD-SA which utilizes the SLD criterion. The study applies SRS, LHS and the nine OLHS strategies to MCSs of two synthetic test cases of SWI. The two test cases are the Henry problem and a two-dimensional radial representation of SWI in a circular island. The comparison demonstrates that the CLD-ESE strategy is the most efficient among the evaluated strategies. This paper also demonstrates how the space-filling characteristics of different OLHS designs change with variations in the input arguments of their optimization algorithms.  相似文献   

1D System identification of a 54‐story steel frame building by seismic interferometry          下载免费PDF全文

Mohammadtaghi Rahmani  Maria I. Todorovska 《地震工程与结构动力学》2014,43(4):627-640

A 54‐story steel, perimeter‐frame building in downtown Los Angeles, California, is identified by a wave method using records of the Northridge earthquake of 1994 (M_L = 6.4, R = 32 km). The building is represented as a layered shear beam and a torsional shaft, characterized by the corresponding velocities of vertically propagating waves through the structure. The previously introduced waveform inversion algorithm is applied, which fits in the least squares sense pulses in low‐pass filtered impulse response functions computed at different stories. This paper demonstrates that layered shear beam and torsional shaft models are valid for this building, within bands that include the first five modes of vibration for each of the North–South (NS), East–West (EW), and torsional responses (0–1.7 Hz for NS and EW, and 0–3.5 Hz for the torsional response). The observed pulse travel time from ground floor to penthouse level is τ ≈1.5 s for NS and EW and τ ≈ 0.9 s for the torsional responses. The identified equivalent uniform shear beam wave velocities are β_eq ≈ 140 m/s for NS and EW responses, and 260 m/s for torsion, and the apparent Q ≈ 25 for the NS and torsional, and ≈14 for the EW response. Across the layers, the wave velocity varied 90–170 m/s for the NS, 80–180 m/s for the EW, and 170–350 m/s for the torsional responses. The identification method is intended for use in structural health monitoring. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

泉州盆地地震效应特征的一维等效线性和二维非线性分析的比较     

董菲蕃  陈国兴  金丹丹 《地震学刊》2014,(2):154-160

采用等效线性动粘弹性模型描述土的动力非线性特性,基于一维等效线性波传法,对泉州盆地地震效应进行了分析;同时,采用修正Martin-Seed-Davidenkov动粘弹塑性模型描述土的动力非线性特性,对泉州盆地非线性地震效应进行了大尺度二维精细化有限元分析,研究了地形地貌和土层横向不均匀性对地震效应的影响。将两种分析结果进行对比,结果表明：①随着基岩输入地震动强度增大,地表峰值加速度PGA放大效应总体呈现减小趋势,中震与小震、大震与小震的地表PGA放大系数之比依次为0.83~0.99、0.72~0.97;②该盆地Ⅲ类场地处,基岩、地表起伏不大,且土层横向分布较均匀,两种方法计算得到的地震效应特征类似;基岩或地表起伏剧烈、土层横向分布明显不均匀的Ⅱ类场地上,二维非线性分析给出的地表PGA放大系数明显大于一维等效线性结果,两种方法得到的地表加速度反应谱及PGA随土层深度的变化特征存在显著差异,二维非线性分析给出的地表加速度反应谱大多呈现双峰甚至多峰现象,且PGA在土层特定深度处存在聚集效应,使PGA随土层深度的变化呈现非单调性。  相似文献   

基于GPU加速的Boussinesq类波浪传播变形数值模型     

孙家文  朱桐  房克照  刘忠波 《海洋工程》2020,38(2):111-119

Boussinesq波浪模型是一类相位解析模型,在时域内求解需要较高的空间和时间分辨率以保证计算精度。为提高计算效率,有必要针对该类模型开展并行算法的研究。与传统的中央处理器(CPU)相比,图形处理器(GPU)有大量的运算器,可显著提高计算效率。基于统一计算设备架构CUDA C语言和图形处理器,实现了Boussinesq模型的并行运算。将本模型的计算结果同CPU数值模拟结果和解析解相比较,发现得到的结果基本一致。同时也比较了CPU端与GPU端的计算效率,结果表明,GPU数值模型的计算效率有明显提升,并且伴随数值网格的增多,提升效果更为明显。  相似文献