| 衰减夹层GPR模拟的时频域全局反射误差 |
| |
| 引用本文： | 张彬, 戴前伟, 尹小波, 冯德山. 2017. 衰减夹层GPR模拟的时频域全局反射误差. 地球物理学报, 60(3): 1168-1178, doi: 10.6038/cjg20170327 |
| |
| 作者姓名： | 张彬 戴前伟 尹小波 冯德山 |
| |
| 作者单位： | 1. 中南大学地球科学与信息物理学院, 长沙 410083; 2. 有色金属成矿预测与地质环境监测教育部重点实验室, 长沙 410083; 3. 中南大学资源勘查与环境地质研究院, 长沙 410083 |
| |
| 基金项目： | 国家自然科学基金项目（41374118，41574116）；中国铁路总公司科技研究开发计划项目（2014G005-B）；湖南省交通科技计划项目（201423）；中南大学创新驱动项目（2015CX008）；湖南省住房和城乡建设厅科技计划项目（BZ201408，BZ201411）；中南大学博士后基金联合资助. |
| |
| 摘 要： | 全局反射误差分析是深入研究探地雷达（GPR）吸收边界条件吸收效率的有力工具.基于常规完全匹配层（PML）的标准交错网格有限差分算法必须满足严格的CFL条件限制，即在单位时间步长内，不容许电磁波传播的距离超过单元网格尺寸.为了提高主区域所有网格节点的计算效率，并有效地吸收传播后期出现的低频隐失波，提出基于非分裂递归卷积完全匹配层（UCPML）的旋转交错网格（RSG）GPR正演算法.该算法采用不同的网格交错策略，并在边界条件中引入了吸收低频隐失波的自由可变因子，使得该算法允许选取较大的时间步长，提高了计算效率，并且实现了对低频隐失波的高效吸收.本文首先给出了RSG差分格式，推导了基于UCPML的RSG差分更新方程，探讨了数值色散的稳定性条件，然后以绕射现象严重的衰减夹层数值模拟为例，分别从波场快照、单道波记录、时间域反射误差（TDRE）、频率域反射误差（FDRE）四个方面分析了UCPML与常规PML的全局反射误差，说明了基于UCPML和RSG的GPR正演算法能更有效地吸收低频隐失波.
|
| 关 键 词： | 探地雷达 数值模拟 衰减夹层 时间域反射误差 频率域反射误差 全局反射误差 |
| 收稿时间： | 2015-03-26 |
| 修稿时间： | 2016-12-06 |
Global reflection errors in the time-frequency domain for GPR simulation in an attenuation interlayer |
| |
| ZHANG Bin, DAI Qian-Wei, YIN Xiao-Bo, FENG De-Shan. 2017. Global reflection errors in the time-frequency domain for GPR simulation in an attenuation interlayer. Chinese Journal of Geophysics (in Chinese), 60(3): 1168-1178, doi: 10.6038/cjg20170327 |
| |
| Authors: | ZHANG Bin DAI Qian-Wei YIN Xiao-Bo FENG De-Shan |
| |
| Affiliation: | 1. School of Geosciences and Info-Physics, Central South University, Changsha 410083, China; 2. Key Laboratory of Metallogenic Prediction of Nonferrous Metals and Geological Environment Monitoring, Ministry of Education, Changsha 410083, China; 3. Institute of Resource Exploration and Environmental Geology, Central South University, Changsha 410083, China |
| |
| Abstract: | Global reflection error analysis is an effective way to compare absorbing availability among different absorbing boundary conditions (ABC) in Ground Penetrating Radar (GPR), in particular the reflection errors in time-frequency domain. When solving the partial differential equations of Maxwell numerically with the method of finite differences (FD), the standard staggered grid based on the regular perfectly matched layer (PML) must meet the strict limit of Courant-Friedrichs-Lewy (CFL) conditions. In order to improve the computation efficiency and absorb efficiency in the later stage, this work studies how the un-split convolution perfectly matched layer (UCPML) and the rotated staggered grid (RSG) are combined to accomplish the efficient absorption of evanescent modes, as well as to improve the computational efficiency without numerical dispersion.With different strategies of the space staggered grid, the RSG scheme using the rotated staggered finite-difference operators can be used to calculate the difference of field components and dielectric parameters with the linear combination value across the diagonal coordinate axes, with a looser limit of CFL conditions. In this method, a larger time step can be set, which eventually leads to a higher computational efficiency. And then the UCPML is implemented in the boundary region, and in addition to the decay parameter, the other two real variable factors are introduced, which is exclusively for evanescent modes at grazing incidence. Then a numerical modeling example with an attenuation interlayer model embedded is presented, which contains a serious diffractions. Finally, a detailed comparison of UCPML and regular PML is made in the aspects of full-time snapshots, single scan signal, and reflection errors in the time-frequency domain..Both the full-time snapshots and single scan signal show that the amplitude and energy of reflection waves with regular PML are larger than the UCPML, which indicates that UCPML is more efficient in absorbing evanescent modes. The Time Domain Reflection Errors (TDRE) results demonstrate that the reflection errors of UCPML are much less than the regular PML at most of the electromagnetic wave full-time propagation, and the FDRE results show that UCPML has an obvious advantage over the regular PML in absorbing low frequency components of evanescent modes..The combination of RSG and UCPML proposed in this work can easily to be programmed and performed, which can improve the computational efficiency and absorb efficiency on long-term evanescent modes at grazing incidence. |
| |
| Keywords: | Ground Penetrating Radar (GPR) Numerical modeling Attenuation interlayer Time Domain Reflection Errors (TDRE) Frequency Domain Reflection Errors (FDRE) Global reflection errors |
| 本文献已被 CNKI 等数据库收录！ |
| 点击此处可从《地球物理学报》浏览原始摘要信息 |
|
点击此处可从《地球物理学报》下载免费的PDF全文 |
|
