3-D Inversion of Z -Axis Tipper Electromagnetic Data Using Finite-Element Method With Unstructured Tetrahedral Grids

Xiaoyue Cao,Xin Huang,Changchun Yin,Liangjun Yan,Yongqi Han
DOI: https://doi.org/10.1109/tgrs.2021.3077510
IF: 8.2
2022-01-01
IEEE Transactions on Geoscience and Remote Sensing
Abstract:<span class="mjpage"><svg xmlns:xlink="http://www.w3.org/1999/xlink" width="1.68ex" height="2.176ex" style="vertical-align: -0.338ex;" viewBox="0 -791.3 723.5 936.9" role="img" focusable="false" xmlns="http://www.w3.org/2000/svg"><g stroke="currentColor" fill="currentColor" stroke-width="0" transform="matrix(1 0 0 -1 0 0)"> <use xlink:href="#MJMATHI-5A" x="0" y="0"></use></g></svg></span> -axis tipper electromagnetic (ZTEM) technique is an airborne electromagnetic method that detects anomalies in the deep earth that are induced by natural sources. Conventional ZTEM forward modeling is generally conducted using structured grids that have limited accuracy and cannot be used to invert complex underground structures and topography. However, unstructured grids can accurately model complex underground structures with fewer cells. Therefore, to effectively model and recover the complex underground structures, we developed a 3-D framework with unstructured tetrahedral grids for ZTEM forward modeling and inversions. The forward problem was formulated using the finite-element method with unstructured tetrahedral grids. To solve the inverse problem, we used a limited-memory quasi-Newton algorithm (L-BFGS) with a parallel direct solver for optimization in order to avoid the explicit calculation of the Hessian matrix and save the memory and computational time. To validate our forward and inversion algorithms, we conducted numerical experiments on three synthetic models and inverted a survey dataset acquired for a tunnel project in Tibet, China. The experimental results demonstrate the effectiveness of our unstructured finite-element and L-BFGS method for modeling and inverting ZTEM data.<svg xmlns="http://www.w3.org/2000/svg" style="display: none;"><defs id="MathJax_SVG_glyphs"><path stroke-width="1" id="MJMATHI-5A" d="M58 8Q58 23 64 35Q64 36 329 334T596 635L586 637Q575 637 512 637H500H476Q442 637 420 635T365 624T311 598T266 548T228 469Q227 466 226 463T224 458T223 453T222 450L221 448Q218 443 202 443Q185 443 182 453L214 561Q228 606 241 651Q249 679 253 681Q256 683 487 683H718Q723 678 723 675Q723 673 717 649Q189 54 188 52L185 49H274Q369 50 377 51Q452 60 500 100T579 247Q587 272 590 277T603 282H607Q628 282 628 271Q547 5 541 2Q538 0 300 0H124Q58 0 58 8Z"></path></defs></svg>
imaging science & photographic technology,remote sensing,engineering, electrical & electronic,geochemistry & geophysics
What problem does this paper attempt to address?