FarquharsonColin G.3,
刘长胜1,2,,
1. 吉林大学仪器科学与电气工程学院, 长春 130061
2. 吉林大学地球信息探测仪器教育部重点实验室, 长春 130061
3. Department of Earth Sciences, Memorial University of Newfoundland, St. John's, Canada A1B3S7
基金项目: 国家自然科学基金面上项目（42074152）以及国家自然科学基金国家重大科研仪器研制项目（41827803），国家留学基金委（201706170173）联合资助
详细信息
作者简介: 张铭, 女, 1989年生, 博士研究生, 主要从事电磁法正反演方法研究.E-mail:zhangm_sunflower@163.com
通讯作者: 刘长胜, 男, 1979年生, 教授, 主要从事电磁探测技术及仪器的研究.E-mail:liuchangsheng@jlu.edu.cn
中图分类号: P631收稿日期:2020-02-14
修回日期:2020-07-28
上线日期:2021-01-10
2.5-D forward modeling of the frequency-domain ground-airborne electromagnetic response in areas with topographic relief
ZHANG Ming1,2,3,,Farquharson Colin G.3,
LIU ChangSheng1,2,,
1. College of Instrumentation and Electrical Engineering, Jilin University, Changchun 130061, China
2. Key Lab of Geo-Exploration Instrumentation, Ministry of Education, Changchun 130061, China
3. Department of Earth Sciences, Memorial University of Newfoundland, St. John's, A1 B3 S7, Canada
More Information
Corresponding author: LIU ChangSheng,E-mail:liuchangsheng@jlu.edu.cn
MSC: P631--> Received Date: 14 February 2020
Revised Date: 28 July 2020
Available Online: 10 January 2021
摘要
摘要:频率域地空电磁探测方法是指在地面布设人工场源，在空中测量电磁场的一种高效的地球物理勘探技术.该方法具有大范围、高通过性、快速测量的优势，尤其适合崎岖山地、沙漠、沼泽、海陆交互带等复杂地貌区域的资源勘查.但是这些地区的地形起伏通常较大，因此分析地形对地空电磁响应的影响具有重要意义.本文利用有限元法对频率域地空电磁响应进行了正演计算，分析了起伏地表条件下的频率域地空电磁响应特征.首先利用傅里叶变换将2.5维问题转化成二维问题，利用伽辽金加权余量法推导了相应的离散有限元方程组.采用任意四边形单元对区域进行不均匀网格剖分，源和异常体附近网格加密处理，保证计算精度，远离目标区域网格逐渐稀疏，模拟无穷远边界，降低对计算资源的要求.在单元内进行插值，将有限元方程组变换为线性方程组，采用总场算法，利用具有一定面积的伪δ函数表达源电流分布，源项近似为分布在以电偶极源为中心的25个节点上.通过求解线性方程组得到波数域电磁响应，再对波数域电磁场响应进行反傅里叶变换从而获得空间域2.5维频率域电磁场值.通过对比2.5维正演结果与均匀半空间解析解，验证了本文算法的精度，同时本文还对地空电磁场与地面电磁场的响应特性进行了对比.
关键词: 地空电磁法/
2.5维/
正演/
起伏地表/
总场算法
Abstract: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.
Key words:Ground-airborne electromagnetic method/
2.5D/
Forward modeling/
Topographic surface/
Total field algorithm
PDF全文下载地址:
http://www.geophy.cn/data/article/export-pdf?id=5fedbacabfedb53bd868238b
