Abstract
Three-dimensional (3-D) electromagnetic (EM) modeling algorithm has been developed using finite element method (FEM) to acquire more efficient interpretation techniques of EM data. When FEM based on nodal elements is applied to EM problem, spurious solutions, so called 'vector parasite', are occurred due to the discontinuity of normal electric fields and may lead the completely erroneous results. Among the methods curing the spurious problem, this study adopts vector element of which basis function has the amplitude and direction. To reduce computational cost and required core memory, complex bi-conjugate gradient (CBCG) method is applied to solving complex symmetric matrix of FEM and point Jacobi method is used to accelerate convergence rate. To verify the developed 3-D EM modeling algorithm, its electric and magnetic field for a layered-earth model are compared with those of layered-earth solution. As we expected, the vector based FEM developed in this study does not cause ny vector parasite problem, while conventional nodal based FEM causes lots of errors due to the discontinuity of field variables. For testing the applicability to high frequencies 100 MHz is used as an operating frequency for the layer structure. Modeled fields calculated from developed code are also well matched with the layered-earth ones for a model with dielectric anomaly as well as conductive anomaly. In a vertical electric dipole source case, however, the discontinuity of field variables causes the conventional nodal based FEM to include a lot of errors due to the vector parasite. Even for the case, the vector based FEM gave almost the same results as the layered-earth solution. The magnetic fields induced by a dielectric anomaly at high frequencies show unique behaviors different from those by a conductive anomaly. Since our 3-D EM modeling code can reflect the effect from a dielectric anomaly as well as a conductive anomaly, it may be a groundwork not only to apply high frequency EM method to the field survey but also to analyze the fold data obtained by high frequency EM method.
유한요소법을 이용한 전자기장의 3차원 모델링은 전자기장의 연속조건을 수치해가 만족하지 못함으로 인해서 발생하는 벡터 기생해(vector parasite)의 문제점을 가지고 있다. 이 연구에서는 벡터 기생해로 인한 오차를 줄이기 위해, 기저함수가 크기와 방향을 가지는 벡터요소를 도입하였다. 유한요소 행렬식은 complex BCG법을 적용하여 계산시간과 기억용량을 줄이고자 하였으며, 반복적인 해의 수렴속도 향상을 위해서 Point Jacobi법을 적용하였다. 개발된 알고리듬을 수직 전기 쌍극자 송신원을 이용한 층서구조 모형에 적용하여 이를 층서구조의 해와 비교함으로써 수치 모델링 알고리듬의 타당성을 검증하였으며, 이 과정에서 기존의 유한요소법에서 발생하는 벡터 기생해의 문제점이 벡터요소를 이용하는 경우에는 나타나지 않는 것을 확인하였다. 개발된 3차원 전자탐사 모델링 기법의 고주파수 영역으로의 적용성을 고찰하기 위하여, 100MHz의 수직 자기 쌍극자 송신원을 이용한 모델링을 유전율 이상층이 존재하는 층서구조 모형에 적용하여, 이를 층서구조 해와 비교하여 알고리듬의 타당성을 확인하였다. 검증된 3차원 전자탐사 모델링 기법을 유전율 이상체에 적용하여 이상체 주변에서의 전기장의 반응을 공간적으로 살펴보았다 이 연구에서 개발된 벡터요소를 사용한 3차원 고주파 전자탐사 모델링 기법은 기존의 전기전도도 이상체 뿐만 아니라 유전율 이상체에 대한 모델링을 가능하게 하여, 고주파 전자탐사법의 새로운 적용 및 해석의 기반을 제공할 수 있을 것으로 기대된다.