Abstract
Fractional Laplacian viscoacoustic/viscoelastic wave equations offer separate controls over amplitude loss and phase dispersion, and have been used in Q-compensated reverse-time migration and full waveform inversion. Previously, the spatially varying-order fractional Laplacians have been solved with the global Fourier pseudo-spectral method by representing the spatially varying order with an average value, which introduces numerical errors into simulations. To reduce the errors, we propose a local pseudo-spectral method, which uses a large number of block-variable values instead of just one to represent the spatially varying order. A numerical implementation scheme for parallel computing, domain decomposition, has been adopted to take advantage of the local pseudo-spectral method, which improves both numerical accuracy and computing efficiency. A tapering internal boundary condition is used to reduce the Fourier artifacts caused by wavefield truncation at subdomain boundaries. An overlap-add communication scheme bewteen subdomains is applied for reducing the additional cost associated with boundary padding and for interpolating the wavefields from different subdomains within the overlapping boundaries. Numerical examples verify the effectiveness of the domain decomposition strategy in improving the accuracy of solving fractional Laplacians in the viscoelastic wave equations.
| Original language | English (US) |
|---|---|
| Pages | 3943-3947 |
| Number of pages | 5 |
| DOIs | |
| State | Published - Jan 1 2019 |
| Event | 88th Society of Exploration Geophysicists International Exposition and Annual Meeting, SEG 2018 - Anaheim, United States Duration: Oct 14 2018 → Oct 19 2018 |
Other
| Other | 88th Society of Exploration Geophysicists International Exposition and Annual Meeting, SEG 2018 |
|---|---|
| Country | United States |
| City | Anaheim |
| Period | 10/14/18 → 10/19/18 |
Fingerprint
All Science Journal Classification (ASJC) codes
- Geophysics
Cite this
}
Solving fractional Laplacian viscoelastic wave equations using domain decomposition. / Xue, Zhiguang; Baek, Hyoungsu; Zhang, Houzhu; Zhao, Yang; Zhu, Tieyuan; Fomel, Sergey.
2019. 3943-3947 Paper presented at 88th Society of Exploration Geophysicists International Exposition and Annual Meeting, SEG 2018, Anaheim, United States.Research output: Contribution to conference › Paper
TY - CONF
T1 - Solving fractional Laplacian viscoelastic wave equations using domain decomposition
AU - Xue, Zhiguang
AU - Baek, Hyoungsu
AU - Zhang, Houzhu
AU - Zhao, Yang
AU - Zhu, Tieyuan
AU - Fomel, Sergey
PY - 2019/1/1
Y1 - 2019/1/1
N2 - Fractional Laplacian viscoacoustic/viscoelastic wave equations offer separate controls over amplitude loss and phase dispersion, and have been used in Q-compensated reverse-time migration and full waveform inversion. Previously, the spatially varying-order fractional Laplacians have been solved with the global Fourier pseudo-spectral method by representing the spatially varying order with an average value, which introduces numerical errors into simulations. To reduce the errors, we propose a local pseudo-spectral method, which uses a large number of block-variable values instead of just one to represent the spatially varying order. A numerical implementation scheme for parallel computing, domain decomposition, has been adopted to take advantage of the local pseudo-spectral method, which improves both numerical accuracy and computing efficiency. A tapering internal boundary condition is used to reduce the Fourier artifacts caused by wavefield truncation at subdomain boundaries. An overlap-add communication scheme bewteen subdomains is applied for reducing the additional cost associated with boundary padding and for interpolating the wavefields from different subdomains within the overlapping boundaries. Numerical examples verify the effectiveness of the domain decomposition strategy in improving the accuracy of solving fractional Laplacians in the viscoelastic wave equations.
AB - Fractional Laplacian viscoacoustic/viscoelastic wave equations offer separate controls over amplitude loss and phase dispersion, and have been used in Q-compensated reverse-time migration and full waveform inversion. Previously, the spatially varying-order fractional Laplacians have been solved with the global Fourier pseudo-spectral method by representing the spatially varying order with an average value, which introduces numerical errors into simulations. To reduce the errors, we propose a local pseudo-spectral method, which uses a large number of block-variable values instead of just one to represent the spatially varying order. A numerical implementation scheme for parallel computing, domain decomposition, has been adopted to take advantage of the local pseudo-spectral method, which improves both numerical accuracy and computing efficiency. A tapering internal boundary condition is used to reduce the Fourier artifacts caused by wavefield truncation at subdomain boundaries. An overlap-add communication scheme bewteen subdomains is applied for reducing the additional cost associated with boundary padding and for interpolating the wavefields from different subdomains within the overlapping boundaries. Numerical examples verify the effectiveness of the domain decomposition strategy in improving the accuracy of solving fractional Laplacians in the viscoelastic wave equations.
UR - http://www.scopus.com/inward/record.url?scp=85059364425&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85059364425&partnerID=8YFLogxK
U2 - 10.1190/segam2018-2998547.1
DO - 10.1190/segam2018-2998547.1
M3 - Paper
AN - SCOPUS:85059364425
SP - 3943
EP - 3947
ER -