GPU acceleration of time gating based reverse time migration using the pseudospectral time-domain algorithm

Xie, Jiangang; Guo, Zichao; Liu, Hai; Liu, Qing

doi:10.1016/j.cageo.2018.05.006

Cited by 13 publications

(5 citation statements)

References 23 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The transmission signal was a Gaussian envelope pulse containing two sinusoidal cycles with a 1 MHz centre frequency. The full-wave simulations were implemented with the pseudospectrum time domain (PSTD) method [37]. The PSTD method assumes that all meshes have square shapes with the same sizes in the computational region.…”

Section: B Simuationmentioning

confidence: 99%

See 1 more Smart Citation

Ray Theory-Based Transcranial Phase Correction for Intracranial Imaging: A Phantom Study

Jiang

et al. 2019

IEEE Access

View full text Add to dashboard Cite

Ultrasonic imaging provides a non-invasive way to diagnose brain disease. However, due to imaging degradation effects from the phase-aberration and reverberation, it is still challenging to achieve an accurate transcranial imaging. The objective of this work is to improve the quality of transcranial imaging. To this end, a ray theory based transcranial phase correction method was proposed to correct the phase abberation induced by cranial bones. With the pre-knowledge of the shape and longitudinal velocity of the cranial bones, the corrected phases are derived by solving eikonal equation (ray-theory). The Ideal Synthetic Aperture Focusing Technique (I-SAFT) is applied for signal acquisitions in simulations and in-vitro phantom experiment with one-element transmitting and all-element receiving method. Dynamic focusing is achieved at each imaging position with I-SAFT and the transcranial imaging distortion is modified with the proposed phase correction method. Simulations and experiments show that the imaging distortion of target circular phantoms was corrected, and the imaging quality is improved after the phase correction. With the proposed method, the average error of the central position of target phantoms decreases from 1.98 mm to 0.21 mm, the eccentricity of fitted ellipse averagely decreases from 0.63 to 0.19, and the average maximum luminance contrast of phantoms improves from 37.36dB to 42.41dB. It is illustrated that the proposed ray-theory based phase correction method might be useful for intracranial imaging.

show abstract

Section: B Simuationmentioning

confidence: 99%

“…The mesh size was assigned as dx = 0.11 mm, approximately one-thirteenth the wavelength λ = 1.48 mm in water and one-third the element spacing of the phase array. The time step was arranged to satisfy the Courant-Friedrich-Levy condition as 0.1 [37]. Two simulations were implemented altogether.…”

Section: B Simuationmentioning

confidence: 99%

Ray Theory-Based Transcranial Phase Correction for Intracranial Imaging: A Phantom Study

Jiang

et al. 2019

IEEE Access

View full text Add to dashboard Cite

show abstract

“…The frugal requirements of the method made it popular in early applications to seismic modeling in the 1980s, given the limited computer memories available at the time (e.g., Kosloff and Baysal, 1982;Fornberg, 1987Fornberg, , 1988Etgen and Dellinger, 1989;Daudt et al, 1989). Recent applications have focused on complex earth models and parallel implementations (see, for instance, Klin et al, 2010;Peng and Cheng, 2016;Xie et al, 2016Xie et al, , 2018. PSTD applications in geophysics are typically defined on unbounded domains or half spaces, thus requiring effective numerical methods to avoid reflections from the computational boundaries of the domain under study.…”

Section: Introductionmentioning

confidence: 99%

Calibration of absorbing boundary layers for geoacoustic wave modeling in pseudo-spectral time-domain methods

Spa,

Rojas,

de la Puente

2023

Geosci. Model Dev.

View full text Add to dashboard Cite

Abstract. This paper develops a calibration methodology of the artificial absorbing techniques typically used by Fourier pseudo-spectral time-domain (PSTD) methods for geoacoustic wave simulations. Specifically, we consider the damped wave equation (DWE) that results from adding a dissipation term to the original wave equation, the sponge boundary layer (SBL) that applies an exponentially decaying factor directly to the wavefields, and finally, a classical split formulation of the perfectly matched layer (PML). These three techniques belong to the same family of absorbing boundary layers (ABLs), where outgoing waves and edge reflections are progressively damped across a grid zone of NABL consecutive layers. The ABLs used are compatible with a pure Fourier formulation of PSTD. We first characterize the three ABLs with respect to multiple sets of NABL and their respective absorption parameters for homogeneous media. Next, we validate our findings in heterogeneous media of increasing complexity, starting with a layered medium and finishing with the SEG/EAGE 3D salt model. Finally, we algorithmically compare the three PSTD–ABL methods in terms of memory demands and computational cost. An interesting result is that PML, despite outperforming the absorption of the other two ABLs for a given NABL value, requires approximately double the computational time. The parameter configurations reported in this paper can be readily used for PSTD simulations in other test cases, and the ABL calibration methodology can be applied to other wave propagation schemes.

show abstract

“…The implementation of high-efficiency FDTD methods with GPUs is currently a subject of exploration (e.g., Dziekonski, Lamecki and Mrozowski 2011;Zygiridis, Kantartzis and Tsiboukis 2014;Xie et al 2018). Like the FDTD method, the symplectic Euler method is also suitable for parallel computing and thus benefits from advances in parallel and multi-core computing.…”

Section: Introductionmentioning

confidence: 99%

Modelling ground‐penetrating radar wave propagation using graphics processor unit parallel implementation of the symplectic Euler method

Fang

Lei

Zhang

et al. 2019

Near Surface Geophysics

View full text Add to dashboard Cite

Inversion of ground‐penetrating radar signals requires accurate and efficient forward modelling. The symplectic Euler method promises good results when simulating ground‐penetrating radar wave propagation in substructures, but its computational efficiency is limited by the same Courant–Friedrichs–Lewy stability condition as the finite‐difference time‐domain method. A two‐dimensional graphics processor unit–accelerated parallel symplectic Euler algorithm is used to simulate ground‐penetrating radar wave propagation. We compared the reflection waveforms as well as the simulation time of the complex underground structure models simulated by the parallel symplectic Euler method with traditional finite‐difference time‐domain method. Results show that the parallel symplectic Euler algorithm achieves the same level of accuracy as the standard finite‐difference time‐domain method. Moreover, it significantly improves the computational efficiency, as the calculation speed is improved by more than 21 times. We verify the performance of the proposed algorithm through a map of the single‐track radar data for a three‐layered pavement model and a simulation wiggle map for a structural damage pavement model. This provides a theoretical basis for accurately interpreting ground‐penetrating radar detection data and efficient forward modelling for the next step of inversion imaging.

show abstract

GPU acceleration of time gating based reverse time migration using the pseudospectral time-domain algorithm

Cited by 13 publications

References 23 publications

Ray Theory-Based Transcranial Phase Correction for Intracranial Imaging: A Phantom Study

Ray Theory-Based Transcranial Phase Correction for Intracranial Imaging: A Phantom Study

Calibration of absorbing boundary layers for geoacoustic wave modeling in pseudo-spectral time-domain methods

Modelling ground‐penetrating radar wave propagation using graphics processor unit parallel implementation of the symplectic Euler method

Contact Info

Product

Resources

About