共查询到19条相似文献，搜索用时 828 毫秒
1.
2.
给出了在非均匀横向各向同性(TI)介质情况下,四阶时间精度、高阶空间精度的一阶速度-应力P-SV波的波动方程交错网格有限差分解法.首先根据一阶速度(应力)波动方程把速度(应力)对时间的一阶和三阶导数转换为应力(速度)对空间的导数,从而在使用四阶时间精度有限差分格式计算某一时刻的波场时只需要前面两个时间步的波场值;然后在空间上采用高阶有限差分格式以提高数值模拟的精度.数值模拟结果和实测垂直地震剖面(VSP)记录符合得很好,说明该方法是可行的. 相似文献
3.
如何更精确高效地模拟大型、复杂系统内域的波动是发展和完善近场波动数值模拟技术的一个重要研究课题。内域波动的数值模拟通常采用计算量较小的显式方法,但现有的时空解耦显式有限元方法的精度只有二阶;低精度不仅影响数值模拟的精度,而且制约着计算效率的提高。鉴于此,本文旨在探索内域波动数值模拟具有更高精度且稳定的显式方法。作者发展了内域波动数值模拟的现有显式有限元解耦技术,提出了一种高精度且稳定的显式数值模拟方法。 相似文献
4.
大型结构地震反应值模拟中的波动输入 总被引：11，自引：6，他引：11
解耦的时域有限元数值模拟技术在考虑土－结构相互作用的大型结构地震反应分析中的应用日趋成熟。为完善其中的波动输入技术，本提出 了一种波动输入时步数值模拟的简便方法，替换对自由场的频域计算，从而使结构反应分析可以完全通过时步数值模拟实现。同时，通过数值实现，对这一实施方案的有效性进行了检验。 相似文献
5.
就大型近场波动的高效数值模拟而言,稳定实现高阶人工边界是一个尚未圆满解决的问题.本文针对使用多次透射公式的SH波动集中质量有限元模拟,依据GKS定理的群速度解释,进一步阐明了人工边界与内域离散格式耦合所导致高频失稳的机理,即两者支持群速度指向内域的外行高频平面谐波,波动能量自发地从人工边界进入內域,从而导致失稳,而这类谐波是由集中质量有限元离散引入的.本文提出了消除此种耦合失稳的一种方法:通过修改有限元刚度阵来改变内域离散格式,并保证修改格式的精度不低于原有格式的精度.理论分析和数值实验表明此法能稳定实现透射边界.本文研究结果具有推广应用前景. 相似文献
6.
《地震工程与工程振动》2017,(2)
本文依据波速有限物理概念,从一维波动方程柯西初值问题的解析解出发,推广应用一种在空间域上采用Lagrange多项式内插、时间域内精确积分的显式方法;并基于包含7个节点的3层网格离散模型,通过权衡精度和稳定性的要求,构建出时空离散精度皆为6阶的稳定递推格式,且具有显式和时空解耦特性。最后,通过算例进行与本文同源的低阶格式(2阶显式格式、4阶显式格式)的对比分析,验证该6阶显式格式的精度和稳定性等理论结果,指出高阶公式对提高计算效率的价值。 相似文献
7.
局部解耦的时域波分析方法 总被引：8，自引：1，他引：7
从模拟散射波的观点出发,推导了波动问题局部解耦的人工边界条件,与有限元内点的时域显式方法相结合建立了波动分析的局部解耦的直接分析方法,文中还从波动分解的概念,建立了外源问题（如地震作用）输入方法。 相似文献
8.
结构动力学方程的显式积分格式 总被引：1，自引：0，他引：1
本文从空间解耦有限元常微分方程组出发，探讨了结构动力学方程的高精度显式积分格式。通过被积函数的拉格朗日多项式内插和分部积分导出了波动数值模拟的一组显式时步积分公式。这组公式是时间和空间解耦的，即波场内任一离散节点在任一时刻的波动数据可以用这组公式依据该节点及其邻近节点在该时刻之前的n＋1个时刻的波动数据显式地算出（n为非负整数），阐明了这组公式的如下特点：第一，其截断误差的量级不超过0（Δt^n＋3），Δt为时间步距。第二，它不仅可用于线性波动的数值模拟，而且可用于本构方程具有强非线性情形。第三，这组公式也可推广应用于一系列数学物理暂态问题的数值求解。针对一个简单的时不变系统初步分析了此组积分格式的稳定性。但是，对其稳定性尚需作进一步研究。 相似文献
9.
一阶弹性波方程错网格高阶差分解法 总被引：26，自引：6，他引：20
提高计算精度和运算效率是所有波场正演方法所追求的目标,本文通过将速度（应力）对时间的奇数阶高阶转化为应力（速度）对空间的导数,运用时间和空间差分精度均可达任意阶的高阶差分法,通过交错网格技术,对一阶速度－应力弹性波方程进行了数值求解;波场快照以及实际模型的正演结果表明,这种求解一阶弹性小听高阶差分解法,和常规的差分法相比网格频散显著减小,精度明显提高,而且可以取较大的空间步长,提高计算效率。 相似文献
10.
提高计算精度和运算效率是所有波场正演方法所追求的目标,本文通过将速度 （应力）对时间的奇数阶高阶寻数转化为应力（速度）对空间的导数,运用时间和空间差分精度 均可达任意阶的高阶差分法,通过交错网格技术,对一阶速度－应力弹性波方程进行了数值求 解．波场快照以及实际模型的正演结果表明,这种求解一阶弹性波方程的高阶差分解法,和 常规的差分法相比网格频散显著减小,精度明显提高,而且可以取较大的空间步长,提高计算 效率。 相似文献
11.
A hybrid method combining finite element and 4th-order finite difference techniques is developed to model SH and P-SV seismic wave propagation in a 2D elastic medium with irregular surface topography. Both the classic staggered grid finite difference scheme and the partially staggered grid scheme are tested. The accuracy of the hybrid method is studied by comparison with a semi-analytical and another numerical method. Subsequently, to study the amplification, numerical simulations of seismic wave propagation in a series of hills are carried out and compared with the single-hill case. Depending on the position of the source in relation to the topography, the ratio between the heights and lengths of the hills or the ratio between the lengths of the hills and the wavelength, the presence of several hills as opposed to a single one can increase the amplification effect due to topography. This study highlights the fact that, when evaluating topographic site effects, surrounding topography must be taken into account in addition to local topography. 相似文献
12.
13.
基于应力、速度混合变量弹性波方程及任意四边形网格差分算子,给出了交错计算应力及速度的非规则网格弹性波应力一速度差分法该方法融合了有限元法能适应复杂形状边界及差分法无需计算刚度阵的特点,具有较高的计算精度,所需计算机存储空间较少,计算效率也很高．基于积分平衡方程引入了任意形状自由表面的边界条件,且通过局部滤波改善了自由表面边界条件的稳定性,使得该方法可应用于考虑地表形状影响的地震波数值模拟 相似文献
14.
A central difference method with low numerical dispersion for solving the scalar wave equation 总被引：4，自引：0，他引：4
In this paper, we propose a nearly‐analytic central difference method, which is an improved version of the central difference method. The new method is fourth‐order accurate with respect to both space and time but uses only three grid points in spatial directions. The stability criteria and numerical dispersion for the new scheme are analysed in detail. We also apply the nearly‐analytic central difference method to 1D and 2D cases to compute synthetic seismograms. For comparison, the fourth‐order Lax‐Wendroff correction scheme and the fourth‐order staggered‐grid finite‐difference method are used to model acoustic wavefields. Numerical results indicate that the nearly‐analytic central difference method can be used to solve large‐scale problems because it effectively suppresses numerical dispersion caused by discretizing the scalar wave equation when too coarse grids are used. Meanwhile, numerical results show that the minimum sampling rate of the nearly‐analytic central difference method is about 2.5 points per minimal wavelength for eliminating numerical dispersion, resulting that the nearly‐analytic central difference method can save greatly both computational costs and storage space as contrasted to other high‐order finite‐difference methods such as the fourth‐order Lax‐Wendroff correction scheme and the fourth‐order staggered‐grid finite‐difference method. 相似文献
15.
Optimal implicit staggered‐grid finite‐difference schemes based on the sampling approximation method for seismic modelling 
下载免费PDF全文
下载免费PDF全文 We propose new implicit staggered‐grid finite‐difference schemes with optimal coefficients based on the sampling approximation method to improve the numerical solution accuracy for seismic modelling. We first derive the optimized implicit staggered‐grid finite‐difference coefficients of arbitrary even‐order accuracy for the first‐order spatial derivatives using the plane‐wave theory and the direct sampling approximation method. Then, the implicit staggered‐grid finite‐difference coefficients based on sampling approximation, which can widen the range of wavenumber with great accuracy, are used to solve the first‐order spatial derivatives. By comparing the numerical dispersion of the implicit staggered‐grid finite‐difference schemes based on sampling approximation, Taylor series expansion, and least squares, we find that the optimal implicit staggered‐grid finite‐difference scheme based on sampling approximation achieves greater precision than that based on Taylor series expansion over a wider range of wavenumbers, although it has similar accuracy to that based on least squares. Finally, we apply the implicit staggered‐grid finite difference based on sampling approximation to numerical modelling. The modelling results demonstrate that the new optimal method can efficiently suppress numerical dispersion and lead to greater accuracy compared with the implicit staggered‐grid finite difference based on Taylor series expansion. In addition, the results also indicate the computational cost of the implicit staggered‐grid finite difference based on sampling approximation is almost the same as the implicit staggered‐grid finite difference based on Taylor series expansion. 相似文献
16.
在地震波场数值模拟中, 交错网格有限差分技术得到了广泛的应用, 但是在弹性模量变化较大时, 通常会因插值而导致模拟误差增大. 旋转交错网格可以很好地克服这个缺点, 因而适合于各向异性介质正演模拟. 但是对于同样大小的网格单元, 旋转交错网格需要的步长比常规交错网格要大, 这会使梯度和散度算子的误差增大因而更易产生空间数值频散. 针对这些问题, 本文提出了旋转交错网格与紧致有限差分相结合的方法, 并基于模拟退火算法进行全局优化, 压制数值频散, 拓宽波数范围. 数值模拟结果表明, 此方法可以有效地压制数值频散, 且具有较高的模拟精度. 相似文献
17.
以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征. 相似文献
18.
本文将DRP/opt MacCormack有限差分格式用于模拟二维各向异性介质中的地震波传播.DRP/opt MacCormack是一种同位网格下的差分格式,避免了传统的交错网格在计算各向异性问题时由于变量插值而导致的误差.而且相对于低阶同位网格差分格式,它具有低色散、低耗散的优点.此格式将中心差分算子分成前向和后向两个空间单边差分,然后在4-6步Runge-Kutta时间积分中使用单边差分组合.在具有垂直对称轴的横向各向同性（VTI）模型下,通过对比DRP/opt MacCormack有限差分和谱元方法的模拟结果,验证了前者具有很高的精度和稳定性.由于实际地质条件下TI介质的对称轴通常是倾斜的(TTI),本文在二维三分量框架下模拟TTI介质中的地震波场.结果显示横波分裂和切平面/反平面运动耦合的特征.数值实验表明DRP/opt MacCormack是一种有效的研究各向异性介质中地震波传播规律的差分格式. 相似文献
19.
This paper presents a Lebedev finite difference scheme on staggered grids for the numerical simulation of wave propagation in an arbitrary 3D anisotropic elastic media. The main concept of the scheme is the definition of all the components of each tensor (vector) appearing in the elastic wave equation at the corresponding grid points, i.e., all of the stresses are stored in one set of nodes while all of the velocity components are stored in another. Meanwhile, the derivatives with respect to the spatial directions are approximated to the second order on two‐point stencils. The second‐order scheme is presented for the sake of simplicity and it is easy to expand to a higher order. Another approach, widely‐known as the rotated staggered grid scheme, is based on the same concept; therefore, this paper contains a detailed comparative analysis of the two schemes. It is shown that the dispersion condition of the Lebedev scheme is less restrictive than that of the rotated staggered grid scheme, while the stability criteria lead to approximately equal time stepping for the two approaches. The main advantage of the proposed scheme is its reduced computational memory requirements. Due to a less restrictive dispersion condition and the way the media parameters are stored, the Lebedev scheme requires only one‐third to two‐thirds of the computer memory required by the rotated staggered grid scheme. At the same time, the number of floating point operations performed by the Lebedev scheme is higher than that for the rotated staggered grid scheme. 相似文献