二维起伏地表条件下坐标变换法直流电场数值模拟

J4 ›› 2009, Vol. 39 ›› Issue (3): 528-534.

二维起伏地表条件下坐标变换法直流电场数值模拟

孙章庆^1|2|孙建国^1|2|张东良^1|2

1.吉林大学地球探测科学与技术学院|长春 130026;2.国土资源部应用地球物理综合解释理论开放实验室|长春 130026

收稿日期:2008-11-26 出版日期:2009-05-26 发布日期:2009-05-26
作者简介:孙章庆(1982-)|男|土家族|重庆黔江人|博士研究生|主要从事地震波传播理论与成像技术、电磁法正演模拟、计算地球物理等方面的研究|Tel：0431-88502592|E-mail:sun_zhangq@163.com
基金资助:
国家高新技术研究发展计划“863”探索导向类课题(2006AA06Z109)

2D DC Electric Field Numerical Modeling Including Surface Topography Using Coordinate Transformation Method

SUN Zhang-qing^{1, 2}, SUN Jian-guo^1,2, ZHANG Dong-liang^{1, 2}

1.College of GeoExploration Science and Technology, Jilin University, Changchun 130026,China;2.Laboratory for Integrated Geophysical Interpretation Theory,Ministry for Land and Resources, Changchun 130026,China

Received:2008-11-26 Online:2009-05-26 Published:2009-05-26

摘要/Abstract

摘要：

针对有限差分方法处理不规则边界难的问题，引入曲化平的思想，实现了一种新的基于坐标变换法的起伏地表条件下直流电场数值模拟方法。通过坐标变换，将物理空间的直流电场满足的基本方程和边界条件变换到计算空间；在计算空间内，利用有限差分方法对变换后的方程和边界条件进行数值分析；最后，将计算结果由计算空间映射回物理空间。精度分析及计算实例表明：本文提出的算法精度高，计算误差主要集中在源附近，且当源离地表较近时平均相对误差在1.0%以下；计算结果满足起伏地表条件下电位及视电阻率的分布规律。

关键词: 起伏地表, 直流电场, 坐标变换, 有限差分

Abstract:

It is difficult to deal with the irregular boundary using finite-difference method in the DC electric field numerical modeling including surface topography. To solve this problem, we introduce a coordinate transformation idea that transforms curved surface into horizontal topography and realize a new finite-difference method basing on the transform. We firstly transform the basic equations and boundary conditions of DC electric field from physical space to computation space, and then resolve the DC fields in computation space using finitedifference method. Finally, we retransform the computation results from computation space into physical space by coordinate mapping. Accuracy analysis and numerical tests show that the computation error congregates around the source, and the average relative error is less than 1% if the source is close enough to the irregular surface. Therefore, the new method has high accuracy, and the computation results can reveal the distribution characteristics of DC electric field and apparent resistivity including surface topography.

Key words: surface topography, DC electric field, coordinate transformation, finite-difference

中图分类号:

P631.3

孙章庆, 孙建国, 张东良. 二维起伏地表条件下坐标变换法直流电场数值模拟[J]. J4, 2009, 39(3): 528-534.

SUN Zhang-qiang, SUN Jian-guo, ZHANG Dong-liang. 2D DC Electric Field Numerical Modeling Including Surface Topography Using Coordinate Transformation Method[J]. J4, 2009, 39(3): 528-534.

参考文献

相关文章 15

[1]	刘明忱, 孙建国, 韩复兴, 孙章庆, 孙辉, 刘志强. 基于自适应加权广义逆矢量方向滤波估计地震同相轴倾角[J]. 吉林大学学报(地球科学版), 2018, 48(3): 881-889.
[2]	李建平, 翁爱华, 李世文, 李大俊, 李斯睿, 杨悦, 唐裕, 张艳辉. 基于球坐标系下有限差分的地磁测深三维正演[J]. 吉林大学学报(地球科学版), 2018, 48(2): 411-419.
[3]	李大俊, 翁爱华, 杨悦, 李斯睿, 李建平, 李世文. 地-井瞬变电磁三维交错网格有限差分正演及响应特性[J]. 吉林大学学报(地球科学版), 2017, 47(5): 1552-1561.
[4]	杨海燕, 岳建华, 徐正玉, 张华, 姜志海. 覆盖层影响下典型地-井模型瞬变电磁法正演[J]. 吉林大学学报(地球科学版), 2016, 46(5): 1527-1537.
[5]	李庆洋, 李振春, 黄建平, 李娜, 苏在荣. 基于贴体全交错网格的起伏地表正演模拟影响因素[J]. 吉林大学学报(地球科学版), 2016, 46(3): 920-929.
[6]	贲放, 刘云鹤, 黄威, 徐驰. 各向异性介质中的浅海海洋可控源电磁响应特征[J]. 吉林大学学报(地球科学版), 2016, 46(2): 581-593.
[7]	王卫平, 曾昭发, 李静, 吴成平. 频率域航空电磁法地形影响和校正方法[J]. 吉林大学学报(地球科学版), 2015, 45(3): 941-951.
[8]	郭振波，李振春. 起伏地表条件下频率-空间域声波介质正演模拟[J]. 吉林大学学报(地球科学版), 2014, 44(2): 683-693.
[9]	巩向博，吕庆田，韩立国，谭尘青. 起伏地表地震波场角度域照明分析[J]. 吉林大学学报(地球科学版), 2013, 43(2): 610-615.
[10]	周晓华, 陈祖斌, 曾晓献, 焦健. 交错网格有限差分法模拟微动信号[J]. J4, 2012, 42(3): 852-857.
[11]	周辉, 谢春临, 王尚旭, 李国发. 复杂地表地震资料叠前深度偏移成像[J]. J4, 2012, 42(1): 262-268.
[12]	孟庆生, 樊玉清, 张珂, 张盟. 高阶有限差分法管波传播数值模拟[J]. J4, 2011, 41(1): 292-298.
[13]	刘四新, 周俊峰, 吴俊军, 曾昭发, 万洪祥. 金属矿钻孔雷达探测的数值模拟[J]. J4, 2010, 40(6): 1479-1484.
[14]	孙章庆, 孙建国, 张东良. 2.5维起伏地表条件下坐标变换法直流电场数值模拟[J]. J4, 2010, 40(2): 425-431.
[15]	杨昊, 孙建国, 韩复兴, 马淑芳. 基于完全三叉树堆排序的波前扩展有限差分地震波走时快速算法[J]. J4, 2010, 40(1): 188-194.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed