The correct interpretation of time-domain or transient electromagnetic (TEM) data relies on effective inversion techniques. Herein, we have developed an effective 3D-TEM inversion method based on the vector finite-element method. The domain was discretized with unstructured tetrahedra to simulate complex geologic structures. We develop a second-order backward-Euler scheme, rigorously with nonuniform time-step sizes, to approximate the time derivative of the electric field. The waveform was directly simulated and considered in the forward modeling and inversion. In a 3D inversion, the most time-consuming element of the workflow is solving a large system of linear equations. To alleviate this problem, we adopt a parallel direct solver for solving a large system of equations involved in forward modeling and sensitivity calculations. We use the Gauss-Newton approach to solve the inverse problem and formulate the least-squares problem to obtain model updates. We avoided the explicit sensitivity calculation by solving the matrix-vector multiplication between the sensitivity matrix (or its transpose) and vector. The developed inversion algorithm is validated using several realistic synthetic models. Finally, we apply the method to practical TEM data with complex topography for coal mine goaf detection. We thus obtain a detailed 3D resistivity model that revealed the layering structure and some localized anomalies related to the coal mine goaf.
Read full abstract