دورية أكاديمية

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

التفاصيل البيبلوغرافية
العنوان: 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
المساهمون: Zhang, W (reprint author), GeoTomo LLC, Houston, TX USA., Univ Sci & Technol China, Sch Earth & Space Sci, Hefei 230026, Anhui, Peoples R China., Peking Univ, Sch Earth & Space Sci, Beijing 100871, Peoples R China., GeoTomo LLC, Houston, TX USA.
المصدر: SCI ; EI
بيانات النشر: geophysical journal international
سنة النشر: 2012
المجموعة: Peking University Institutional Repository (PKU IR) / 北京大学机构知识库
مصطلحات موضوعية: 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
الوصف: 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
نوع الوثيقة: journal/newspaper
اللغة: English
تدمد: 0956-540X
العلاقة: GEOPHYSICAL JOURNAL INTERNATIONAL.2012,190,(1),358-378.; 874283; http://hdl.handle.net/20.500.11897/311150Test; WOS:000305080900026
DOI: 10.1111/j.1365-246X.2012.05472.x
الإتاحة: https://doi.org/20.500.11897/311150Test
https://doi.org/10.1111/j.1365-246X.2012.05472.xTest
https://hdl.handle.net/20.500.11897/311150Test
رقم الانضمام: edsbas.81DA0F4
قاعدة البيانات: BASE
الوصف
تدمد:0956540X
DOI:10.1111/j.1365-246X.2012.05472.x