TY - JOUR
T1 - Accelerated 2.5-D inversion of airborne transient electromagnetic data using reduced 3-D meshing
AU - Engebretsen, Kim Wann
AU - Zhang, Bo
AU - Fiandaca, Gianluca
AU - Madsen, Line Meldgaard
AU - Auken, Esben
AU - Christiansen, Anders Vest
N1 - Publisher Copyright:
© 2022 The Author(s) 2022. Published by Oxford University Press on behalf of The Royal Astronomical Society.
PY - 2022/7/1
Y1 - 2022/7/1
N2 - Airborne systems collecting transient electromagnetic data are able to gather large amounts of data over large areas in a very short time. These data are most often interpreted through 1-D inversions, due to the availability of robust, fast and efficient codes. However, in areas where the subsurface contains complex structures or large conductivity contrasts, 1-D inversions may introduce artefacts into the models, which may prevent correct interpretation of the results. In these cases, 2-D or 3-D inversion should be used. Here, we present a 2.5-D inversion code using 3-D forward modelling combined with a 2-D model grid. A 2.5-D inversion is useful where the flight lines are spaced far apart, in which case a 3-D inversion would not add value in relation to the added computational cost and complexity. By exploiting the symmetry of the transmitter and receiver system we are able to perform forward calculations on a reduced 3-D mesh using only half the domain transecting the centre of the transmitter and receiver system. The forward responses and sensitivities from the reduced 3-D mesh are projected onto a structured 2-D model grid following the flight direction. The difference in forward calculations is within 1.4 per cent using the reduced mesh compared to a full 3-D solution. The inversion code is tested on a synthetic example constructed with complex geology and high conductivity contrasts and the results are compared to a 1-D inversion. We find that the 2.5-D inversion recovers both the conductivity values and shape of the true model with a significantly higher accuracy than the 1-D inversion. Finally, the results are supported by a field case using airborne TEM data from the island of Mayotte. The inverted flight line consisted of 418 soundings, and the inversion spent an average of 6750 s per iteration, converging in 16 iterations with a peak memory usage of 97 GB, using 18 logical processors. In general, the total time of the 2-D inversions compared to a full 3-D inversion is reduced by a factor of 2.5 while the memory consumption was reduced by a factor of 2, reflecting the half-mesh approach.
AB - Airborne systems collecting transient electromagnetic data are able to gather large amounts of data over large areas in a very short time. These data are most often interpreted through 1-D inversions, due to the availability of robust, fast and efficient codes. However, in areas where the subsurface contains complex structures or large conductivity contrasts, 1-D inversions may introduce artefacts into the models, which may prevent correct interpretation of the results. In these cases, 2-D or 3-D inversion should be used. Here, we present a 2.5-D inversion code using 3-D forward modelling combined with a 2-D model grid. A 2.5-D inversion is useful where the flight lines are spaced far apart, in which case a 3-D inversion would not add value in relation to the added computational cost and complexity. By exploiting the symmetry of the transmitter and receiver system we are able to perform forward calculations on a reduced 3-D mesh using only half the domain transecting the centre of the transmitter and receiver system. The forward responses and sensitivities from the reduced 3-D mesh are projected onto a structured 2-D model grid following the flight direction. The difference in forward calculations is within 1.4 per cent using the reduced mesh compared to a full 3-D solution. The inversion code is tested on a synthetic example constructed with complex geology and high conductivity contrasts and the results are compared to a 1-D inversion. We find that the 2.5-D inversion recovers both the conductivity values and shape of the true model with a significantly higher accuracy than the 1-D inversion. Finally, the results are supported by a field case using airborne TEM data from the island of Mayotte. The inverted flight line consisted of 418 soundings, and the inversion spent an average of 6750 s per iteration, converging in 16 iterations with a peak memory usage of 97 GB, using 18 logical processors. In general, the total time of the 2-D inversions compared to a full 3-D inversion is reduced by a factor of 2.5 while the memory consumption was reduced by a factor of 2, reflecting the half-mesh approach.
KW - Controlled source electromagnetics (CSEM)
KW - Electrical properties
KW - Hydrogeophysics
KW - Inverse theory
KW - Numerical modelling
UR - http://www.scopus.com/inward/record.url?scp=85128492649&partnerID=8YFLogxK
U2 - 10.1093/gji/ggac077
DO - 10.1093/gji/ggac077
M3 - Article
AN - SCOPUS:85128492649
SN - 0956-540X
VL - 230
SP - 643
EP - 653
JO - Geophysical Journal International
JF - Geophysical Journal International
IS - 1
ER -