共查询到20条相似文献，搜索用时 22 毫秒
1.
Application of the nearly perfectly matched layer for seismic wave propagation in 2D homogeneous isotropic media 总被引：1，自引：0，他引：1
Jingyi Chen 《Geophysical Prospecting》2011,59(4):662-672
Numerical modelling plays an important role in helping us understand the characteristics of seismic wave propagation. The presence of spurious reflections from the boundaries of the truncated computational domain is a prominent problem in finite difference computations. The nearly perfectly matched layer has been proven to be a very effective boundary condition to absorb outgoing waves in both electromagnetic and acoustic media. In this paper, the nearly perfectly matched layer technique is applied to elastic isotropic media to further test the method's absorbing ability. The staggered‐grid finite‐difference method (fourth‐order accuracy in space and second‐order accuracy in time) is used in the numerical simulation of seismic wave propagation in 2D Cartesian coordinates. In the numerical tests, numerical comparisons between the nearly perfectly matched layer and the convolutional perfectly matched layer, which is considered the best absorbing layer boundary condition, is also provided. Three numerical experiments demonstrate that the nearly perfectly matched layer has a similar performance to the convolutional perfectly matched layer and can be a valuable alternative to other absorbing layer boundary conditions. 相似文献
2.
We derive a governing second-order acoustic wave equation in the time domain with a perfectly matched layer absorbing boundary condition for general inhomogeneous media. Besides, a new scheme to solve the perfectly matched layer equation for absorbing reflections from the model boundaries based on the rapid expansion method is proposed. The suggested scheme can be easily applied to a wide class of wave equations and numerical methods for seismic modelling. The absorbing boundary condition method is formulated based on the split perfectly matched layer method and we employ the rapid expansion method to solve the derived new perfectly matched layer equation. The use of the rapid expansion method allows us to extrapolate wavefields with a time step larger than the ones commonly used by traditional finite-difference schemes in a stable way and free of dispersion noise. Furthermore, in order to demonstrate the efficiency and applicability of the proposed perfectly matched layer scheme, numerical modelling examples are also presented. The numerical results obtained with the put forward perfectly matched layer scheme are compared with results from traditional attenuation absorbing boundary conditions and enlarged models as well. The analysis of the numerical results indicates that the proposed perfectly matched layer scheme is significantly effective and more efficient in absorbing spurious reflections from the model boundaries. 相似文献
3.
Mesh-free discretization, flexibly distributing nodes without computationally expensive meshing process, is able to deal with staircase problem, oversampling and undersampling problems and saves plenty of nodes through distributing nodes suitably with respect to irregular boundaries and model parameters. However, the time-domain mesh-free discretization usually exhibits poorer stability than that in regular grid discretization. In order to reach unconditional stability and easy implementation in parallel computing, we develop the frequency-domain finite-difference method in a mesh-free discretization, incorporated with two perfectly matched layer boundary conditions. Furthermore, to maintain the flexibility of mesh-free discretization, the nodes are still irregularly distributed in the absorbing zone, which complicates the situation of artificial boundary reflections. In this paper, we implement frequency-domain acoustic wave modelling in a mesh-free system. First, we present the perfectly matched layer boundary condition to suppress spurious reflections. Moreover, we develop the complex frequency shifted–perfectly matched layer boundary condition to improve the attenuation of grazing waves. In addition, we employ the radial-basis-function-generated finite difference method in the mesh-free discretization to calculate spatial derivatives. The numerical experiment on a rectangle homogeneous model shows the effectiveness of the perfectly matched layer boundary condition and the complex frequency shifted–perfectly matched layer boundary condition, and the latter one is better than the former one when absorbing large angle incident waves. The experiment on the Marmousi model suggests that the complex frequency shifted–perfectly matched layer boundary condition works well for complicated models. 相似文献
4.
模拟弹性波在孔隙介质中传播，需要稳定有效的吸收边界来消除或尽可能的减小由人工边界引起的虚假反射. 本文在前人工作基础上，首次建立了弹性孔隙介质情况下完全匹配层吸收边界的高阶速度-应力交错网格有限差分算法，并详细讨论了完全匹配层的构建及其有限差分算法实现. 首先，本文通过均匀孔隙模型的数值解与解析解的对比，验证所提出的数值方法的正确性；然后，本文考察了完全匹配层对不同入射角度入射波和自由表面上的瑞利波的吸收性能，将完全匹配层与廖氏和阻尼吸收边界进行了对比，研究了这三种吸收边界在不同吸收厚度情况下对弹性波吸收能力. 数值结果表明，在孔隙介质中，完全匹配层作为吸收边界能十分有效地吸收衰减外行波，无论对体波还是面波，是一种高效边界吸收算法. 相似文献
5.
6.
The complex‐valued first‐arrival traveltime can be used to describe the properties of both velocity and attenuation as seismic waves propagate in attenuative elastic media. The real part of the complex‐valued traveltime corresponds to phase arrival and the imaginary part is associated with the amplitude decay due to energy absorption. The eikonal equation for attenuative vertical transversely isotropic media discretized with rectangular grids has been proven effective and precise to calculate the complex‐valued traveltime, but less accurate and efficient for irregular models. By using the perturbation method, the complex‐valued eikonal equation can be decomposed into two real‐valued equations, namely the zeroth‐ and first‐order traveltime governing equations. Here, we first present the topography‐dependent zeroth‐ and first‐order governing equations for attenuative VTI media, which are obtained by using the coordinate transformation from the Cartesian coordinates to the curvilinear coordinates. Then, we apply the Lax–Friedrichs sweeping method for solving the topography‐dependent traveltime governing equations in order to approximate the viscosity solutions, namely the real and imaginary parts of the complex‐valued traveltime. Several numerical tests demonstrate that the proposed scheme is efficient and accurate in calculating the complex‐valued P‐wave first‐arrival traveltime in attenuative VTI media with an irregular surface. 相似文献
7.
Seismic anisotropy has an important influence on seismic data processing and interpretation. Although the frequency-domain seismic wavefield simulation has a problem of solving the large scale linear sparse matrix due to the computational limitations, it has some advantages over the time-domain seismic wavefield simulation including efficient inversion using only a limited number of frequency components and easy implementation of multiple sources. To accurately simulate seismic wave propagation in the frequency domain, we also need to choose the absorbing boundary conditions to absorb artificial reflections from edges of the model as we do in the time domain. Compared with the classical boundary conditions including the perfectly matched layer and complex frequency-shifted perfectly matched layer, the complex frequency-shifted multi-axial perfectly matched layer has been proven to effectively suppress the unwanted reflections at grazing incidence and solve the instability problem in the time-domain seismic numerical modelling in anisotropic elastic media. In this paper, we propose to extend the complex frequency-shifted multi-axial perfectly matched layer absorbing boundary condition to the frequency-domain seismic wavefield simulation in anisotropic elastic media. To test the validity of our proposed algorithm, we compare the results (snapshots and seismograms) of the frequency-domain seismic wavefield simulation with those of the time-domain modelling. The model studies indicate that the complex frequency-shifted multi-axial perfectly matched layer absorbing boundary condition is stable in the frequency-domain seismic wavefield simulation in anisotropic media, and provides better absorbing performance than the complex frequency-shifted perfectly matched layer boundary condition. 相似文献
8.
根据具有垂直对称轴的横向各向同性（VTI）介质中的一阶准P波方程,导出了该方程在交错网格中逆时延拓的高阶有限差分格式,给出了其稳定性条件,采用完全匹配层吸收边界条件解决边界反射问题,分别应用下行波最大能量法和归一化互相关成像条件, 实现了VTI介质中准P波方程的叠前逆时深度偏移．各向异性Marmousi模型的试算结果表明,VTI介质准P波方程叠前逆时深度偏移算法不受地下构造倾角和介质横向速度变化的限制,对复杂模型具有良好的成像能力;应用归一化互相关成像条件能得到更好的成像效果．对比该模型的各向异性和各向同性逆时偏移剖面表明,在各向异性地区采集的纵波数据用各向异性偏移算法理论上能得到更好的成像结果． 相似文献
9.
Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations 总被引：2，自引：0，他引：2
The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second-order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite-element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media. 相似文献
10.
The multi-axial perfectly matched layer (M-PML) is a material boundary condition for wave propagation problems in unbounded domains. It is obtained by extending the formulation of the split-field perfectly matched layer to a more general absorptive medium, for which damping profiles are specified along all dimensions of the problem. Under the hypothesis of small damping, it has been demonstrated that the stability of the system of partial differential equations of the M-PML can be related to the ratio of the damping profiles, and stable M-PML terminations for isotropic and orthotropic elastic media have been constructed. In the present work, we use the Routh–Horwitz determinants to demonstrate that the conclusions regarding the stability of M-PML for isotropic media for small damping are in fact valid for the more general case of damping coefficients of any (positive) value. The effectiveness of the M-PML is demonstrated by constructing stable terminations for the abovementioned media. The stability analysis is presented for 2-D in-plane (P-SV) wave propagation in elastic isotropic continua. 相似文献
11.
Multi‐block finite‐difference method for 3D elastodynamic simulations in anisotropic subhorizontally layered media 
下载免费PDF全文
下载免费PDF全文 Prediction of elastic full wavefields is required for reverse time migration, full waveform inversion, borehole seismology, seismic modelling, etc. We propose a novel algorithm to solve the Navier wave equation, which is based on multi‐block methodology for high‐order finite‐difference schemes on curvilinear grids. In the current implementation, the blocks are subhorizontal layers. Smooth anisotropic heterogeneous media in each layer can have strong discontinuities at the interfaces. A curvilinear adaptive hexahedral grid in blocks is generated by mapping the original 3D physical domain onto a parametric cube with horizontal layers and interfaces. These interfaces correspond to the main curvilinear physical contrast interfaces of a subhorizontally layered formation. The top boundary of the parametric cube handles the land surface with smooth topography. Free‐surface and solid–solid transmission boundary conditions at interfaces are approximated with the second‐order accuracy. Smooth media in the layers are approximated up to sixth‐order spatial schemes. All expected properties of the developed algorithm are demonstrated in numerical tests using corresponding parallel message passing interface code. 相似文献
12.
13.
The scaled boundary finite‐element method is extended to simulate time‐harmonic responses of non‐homogeneous unbounded domains with the elasticity modulus and mass density varying as power functions of spatial coordinates. The unbounded domains and the elasticity matrices are transformed to the scaled boundary coordinates. The scaled boundary finite‐element equation in displacement amplitudes are derived directly from the governing equations of elastodynamics. To enforce the radiation condition at infinity, an asymptotic expansion of the dynamic‐stiffness matrix for high frequency is developed. The dynamic‐stiffness matrix at lower frequency is obtained by numerical integration of ordinary differential equations. Only the boundary is discretized yielding a reduction of the spatial dimension by one. No fundamental solution is required. Material anisotropy is modelled without additional efforts. Examples of two‐ and three‐dimensional non‐homogeneous isotropic and transversely isotropic unbounded domains are presented. The results demonstrate the accuracy and simplicity of the scaled boundary finite‐element method. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
14.
瑞雷面波经常被用来反演地表浅层横波速度，受到越来越广泛的关注。对瑞雷波的研究一般都基于完全弹性介质，而实际地层更接近黏弹性介质，对黏弹性介质中的瑞雷面波进行模拟更具实际意义。本文采用广义标准线性体模型来描述黏弹性介质，并采用交错网格有限差分法对考虑水平自由表面的黏弹介质进行正演模拟，再与弹性介质中的结果进行对比分析。首先采用非线性最优化算法根据期望常数品质因子直接求取松弛时间来拟合常Q模型，并给出广义标准线性固体的具体算例，实施自由表面条件时采用声学-弹性边界近似法，通过剪切模量不变来考虑自由表面上、下横向应力保持连续的条件。对于非自由表面，采用非分裂的多轴卷积完全匹配层来吸收波场。然后对几种典型的数值模型进行正演模拟计算，数值解与解析解的对比验证了本文方法的准确性与有效性，正演结果的对比表明波场尤其是面波频散会受黏弹性影响，因此有必要在面波勘探中考虑黏弹性因素。 相似文献
15.
弹性波正演模拟中改进的非分裂式PML实现方法(英文) 总被引：6，自引：3，他引：3
在弹性波有限差分正演模拟中,吸收边界条件常用来吸收截断边界处引入的不期望边界反射,其中完全匹配层(PML)吸收边界条件被认为是目前最理想的吸收边界条件。但是PML吸收边界条件的传统实现却存在着很大不足:全局分裂式PML吸收边界条件实现简单但是需要占用太多内存;局部分裂式PML吸收边界条件需要考虑多个边界和角点区域,编程实现非常复杂;非分裂式PML吸收边界条件由于涉及卷积运算,计算量很大。本文基于非分裂式PML吸收边界条件,结合复频移伸展函数,提出了一种新的数值实现方法,其计算方程简单、占用内存小、编程实现容易,是对PML介质理论数值实现的改进和完善。 相似文献
16.
We developed a reverse‐time migration scheme that can image regions with rugged topography without requiring any approximations by adopting an irregular, unstructured‐grid modelling scheme. This grid, which can accurately describe surface topography and interfaces between high‐velocity‐contrast regions, is generated by Delaunay triangulation combined with the centroidal Voronoi tessellation method. The grid sizes vary according to the migration velocities, resulting in significant reduction of the number of discretized nodes compared with the number of nodes in the conventional regular‐grid scheme, particularly in the case wherein high near‐surface velocities exist. Moreover, the time sampling rate can be reduced substantially. The grid method, together with the irregular perfectly matched layer absorbing boundary condition, enables the proposed scheme to image regions of interest using curved artificial boundaries with fewer discretized nodes. We tested the proposed scheme using the 2D SEG Foothill synthetic dataset. 相似文献
17.
John P. Wolf 《地震工程与结构动力学》2002,31(1):15-32
The scaled boundary finite‐element method is a powerful semi‐analytical computational procedure to calculate the dynamic stiffness of the unbounded soil at the structure–soil interface. This permits the analysis of dynamic soil–structure interaction using the substructure method. The response in the neighbouring soil can also be determined analytically. The method is extended to calculate numerically the response throughout the unbounded soil including the far field. The three‐dimensional vector‐wave equation of elasto‐dynamics is addressed. The radiation condition at infinity is satisfied exactly. By solving an eigenvalue problem, the high‐frequency limit of the dynamic stiffness is constructed to be positive definite. However, a direct determination using impedances is also possible. Solving two first‐order ordinary differential equations numerically permits the radiation condition and the boundary condition of the structure–soil interface to be satisfied sequentially, leading to the displacements in the unbounded soil. A generalization to viscoelastic material using the correspondence principle is straightforward. Alternatively, the displacements can also be calculated analytically in the far field. Good agreement of displacements along the free surface and below a prism foundation embedded in a half‐space with the results of the boundary‐element method is observed. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
18.
在进行地震波模拟计算的过程中用有限的计算区域模拟地下无限空间,需要进行边界截断.为了在边界处不产生虚假反射影响模拟结果,需要引入吸收边界条件.本文采用的近似完全匹配层是一种新型非分裂完全匹配层,计算效率较高.同时相比于其他非分裂完全匹配层,其还具有不改变方程的形式、易于实现等优势.但是当入射波角度较大,边界吸收效果变弱,且残留在边界中的能量使近似完全匹配层变得极其不稳定.多轴复频移近似完全匹配层的提出就是为了改善对大角度入射波的吸收并且提高边界的稳定性.通过实验模拟和矩阵特征值灵敏度来研究多轴复频移近似完全匹配层的吸收效果及稳定性.结果表明该方法不仅能够吸收掠入波,而且对常规入射波的吸收也得到提升,同时拥有更好的稳定性. 相似文献
19.
传统的完全匹配层技术是一种能够较为有效地消除边界反射的边界条件,但是当表层为泊松比较高的自由表面时,该技术可能会产生不稳定的现象.针对传统的完全匹配层技术固有的不稳定和掠射情况下吸收效果不佳等缺陷,发展了多轴完全匹配层、卷积完全匹配层以及将两者结合的多轴卷积完全匹配层等3种边界条件.本文介绍了水平自由表面的不同处理方法以及传统、多轴、卷积和多轴卷积等4种完全匹配层条件的原理,通过二维半无限空间模型的交错网格有限差分正演模拟对比,分析了几种自由边界实施方法在这几种完全匹配层条件下的稳定性,并通过提取单道波形与解析解进行对比,定性分析了水平自由表面几种不同处理方法的准确性以及各自的适用条件. 结果表明,泊松比和水平自由表面实施方法对波场模拟效果及其稳定性有重要影响. 相似文献
20.
以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征. 相似文献