论文部分内容阅读
近年来,随着数字处理技术和电子技术的飞速发展,探地雷达(GPR)的实际应用范围迅速扩大,现已覆盖考古、矿产资源勘探、水文地质调查、岩土勘查、无损检测、工程建筑物结构调查、军事等众多领域,解决了很多工程实际问题,成为浅层勘探的有力工具.而探地雷达的理论研究与实际的应用相比,具有明显的滞后性.但是解释人员要达到精确地对探地雷达实际资料的进行解析,必须事先了解地质体的雷达反射剖面的特征,所以作为反演与解释基础的复杂地电模型的探地雷达正演模拟技术,就成了探地雷达理论研究的主要内容之一.本文以麦克斯韦两个旋度方程为基本出发点,运用K.S.Yee的空间网格模型理论和时域有限差分法的基本原理,推导出二维空间的探地雷达正演方程组,并详细地分析了差分格式中半空间步长与半时间步长的实现方法,及其雷达波电场与磁场分量在计算机上相互关系的C程序实现.然后讨论了数值频散关系及其产生原因,通过同时考虑时域有限差分法及Yee氏网格的特点,推导出了符合探地雷达实际传播规律的理想频散关系,作者自制了探地雷达正演程序,并分别计算了Mur超吸收边界条件及无边界条件下的雷达地电模型,通过对比可知,超吸收边界条件可利用,大大地减少截断边界处的干扰波,达到用有限区域达到在无限空间传播的效果.最后作者利用自制程序,对“V”字形和同一斜面上的五个圆的两个典型的探地雷达地电模型进行了正演模拟,得到了正演剖面图,消除了边界反射后的雷达剖面能很好地指导工作人员对雷达实测剖面的地质解释,同时使正演研究更符合实际的地质情况.
In recent years, with the rapid development of digital processing and electronic technologies, the scope of GPR applications has expanded rapidly. Now it has covered archeology, exploration of mineral resources, hydrogeological investigation, geotechnical investigation, non-destructive testing, engineering construction Physical structure investigation, military and many other fields to solve many practical engineering problems and become a powerful tool for shallow exploration.The theory of ground penetrating radar has obvious hysteresis compared with the actual application.However, the interpreter should achieve the accurate To analyze the actual data of GPR, we must know the characteristics of the radar reflection profile of the geological body in advance. Therefore, the GPR forward modeling technology of complex geoelectric model, which is the basis of inversion and interpretation, becomes a ground penetrating radar Theoretical study of one of the main contents of this paper Maxwell two curl equations as the basic starting point, the use of KSYee’s space grid model theory and the basic principle of time-domain finite difference method, derived two-dimensional space GPR forward Equations, and detailed analysis of the half-space step and half-time steps in differential format, and its radar wave electric field and magnetic field components And the C program which is related to each other on the computer.And then discuss the relationship between numerical dispersion and its causes, and by considering the characteristics of time-domain finite difference method and Yee’s grid, the ideal that meets the actual propagation law of GPR Dispersion relationship, the author made a GPR forward program, and calculated the Mur superabsorbent boundary conditions and without the radar geoelectric model under the boundary conditions, by contrast we can see that the conditions of the superabsorbent can be used to greatly reduce the cut-off boundary The interference wave can reach the effect of spreading in the infinite space with a limited area.Finally, using the self-made program, the authors performed two typical ground penetrating radar geoelectric models with “V” shape and five circles on the same slope The simulation results show that the radar profile after eliminating the boundary reflection can well guide the staff to interpret the geological interpretation of the radar profile and at the same time make the forward modeling more in line with the actual geological conditions.