Three-dimensional elastic wave numerical modelling in the presence of   surface topography by a collocated-grid finite-difference method on   curvilinear grids

doi:10.1111/j.1365-246X.2012.05472.x

CORC > 北京大学 > 地球与空间科学学院

	Three-dimensional elastic wave numerical modelling in the presence of surface topography by a collocated-grid finite-difference method on curvilinear grids
	Zhang, Wei ; Zhang, Zhenguo ; Chen, Xiaofei
刊名	geophysical journal international
	2012
关键词	Computational seismology Numerical solutions Wave propagation ADAPTIVE MESH REFINEMENT PLANAR FREE-SURFACE DISCONTINUOUS GALERKIN METHOD SPECTRAL ELEMENT METHOD IRREGULAR FREE-SURFACE HIGH-ORDER PSEUDOSPECTRAL METHOD VELOCITY-STRESS HETEROGENEOUS MEDIA SEISMIC MOTION
DOI	10.1111/j.1365-246X.2012.05472.x
英文摘要	Surface topography has been considered a difficult task for seismic wave numerical modelling by the finite-difference method (FDM) because the most popular staggered finite-difference scheme requires a rectilinear grid. Even though there are numerous collocated grid schemes in other computational fields that could be used to solve the first-order hyperbolic equations, the lack of a stable free-surface boundary condition implementation for curvilinear grids also obstructs the adoption of curvilinear grids in seismic wave FDM modelling. In this study, we use generalized curvilinear grids that can fit the surface topography to discretize the computational domain and describe the implementation of a collocated grid finite-difference scheme, a higher order MacCormack scheme, to solve the first-order hyperbolic velocity-stress equations on the curvilinear grid. To achieve a sufficiently accurate and stable free-surface boundary condition implementation on the curvilinear grids, we propose the traction image method that antisymmetrically images the traction components instead of the stress components to the ghost points above the free surface. Since the velocity derivatives at the free surface are provided by the free-surface condition, we use a compact scheme to compute the velocity derivatives near the free surface and avoid the use of velocity values on the ghost points. Numerical tests verify that using the curvilinear grid, the collocated finite-difference scheme and the traction image technique can simulate seismic wave propagation in the presence of surface topography with sufficient accuracy.; Geochemistry & Geophysics; SCI(E); EI; 15; ARTICLE; 1; 358-378; 190
语种	英语
内容类型	期刊论文
源URL	[http://ir.pku.edu.cn/handle/20.500.11897/311150]
专题	地球与空间科学学院
推荐引用方式 GB/T 7714	Zhang, Wei,Zhang, Zhenguo,Chen, Xiaofei. Three-dimensional elastic wave numerical modelling in the presence of surface topography by a collocated-grid finite-difference method on curvilinear grids[J]. geophysical journal international,2012.
APA	Zhang, Wei,Zhang, Zhenguo,&Chen, Xiaofei.(2012).Three-dimensional elastic wave numerical modelling in the presence of surface topography by a collocated-grid finite-difference method on curvilinear grids.geophysical journal international.
MLA	Zhang, Wei,et al."Three-dimensional elastic wave numerical modelling in the presence of surface topography by a collocated-grid finite-difference method on curvilinear grids".geophysical journal international (2012).