王敏玲^1,2,,
王洪华^1,2,,,
张智^1,2,
徐志锋^1,2
1. 桂林理工大学地球科学学院, 广西桂林 541004
2. 广西隐伏金属矿产勘查重点实验室, 广西桂林 541004

基金项目: 国家自然科学基金项目（41604039，41604102，41574078，41274070），广西自然科学基金项目（2016GXNSFBA380215，2016GXNSFBA380082，2015GXNSFAA139238，2014GXNSFAA118305），广西中青年教师基础能力提升项目（KY2016YB199）资助


详细信息

作者简介: 王敏玲, 女, 1988年生, 博士, 讲师, 主要从事地球物理成像理论及应用研究.E-mail:wangml@glut.edu.cn

通讯作者: 王洪华, 男, 1986年生, 博士, 讲师, 主要从事探地雷达的正演模拟与反演成像研究.E-mail:wanghonghua5@163.com

中图分类号: P631

收稿日期:2017-04-28
修回日期:2018-07-04
上线日期:2018-08-05

Reverse time migration in Ground Penetrating Radar based on the excitation amplitude imaging condition

WANG MinLing^1,2,,
WANG HongHua^1,2,,,
ZHANG Zhi^1,2,
XU ZhiFeng^1,2
1. College of Earth Sciences, Guilin University of Technology, Guangxi Guilin 541004, China
2. Guangxi Key Laboratory of Hidden Metallic Ore Deposits Exploration, Guangxi Guilin 541004, China


More Information

Corresponding author: WANG HongHua,E-mail:wanghonghua5@163.com

MSC: P631

--> Received Date: 28 April 2017
Revised Date: 04 July 2018
Available Online: 05 August 2018


摘要
摘要:针对基于互相关成像条件的探地雷达（GPR）逆时偏移计算效率低、存储量大及易产生低频假象的不足，本文将激发振幅成像条件应用于GPR逆时偏移成像中.通过在源点电磁波场正向传播过程计算每个网格点的能量密度，并保存最大能量密度的时刻和相应的电磁波场值；在接收点电磁波场逆向传播过程提取每个网格点最大能量密度时刻及对应的电磁波场值，并利用保存的最大能量源点电磁波场及走时做归一化，从而获得了依赖反射系数成像剖面，避免了源点正向传播电磁波场的存储和重建.此外，为了提高电磁波场的模拟精度，采用了基于三角形剖分的时间域有限单元法（FETD）计算电磁波正向和逆向传播过程.最后通过模型试算表明：激发振幅成像条件相比于归一化互相关成像条件，成像结果低频噪声更弱，空间分辨率更高，计算效率提高了近2倍.
关键词: 探地雷达/
逆时偏移/
激发振幅成像条件/
能量密度/
归一化互相关成像条件

Abstract:In order to overcome the shortcoming of lower calculation efficiency, larger storage and low frequency illusion, which can be easily produced in Ground Penetrating Radar (GPR) reverse time migration (RTM) based on normalized correlation imaging condition, a stable excitation amplitude imaging condition for GPR reverse time migration is proposed. In the propagation of the source electromagnetic wave filed extrapolation along a positive time axis, the energy density for total grids are computed at each time step, and the travel time as well as electromagnetic wavefiled values corresponding to the maximum energy density are stored to extract and normalize the receiver electromagnetic wave field, the imaging profile which rely on reflectance can be obtained. Compared to the normalized cross-correlation (NC) imaging condition, this strategy will save a great amount of memory resource and significantly improve the computational efficiency for RTM. In order to improve the simulation accuracy of electromagnetic wavefield, the finite element time domain (FETD) method based on triangulation is used to calculate the forward and reverse propagation process of electromagnetic wave. At last the feasibility and accuracy of the proposed GPR RTM algorithm are validated by numerical tests of inclined model and fluctuation interface model by compared with the results of RTM based on the NC imaging condition. The compared results demonstrated that the RTM based on the excitation imaging condition has a higher spatial resolution and can produce less low-frequency artifacts compared with the RTM results based on NC imaging condition. Moreover, the calculation efficiency of the excitation imaging condition has improved nearly 2 times because computer storage capacity due to that it needs not to store source electromagnetic field.
Key words:Ground Penetrating Radar (GPR)/
Reverse Time Migration (RTM)/
The excitation amplitude imaging condition/
Energy density/
The normalized correlation imaging condition



PDF全文下载地址:

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