Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Earth Sci.

Sec. Solid Earth Geophysics

This article is part of the Research TopicThe State-of-Art Techniques of Seismic Imaging for the Deep and Ultra-deep Hydrocarbon Reservoirs - Volume IIIView all 8 articles

Elastic wave equation numerical simulation with time-space domain dispersion-relation-based staggered-grid finite-difference method

Provisionally accepted
Xuewen  ShiXuewen Shi1,2Chang  WangChang Wang1,2,3Wei  LiuWei Liu4,5*Dongjun  ZhangDongjun Zhang1,2Hu  ZiDuoHu ZiDuo5,6Yanwen  FengYanwen Feng1,2Yang  GaoYang Gao1,2Zhiwen  ZhouZhiwen Zhou1,2
  • 1Shale Gas Research Institute, PetroChina Southwest Oil & Gas Field Company, Chengdu, China
  • 2Sichuan Provincial Key Laboratory of Shale Gas Evaluation and Exploitation, Chengdu, China
  • 3BGP Inc., China National Petroleum Corporation, Zhuozhou, China
  • 4Research Institute of Petroleum Exploration & Development-Northwest, PetroChina, Lanzhou, China
  • 5PetroChina Key Laboratory of Reservoir Description, Lanzhou, China
  • 6PetroChina Research Institute of Petroleum Exploration & Development-Northwest, Lanzhou, China

The final, formatted version of the article will be published soon.

Compared to the conventional high-order staggered-grid finite-difference method (C-SFD), the 3 time-space domain dispersion-relation-based high-order staggered-grid finite-difference method 4 (TS-SFD) is able to suppress the numerical dispersion more efficiently to achieve higher modelling 5 accuracy when applied to numerically solving the velocity-stress acoustic equation. The enhanced 6 modelling accuracy of TS-SFD is attributed to the difference coefficients calculating approach 7 based on the time-space domain dispersion-relation, which makes the difference coefficients 8 change adaptively with the propagation velocity of seismic wave in the medium. However, when 9 numerical simulation of the velocity-stress elastic wave equation is conducted with TS-SFD, if the 10 difference coefficients are calculated based on the time-space domain P-wave dispersion-relation, 11 the modelling accuracy of P-wave is high while the modelling accuracy of S-wave is low, and 12 vice versa. In order to address the limitation of TS-SFD in ensuring high modelling accuracy 13 for both Pand S-wave, in this paper we propose a novel strategy to conduct elastic wave 14 numerical simulation with TS-SFD, in which the decoupled Pand S-wave elastic wave equation 15 are numerically solved with TS-SFD, and the difference coefficients are calculated based on the 16 time-space domain Pand S-wave dispersion-relation respectively. Numerical dispersion analysis 17 and numerical simulation examples show that the simulation strategy proposed in this paper 18 can ensure high simulation accuracy for both Pand S-wave. Stability analysis shows that this 19 simulation strategy can effectively improve the stability of elastic wave numerical simulation with 20 TS-SFD. In addition, the simulation strategy has the added advantage of automatically separating 21 the Pand S-wave, which provides a solid foundation for the subsequent analysis of Pand 22 S-wave propagation characteristics in elastic media and elastic wave reverse time migration.

Keywords: Elastic wave, numerical modelling, numerical dispersion, Time-space domain, Difference coefficient calculation

Received: 29 Jul 2025; Accepted: 30 Nov 2025.

Copyright: © 2025 Shi, Wang, Liu, Zhang, ZiDuo, Feng, Gao and Zhou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

* Correspondence: Wei Liu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.