摘要
本文给出一种用有限差分法(FD)计算二维(2D)介质大地电磁感应的方法。该方法在计算区域的两侧及上部用了J. T. Weaver等人多次改进后的边条件。为了使用叠代法求解离散化得到的高阶线性方程组以节省消去法所必需的内存和计算时间,提出了地球表面边条件的一种高精度差分格式。本文的方法可以计算信号周期为10^(-1)——10~3秒的电磁感应响应。
A finite difference algorithm to calculate two dimensional (2D) geoelectromagnetic induction responses is given. Boundary conditions improved for several times by J. T. Weaver et al. are employed at both sides and upper edge of the grid. In order to use iterative method to solve the large linear system obtained by discretization for the purpose of reducing computer storage and CPU time required by the use of elimination, higher-order difference schemes of boundary conditions at the surface of the earth is proposed. Ceoelectromagnetic induction responses with source periods ranging from 0.1 sec. to 1000 sec. could be calculated by this approach.
出处
《北京大学学报（自然科学版）》
CAS
CSCD
北大核心
1992年第3期352-360,共9页
Acta Scientiarum Naturalium Universitatis Pekinensis
关键词
大地
电磁感应
二维
导体
叠代法
Geoelectromagnetic induction
Two-dimensional conductivities
Finite difference modelling
Iterative method
