Mathematical Problems in Engineering | Volume 2017 ,2017-05-15 |
Finite Element Method for Modeling 3D Resistivity Sounding on Anisotropic Geoelectric Media | |
Research Article | |
Tao Song ^{1} , ^{2} Yun Liu ^{1} Yun Wang ^{1} | |
Show affiliations | |
DOI：10.1155/2017/8027616 | |
Received 2017-01-11, accepted for publication 2017-04-09, Published 2017-04-09 | |
摘要
A 3D DC finite element method forward program is developed in this paper for anisotropic geoelectric media. Both total and secondary field approaches have been implemented. In this paper, we focused on the structured grid scheme. The modeling shows that the symmetry of the structure grid determines the symmetry of the response potential around the point source for an anisotropic half-space. Through numerical modeling with three kinds of coarse meshes by the total field approach and secondary field approach separately, a higher accuracy can be achieved via the secondary field approach. And the relative error between the numerical solution and the analytical solution is less than 2%. Modeling results contrasting with previous scholars also verify the correctness of the algorithm. Then, the numerical results of 3D models with anisotropic properties were presented and compared to models with isotropic properties. These results clearly illustrated the strong effect of anisotropy and the problems in interpretation if anisotropy was not properly addressed. This work will establish the foundation for our future effort of building a 3D inversion program with arbitrary anisotropy.
授权许可
Copyright © 2017 Tao Song et al. 2017
This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
图表
The global region subdivision with coarse area. (a) 2D global mesh design; (b) 3D global mesh design.
The global region subdivision with coarse area. (a) 2D global mesh design; (b) 3D global mesh design.
Subdivision of the region. (a) Five-tetrahedron scheme A; (b) five-tetrahedron scheme B; (c) 3D global mesh with scheme A; (d) 3D global mesh with scheme B; (e) 3D global mesh with schemes A and B.
Subdivision of the region. (a) Five-tetrahedron scheme A; (b) five-tetrahedron scheme B; (c) 3D global mesh with scheme A; (d) 3D global mesh with scheme B; (e) 3D global mesh with schemes A and B.
Subdivision of the region. (a) Five-tetrahedron scheme A; (b) five-tetrahedron scheme B; (c) 3D global mesh with scheme A; (d) 3D global mesh with scheme B; (e) 3D global mesh with schemes A and B.
Subdivision of the region. (a) Five-tetrahedron scheme A; (b) five-tetrahedron scheme B; (c) 3D global mesh with scheme A; (d) 3D global mesh with scheme B; (e) 3D global mesh with schemes A and B.
Subdivision of the region. (a) Five-tetrahedron scheme A; (b) five-tetrahedron scheme B; (c) 3D global mesh with scheme A; (d) 3D global mesh with scheme B; (e) 3D global mesh with schemes A and B.
Forward modeling accuracy in percentage for three kinds of subdivision relative to analytical solutions by total potential approach: (a) by grid shown in Figure 2(c), (b) by grid shown in Figure 2(d), and (c) by grid shown in Figure 2(e).
Forward modeling accuracy in percentage for three kinds of subdivision relative to analytical solutions by total potential approach: (a) by grid shown in Figure 2(c), (b) by grid shown in Figure 2(d), and (c) by grid shown in Figure 2(e).
Forward modeling accuracy in percentage for three kinds of subdivision relative to analytical solutions by total potential approach: (a) by grid shown in Figure 2(c), (b) by grid shown in Figure 2(d), and (c) by grid shown in Figure 2(e).
Forward modeling accuracy in percentage for three kinds of coarse meshes relative to analytical solutions by total potential approach: (a) by Mesh 1, (b) by Mesh 2, and (c) by Mesh 3.
Forward modeling accuracy in percentage for three kinds of coarse meshes relative to analytical solutions by total potential approach: (a) by Mesh 1, (b) by Mesh 2, and (c) by Mesh 3.
Forward modeling accuracy in percentage for three kinds of coarse meshes relative to analytical solutions by total potential approach: (a) by Mesh 1, (b) by Mesh 2, and (c) by Mesh 3.
Forward modeling accuracy in percentage relative to analytical solutions with three kinds of sparse meshes by secondary potential approach: (a) by Mesh 1, (b) by Mesh 2, and (c) by Mesh 3.
Forward modeling accuracy in percentage relative to analytical solutions with three kinds of sparse meshes by secondary potential approach: (a) by Mesh 1, (b) by Mesh 2, and (c) by Mesh 3.
Forward modeling accuracy in percentage relative to analytical solutions with three kinds of sparse meshes by secondary potential approach: (a) by Mesh 1, (b) by Mesh 2, and (c) by Mesh 3.
A two-layer model.
A two-layer model with azimuthal anisotropy.
Apparent resistivity of pole-pole array along x-direction for model shown in Figure 7.
A two-layer model consisting of a TTI overburden layer with variable dip over an isotropic half-space.
Apparent resistivities along x- and y-direction for pole-pole array with changing dipping angles.
A 2D anisotropic model.
Dipole-dipole array apparent resistivity of model shown in Figure 11. The resistivity of the isotropic model is described as follows: the resistivity of the half-space is 20 Ω⋅m and the resistivity of the anomaly is 1 Ω⋅m. (a) Dipole spacing r=4 m. (b) Dipole spacing r=20 m.
Dipole-dipole array apparent resistivity of model shown in Figure 11. The resistivity of the isotropic model is described as follows: the resistivity of the half-space is 20 Ω⋅m and the resistivity of the anomaly is 1 Ω⋅m. (a) Dipole spacing r=4 m. (b) Dipole spacing r=20 m.
Comparison of 3D and 2D FEM numerical solutions of the model shown in Figure 11. (a) Apparent resistivity pseudosections along the central line by 3D FEM. (b) Apparent resistivity pseudosections by 2D FEM.
Comparison of 3D and 2D FEM numerical solutions of the model shown in Figure 11. (a) Apparent resistivity pseudosections along the central line by 3D FEM. (b) Apparent resistivity pseudosections by 2D FEM.
(a) A 3D anisotropic cube in a homogeneous isotropic half-space. The principal resistivity of the cube is given by ρx/ρy/ρz = 10/10/100 Ω⋅m and the background resistivity is 100 (Ω⋅m). (b) 21 survey lines were deployed along both x- and y-direction and the solid lines outline the boundary of the buried cube.
(a) A 3D anisotropic cube in a homogeneous isotropic half-space. The principal resistivity of the cube is given by ρx/ρy/ρz = 10/10/100 Ω⋅m and the background resistivity is 100 (Ω⋅m). (b) 21 survey lines were deployed along both x- and y-direction and the solid lines outline the boundary of the buried cube.
Apparent resistivity pseudosection of the model shown in Figure 14. (a) Apparent resistivity pseudosection of line X11 for the isotropic model. (b) Apparent resistivity pseudosection of line X07 for the isotropic model. (c) Apparent resistivity pseudosection of line X11 for the anisotropic model. (d) Apparent resistivity pseudosection of line X07 for the anisotropic model.
Apparent resistivity pseudosection of the model shown in Figure 14. (a) Apparent resistivity pseudosection of line X11 for the isotropic model. (b) Apparent resistivity pseudosection of line X07 for the isotropic model. (c) Apparent resistivity pseudosection of line X11 for the anisotropic model. (d) Apparent resistivity pseudosection of line X07 for the anisotropic model.
Apparent resistivity pseudosection of the model shown in Figure 14. (a) Apparent resistivity pseudosection of line X11 for the isotropic model. (b) Apparent resistivity pseudosection of line X07 for the isotropic model. (c) Apparent resistivity pseudosection of line X11 for the anisotropic model. (d) Apparent resistivity pseudosection of line X07 for the anisotropic model.
Apparent resistivity pseudosection of the model shown in Figure 14. (a) Apparent resistivity pseudosection of line X11 for the isotropic model. (b) Apparent resistivity pseudosection of line X07 for the isotropic model. (c) Apparent resistivity pseudosection of line X11 for the anisotropic model. (d) Apparent resistivity pseudosection of line X07 for the anisotropic model.
Pseudosection of the apparent resistivity for dipole-dipole array when the dipping angle α=45°. (a) Pseudosection of apparent resistivity at X11. (b) Pseudosection of apparent resistivity at X07. (c) Pseudosection of apparent resistivity at Y11. (d) Pseudosection of apparent resistivity Y07.
Pseudosection of the apparent resistivity for dipole-dipole array when the dipping angle α=45°. (a) Pseudosection of apparent resistivity at X11. (b) Pseudosection of apparent resistivity at X07. (c) Pseudosection of apparent resistivity at Y11. (d) Pseudosection of apparent resistivity Y07.
Pseudosection of the apparent resistivity for dipole-dipole array when the dipping angle α=45°. (a) Pseudosection of apparent resistivity at X11. (b) Pseudosection of apparent resistivity at X07. (c) Pseudosection of apparent resistivity at Y11. (d) Pseudosection of apparent resistivity Y07.
Pseudosection of the apparent resistivity for dipole-dipole array when the dipping angle α=45°. (a) Pseudosection of apparent resistivity at X11. (b) Pseudosection of apparent resistivity at X07. (c) Pseudosection of apparent resistivity at Y11. (d) Pseudosection of apparent resistivity Y07.
Pseudosection of the apparent resistivity for dipole-dipole array when the dipping angle α=90°. (a) Pseudosection of apparent resistivity at X11. (b) Pseudosection of apparent resistivity at X07. (c) Pseudosection of apparent resistivity at Y11. (d) Pseudosection of apparent resistivity at Y07.
Pseudosection of the apparent resistivity for dipole-dipole array when the dipping angle α=90°. (a) Pseudosection of apparent resistivity at X11. (b) Pseudosection of apparent resistivity at X07. (c) Pseudosection of apparent resistivity at Y11. (d) Pseudosection of apparent resistivity at Y07.
Pseudosection of the apparent resistivity for dipole-dipole array when the dipping angle α=90°. (a) Pseudosection of apparent resistivity at X11. (b) Pseudosection of apparent resistivity at X07. (c) Pseudosection of apparent resistivity at Y11. (d) Pseudosection of apparent resistivity at Y07.
Pseudosection of the apparent resistivity for dipole-dipole array when the dipping angle α=90°. (a) Pseudosection of apparent resistivity at X11. (b) Pseudosection of apparent resistivity at X07. (c) Pseudosection of apparent resistivity at Y11. (d) Pseudosection of apparent resistivity at Y07.
通讯作者
Tao Song.State Key Laboratory of Ore Deposit Geochemistry, Institute of Geochemistry, Chinese Academy of Sciences, Guiyang 550081, China, cas.cn;University of Chinese Academy of Sciences, Beijing 100049, China, ucas.ac.cn.songtao@mail.gyig.ac.cn
推荐引用方式
Tao Song,Yun Liu,Yun Wang. Finite Element Method for Modeling 3D Resistivity Sounding on Anisotropic Geoelectric Media. Mathematical Problems in Engineering ,Vol.2017(2017)
您觉得这篇文章对您有帮助吗？
分享和收藏
参考文献
[1] | D. F. Pridmore, G. W. Hohmann, S. H. Ward, W. R. Sill. et al.(1981). An investigation of finite-element modeling for electrical and electromagnetic data in three dimensions. Geophysics.46(7):1009-1024. DOI: 10.1190/1.1441239. |
[2] | Y. Sasaki. (1994). 3-D resistivity inversion using the finite-element method. Geophysics.59(11):1839-1848. DOI: 10.1190/1.1441239. |
[3] | C. Yin. (2000). Geoelectrical inversion for a one-dimensional anisotropic model and inherent non-uniqueness. Geophysical Journal International.140(1):11-23. DOI: 10.1190/1.1441239. |
[4] | Z. Bing, S. A. Greenhalgh. (2001). Finite element three-dimensional direct current resistivity modelling: accuracy and efficiency considerations. Geophysical Journal International.145(3):679-688. DOI: 10.1190/1.1441239. |
[5] | Y. Li, K. Spitzer. (2005). Finite element resistivity modelling for three-dimensional structures with arbitrary anisotropy. Physics of the Earth and Planetary Interiors.150(1-3):15-27. DOI: 10.1190/1.1441239. |
[6] | W. Wang. (2013). Study of 3D arbitrary anisotropic resistivity modeling and interpretation using unstructured finite element methods [Ph.D. thesis]. DOI: 10.1190/1.1441239. |
[7] | C. Yin, P. Weidelt. (1999). Geoelectrical fields in a layered earth with arbitrary anisotropy. Geophysics.64(2):426-434. DOI: 10.1190/1.1441239. |
[8] | A. Dey, H. F. Morrison. (1979). Resistivity modelling for arbitrarily shaped two-dimensional structures. Geophysical Prospecting.27(1):106-136. DOI: 10.1190/1.1441239. |
[9] | S. A. Greenhalgh, B. Zhou, M. Greenhalgh, L. Marescot. et al.(2009). Explicit expressions for the Fréchet derivatives in 3D anisotropic resistivity inversion. Geophysics.74(3):F31-F43. DOI: 10.1190/1.1441239. |
[10] | P. Li, N. F. Uren. (1997). Analytical solution for the point source potential in an anisotropic 3-D half-space. I. Two-horizontal-layer case. Mathematical and Computer Modelling.26(5):9-27. DOI: 10.1190/1.1441239. |
[11] | T. Lowry, M. B. Allen, P. N. SHive. (1989). Singularity removal: a refinement of resistivity modeling techniques. Geophysics.54(3):700-707. DOI: 10.1190/1.1441239. |
[12] | B. Zhou, M. Greenhalgh, S. A. Greenhalgh. (2009). 2.5-D/3-D resistivity modelling in anisotropic media using Gaussian quadrature grids. Geophysical Journal International.176(1):63-80. DOI: 10.1190/1.1441239. |
[13] | K. Spitzer. (1995). A 3-D finite-difference algorithm for DC resistivity modelling using conjugate gradient methods. Geophysical Journal International.123(3):903-914. DOI: 10.1190/1.1441239. |
浏览 | 143次 |
下载全文 | 35次 |
评分次数 | 0次 |
用户评分 | 0.0分 |
分享 | 0次 |