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

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]. 吉林大学学报(地球科学版), 2025, 55(2): 627-645.
[2]	嵇艳鞠, 邓昌伟, 王宇航, 刘航, 吴琼. 基于时域有限差分法的电性源感应-极化效应三维数值模拟[J]. 吉林大学学报(地球科学版), 2024, 54(4): 1373-1382.
[3]	张志佳, 孙章庆, 王瑞湖, 韩复兴, 陆济璞, 刘亚光, 岑文攀, 高正辉. 复杂地表起伏多重变加密网格有限差分法波动方程数值模拟[J]. 吉林大学学报(地球科学版), 2023, 53(2): 598-.
[4]	刘志强, 黄磊, 李钢柱, 许磊, 李文宝. 崎岖海底对下伏水平地层反射波特征的影响[J]. 吉林大学学报(地球科学版), 2023, 53(1): 274-282.
[5]	韩丽, 胥铁潇, 金胜昔, 董文宇, 嵇艳鞠. 基于双电层模型的频域震电响应数值模拟 [J]. 吉林大学学报(地球科学版), 2022, 52(3): 737-743.
[6]	韩复兴, 王若雯, 孙章庆, 高正辉. 地震声波数值模拟中人工边界条件的差别与组合[J]. 吉林大学学报(地球科学版), 2022, 52(1): 261-.
[7]	岳崇旺, 王祝文. Kelvin-Voigt黏弹性介质井孔声场有限差分数值模拟与波形特征[J]. 吉林大学学报(地球科学版), 2021, 51(4): 1268-1275.
[8]	邵广周, 赵凯鹏, 吴华. 基于MPI的面波有限差分正演模拟[J]. 吉林大学学报(地球科学版), 2020, 50(1): 294-303.
[9]	刘明忱, 孙建国, 韩复兴, 孙章庆, 孙辉, 刘志强. 基于自适应加权广义逆矢量方向滤波估计地震同相轴倾角[J]. 吉林大学学报(地球科学版), 2018, 48(3): 881-889.
[10]	李建平, 翁爱华, 李世文, 李大俊, 李斯睿, 杨悦, 唐裕, 张艳辉. 基于球坐标系下有限差分的地磁测深三维正演[J]. 吉林大学学报(地球科学版), 2018, 48(2): 411-419.
[11]	李大俊, 翁爱华, 杨悦, 李斯睿, 李建平, 李世文. 地-井瞬变电磁三维交错网格有限差分正演及响应特性[J]. 吉林大学学报(地球科学版), 2017, 47(5): 1552-1561.
[12]	杨海燕, 岳建华, 徐正玉, 张华, 姜志海. 覆盖层影响下典型地-井模型瞬变电磁法正演[J]. 吉林大学学报(地球科学版), 2016, 46(5): 1527-1537.
[13]	李庆洋, 李振春, 黄建平, 李娜, 苏在荣. 基于贴体全交错网格的起伏地表正演模拟影响因素[J]. 吉林大学学报(地球科学版), 2016, 46(3): 920-929.
[14]	贲放, 刘云鹤, 黄威, 徐驰. 各向异性介质中的浅海海洋可控源电磁响应特征[J]. 吉林大学学报(地球科学版), 2016, 46(2): 581-593.
[15]	王卫平, 曾昭发, 李静, 吴成平. 频率域航空电磁法地形影响和校正方法[J]. 吉林大学学报(地球科学版), 2015, 45(3): 941-951.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed