-
采用矢量有限元法实现了三维电磁扩散场数值模拟, 并成功将其应用在大地电磁的正演研究中. 为灵活精确地拟合起伏地形和地下不规则构造, 采用由不规则四面体单元组成的非结构化网格, 可根据模型设计的需要调整网格的大小. 引入了基于二次场理论, 将解析的一次场从总场中扣除, 直接计算二次场, 使得误差仅局限于相对较小的二次场, 以提高总场计算精度. 常规的节点有限元法不满足电性分界面上法向电场不连续和无源区单元内电流密度无散, 违反麦克斯韦方程组. 为克服节点有限元法的弊端, 使用矢量有限元法求解基于二次电场的偏微分方程. 另外, 在算法设计中, 考虑了磁导率参数的变化, 可以模拟磁导率不均匀的模型. 通过与COMMEMI模型已发表的结果对比, 证明了本文算法的正确性和精确性. 为突显非结构网格优势, 计算了椭球异常体模型和任意地形模型的MT响应, 并详细讨论了地形和磁化效应对三维数值模拟结果的影响.We present a newly developed algorithm for three-dimensional (3D) magnetotelluric model based on the vector finite element method. In this paper, unstructured grids which are composed of irregular tetrahedrons are used in our finite element model, which can be refined locally and adaptively according to the complex geometry of computational domain or subsurface structure. For obtaining more accurate solutions, secondary field rather than total field is numerically computed, which makes the errors limited to relatively small secondary field. Traditional node-based finite element method does not satisfy the condition that the normal electrical field is discontinuous at the interface of electric separatrix and the electrical current density is divergence-free throughout the regions without source, which obviously violates Maxwell equations. In order to overcome these drawbacks of node-based finite element, vector finite element method is employed to solve the secondary field-based partial differential equation. Moreover, in this study, the heterogeneous permeability is taken into consideration in our algorithm, as a consequence, which can deal with heterogeneous permeability model. The accuracy of our approach is verified by comparing with previously published numerical simulations of a COMMEMI-3D model. The advantages of our approach are also illustrated by the numerical simulations of model with arbitrary topography and complicated anomalous body. In addition, the simulation results of the irregular anomaly and the complex model of arbitrary terrain are both ideal results. It proves that the effect of 3D terrain is more serious and complex than that of two-dimensional (2D) terrain, and dealing with 3D models or data with 2D algorithm may bring in large errors. In the area where the magnetic permeability is abnormal, the magnetic permeability has an important influence on the numerical simulation results, and the magnetic permeability must be treated as an independent parameter in the magnetotellurics survey.
-
Keywords:
- electromagnetic diffusion /
- vector finite element method /
- numerical simulation /
- magnetization /
- magnetotellurics
[1] Hohmann G W 1983 Surv. Geophys. 6 27Google Scholar
[2] Zhdanov M S, Lee S K, Yoshioka K 2006 Geophysics 71 G333Google Scholar
[3] 陈桂波, 汪宏年, 姚敬金, 韩子夜 2009 物理学报 58 3848Google Scholar
Chen G B, Wang H N, Yao J J, Han Z Y 2009 Acta Phys. Sin. 58 3848Google Scholar
[4] 谭捍东, 余钦范, John B 2003 地球物理学报 46 705Google Scholar
Tang H D, Yu Q F, John B 2003 Chin. J. Geophys. 46 705Google Scholar
[5] Streich R 2009 Geophysics 74 F95Google Scholar
[6] Smith J T 1996 Geophysics 61 1308Google Scholar
[7] Wang W, Wu X P, Spitzer K 2013 Geophys. J. Int. 193 734Google Scholar
[8] Coggon J H 1971 Geophysics 36 132Google Scholar
[9] Fox R C, Hohmann G W, Killpack T J 1980 Geophysics 45 75Google Scholar
[10] Wannamaker P E, Stodt J A, Rijo L 1987 Geophys. J. Int. 88 277Google Scholar
[11] 李勇, 吴小平, 林品荣 2013 地球物理学进展 27 2452
Li Y, Wu X P, Lin P R 2013 Prog. Geophys. 27 2452
[12] Liu C, Ren Z, Tang J 2008 Appl. Geophys. 5 170Google Scholar
[13] Key K, Weiss C 2006 Geophysics 71 G291Google Scholar
[14] Franke A, Börner R U, Spitzer K 2007 Geophys. J. Int. 171 71Google Scholar
[15] Jin J 2002 The Finite Element Method in Electromagnetics (New York: John Wiley & Sons)
[16] 柳建新, 童孝忠, 郭荣文 2012 大地电磁测深法勘探 (北京: 科学出版社) 第48−50页
Liu J X, Tong X Z, Guo R W 2012 The Magnetotelluric Sounding Method (Beijing: Science Press) pp48−50
[17] Unsworth M J 2010 Surv. Geophys. 31 137Google Scholar
[18] Rao C K, Jones A G, Moorkamp M 2014 Geophys. J. Int. 198 737Google Scholar
[19] Juanatey M, Hubert J, Tryggvason A 2013 Geophys. Prospect. 61 200Google Scholar
[20] Egger A E, Jonathan M G, Glen J M G 2014 Geol. Soc. Am. Bull. 126 523Google Scholar
[21] Abdelzaher M, Nishijima J, Saibi H 2012 Pure Appl. Geophys. 169 1679Google Scholar
[22] 徐世浙 1994 地球物理中的有限单元法 (北京: 科学出版社) 第220−224页
Xu S Z 1994 FEM in Geophysics (Beijing: Science Press) pp220−224 (in Chinese)
[23] Farquharson C G, Miensopust M P 2011 J, Appl. Geophys. 75 699Google Scholar
[24] Zhdanov M S, Varentsov I M, Weaver J T 1997 J. Appl. Geophys. 37 133Google Scholar
[25] Mitsuhata Y, Uchida T 2004 Geophysics 69 108Google Scholar
[26] Ren Z, Kalscheuer, Greenhalgh S 2014 Geophysics 79 E255Google Scholar
[27] Mukherjee S, Everett M E 2011 Geophysics 76 F215Google Scholar
-
-
[1] Hohmann G W 1983 Surv. Geophys. 6 27Google Scholar
[2] Zhdanov M S, Lee S K, Yoshioka K 2006 Geophysics 71 G333Google Scholar
[3] 陈桂波, 汪宏年, 姚敬金, 韩子夜 2009 物理学报 58 3848Google Scholar
Chen G B, Wang H N, Yao J J, Han Z Y 2009 Acta Phys. Sin. 58 3848Google Scholar
[4] 谭捍东, 余钦范, John B 2003 地球物理学报 46 705Google Scholar
Tang H D, Yu Q F, John B 2003 Chin. J. Geophys. 46 705Google Scholar
[5] Streich R 2009 Geophysics 74 F95Google Scholar
[6] Smith J T 1996 Geophysics 61 1308Google Scholar
[7] Wang W, Wu X P, Spitzer K 2013 Geophys. J. Int. 193 734Google Scholar
[8] Coggon J H 1971 Geophysics 36 132Google Scholar
[9] Fox R C, Hohmann G W, Killpack T J 1980 Geophysics 45 75Google Scholar
[10] Wannamaker P E, Stodt J A, Rijo L 1987 Geophys. J. Int. 88 277Google Scholar
[11] 李勇, 吴小平, 林品荣 2013 地球物理学进展 27 2452
Li Y, Wu X P, Lin P R 2013 Prog. Geophys. 27 2452
[12] Liu C, Ren Z, Tang J 2008 Appl. Geophys. 5 170Google Scholar
[13] Key K, Weiss C 2006 Geophysics 71 G291Google Scholar
[14] Franke A, Börner R U, Spitzer K 2007 Geophys. J. Int. 171 71Google Scholar
[15] Jin J 2002 The Finite Element Method in Electromagnetics (New York: John Wiley & Sons)
[16] 柳建新, 童孝忠, 郭荣文 2012 大地电磁测深法勘探 (北京: 科学出版社) 第48−50页
Liu J X, Tong X Z, Guo R W 2012 The Magnetotelluric Sounding Method (Beijing: Science Press) pp48−50
[17] Unsworth M J 2010 Surv. Geophys. 31 137Google Scholar
[18] Rao C K, Jones A G, Moorkamp M 2014 Geophys. J. Int. 198 737Google Scholar
[19] Juanatey M, Hubert J, Tryggvason A 2013 Geophys. Prospect. 61 200Google Scholar
[20] Egger A E, Jonathan M G, Glen J M G 2014 Geol. Soc. Am. Bull. 126 523Google Scholar
[21] Abdelzaher M, Nishijima J, Saibi H 2012 Pure Appl. Geophys. 169 1679Google Scholar
[22] 徐世浙 1994 地球物理中的有限单元法 (北京: 科学出版社) 第220−224页
Xu S Z 1994 FEM in Geophysics (Beijing: Science Press) pp220−224 (in Chinese)
[23] Farquharson C G, Miensopust M P 2011 J, Appl. Geophys. 75 699Google Scholar
[24] Zhdanov M S, Varentsov I M, Weaver J T 1997 J. Appl. Geophys. 37 133Google Scholar
[25] Mitsuhata Y, Uchida T 2004 Geophysics 69 108Google Scholar
[26] Ren Z, Kalscheuer, Greenhalgh S 2014 Geophysics 79 E255Google Scholar
[27] Mukherjee S, Everett M E 2011 Geophysics 76 F215Google Scholar
计量
- 文章访问数: 8765
- PDF下载量: 71
- 被引次数: 0