2.5维起伏地表条件下频率域地空电磁正演模拟

频率域地空电磁探测方法是指在地面布设人工场源,在空中测量电磁场的一种高效的地球物理勘探技术.该方法具有大范围、高通过性、快速测量的优势,尤其适合崎岖山地、沙漠、沼泽、海陆交互带等复杂地貌区域的资源勘查.但是这些地区的地形起伏通常较大,因此分析地形对地空电磁响应的影响具有重要意义.本文利用有限元法对频率域地空电磁响应进行了正演计算,分析了起伏地表条件下的频率域地空电磁响应特征.首先利用傅里叶变换将2.5维问题转化成二维问题,利用伽辽金加权余量法推导了相应的离散有限元方程组.采用任意四边形单元对区域进行不均匀网格剖分,源和异常体附近网格加密处理,保证计算精度,远离目标区域网格逐渐稀疏,模拟无穷远边界,降低对计算资源的要求.在单元内进行插值,将有限元方程组变换为线性方程组,采用总场算法,利用具有一定面积的伪δ函数表达源电流分布,源项近似为分布在以电偶极源为中心的25个节点上.通过求解线性方程组得到波数域电磁响应,再对波数域电磁场响应进行反傅里叶变换从而获得空间域2.5维频率域电磁场值.通过对比2.5维正演结果与均匀半空间解析解,验证了本文算法的精度,同时本文还对地空电磁场与地面电磁场的响应特性进行了对比.

2.5-D forward modeling of the frequency-domain ground-airborne electromagnetic response in areas with topographic relief

Ground-airborne frequency-domain electromagnetic method is an effective and efficient geophysical prospecting method. It usually places the controlled source on the ground and measure the electromagnetic field in the air. This method can achieve fast measurement in large scale, is especially suitable for resource exploration in rugged mountains, deserts, swamps, sea-land interaction zones and other areas with complex topography relief. However, in these areas, the topography relief is serious, so the study of the topographic effect on ground-airborne electromagnetic response is meaningful. In this paper, the 2.5D frequency domain ground-airborne electromagnetic response of topography earth model was calculated based on the finite element method and the response characteristic is analyzed. Fourier transform was applied to the Maxwell's equations to transform a 2.5D problem into a 2D problem, and then the discrete finite element equation was deduced using the Galerkin weighted residual method. The surveyed region is subdivided into many arbitrary quadrilateral elements with different sizes, the ones region near the source and the target have smaller size to ensure the accuracy of numerical calculation, the ones away from the target region have bigger size to simulate infinite boundary and reduce the requirements for computing resource. Interpolation is conducted in the element to transform the finite element equations to the linear system of equations. Using the total field algorithm, the pseudo-delta function with a certain area is employed to approximate the distribution of the source current, the source terms are approximately distributed at 25 nodes centered on the electric dipole source. By solving the equation, the solution of electromagnetic field in wavenumber domain was obtained, then the solution of electromagnetic field in spatial domain can be calculated by the Invers Fourier Transformation. The accuracy of our algorithm was verified by the comparison of our 2.5D forward modeling results and the analytical results. We also compared the electromagnetic field response characteristics measured in the air and on the ground.