Original Article

Decomposition of Displacement Field into the Irrotational and Solenoidal Component Using Fast Fourier Transform

Abstract

Purpose: A new code based on Helmholtz decomposition is presented to separate longitudinal (pressure) and transverse (shear) components of a mixed wave field. This algorithm will help isolate shear or pressure components of an elastic wave to further concentrate on each specific wave and its physical characteristics, particularly in medical imaging instrument development and image processing techniques.

Materials and Methods: Using the combination of Fourier transform and Helmholtz decomposition, first, the mathematical basis of the work is prepared. After reaching a usable formula, this basis is embedded in the Code written in MATLAB program. Then, various test data containing shear and pressure waves were created and fed to the Code to evaluate its ability to decompose the displacements into the shear and pressure waves.

Results: This new algorithm successfully isolated the transverse and longitudinal wavefront of the mixed wavefield. The Code demonstrated 100% accuracy for separating the shear wave and more than 99% for the pressure wave. Moreover, the background noise was kept under 0.03% in every step.

Conclusion: The results show that using Helmholtz decomposition in Fourier space on 3D data can help decompose a displacement field into its irrotational and solenoidal components with high accuracy. A weak dependency on wave thickness and contrast was observed, but the algorithm's accuracy never fell below 99%.

1- Wolfgang Sprössig, "On Helmholtz decompositions and their generalizations—an overview." Mathematical Methods in the Applied Sciences, Vol. 33 (No. 4), pp. 374-83, (2010).
2- Jia Yan and Paul Sava, "Isotropic angle-domain elastic reverse-time migration." Geophysics, Vol. 73 (No. 6), pp. S229-S39, (2008).
3- ROBERT Sun, "Separating P-and S-waves in a prestack 2-dimensional elastic seismogram." in 61st EAGE Conference and Exhibition, (1999): European Association of Geoscientists & Engineers, pp. cp-132-00147.
4- Robert Sun, Jinder Chow, and Kuang-Jung Chen, "Phase correction in separating P-and S-waves in elastic data." Geophysics, Vol. 66 (No. 5), pp. 1515-18, (2001).
5- Robert Sun, George A McMechan, Hsu-Hong Hsiao, and Jinder Chow, "Separating P-and S-waves in prestack 3D elastic seismograms using divergence and curl." Geophysics, Vol. 69 (No. 1), pp. 286-97, (2004).
6- Robert Sun, George A McMechan, Chen-Shao Lee, Jinder Chow, and Chen-Hong Chen, "Prestack scalar reverse-time depth migration of 3D elastic seismic data." Geophysics, Vol. 71 (No. 5), pp. S199-S207, (2006).
7- Xiang Xiao and W Scott Leaney, "Local vertical seismic profiling (VSP) elastic reverse-time migration and migration resolution: Salt-flank imaging with transmitted P-to-S waves." Geophysics, Vol. 75 (No. 2), pp. S35-S49, (2010).
8- Qunshan Zhang and George A McMechan, "2D and 3D elastic wavefield vector decomposition in the wavenumber domain for VTI media." Geophysics, Vol. 75 (No. 3), pp. D13-D26, (2010).
9- Jiubing Cheng and Sergey Fomel, "Fast algorithms for elastic-wave-mode separation and vector decomposition using low-rank approximation for anisotropic media." Geophysics, Vol. 79 (No. 4), pp. C97-C110, (2014).
10- Hejun Zhu, "Elastic wavefield separation based on the Helmholtz decompositionElastic wavefield separation." Geophysics, Vol. 82 (No. 2), pp. S173-S83, (2017).
11- Jidong Yang, Hejun Zhu, Wenlong Wang, Yang Zhao, and Houzhu Zhang, "Isotropic elastic reverse time migration using the phase-and amplitude-corrected vector P-and S-wavefields." Geophysics, Vol. 83 (No. 6), pp. S489-S503, (2018).
12- Yang Zhao, Houzhu Zhang, Jidong Yang, and Tong Fei, "Reducing artifacts of elastic reverse time migration by the deprimary technique." Geophysics, Vol. 83 (No. 6), pp. S569-S77, (2018).
13- Bingluo Gu, Zhiming Ren, Qingqing Li, Jianguang Han, and Zhenchun Li, "An application of vector wavefield decomposition to 3D elastic reverse time migration and field data test." Computers & Geosciences, Vol. 131pp. 112-31, (2019).
14- Harsh Bhatia, Gregory Norgard, Valerio Pascucci, and Peer-Timo Bremer, "The Helmholtz-Hodge decomposition—a survey." IEEE Transactions on visualization and computer graphics, Vol. 19 (No. 8), pp. 1386-404, (2012).
15- Cherif Amrouche, Christine Bernardi, Monique Dauge, and Vivette Girault, "Vector potentials in three‐dimensional non‐smooth domains." Mathematical Methods in the Applied Sciences, Vol. 21 (No. 9), pp. 823-64, (1998).
Files
IssueVol 10 No 4 (2023) QRcode
SectionOriginal Article(s)
DOI https://doi.org/10.18502/fbt.v10i4.13729
Keywords
Helmholtz Decomposition Fourier Transform Vector Field Longitudinal Waves Transverse Waves

Rights and permissions
Creative Commons License This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.
How to Cite
1.
Bahrami Gorji R, Mohammadi M, Makkiabadi B. Decomposition of Displacement Field into the Irrotational and Solenoidal Component Using Fast Fourier Transform. Frontiers Biomed Technol. 2023;10(4):465-479.