Second-generation wavelet fi nite element based on the lifting scheme for GPR simulation*
Feng De-Shan 1,2, Zhang Hua 1,2, and Wang Xun 1,2
1. School of Geosciences and Info-Physics, Central South University, Changsha 410083, China.
2. Key Laboratory of Metallogenic Prediction of Nonferrous Metals, Ministry of Education, Changsha 410083, China.
Abstract:
Ground-penetrating radar (GPR) is a highly efficient, fast and non-destructive exploration method for shallow surfaces. High-precision numerical simulation method is employed to improve the interpretation precision of detection. Second-generation wavelet finite element is introduced into the forward modeling of the GPR. As the finite element basis function, the second-generation wavelet scaling function constructed by the scheme is characterized as having multiple scales and resolutions. The function can change the analytical scale arbitrarily according to actual needs. We can adopt a small analysis scale at a large gradient to improve the precision of analysis while adopting a large analytical scale at a small gradient to improve the efficiency of analysis. This approach is beneficial to capture the local mutation characteristics of the solution and improve the resolution without changing mesh subdivision to realize the efficient solution of the forward GPR problem. The algorithm is applied to the numerical simulation of line current radiation source and tunnel non-dense lining model with analytical solutions. Result show that the solution results of the secondgeneration wavelet finite element are in agreement with the analytical solutions and the conventional finite element solutions, thereby verifying the accuracy of the second-generation wavelet finite element algorithm. Furthermore, the second-generation wavelet finite element algorithm can change the analysis scale arbitrarily according to the actual problem without subdividing grids again. The adaptive algorithm is superior to traditional scheme in grid refi nement and basis function order increase, which makes this algorithm suitable for solving complex GPR forward-modeling problems with large gradient and singularity.