共查询到20条相似文献，搜索用时 15 毫秒
1.
New alternative formulations of exact boundary conditions for arbitrary three-dimensional (3D) free-surface topographies on seismic media have been derived. They are shown to be equivalent to previously published formulations, thereby verifying the validity of each set of formulations. The top of a curved grid represents the free-surface topography while the interior of the grid represents the physical medium. We assume the velocity–stress version of the viscoelastic wave equations to be valid in this grid before transforming the equations to a rectangular grid. In order to perform the numerical discretization we apply the latter version of the equations for seismic wave propagation simulation in the medium. The numerical discretization of the free-surface topography boundary conditions by second-order finite differences (FDs) is shown, as well as the spatially unconditional stability of the resulting system of equations. The FD order is increased by two for each point away from the free surface up to eight, which is the order used in the interior. We use staggered grids in both space and time and the second-order leap-frog and Crank– Nicholson methods for wavefield time propagation. An application using parameters typical of teleseismic earthquakes and explosions is presented using a 200 × 100 km2 area of real topography from southwestern Norway over a homogeneous medium. A dipping plane wave simulates a teleseismic P-wave incident on the surface topography. Results show conversion from P- to Rg- (short period fundamental mode Rayleigh) waves in the steepest and/or roughest topography, as well as attenuated waves in valleys and fjords. The codes are parallelized for simulation on fast supercomputers and PC-clusters to model high frequencies and/or large areas. 相似文献
2.
We have pursued two-dimensional (2D) finite-difference (FD) modelling of seismic scattering from free-surface topography. Exact free-surface boundary conditions for the particle velocities have been derived for arbitrary 2D topographies. The boundary conditions are combined with a velocity–stress formulation of the full viscoelastic wave equations. A curved grid represents the physical medium and its upper boundary represents the free-surface topography. The wave equations are numerically discretized by an eighth-order FD method on a staggered grid in space, and a leap-frog technique and the Crank–Nicholson method in time.
In order to demonstrate the capabilities of the surface topography modelling technique, we simulate incident point sources with a sinusoidal topography in seismic media of increasing complexities. We present results using parameters typical of exploration surveys with topography and heterogeneous media. Topography on homogeneous media is shown to generate significant scattering. We show additional effects of layering in the medium, with and without randomization, using a von Kármán realization of apparent anisotropy. Synthetic snapshots and seismograms indicate that prominent surface topography can cause back-scattering, wave conversions and complex wave patterns which are usually discussed in terms of inter-crust heterogeneities. 相似文献
In order to demonstrate the capabilities of the surface topography modelling technique, we simulate incident point sources with a sinusoidal topography in seismic media of increasing complexities. We present results using parameters typical of exploration surveys with topography and heterogeneous media. Topography on homogeneous media is shown to generate significant scattering. We show additional effects of layering in the medium, with and without randomization, using a von Kármán realization of apparent anisotropy. Synthetic snapshots and seismograms indicate that prominent surface topography can cause back-scattering, wave conversions and complex wave patterns which are usually discussed in terms of inter-crust heterogeneities. 相似文献
3.
4.
A 2D numerical finite-difference algorithm accounting for surface topography is presented. Higher-order, dispersion-bounded, cost-optimized finite-difference operators are used in the interior of the numerical grid, while non-reflecting absorbing boundary conditions are used along the edges. Transformation from a curved to a rectangular grid achieves the modelling of the surface topography. We use free-surface boundary conditions along the surface. In order to obtain complete modelling of the effects of wave propagation, it is important to account for the surface topography, otherwise near-surface effects, such as scattering, are not modelled adequately. Even if other properties of the medium, for instance randomization, can improve numerical simulations, inclusion of the surface topography makes them more realistic. 相似文献
5.
本文以基于改进BISQ模型的二维双相各向同性介质一阶速度-应力方程为基础，推导出了曲线坐标系下对应的方程，然后采用低频散、低耗散的同位网格MacCormack有限差分法来离散方程，并采用紧致的单边MacCormack差分格式结合牵引力镜像法来施加自由地表边界条件，实现了地震波场数值模拟.曲线网格有限差分法采用贴体网格来描述自由表面，地表的网格线紧贴地形，避免了台阶近似造成的数值散射.数值模拟结果表明，在双相介质起伏自由地表和分界面处，各类波型复杂的反射透射规律可以清晰展现，曲线网格有限差分法可以精确地解决地震波在含起伏地表的双相各向同性介质中的传播问题. 相似文献
6.
In this paper, we deduced the corresponding first-order velocity–stress equation for curvilinear coordinates from the first-order velocity–stress equation based on the modified Biot/squirt model for a two-dimensional two-phase medium. The equations are then numerically solved by an optimized high-order non-staggered finite difference scheme, that is, the dispersion relation preserving/optimization MacCormack scheme. To implement undulating free-surface topography, we derive an analytical relationship between the derivatives of the particle velocity components and use the compact finite-difference scheme plus a traction-image method. In the undulating free surface and the undulating subsurface interface of two-phase medium, the complex reflected wave and transmitted wave can be clearly recognized in the numerical simulation results. The simulation results show that the curvilinear-grid finite-difference method, which uses a body-conforming grid to describe the undulating surface, can accurately reduce the numerical scattering effect of seismic wave propagation caused by the use of ladder-shaped grid to fit the surfaces when undulating topography is present in a two-phase isotropic medium. 相似文献
7.
Seismic wavefield modeling in media with fluid-filled fractures and surface topography 总被引：3，自引：0，他引：3
We present a finite difference (FD) method for the simulation of seismic wave fields in fractured medium with an irregular (non-flat) free surface which is beneficial for interpreting exploration data acquired in mountainous regions. Fractures are introduced through the Coates-Schoenberg approach into the FD scheme which leads to local anisotropic properties of the media where fractures are embedded. To implement surface topography, we take advantage of the boundary-conforming grid and map a rectangular grid onto a curved one. We use a stable and explicit second-order accurate finite difference scheme to discretize the elastic wave equations (in a curvilinear coordinate system) in a 2D heterogeneous transversely isotropic medium with a horizontal axis of symmetry (HTI). Efficiency tests performed by different numerical experiments clearly illustrate the influence of an irregular free surface on seismic wave propagation in fractured media which may be significant to mountain seismic exploration. The tests also illustrate that the scattered waves induced by the tips of the fracture are re-scattered by the features of the free surface topography. The scattered waves provoked by the topography are re-scattered by the fractures, especially Rayleigh wave scattering whose amplitudes are much larger than others and making it very difficult to identify effective information from the fractures. 相似文献
8.
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. 相似文献
9.
Application of a perfectly matched layer in seismic wavefield simulation with an irregular free surface 总被引：2，自引：0，他引：2 
下载免费PDF全文
下载免费PDF全文 Haiqiang Lan Jingyi Chen Zhongjie Zhang Youshan Liu Jianguo Zhao Ruiqi Shi 《Geophysical Prospecting》2016,64(1):112-128
Recently, an effective and powerful approach for simulating seismic wave propagation in elastic media with an irregular free surface was proposed. However, in previous studies, researchers used the periodic condition and/or sponge boundary condition to attenuate artificial reflections at boundaries of a computational domain. As demonstrated in many literatures, either the periodic condition or sponge boundary condition is simple but much less effective than the well‐known perfectly matched layer boundary condition. In view of this, we intend to introduce a perfectly matched layer to simulate seismic wavefields in unbounded models with an irregular free surface. We first incorporate a perfectly matched layer into wave equations formulated in a frequency domain in Cartesian coordinates. We then transform them back into a time domain through inverse Fourier transformation. Afterwards, we use a boundary‐conforming grid and map a rectangular grid onto a curved one, which allows us to transform the equations and free surface boundary conditions from Cartesian coordinates to curvilinear coordinates. As numerical examples show, if free surface boundary conditions are imposed at the top border of a model, then it should also be incorporated into the perfectly matched layer imposed at the top‐left and top‐ right corners of a 2D model where the free surface boundary conditions and perfectly matched layer encounter; otherwise, reflections will occur at the intersections of the free surface and the perfectly matched layer, which is confirmed in this paper. So, by replacing normal second derivatives in wave equations in curvilinear coordinates with free surface boundary conditions, we successfully implement the free surface boundary conditions into the perfectly matched layer at the top‐left and top‐right corners of a 2D model at the surface. A number of numerical examples show that the perfectly matched layer constructed in this study is effective in simulating wave propagation in unbounded media and the algorithm for implementation of the perfectly matched layer and free surface boundary conditions is stable for long‐time wavefield simulation on models with an irregular free surface. 相似文献
10.
基于应力、速度混合变量弹性波方程及任意四边形网格差分算子,给出了交错计算应力及速度的非规则网格弹性波应力一速度差分法该方法融合了有限元法能适应复杂形状边界及差分法无需计算刚度阵的特点,具有较高的计算精度,所需计算机存储空间较少,计算效率也很高．基于积分平衡方程引入了任意形状自由表面的边界条件,且通过局部滤波改善了自由表面边界条件的稳定性,使得该方法可应用于考虑地表形状影响的地震波数值模拟 相似文献
11.
We present an original implementation of the free-surface boundary condition in a mesh-free finite-difference method for simulating elastic wave propagation in the frequency domain. For elastic wave modelling in the frequency domain, the treatment of free surfaces is a key issue which requires special consideration. In the present study, the free-surface boundary condition is directly implemented at node positions located on the free-surface. Flexible nature of the mesh-free method for nodal distribution enables us to introduce topography into numerical models in an efficient manner. We investigate the accuracy of the proposed implementation by comparing numerical results with an analytical solution. The results show that the proposed method can calculate surface wave propagation even for an inclined free surface with substantial accuracy. Next, we calculate surface wave propagation in a model with a topographic surface using our method, and compare the numerical result with that using the finite-element method. The comparison shows the excellent agreement with each other. Finally, we apply our method to the SEG foothill model to investigate the effectiveness of the proposed method. Since the mesh-free method has high flexibility of nodal distribution, the proposed implementation would deal with models of topographic surface with sufficient accuracy and efficiency. 相似文献
12.
本文对基于弹性波动理论的多波多分量高斯束偏移进行了完整且详细的分析和公式推导,实现了3D空间多分量(矢量)波场的直接成像.由于当前多数基于弹性波动方程的偏移方法只是假设应力边界条件为自由地表边界条件,这种假设不符合垂直地震剖面(VSP)和海底电缆(OBC)等地震数据.为此本文详细分析了实际应用中常见的三种弹性各向同性介质模型的应力边界条件:自由空间、海底和自由地表模型.在上行传播假设情况下,获得了应力边界条件与位移边界条件的关系式.在此基础上,准确推导了3D多波多分量高斯束波场延拓和偏移成像公式,并在偏移过程中实现了完整的多波型自动分离.由于常规的互相关成像条件不适用于矢量波场成像,本文引用了散度/旋度互相关成像条件.通过约定PS转换波的正向旋转方向解决了3D空间PS成像极性翻转问题.利用2D和3D模型数据偏移成像验证了我们所提出的多波多分量高斯束偏移方法的可行性. 相似文献
13.
A finite-volume approach for 2D magnetotellurics modeling with arbitrary topographies 总被引：2，自引：0，他引：2
A novel finite-volume approach for complicated 2D magnetotellurics (MT) problems with arbitrarily surface topography is presented. An edge-surface integral balance equation is derived by employing a conservation law on the generalized 2D MT boundary value problem. A triangular grid is used to discretize the 2D conductivity model so that we can deal with arbitrarily complex cases with surface topography. The node-centered finite-volume algorithm is used to derive the final system of linear equations on a dual mesh of the triangular grid, which is solved by a robust direct solver. Three synthetic models verify the accuracy of the presented finite-volume algorithm and its capability of dealing with surface topography. 相似文献
14.
Christopher Juhlin 《Geophysical Prospecting》1995,43(6):843-858
The velocity-stress formulation for propagation of elastic seismic waves through 2D heterogeneous transversely isotropic media of arbitrary orientation is presented. The equations are recast into a finite-difference scheme and solved numerically using fourth-order spatial operators and a second-order temporal operator on a staggered grid. Absorbing, free-surface and symmetry boundary conditions have been implemented. Test cases compare well with other published solutions. Synthetic seismograms are calculated over two idealized models: (i) vertical fractures in granite with a dolerite sill reflector and (ii) a dipping anisotropic shale. Comparisons with the isotropic counterparts show significant differences which may have to be accounted for in seismic processing in the future. 相似文献
15.
本文基于Kjartansson常Q模型理论，推导了常Q衰减介质中黏声波和黏弹性波的速度-应力方程，并采用基于二项式窗函数的优化交错网格有限差分方法进行了数值模拟，同时引入不分裂的复频移卷积完全匹配层(CPML)吸收边界条件，以消除边界反射.使用基于自适应时间步长记忆方法的中心差分近似时间分数阶导数，与常用的短时记忆方法相比，提高了波动方程的离散化精度和计算效率.通过对比均匀模型下声波的数值解与解析解，验证了算法的精确性，并进一步分析了不同品质因子下地震波的频散及衰减特征.对BP盐丘模型的数值模拟结果可以较好地反映本文数值方法对复杂介质的适应性及频散压制效果. 相似文献
16.
传统的完全匹配层技术是一种能够较为有效地消除边界反射的边界条件,但是当表层为泊松比较高的自由表面时,该技术可能会产生不稳定的现象.针对传统的完全匹配层技术固有的不稳定和掠射情况下吸收效果不佳等缺陷,发展了多轴完全匹配层、卷积完全匹配层以及将两者结合的多轴卷积完全匹配层等3种边界条件.本文介绍了水平自由表面的不同处理方法以及传统、多轴、卷积和多轴卷积等4种完全匹配层条件的原理,通过二维半无限空间模型的交错网格有限差分正演模拟对比,分析了几种自由边界实施方法在这几种完全匹配层条件下的稳定性,并通过提取单道波形与解析解进行对比,定性分析了水平自由表面几种不同处理方法的准确性以及各自的适用条件. 结果表明,泊松比和水平自由表面实施方法对波场模拟效果及其稳定性有重要影响. 相似文献
17.
During seismic wave propagation on a free surface, a strong material contrast boundary develops in response to interference by P- and S- waves to create a surfacewave phenomenon. To accurately determine the effects of this interface on surface-wave propagation, the boundary conditions must be accurately modeled. In this paper, we present a numerical approach based on the dynamic poroelasticity for a space–time-domain staggeredgrid finite-difference simulation in porous media that contain a free-surface boundary. We propose a generalized stess mirror formulation of the free-surface boundary for solids and fluids in porous media for the grid mesh on which lays the free-surface plane. Its analog is that used for elastic media, which is suitable for precise and stable Rayleigh-type surface-wave modeling. The results of our analysis of first kind of Rayleigh (R1) waves obtained by this model demonstrate that the discretization of the mesh in a similar way to that for elastic media can realize stable numerical solutions with acceptable precision. We present numerical examples demonstrating the efficiency and accuracy of our proposed method. 相似文献
18.
We recently proposed an efficient hybrid scheme to absorb boundary reflections for acoustic wave modelling that could attain nearly perfect absorptions. This scheme uses weighted averaging of wavefields in a transition area, between the inner area and the model boundaries. In this paper we report on the extension of this scheme to 2D elastic wave modelling with displacement‐stress formulations on staggered grids using explicit finite‐difference, pseudo‐implicit finite‐difference and pseudo‐spectral methods. Numerical modelling results of elastic wave equations with hybrid absorbing boundary conditions show great improvement for modelling stability and significant absorption for boundary reflections, compared with the conventional Higdon absorbing boundary conditions, demonstrating the effectiveness of this scheme for elastic wave modelling. The modelling results also show that the hybrid scheme works well in 2D rotated staggered‐grid modelling for isotropic medium, 2D staggered‐grid modelling for vertically transversely isotropic medium and 2D rotated staggered‐grid modelling for tilted transversely isotropic medium. 相似文献
19.
本文基于显式有限元法研究了地震波垂直入射时三维复杂山谷地形对地震地面运动的影响,在数值分析中应用了三维化二维的解法和黏弹性人工边界的处理方法,实现了地震波垂直输入下三维复杂场地地震动数值模拟,并验证了该方法的合理性.以四川桃坪地区一山谷地形作为研究对象,基于地表高程数据分别建立了二维和三维场地模型,对比研究表明:在复杂地形情况下考虑二、三维模型时具有明显差异,三维模型能更真实地反映地形变化对地震动的影响,复杂地形条件下有必要考虑三维实际场地模型.本文对边界自由场的处理方法也可用于处理三维复杂场地地震动斜入射问题,为三维复杂地形场地地震效应研究提供参考. 相似文献
20.
研究了三维弱各向异性近似下,利用伪P波(伪纵波)模拟弹性波场P分量在倾斜对称轴的横向各向同性(TTI)介质中的传播过程,并对比了分别基于弹性Hooke定律、弹性波投影和运动学色散方程所建立的三种二阶差分伪P波方程的正演特点.目前这些伪P波方程数值计算主要采用规则网格差分,但是规则网格在TTI模拟中有低效率、低精度以及不稳定的缺点.为了提高计算的精度,本文构建出相应方程的交错网格有限差分格式.通过对比伪P波方程在三维TTI介质中不同的数值模拟的表达形式,本文认为基于色散方程所建立的伪P波方程在模拟弹性波中P波传播的过程中具有最小的噪声.本文分析不同的各向同性对称轴空间角度的频散特征,并引入适当的横波速度维持计算的稳定.二维模型算例表明,本文提出的交错网格正演算法可以得到稳定光滑的伪P波正演波场.使用本文交错网格算法对二维BP TTI模型的逆时偏移也具有较稳定的偏移结果. 相似文献