An efficient regularized inversion approach for self-potential data interpretation of ore exploration using a mix of logarithmic and non-logarithmic model parameters |
| |
| Affiliation: | 1. Departamento de Matemáticas, Universidad de Oviedo, Oviedo, Spain;2. ETSI en Topografía, Geodesia y Cartografía, Universidad Politécnica de Madrid, Madrid, Spain;1. Physics Department, Institut Teknologi Sepuluh Nopember (ITS), Jl. Arif Rahman Hakim, Surabaya 60111, Indonesia;2. Sidoarjo Mud Flow Mitigation Agency (BPLS), Indonesia;3. Geophysical Engineering Department, Institut Teknologi Sepuluh Nopember (ITS), Jl. Arif Rahman Hakim, Surabaya 60111, Indonesia;4. University of Lisbon — DEGGE-IDL, Campo Grande Ed. C8, 1749-016, Lisboa, Portugal;1. School of Information and Communication Technology, Griffith University, Nathan, Brisbane, QLD 4111, Australia;2. Griffith College, Mt Gravatt, Brisbane, QLD 4122, Australia;1. ETSI en Topografía, Geodesia y Cartografía, Universidad Politécnica de Madrid, Madrid, Spain;2. Departamento de Matemáticas, Universidad de Oviedo, Oviedo, Spain;3. Laboratoire GET (Université de Toulouse, CNRS, IRD, CNES), Bureau Gravimétrique International (BGI), Toulouse, France;4. Université de Toulouse, Mines Albi, CNRS, Centre RAPSODEE, Albi cedex, France;1. Department of Physics, Institut Teknologi Sepuluh Nopember, Kampus ITS Sukolilo Surabaya-60111, Indonesia;2. Department of Geophysical Engineering, Institut Teknologi Sepuluh Nopember, Kampus ITS Sukolilo Surabaya-60111, Indonesia |
| |
| Abstract: | A very fast and efficient approach to self-potential (SP) data inversion for ore exploration has been developed. This approach is based on Tikhonov regularization and the conjugate gradient method, and simultaneously inverts for the depth (z), electric dipole moment (k), and angle of polarization (θ) of a buried anomalous body from SP data measured along a profile. This inversion algorithm works iteratively, and solves for z and k in the logarithmic-space (log(z) and log(k)), and solves for θ in the linear-space (non-logarithmic). It is found that the original inversion formulation that uses the model parameters themselves (z, k and θ) is unstable and divergent. It is also found that the inversion formulation that uses the logarithm of the model parameters (log(z), log(k) and log(θ)) is unstable and divergent. Rather, the new inversion scheme that is based on the aforementioned mixed log-linear combination of the model parameters (log(z), log(k), and θ) overcomes and eliminates the mentioned instability and divergence problems. The sensitivity analysis and numerical experiments investigated have indicated that the new approach has a far better and far more optimized minimization search direction. This proposed technique fits the observed data by some geometrically simple body in the restricted class of vertical cylinder, horizontal cylinder, and sphere models. The applicability of the algorithm has been demonstrated on various reliable synthetic data sets with and without noise. The algorithm has been carefully and successfully applied to six real data examples, with ore bodies buried in different complex geologic settings and at various depths in the subsurface. The method is shown to be highly applicable for mineral exploration, and is of particular value in cases where the SP observed data is due to ore body embedded in the subsurface. On average, it took about 40 s of computation (not CPU) time on a 1 GHz PC. |
| |
| Keywords: | |
| 本文献已被 ScienceDirect 等数据库收录！ |
|
