2D electromagnetic simulation for ground penetrating radar with a topographic ground surface by the curvilinear collocated-grid finite-difference method combined with equivalent field method

Zhang，Heng1,4; Sun，Yao Chong4; Ren，Hengxin2,3,4; Ma，Bowen5; Zhang，Wei2,3,4; Huang，Qinghua5; Chen，Xiaofei2,3,4

2022-11-01

10.1016/j.jappgeo.2022.104812

JOURNAL OF APPLIED GEOPHYSICS   Impact Factor & Quartile

0926-9851

1879-1859

In this paper, we extend the curvilinear collocated-grid finite-difference method, which has been developed in the study of seismology, to the 2D electromagnetic simulation of ground penetrating radar (GPR). The method combines curvilinear coordinate and collocated grid, so it can better describe the geometry of the irregular interfaces than traditional finite difference method and is more efficient than finite element method which is suitable for complex geometries. Due to the differences in interfacial condition between elastic waves and electromagnetic waves, a novel scheme of equivalent electromagnetic fields is introduced to explicitly ensure discontinuous interfacial boundary conditions with collocated grid. Furthermore, this approach eliminates the need for orthogonality of the grid at the interface, which reduces the complexity of mesh generation. The proposed method is verified for a series of ground penetrating radar application models, by comparing synthetic waveforms with independent reference solutions. The perfect consistency of the comparisons indicates the accuracy of our new method in the simulation of GPR. Compared with the second-order gprMax, the proposed method is proved to be more efficient in terms of computing time and memory. Combination of collocated grid and the novel equivalent field method leads to reliable and accurate solutions of EM fields even on the ground surface.

Curvilinear collocated grid finite difference method Electromagnetic modelling Equivalent field method Ground penetrating radar Topographic ground surface

SCI

English

Corresponding

National Natural Science Foundation of China[41974081];National Natural Science Foundation of China[42022027];

Geology ; Mining & Mineral Processing

Geosciences, Multidisciplinary ; Mining & Mineral Processing

WOS:000862589800004

ELSEVIER

GEOSCIENCES

2-s2.0-85138454094

Scopus

1.School of Astronautics,Harbin Institute of Technology,Harbin,Heilongjiang,150001,China
2.Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou),Guangzhou,Guangdong,511458,China
3.Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology,Shenzhen,Guangdong,518000,China
4.Department of Earth and Space Sciences,Southern University of Science and Technology,Shenzhen,Guangdong,518000,China
5.Department of Geophysics,School of Earth and Space Sciences,Peking University,Beijing,100871,China

Zhang，Heng,Sun，Yao Chong,Ren，Hengxin,et al. 2D electromagnetic simulation for ground penetrating radar with a topographic ground surface by the curvilinear collocated-grid finite-difference method combined with equivalent field method[J]. JOURNAL OF APPLIED GEOPHYSICS,2022,206.

Zhang，Heng.,Sun，Yao Chong.,Ren，Hengxin.,Ma，Bowen.,Zhang，Wei.,...&Chen，Xiaofei.(2022).2D electromagnetic simulation for ground penetrating radar with a topographic ground surface by the curvilinear collocated-grid finite-difference method combined with equivalent field method.JOURNAL OF APPLIED GEOPHYSICS,206.

Zhang，Heng,et al."2D electromagnetic simulation for ground penetrating radar with a topographic ground surface by the curvilinear collocated-grid finite-difference method combined with equivalent field method".JOURNAL OF APPLIED GEOPHYSICS 206(2022).