Solving fractional Laplacian visco-acoustic wave equations on complex-geometry domains using Grunwald-formula based radial basis collocation method
Solving fractional Laplacian visco-acoustic wave equations on complex-geometry domains using Grunwald-formula based radial basis collocation method
复制标题
使用基于 Grunwald 公式的径向基配置方法求解复杂几何域上的分数拉普拉斯粘声波方程
DOI：
10.1016/j.camwa.2019.10.013
复制
发表时间：
2020
影响因子：
2.9
通讯作者：
Guofei Pang
中科院分区：
文献类型：
--
作者：
Yiran Xu;Jingye Li;Xiaohong Chen;Guofei Pang
We propose a radial basis function collocation method (RBF method) to solve fractional Laplacian visco-acoustic wave equation for the Earth media having heterogeneous velocity model and complex geometry. Unlike the fractional Laplacian wave equation proposed in Zhu and Harris (2014), the wave equation we consider has a different definition for the fractional Laplacian. Specifically, spectral and Riesz fractional Laplacians are considered in Zhu and Harris (2014) and the present paper, respectively. Accordingly, the Fourier pseudospectral method (FPS method) and the RBF method are employed to solve the spectral and the Riesz fractional Laplacian wave equations. The two wave equations are observed to produce obviously different wavefields. We demonstrate the validity and flexibility of the proposed RBF method by considering five benchmarks of seismic forward modeling: (1) two-dimensional Earth media with four types of velocity models (homogeneous, two-layer, homogeneous but complex-geometry, and heterogeneous models) and (2) a spherical medium with homogeneous velocity model. We make a three-way comparison among numerical solutions to the Riesz fractional Laplacian, the spectral fractional Laplacian, and the integer-order visco-acoustic wave equations, and observe that when wave attenuation is weak the Riesz wave equation yields more similar wavefield to that of the integer-order wave equation than the spectral wave equation does. Furthermore, uniform and quasi-uniform layouts for collocation points of the RBF method are considered, and the latter layout turns out to be economical since it can preserve the solution accuracy with the minimum number of collocation points. The RBF method is truly mesh-free and dimension-free and can easily handle high-dimensional, irregular domains. Additionally, the method is easier to implement than element-based methods, such as finite element and spectral element methods, for discretizing the Riesz fractional Laplacian.
