A finite element method is developed for simulating frequency domain electromagnetic responses due to a dipole source in the 2-D conductive structures. Computing costs are considerably minimized by reducing the full t...A finite element method is developed for simulating frequency domain electromagnetic responses due to a dipole source in the 2-D conductive structures. Computing costs are considerably minimized by reducing the full three-dimensional problem to a series of two-dimensional problems. This is accomplished by transforming the problem into y-wave number (Ky) domain using Fourier transform and the y-axis is parallel to the structural strike. In the Ky domain, two coupled partial differential equations for magnetic field Hy and electric field Ey are derived. For a specific value of Ky, the coupled equations are solved by the finite element method with isoparametric elements in the x-z plane. Application of the inverse Fourier transform to the Ky, domain provides the electric and magnetic fields in real space. The equations derived can be applied to general complex two-dimensional structures containing either electric or magnetic dipole source in any direction. In the modeling of the electromagnetic measurement, we adopted a pseudo-delta function to distribute the dipole source current and circumvent the problem of singularity at the source point. Moreover, the suggested method used isoparametric finite elements to accommodate the complex subsurface formation. For the large scale linear system derived from the discretization of the Maxwell's equations, several iterative solvers were used and compared to select the optimal one. A quantitative test of accuracy was presented which compared the finite element results with analytic solutions for a dipole source in homogeneous space for different ranges and different wave numbers Ky. to validate the addressed the effects of the distribution range τ of the homogeneous medium. code and check its effectiveness. In addition, we pseudo-delta function on the numerical results in展开更多
Numerical calculation for two integral transforms in 2.5-D transient electromagnetic forward is a difficult and key task, namely, the inverse Fourier transform and the inverse Laplace transform. Some effective algorit...Numerical calculation for two integral transforms in 2.5-D transient electromagnetic forward is a difficult and key task, namely, the inverse Fourier transform and the inverse Laplace transform. Some effective algorithms for them were described. Based on the known algorithms in DC resistivity on wave-number distribution and selection, we proposed a principle on how to choose the least wave-number concerning the central-loop transient electromagnetic method. First, observe the behavior of transformation function curve with regard to wave-number in Fourier domain. In the light of its asymptote, ascertain the coverage scope of wave-number. Compared with analytic solution, the least wave-number in Fourier domain can be derived. Furthermore, the Laplace numerical inversion algorithm which needs only a few Laplace variables in pure real domain was also introduced here. The procedure was applied to forward modeling on transient electromagnetic field of a vertical magnetic dipole over uniform half-space to demonstrate them effectiveness and general applicability.展开更多
We highlighted the flexibility of using unstructured mesh together with the local refinement by a resistivity model with complicated topography. The effect of topography is emphasized. Based on this, we calculated a s...We highlighted the flexibility of using unstructured mesh together with the local refinement by a resistivity model with complicated topography. The effect of topography is emphasized. Based on this, we calculated a specific class of layered models and found that the accuracy is not always satisfactory by utilizing the standard approach. As an improvement, we employed the layered earth as the reference model to calculate the wavenumbers. The comparison demonstrates that the accuracy is considerably improved by using this enhanced approach.展开更多
文摘A finite element method is developed for simulating frequency domain electromagnetic responses due to a dipole source in the 2-D conductive structures. Computing costs are considerably minimized by reducing the full three-dimensional problem to a series of two-dimensional problems. This is accomplished by transforming the problem into y-wave number (Ky) domain using Fourier transform and the y-axis is parallel to the structural strike. In the Ky domain, two coupled partial differential equations for magnetic field Hy and electric field Ey are derived. For a specific value of Ky, the coupled equations are solved by the finite element method with isoparametric elements in the x-z plane. Application of the inverse Fourier transform to the Ky, domain provides the electric and magnetic fields in real space. The equations derived can be applied to general complex two-dimensional structures containing either electric or magnetic dipole source in any direction. In the modeling of the electromagnetic measurement, we adopted a pseudo-delta function to distribute the dipole source current and circumvent the problem of singularity at the source point. Moreover, the suggested method used isoparametric finite elements to accommodate the complex subsurface formation. For the large scale linear system derived from the discretization of the Maxwell's equations, several iterative solvers were used and compared to select the optimal one. A quantitative test of accuracy was presented which compared the finite element results with analytic solutions for a dipole source in homogeneous space for different ranges and different wave numbers Ky. to validate the addressed the effects of the distribution range τ of the homogeneous medium. code and check its effectiveness. In addition, we pseudo-delta function on the numerical results in
基金Project(40344022) supported by the National Natural Science Foundation of China
文摘Numerical calculation for two integral transforms in 2.5-D transient electromagnetic forward is a difficult and key task, namely, the inverse Fourier transform and the inverse Laplace transform. Some effective algorithms for them were described. Based on the known algorithms in DC resistivity on wave-number distribution and selection, we proposed a principle on how to choose the least wave-number concerning the central-loop transient electromagnetic method. First, observe the behavior of transformation function curve with regard to wave-number in Fourier domain. In the light of its asymptote, ascertain the coverage scope of wave-number. Compared with analytic solution, the least wave-number in Fourier domain can be derived. Furthermore, the Laplace numerical inversion algorithm which needs only a few Laplace variables in pure real domain was also introduced here. The procedure was applied to forward modeling on transient electromagnetic field of a vertical magnetic dipole over uniform half-space to demonstrate them effectiveness and general applicability.
基金supported by the National High Technology Research and Development Program of China (863 Program) (No. 2007AA06Z134)the National Natural Science Foundation of China (No. 40874072)
文摘We highlighted the flexibility of using unstructured mesh together with the local refinement by a resistivity model with complicated topography. The effect of topography is emphasized. Based on this, we calculated a specific class of layered models and found that the accuracy is not always satisfactory by utilizing the standard approach. As an improvement, we employed the layered earth as the reference model to calculate the wavenumbers. The comparison demonstrates that the accuracy is considerably improved by using this enhanced approach.