冯德山^1,2,,
王珣^1,2,,
1. 中南大学地球科学与信息物理学院, 长沙 410083
2. 有色金属成矿预测教育部重点实验室, 长沙 410083

基金项目: 国家自然科学基金项目（41574116），中南大学创新驱动项目（2015CX008），教育部新世纪优秀人才支持计划（NCET-12-0551），中南大学教师研究基金（2014JSJJ001），中南大学升华育英人才计划，湖湘青年创新创业平台培养对象项目共同资助


详细信息

作者简介: 冯德山, 男, 1978年生, 汉族, 湖南祁阳人, 博士, 教授, 从事地球物理数据处理与正反演研究.E-mail:fengdeshan@126.com

通讯作者: 王珣, 男, 1990年生, 博士研究生, 主要从事电磁法数值模拟研究.E-mail:wangxun0727@csu.edu.cn

中图分类号: P631

收稿日期:2017-07-30
修回日期:2018-03-12
上线日期:2018-09-05

Adaptive multi-scale second-generation wavelet collocation method numerical simulation of Ground Penetrating Radar

FENG DeShan^1,2,,
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


More Information

Corresponding author: WANG Xun,E-mail:wangxun0727@csu.edu.cn

MSC: P631

--> Received Date: 30 July 2017
Revised Date: 12 March 2018
Available Online: 05 September 2018


摘要
摘要:基于第二代小波变换的提升方案构造了插值小波，将雷达波场函数进行了二维小波变换，得到所有尺度上与计算网格相联系的小波系数和尺度系数.对所有尺度上的小波系数进行分析，根据解的局部性与小波系数阈值的控制，实现网格压缩和配点的自适应调节.保留大于给定阈值的小波系数及对应网格点，令小于给定阈值的小波系数为零，并舍弃其对应网格点.达到光滑区域采用较少的计算网格点，在奇异性较大的区域采用较多的计算网格点的目的.通过对自适应网格进行邻域校正、重构检查等附加修正，推导了场值更新的显式时间迭代方案.最后，以均匀、阶梯与复杂三个典型GPR模型为例，与常规数值计算结果对比表明：自适应小波配点法（AWCM）利用第二代小波的多尺度分解和快速变换的特点，可以使计算网格随着时间步适应解的移动和变化，允许计算资源更有效地使用，具有高压缩率，达到跟踪奇异性的目的，特别适合于探地雷达正演中波传问题的模拟.
关键词: 探地雷达/
自适应小波配点法/
多尺度/
第二代小波/
插值尺度函数

Abstract:In this work, an interpolating wavelet is constructed by the lifting scheme of the second generation wavelet transformation, by employing the two-dimensional wavelet transform into the wave field function of radar, both the wavelet coefficients and scale coefficients that corresponding to the computational grid of different scales can be all calculated. And the wavelet coefficients of different scales are sequentially analyzed according to the locality of solution and the threshold of wavelet coefficient which facilitates the implementation of mesh compression and self-adaptation adjustment to collocation point. Since the wavelet coefficients that larger than a given threshold are reserved, as well as its corresponding grid points, while others that less than the given threshold are set to zero and discarded. Subsequently, the strategy of using fewer computational grid points in smoother region and more computational grid points in singularity region can be implemented and have effect. In addition, an explicit scheme of time iteration for updating the field values is deduced by using an additional corrections techniques to the adaptive mesh, for instance, the adding adjacent zone and the reconstruction check. Finally, the validation of the proposed method is carried out by using three typical simulation tests, which involving a homogeneous, a stair and a complicated model. Compared with the conventional method, the results demonstrate that, the advantage in multi-scale decomposition and fast transform of second generation wavelet, the proposed method can make the computing grid adapted to the movement or change of the solution with time step, and make the computing resources be utilized more efficiently, which characterizes proposed methods high compression rate and tracks the singularity region more accurately. Thus, the proposed adaptive wavelet collocation method (AWCM) is considered as the suitable solution for the simulation of wave-propagation problem in GPR.
Key words:Ground Penetrating Radar/
Adaptive wavelets collocation method/
Multi-scale/
Second-generation wavelet/
Interpolation scale function



PDF全文下载地址:

http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_14685