none 10.1190/geo2023-0324.1 Society of Exploration Geophysicists Society of Exploration Geophysicists 186 171498096 5981 2025121212264280500 10.1190 2025-12-12T17:26:53Z 2025-02-06T09:37:32Z 0 Geophysics 0016-8033 1942-2156 01 05 2025 90 3 Moosoo Won Khalifa University of Sciences and Technology 1 , Earth Sciences, Abu Dhabi, UAE and , Earthquake Research Center, Daejeon, South Korea. 100058280@ku.ac.ae ; moosoowon@kigam.re.kr (corresponding author). Korea Institute of Geoscience and Mineral Resources 1 , Earth Sciences, Abu Dhabi, UAE and , Earthquake Research Center, Daejeon, South Korea. 100058280@ku.ac.ae ; moosoowon@kigam.re.kr (corresponding author). http://orcid.org/0000-0001-7081-9757 Bing Zhou Khalifa University of Sciences and Technology 2 , Earth Sciences, Abu Dhabi, UAE. bing.zhou@ku.ac.ae . http://orcid.org/0000-0002-6813-5265 Xu Liu King Fahd University of Petroleum and Minerals 3 , College of Petroleum Engineering and Geosciences, Dhahran, Saudi Arabia. xu.liu@kfupm.edu.sa . http://orcid.org/0000-0003-2835-2844 Mohamed Jamal Zemerly Khalifa University of Sciences and Technology 4 , Electrical Engineering and Computer Science Department, Abu Dhabi, UAE. jamal.zemerly@ku.ac.ae . http://orcid.org/0000-0003-1845-5946 Mohammad Al-Khaleel Khalifa University of Sciences and Technology 5 , Mathematics Department, Abu Dhabi, UAE and , Mathematics Department, Irbid, Jordan. mohammad.alkhaleel@ku.ac.ae . Yarmouk University 5 , Mathematics Department, Abu Dhabi, UAE and , Mathematics Department, Irbid, Jordan. mohammad.alkhaleel@ku.ac.ae . http://orcid.org/0000-0001-6266-373X Mohamed Kamel Riahi Khalifa University of Sciences and Technology 6 , Mathematics Department, Abu Dhabi, UAE. mohamed.riahi@ku.ac.ae . http://orcid.org/0000-0001-7987-1391 ABSTRACT Accurately modeling seismic wave propagation in complex subsurface geologic structures is crucial for understanding their properties. However, 3D seismic wave modeling can be computationally demanding, requiring a significant amount of computer memory and runtime. To address this issue, a more efficient and accurate approach known as 2.5D modeling can be used when the subsurface geologic structure is two dimensions. We develop a new numerical method for 2.5D time-domain viscoelastic wave modeling characterized by the subdomain Chebyshev differentiation for accurate spatial derivatives, a Taylor-series recursive approach for accurate computations of temporal convolutions, a novel transform of the complex-domain computations into real-domain implementation, and fully parallel computing for high computational efficiency. Comparing our results with 3D analytical and numerical reference solutions, we demonstrate that our method offers satisfactory accuracies, excellent computational efficiencies, and a powerful capability of modeling 3D wavefields in complex 2D heterogeneous viscoacoustic and viscoelastic geologic models. Overall, our findings provide a robust and efficient approach for modeling seismic wave propagation in complex 2D subsurface geologic structures. 01 05 2025 T53 T67 Korea Institute of Geoscience and Mineral Resources http://dx.doi.org/10.13039/501100003700 IP2020-018 College of Petroleum Engineering and Geosciences, King Fahd University of Petroleum and Minerals http://dx.doi.org/10.13039/501100020361 H2FC2311 1 10.1190/crossmark-policy library.seg.org true 2023-05-24 2025-01-10 2025-01-16 2025-02-06 Korea Institute of Geoscience and Mineral Resources http://dx.doi.org/10.13039/501100003700 IP2020-018 College of Petroleum Engineering and Geosciences, King Fahd University of Petroleum and Minerals http://dx.doi.org/10.13039/501100020361 H2FC2311 10.1190/geo2023-0324.1 https://pubs.geoscienceworld.org/geophysics/article/90/3/T53/653607/An-accurate-and-efficient-numerical-method-for-2 https://pubs.geoscienceworld.org/seg/geophysics/article-pdf/90/3/T53/7170083/geo-2023-0324.1.pdf https://pubs.geoscienceworld.org/seg/geophysics/article-pdf/90/3/T53/7170083/geo-2023-0324.1.pdf https://pubs.geoscienceworld.org/geophysics/article/90/3/t53/653607/an-accurate-and-efficient-numerical-method-for-2 Quantitative seismology Aki 2002 10.1190/1.3255610 Aldridge D. F. Preston L. A. , 2009, Dispersion and attenuation for the anelastic velocity-memory-stress system: 79th Annual International Meeting, SEG, Expanded Abstracts, 3584–3589, doi: http://dx.doi.org/10.1190/1.3255610. 10.1190/1.1817065 Aldridge D. F. Symons N. P. , 2002, Finite integro-difference algorithm for 3D an acoustic wave propagation: 72nd Annual International Meeting, SEG, Expanded Abstracts, 1911–1914, doi: http://dx.doi.org/10.1190/1.1817065. 0016-8033 Geophysics Auer 78 6 R235 2013 10.1190/geo2012-0382.1 A critical appraisal of asymptotic 3D-to-2D data transformation in full-waveform seismic crosshole tomography 0016-8033 Geophysics Bai 82 4 WA83 2017 10.1190/geo2016-0596.1 Waveform inversion for attenuation estimation in anisotropic media 0016-8033 Geophysics Bielak 76 6 T131 2011 10.1190/geo2011-0019.1 Memory-efficient displacement-based internal friction for wave propagation simulation 0956-540X Geophysical Journal International Blanc 205 414 2016 10.1093/gji/ggw024 Highly accurate stability-preserving optimization of the Zener viscoelastic model, with application to wave propagation in the presence of strong attenuation 0016-8033 Geophysics Blanch 60 176 1995 10.1190/1.1443744 Modeling of a constant Q: Methodology and algorithm for an efficient and optimally inexpensive viscoelastic technique 0016-8025 Geophysical Prospecting Bleistein 34 686 1986 10.1111/j.1365-2478.1986.tb00488.x Two-and-one-half dimesional in plane wave propagation Numerical methods for ordinary differential equations Butcher 2016 10.1002/9781119121534 Wave fields in real media: Wave propagation in anisotropic, anelastic, porous and electromagnetic media Carcione 2007 0956-540X Geophysical Journal International Carcione 95 597 1988 10.1111/j.1365-246X.1988.tb06706.x Wave propagation simulation in a linear viscoelastic medium 0016-8033 Geophysics Carcione 53 769 1988 10.1190/1.1442512 Viscoacoustic wave propagation simulation in the Earth A time-domain finite-difference method with attenuation by a recursive algorithm Cheng 1996 0016-8033 Geophysics Chichinina 74 5 WB165 2009 10.1190/1.3173806 Attenuation anisotropy in the linear-slip model: Interpretation of physical modeling data 0956-540X Geophysical Journal International Day 78 105 1984 10.1111/j.1365-246X.1984.tb06474.x Numerical simulation of attenuated wavefields using a Padé approximant method 10.1190/segam2014-1240.1 Doyon B. Giroux B. , 2014, Practical aspects of 2.5D frequency-domain finite-difference modelling of viscoelastic waves: 84th Annual International Meeting, SEG, Expanded Abstracts, 3482–3486, doi: http://dx.doi.org/10.1190/segam2014-1240.1. 0016-8033 Geophysics Emmerich 52 1252 1987 10.1190/1.1442386 Incorporation of attenuation into time-domain computations of seismic wave field 0016-8033 Geophysics Ernst 72 5 J53 2007 10.1190/1.2761848 Application of a new 2D time-domain full-waveform inversion scheme to crosshole radar data 0956-540X Geophysical Journal International Fabien-Ouellet 209 1718 2017 10.1093/gji/ggx110 Time domain viscoelastic full waveform inversion 0956-540X Geophysical Journal International Fan 206 1539 2016 10.1093/gji/ggw228 Two-dimensional time-domain finite-difference modeling for viscoelastic seismic wave propagation 0956-540X Geophysical Journal International Fichtner 198 1878 2014 10.1093/gji/ggu228 Models and Fréchet kernels for frequency-(in)dependent Q 0956-540X Geophysical Journal International Forbriger 198 1387 2014 10.1093/gji/ggu199 Line-source simulation for shallow-seismic data — Part 1: Theoretical background 0956-540X Geophysical Journal International Furumura 124 820 1996 10.1111/j.1365-246X.1996.tb05640.x 2.5-D modelling of elastic waves using the pseudospectral method 0016-8033 Geophysics Guo 84 6 T381 2019 10.1190/geo2018-0747.1 Modeling the viscoelastic effects in P-waves with modified viscoacoustic wave propagation 0956-540X Geophysical Journal International Hao 219 1939 2019 The generalized standard-linear-solid model and the corresponding viscoacoustic wave equations revisited 0016-8033 Geophysics Hao 86 4 T239 2021 10.1190/geo2020-0548.1 Nearly constant Q models of the generalized standard linear solid type and the corresponding wave equations 0956-540X Geophysical Journal International Hestholm 139 852 1999 10.1046/j.1365-246x.1999.00994.x Three-dimensional finite difference viscoelastic wave modelling including surface topography 0956-540X Geophysical Journal International Igel 136 559 1999 10.1046/j.1365-246x.1999.00758.x Wave propagation in three-dimensional spherical sections by the Chebyshev spectral method Applied Geophysics Kai 17 561 2021 10.1007/s11770-020-0838-2 Characterization for wave equations in viscoelastic media based on the constant Q property 0018-926X IEEE Transactions on Antennas and Propagation Kelley 44 792 1996 10.1109/8.509882 Piecewise linear recursive convolution for dispersive media using FDTD 0016-8033 Geophysics Komatitsch 72 5 SM155 2007 10.1190/1.2757586 An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation 0016-8033 Geophysics Kosloff 55 734 1990 10.1190/1.1442885 Solution of the equations of dynamic elasticity by a Chebychev spectral method 10.1190/IGC2018-232 Lin B. Zhou B. , 2018, Subdomain Chebyshev spectral method for 2.5-D seismic wave modeling in VTI media: CPS/SEG International Geophysical Conference. Canadian Journal of Exploration Geophysics Liu 29 177 1993 Velocity and attenuation anisotropy caused by microcracks and macrofractures in a multiazimuth reverse VSP 0956-540X Geophysical Journal International Liu 186 1285 2011 10.1111/j.1365-246X.2011.05106.x 2.5-D poroelastic wave modelling in double porosity media 0956-540X Geophysical Journal International Liu 207 1472 2016 10.1093/gji/ggw345 Generalized poroviscoelastic model based on effective Biot theory and its application to borehole guided wave analysis 0018-9375 IEEE Transactions on Electromagnetic Compatibility Luebbers 32 222 1990 10.1109/15.57116 A frequency-dependent finite-difference time-domain formulation for dispersive materials 0956-540X Geophysical Journal International Martin 179 333 2009 10.1111/j.1365-246X.2009.04278.x An unsplit convolutional perfectly matched layer technique improved at grazing incidence for the viscoelastic wave equation 0029-5981 International Journal for Numerical Methods in Engineering Matzen 88 951 2011 10.1002/nme.3205 An efficient finite element time-domain formulation for the elastic second-order wave equation: A non-split complex frequency shifted convolutional PML 0094-8276 Geophysical Research Letters Moczo 32 2005 10.1029/2004GL021598 On the rheological models used for time-domain methods of seismic wave propagation Advances in geophysics Moczo 421 2007 0016-8025 Geophysical Prospecting Novais 53 523 2005 10.1111/j.1365-2478.2005.00488.x 2.5D finite-difference solution of the acoustic wave equation 0016-8033 Geophysics Pica 55 284 1990 10.1190/1.1442836 Nonlinear inversion of seismic reflection data in a laterally invariant medium Archives of Computational Methods in Engineering Pled 29 471 2022 10.1007/s11831-021-09581-y Review and recent developments on the perfectly matched layer (PML) method for the numerical modeling and simulation of elastic wave propagation in unbounded domains 0956-540X Geophysical Journal International Qiao 217 798 2019 10.1093/gji/ggz050 Theory and modelling of constant-Q viscoelastic anisotropic media using fractional derivative 0895-2477 Microwave and Optical Technology Letters Roden 27 334 2000 10.1002/1098-2760(20001205)27:5<334::AID-MOP14>3.0.CO;2-A Convolution PML (CPML): An efficient FDTD implementation of the CFS-PML for arbitrary media 0016-8033 Geophysics Robertsson 61 1921 1996 10.1190/1.1444107 A numerical free-surface condition for elastic/viscoelastic finite-difference modeling in the presence of topography 0956-540X Geophysical Journal International Schäer 198 1405 2014 10.1093/gji/ggu171 Line-source simulation for shallow-seismic data — Part 2: Full-waveform inversion — A synthetic 2-D case study Exploration Geophysics Sinclair 38 225 2007 10.1071/EG07025 2.5D modelling of elastic waves in transversely isotropic media using the spectral element method 0016-8033 Geophysics Song 60 784 1995 10.1190/1.1443817 Frequency- domain acoustic-wave modeling and inversion of crosshole data: Part I — 2.5-D modeling method 0956-540X Geophysical Journal International Takenaka 125 5 1996 10.1111/j.1365-246X.1996.tb00001.x A 2.5-D time-domain elastodynamic equation for plane-wave incidence 0956-540X Geophysical Journal International Tessmer 121 557 1995 10.1111/j.1365-246X.1995.tb05733.x 3-D seismic modelling of general material anisotropy in the presence of the free surface by a Chebyshev spectral method Spectral methods in MATLAB Trefethen 2000 10.1137/1.9780898719598 Journal of Geophysical Research: Solid Earth Wang 127 2022 Propagating seismic waves in VTI attenuating media using fractional viscoelastic wave equation 0016-8033 Geophysics Williamson 60 591 1995 10.1190/1.1443798 A critical review of acoustic wave modeling procedures in 2.5 dimensions Yang 2022 Numerical Fréchet derivatives of the displacement tensor for 2.5-D frequency-domain seismic full-waveform inversion in viscoelastic TTI media 0016-8033 Geophysics Yang 85 6 T315 2020 10.1190/geo2020-0143.1 A new generalized stiffness reduction method for 2D/2.5D frequency-domain seismic wave modeling in viscoelastic anisotropic media 0016-8033 Geophysics Yang 87 2 T85 2021 10.1190/geo2021-0166.1 Frequency-domain seismic data transformation from point source to line source for 2D viscoelastic anisotropic media 0033-4553 Pure and Applied Geophysics Yang 178 2999 2021 10.1007/s00024-021-02775-2 A generalized 2.5-D time-domain seismic wave equation to accommodate various elastic media and boundary conditions 0956-540X Geophysical Journal International Zhao 211 164 2017 10.1093/gji/ggx273 2.5-D frequency-domain viscoelastic wave modelling using finite-element method 0033-4553 Pure and Applied Geophysics Zhou 163 1399 2006 10.1007/s00024-006-0081-7 An adaptive wavenumber sampling strategy for 2.5D seismic-wave modeling in the frequency domain 0956-540X Geophysical Journal International Zhou 188 223 2012 10.1111/j.1365-246X.2011.05246.x Wavenumber sampling strategies for 2.5-D frequency-domain seismic wave modelling in general anisotropic media Computers & Geosciences Zhou 39 18 2012 10.1016/j.cageo.2011.06.005 2.5-D frequency-domain seismic wave modeling in heterogeneous, anisotropic media using a Gaussian quadrature grid technique 0016-8033 Geophysics Zhou 63 1813 1998 10.1190/1.1444475 Composite boundary-valued solution of the 2.5-D Green’s function for arbitrary acoustic media Journal of Applied Mathematics and Physics Zhou 03 358 2015 10.4236/jamp.2015.33047 Subdomain Chebyshev spectral method for 2D and 3D numerical differentiations in a curved coordinate system 10.1190/image2022-3751043.1 Zhou B. Won M. Jin C. , 2022, New numerical schemes for time-domain seismic wave modeling in viscoelastic media: Second International Meeting for Applied Geoscience & Energy, SEG, Expanded Abstracts, 2384–2388, doi: http://dx.doi.org/10.1190/image2022-3751043.1.