2.5D forward modeling and inversion of frequency-domain airborne electromagnetic data

Frequency-domain airborne electromagnetics is a proven geophysical exploration method. Presently, the interpretation is mainly based on resistivity—depth imaging and one-dimensional layered inversion; nevertheless, it is difficult to obtain satisfactory results for two- or three-dimensional complex earth structures using 1D methods. 3D forward modeling and inversion can be used but are hampered by computational limitations because of the large number of data. Thus, we developed a 2.5D frequency-domain airborne electromagnetic forward modeling and inversion algorithm. To eliminate the source singularities in the numerical simulations, we split the fields into primary and secondary fields. The primary fields are calculated using homogeneous or layered models with analytical solutions, and the secondary (scattered) fields are solved by the finite-element method. The linear system of equations is solved by using the large-scale sparse matrix parallel direct solver, which greatly improves the computational efficiency. The inversion algorithm was based on damping least-squares and singular value decomposition and combined the pseudo forward modeling and reciprocity principle to compute the Jacobian matrix. Synthetic and field data were used to test the effectiveness of the proposed method.