基于矢量有限元的大地电磁快速三维正演研究

顾观文, 武晔, 石砚斌

Research on fast three-dimensional forward algorithm of magnetotelluric sounding based on vector finite element

GU Guan-Wen, WU Ye, SHI Yan-Bin

表5 K型层状模型矢量有限元解与解析解对比

Table 5 Comparison between vector finite element solution and analytical solution of K-type layered model

频点 /Hz 视电阻率(ρxy,ρyx)/(Ω·m) 相位φ/(°) 矢量有限元解 解析解 误差/% 矢量有限元解 解析解 误差/% 10000 99.4928 99.99995 0.507154 44.9823 44.99998 0.039291 8000 99.5179 100.0004 0.482464 44.9795 44.99994 0.04542 5000 99.5522 99.9958 0.443615 44.9822 45.00155 0.042987 2000 99.9938 100.3013 0.306591 44.9731 44.9962 0.05134 1000 98.5514 99.11594 0.569571 45.4799 45.46765 0.02695 500 94.8705 94.84442 0.0275 43.8641 44.11071 0.55908 200 111.446 109.5436 1.73666 42.1518 41.35312 1.93136 100 124.19 128.152 3.091671 47.226 45.89309 2.90438 50 115.34 122.0064 5.46397 54.8037 54.506 0.54618 10 55.4147 56.83843 2.504864 64.9559 65.42769 0.721085 5 38.6174 39.2693 1.660063 64.5842 64.96375 0.584241 2 25.6211 25.88921 1.03559 61.7858 62.03061 0.394665 1 20.052 20.20837 0.773778 58.9626 59.13016 0.28337 0.5 16.5801 16.68206 0.611172 56.1309 56.24364 0.200444 0.1 12.6035 12.65607 0.415342 50.8822 50.92886 0.091626 0.05 11.7748 11.82107 0.391403 49.346 49.36433 0.037132 0.01 10.7493 10.77999 0.284731 46.9822 47.06061 0.166611 0.005 10.5349 10.54576 0.102932 46.4109 46.47605 0.140175 0.001 10.2412 10.24057 0.00616 45.6579 45.67163 0.030065 0.0005 10.17 10.16953 0.00462 45.4711 45.47688 0.012716 0.0001 10.0755 10.07547 0.00035 45.2136 45.21445 0.001873