|
|
The application of reverse time migration to GPR data processing |
Jun-Bo GONG1, Hong-Hua WANG1,2, Min-Ling WANG1,2, Ze-Ming LUO1 |
1. School of Geosciences,Guilin University of Technology,Guilin 541004,China 2. Guangxi Key Laboratory of Hidden Metallic Ore Deposits Exploration,Guilin 541004,China |
|
|
Abstract As a shallow geophysical detection technology,Ground Penetrating Radar (GPR) has been widely used in engineering detection.In order to improve the interpretation accuracy of radar data,this paper introduces reverse time migration (RTM) to process and interpret GPR observed data.Firstly,Finite Difference Time Domain (FDTD) equations of electromagnetic wave in two-dimensional TM mode,CFL numerical stability condition and dispersion relations are derived from Maxwell equations.On such a basis,the principle of RTM and the formula of zero time imaging condition are elaborated,and the corresponding program of RTM is compiled.Two typical voids models are established.The simulated data of two typical void GPR models are processed by RTM to verify the validity of GPR RTM in improving imaging resolution.Compared with Kirchhoff migration result,the homing of reflected wave is more accurate and the convergency of diffracted wave is more complete in the RTM image.Finally,RTM is applied to processing the observed data,which is observed on the underground pipeline.The pipeline is accurately delineated,which provides a scientific basis for the follow-up construction and disposal of the road.
|
Received: 06 September 2018
Published: 15 August 2019
|
|
Corresponding Authors:
Hong-Hua WANG
|
|
|
|
|
Sketch map of reverse time migration with common offset GPR data a—simulation of exploding reflector;b—imaging of exploding reflector
|
|
Schematic map of GPR model 1
|
|
Simulated profile of model 1
|
|
Kirchhoff migration and reverse time migration profiles of model 1 a—kirchhoff migration profile;b—reverse time migration profile
|
|
Schematic map of GPR model 2
|
|
Simulated profile of model 2
|
|
Kirchhoff migration and reverse time migration profiles of model 2 a—kirchhoff migration profile;b—reverse time migration profile
|
|
Original GPR data profile of survey 1
|
|
Results of conventional data processing
|
|
Kirchhoff migration and reverse time migration profiles a—kirchhoff migration profile;b—reverse time migration profile
|
[1] |
范占锋, 李天斌, 孟陆波 . 探地雷达在公路隧道超前地质预报中的应用[J]. 物探与化探, 2010,34(1):119-122.
|
[1] |
Fan Z F, Li T B, Meng L B . Advanced geological forecast of application of GPR in road tunnel[J]. Geophysical and Geochemical Exploration, 2010,34(1):119-122.
|
[2] |
王正成, 吴晔 . 探地雷达隧道衬砌质量检测技术[J]. 物探与化探, 2013,32(6):1152-1156.
|
[2] |
Wang Z C, Wu Y . A tentative discussion on tunnel lining quality of GPR detection technology[J]. Geophysical and Geochemical Exploration, 2013,32(6):1152-1156.
|
[3] |
王亮, 王绪本, 李正文 . 探地雷达在金沙遗址考古探测中的应用研究[J]. 物探与化探, 2008,32(4):401-403.
|
[3] |
Wang L, Wang X B, Li Z W . The application of ground penetrating radar to the archeological exploration of Jinsha ruins[J]. Geophysical and Geochemical Exploration, 2008,32(4):401-403.
|
[4] |
袁明德 . 探地雷达探测地下管线的能力[J]. 物探与化探, 2002,26(2):152-155,162.
|
[4] |
Yuan M D . The capacity of the ground-penetrating radar for detecting underground pipelines[J]. Geophysical and Geochemical Exploration, 2002,26(2):152-155,162.
|
[5] |
肖建平, 胡章棚, 柳建新 , 等. 铁路路基病害检测雷达信号中的强干扰谱分析及滤波处理[J]. 物探化探计算技术, 2017,39(2):183-187.
|
[5] |
Xiao J P, Hu Z P, Liu J X , et al. Spectrum analysis and filtering on suppressing intensive noise in train-mounted GPR data[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2017,39(2):183-187.
|
[6] |
曹震峰, 杨世福, 宋世荣 , 等. 探地雷达数据处理方法及其应用[J]. 地质与勘探, 1996,32(1):34-42.
|
[6] |
Cao Z F, Yang S F, Song S R , et al. Ground penetrating radar data processing method and its application[J]. Geology and Exploration, 1996,32(1):34-42.
|
[7] |
戴前伟, 冯德山, 何继善 . Kirchhoff偏移法在探地雷达正演图像处理中的应用[J]. 地球物理学进展, 2005,20(3):849-853.
|
[7] |
Dai Q W, Feng D S, He J S . The application of Kirchhoff's migration method in the image processing of the ground penetrating radar forward simulate[J]. Progress in Geophysics, 2005,20(3):849-853.
|
[8] |
Bitri A, Grandjean G . Frequency-wavenumber modelling and migration of 2D GPR data in moderately heterogeneous dispersive media[J]. Geophysical Prospecting, 1998,46(3):287-301.
|
[9] |
张智, 刘有山, 徐涛 , 等. 弹性波逆时偏移中的稳定激发振幅成像条件[J]. 地球物理学报, 2013,56(10):3523-3533.
|
[9] |
Zhang Z, Liu Y S, Xu T , et al. A stable excitation amplitude imaging condition for reverse time migration in elastic wave equation[J]. Chinese Journal of Geophysics, 2013,56(10):3523-3533.
|
[10] |
Fisher E, Mcmechan G A, Annan A P , et al. Examples of reverse-time migration of single channel,ground-penetrating radar profiles[J]. Geophysics, 1992,57(4):577-586.
|
[11] |
鲁兴林, 钱荣毅 . 地质雷达有限差分逆时偏移方法研究[J]. 地球物理学进展, 2017,32(2):885-890.
|
[11] |
Lu X L, Qian R Y . Ground-penetrating radar finite-difference reverse time migration[J]. Progress in Geophysics, 2017,32(2):885-890.
|
[12] |
周辉, 王兆磊, 韩波 , 等. 同时实现地质雷达数据地形校正和偏移成像的方法[J]. 吉林大学学报:地球科学版, 2004,34(3):459-463.
|
[12] |
Zhou H, Wang Z L, Han B , et al. Terrain correction and migration of GPR profile fulfilled simultaneously using reverse-time migration[J]. Journal of Jiling University:Earth Science Edition, 2004,34(3):459-463.
|
[13] |
傅磊, 刘四新, 刘澜波 , 等. 机载探地雷达数值模拟及逆时偏移成像[J]. 地球物理学报, 2014,57(5):1636-1646.
|
[13] |
Fu L, Liu S X, Liu L B , et al. Airborne ground penetrating radar numerical simulation and reverse time migration[J]. Chinese Journal of Geophysics, 2014,57(5):1636-1646.
|
[14] |
底青云, 许琨, 王妙月 . 衰减雷达波有限元偏移[J] . 地球物理学报, 2000,43(2):257-263.
|
[14] |
Di Q Y, Xu K, Wang M Y . The attenuated radar wave migration with finite element method[J]. Chinese Journal of Geophysics, 2000,43(2):257-263.
|
[15] |
薛桂霞, 邓世坤, 刘秀娟 . 逆时偏移在探地雷达信号处理中的应用[J]. 煤田地质与勘探, 2004,32(1):55-57.
|
[15] |
Xue G X, Deng S K, Liu X J . An application of reverse-time migration in the ground-penetrating radar data processing[J]. Coal Geology & Exploration, 2004,32(1):55-57.
|
[16] |
韩波 . 遗传算法在地质雷达成像中的应用研究[D]. 青岛:中国海洋大学, 2005.
|
[16] |
Han B . The application of genetic algorithm in GPR prestack migration[D]. Qingdao:Ocean University of China, 2005.
|
[17] |
雷文太, 黄春琳, 粟毅 . 非均匀空间采样下的表层穿透雷达逆时偏移成像算法[J]. 工程地球物理学报, 2006,3(3):161-166.
|
[17] |
Lei W T, Huang C L, Su Y . Surface penetrating radar of non-uniform spatial sampling and reverse time migration imagine algorithm[J]. Chinese Journal of Engineering Geophysics, 2006,3(3):161-166.
|
[18] |
石明, 戴前伟, 冯德山 , 等. 高陡起伏地形 GPR 偏移的波场延拓及精度分析[J]. 中南大学学报:自然科学版, 2016,47(7):2448-2455.
|
[18] |
Shi M, Dai Q W, Feng D S , et al. Precision analysis of wave-field extrapolation in migration with steep-dip rough terrains for GPR[J]. Journal of Central South University:Science and Technology, 2016,47(7):2448-2455.
|
[19] |
葛德彪, 闫玉波 . 电磁波时域有限差分方法[M]. 西安: 西安电子科技大学出版社, 2002.
|
[19] |
Ge D B, Yan Y B. Finite difference time domain method of electromagnetic wave[M]. Xi’an: Xidian University Press, 2002.
|
[20] |
梁小强, 杨道学, 张可能 , 等. FDTD 数值模拟在 GPR 管线探测中的应用[J]. 地球物理学进展, 2017,32(4):1803-1807.
|
[20] |
Liang X Q, Yang D X, Zhang K N , et al. Application of FDTD numerical simulation of Ground Penetrating Radar in pipeline detection[J]. Progress in Geophysics, 2017,32(4):1803-1807.
|
[21] |
Berenger J P . A perfectly matched layer for the absorption of electromagnetic waves[J]. Journal of Computational Physics, 1994,114(2):185-200.
|
[22] |
薛桂霞, 刘秀娟, 邓世坤 . 雷达波偏移方法研究[J]. 地球物理学进展, 2004,19(4):903-908.
|
[22] |
Xue G X, Liu X J, Deng S K . Research of migration method of radar wave[J]. Progress in Geophysics, 2004,19(4):903-908.
|
[1] |
XIAO Guan-Hua, ZHANG Wei, CHEN Heng-Chun, ZHUO Wu, WANG Yan-Jun, REN Li-Ying. Application of shallow seismic reflection surveys in the exploration of urban underground space in Jinan[J]. Geophysical and Geochemical Exploration, 2022, 46(1): 96-103. |
[2] |
YONG Fan, LIU Zi-Long, JIANG Zheng-Zhong, LUO Shui-Yu, LIU Jian-Sheng. The key technology of shallow imaging in urban 3D seismic data processing[J]. Geophysical and Geochemical Exploration, 2021, 45(5): 1266-1274. |
|
|
|
|