A Finite Difference Method for Modeling the DC Electrical Potential Field Including Surface Topography

Jianguo Sun,Dongliang Zhang,Zhangqing Sun
DOI: https://doi.org/10.1190/1.3255842
2009-01-01
Abstract:PreviousNext No AccessSEG Technical Program Expanded Abstracts 2009A finite difference method for modeling the DC electrical potential field including surface topographyAuthors: Jianguo SunDongliang ZhangZhangqing SunJianguo SunCollege of Geoexploration Science and Technology, Department of Geophysics, Jilin UniversitySearch for more papers by this author, Dongliang ZhangCollege of Geoexploration Science and Technology, Department of Geophysics, Jilin UniversitySearch for more papers by this author, and Zhangqing SunCollege of Geoexploration Science and Technology, Department of Geophysics, Jilin UniversitySearch for more papers by this authorhttps://doi.org/10.1190/1.3255842 SectionsSupplemental MaterialAboutPDF/ePub ToolsAdd to favoritesDownload CitationsTrack CitationsPermissions ShareFacebookTwitterLinked InRedditEmail Abstract For modeling the direct current (DC) electrical potential field in the earth model with complex surface relief, we develop a finite difference (FD) method that avoids the direct numerical implementation of the surface boundary condition. The key point of the method is to determine the electrical potential at the points in the earth's surface by using the fact that the DC equipotential lines or surfaces are perpendicular to the earth's surface. In the practical implementation, the unit normal vector to the earth's surface is firstly determined at a given surface point and then the normal is extended into the earth until it reaches a point that is a regular grid point located in the nearest vicinity of the surface, or that is located in a regular grid square. Next, the DC potential at the endpoint of the extended normal line is computed by bilinear interpolation or by direct assignment, according to whether the endpoint is located in a regular grid square or coincident with a regular grid point lying in the nearest vicinity of the surface. Finally, the potential value at the endpoint is assigned to the considered surface point. After this assignment step is finished, a linear equation system can be established by using the conventional FD stencils and by using the computed surface potential values. Numerical tests in 2D and 2.5D cases show that the method presented here is simple, flexible, efficient, and accurate enough for solving the modeling problem in DC electrical methods when the earth's surface has a complex geometry.Permalink: https://doi.org/10.1190/1.3255842FiguresReferencesRelatedDetails SEG Technical Program Expanded Abstracts 2009ISSN (print):1052-3812 ISSN (online):1949-4645Copyright: 2009 Pages: 4338 publication data© 2009 Copyright © 2009 Society of Exploration GeophysicistsPublisher:Society of Exploration Geophysicists HistoryPublished Online: 14 Oct 2009 CITATION INFORMATION Jianguo Sun, Dongliang Zhang, and Zhangqing Sun, (2009), "A finite difference method for modeling the DC electrical potential field including surface topography," SEG Technical Program Expanded Abstracts : 664-668. https://doi.org/10.1190/1.3255842 Plain-Language Summary PDF DownloadLoading ...
What problem does this paper attempt to address?