Parallel 3-D simulation of seismic wave propagation in heterogeneous anisotropic media: a grid method approach

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

	Parallel 3-D simulation of seismic wave propagation in heterogeneous anisotropic media: a grid method approach
	Gao, HW; Zhang, JF
刊名	GEOPHYSICAL JOURNAL INTERNATIONAL
	2006-06-01
卷号	165 期号:3 页码:875-888
关键词	numerical techniques seismic anisotropy seismic modelling synthetic seismograms seismic wave propagation topography
ISSN号	0956-540X
DOI	10.1111/j.1365-246X.2006.02946.x
文献子类	Article
英文摘要	This paper presents a parallel numerical technique for modelling wave propagation in 3-D heterogeneous anisotropic media. The scheme is developed by following a so-called 3-D grid method of the elastic-isotropic case. The proposed parallel algorithm needs small data exchanges between subdomains in contrast to that developed based on other numerical techniques; therefore, it is more suitable for a PC-Cluster. The algorithm is implemented on a mesh of mixed tetrahedrons and parallelepipedons, thus providing an accurate description of arbitrary 3-D surface and interface topographies and an easy generation of a non-uniform, unstructured mesh. The unstructured mesh means that the proposed algorithm can reduce the memory requirement by flexibly assigning small grid spacing in regions with low velocities and larger grid spacing in regions with higher velocities. Like the 3-D grid method, the resulting anisotropic scheme naturally satisfies the free-surface boundary conditions of arbitrary surface topography. As a result, the near-surface scattering effects can be more accurately modelled. The proposed scheme can handle a general anisotropy without any interpolations. In this paper, the transversely isotropic medium with a tilted symmetry axis, as typically caused by a system of parallel cracks or fine layers, is discussed in detail. A paraxial absorbing boundary condition in a 3-D general anisotropic case is also proposed. Comparisons with analytical solutions demonstrate the accuracy of the parallel algorithm. Computed 3-D radiation patterns illustrate shear-wave splitting, as predicted by the theory. We show the generality and flexibility of the algorithm by modelling wave propagation in an anisotropic half-space with a hemispherical crater on the surface and in mixed isotropic/anisotropic models with horizontal and inclined interfaces.
WOS关键词	ABSORBING BOUNDARY-CONDITIONS ; SPECTRAL ELEMENT METHOD ; ELASTIC-WAVES ; SOLIDS
WOS研究方向	Geochemistry & Geophysics
语种	英语
出版者	BLACKWELL PUBLISHING
WOS记录号	WOS:000238649000013
内容类型	期刊论文
源URL	[http://ir.iggcas.ac.cn/handle/132A11/66688]
专题	中国科学院地质与地球物理研究所
通讯作者	Zhang, JF
作者单位	Chinese Acad Sci, Inst Geol & Geophys, State Key Lab Lithosphere Evolut, Beijing 100029, Peoples R China
推荐引用方式 GB/T 7714	Gao, HW,Zhang, JF. Parallel 3-D simulation of seismic wave propagation in heterogeneous anisotropic media: a grid method approach[J]. GEOPHYSICAL JOURNAL INTERNATIONAL,2006,165(3):875-888.
APA	Gao, HW,&Zhang, JF.(2006).Parallel 3-D simulation of seismic wave propagation in heterogeneous anisotropic media: a grid method approach.GEOPHYSICAL JOURNAL INTERNATIONAL,165(3),875-888.
MLA	Gao, HW,et al."Parallel 3-D simulation of seismic wave propagation in heterogeneous anisotropic media: a grid method approach".GEOPHYSICAL JOURNAL INTERNATIONAL 165.3(2006):875-888.