摘要
We conducted a study on the numerical calculation and response analysis of a transient electromagnetic field generated by a ground source in geological media. One solution method, the traditional discrete image method, involves complex operation, and its digital filtering algorithm requires a large number of calculations. To solve these problems, we proposed an improved discrete image method, where the following are realized： the real number of the electromagnetic field solution based on the Gaver-Stehfest algorithm for approximate inversion, the exponential approximation of the objective kernel function using the Prony method, the transient electromagnetic field according to discrete image theory, and closed-form solution of the approximate coefficients. To verify the method, we tentatively calculated the transient electromagnetic field in a homogeneous model and compared it with the results obtained from the Hankel transform digital filtering method. The results show that the method has considerable accuracy and good applicability. We then used this method to calculate the transient electromagnetic field generated by a ground magnetic dipole source in a typical geoelectric model and analyzed the horizontal component response of the induced magnetic field obtained from the ＂ground excitation-stratum measurement method. We reached the conclusion that the horizontal component response of a transient field is related to the geoelectric structure, observation time, spatial location, and others. The horizontal component response of the induced magnetic field reflects the eddy current field distribution and its vertical gradient variation. During the detection of abnormal objects, positions with a zero or comparatively large offset were selected for the drill- hole measurements or a comparatively long observation delay was adopted to reduce the influence of the ambient field on the survey results. The discrete image method and forward calculation results in this paper can be used as references for relevant research.
针对地面场源在地层介质中所产生瞬变电磁场的数值计算方法和响应分析问题,本文展开研究.解法方面,传统离散镜像法多采用复数运算、数字滤波等方法计算量大,针对这些问题,提出一种改进的离散镜像方法:基于Gaver-Stehfest概率变换算法将电磁场解式实数化,选用Prony方法对目标核函数进行指数级数逼近,根据离散镜像原理和近似系数闭合式求解瞬变电磁场.通过试算均质模型中瞬变电磁场并对比汉克尔变换的数字滤波法所得结果,证明该方法有效且具有较好的精度和适用性.继而基于该方法计算地表磁偶源在典型地电模型中产生的瞬变电磁场,对"地面激发-地层中测量"方式所得感应磁场水平分量响应进行分析并得出结论瞬变场水平分量响应与地电结构、观测时间、空间位置等因素有关,感应磁场水平分量响应反映出涡流场分布及其垂向梯度变化情况,在探测异常体的工作中应尽量选取零偏移距、较大偏移距位置钻孔或用较大观测延时以减小背景场对勘测结果的影响。文中所用离散镜像方法与正演计算结论可为相关研究工作提供参考依据。
基金
supported by the Young Scientists Fund of the National Natural Science Foundation of China(No.41304082)
the Young Scientists Fund of the Natural Science Foundation of Hebei Province(No.D2014403011)
the Geological survey project of China Geological Survey(No.12120114090201)
