AarhusInv - inversion code


AarhusInv is a high performance modeling and inversion code supporting a variety of geophysical data types, configurations, and source-receiver types. It provides efficient and high precision modelling and inversion of any airborne EM configuration and ground based TEM, GCM, DC, MRS and time domain IP configurations.


Key feautures

  • Layer, smooth and sharp model description

  • Model description in a regular 3D grid with scatted data records

  • Lateral and vertical constraining of any model parameters

  • Prior constraints on any model parameters

  • Joint or constraints inversion of different data types and methods

  • Depth of investigation estimates (DOI)



In general, AarhusInv uses a local 1D model description. The models are laterally and spatially constrained forming pseudo 2D and 3D model spaces. The inversion code is designed to handle data from very big airborne EM surveys and support multi CPU. Besides the resistivity models, the inversion code also calculates a depth of investigation value (DOI) and a data residual for each resistivity model. The inversion setups are: 


  • LCI-setup: The models are laterally constrained along the flight lines forming a 2D model space.

  • SCI-setup: The models are laterally constrained along the flight lines and across the flight lines, resulting in a 3D model space.

  • For DC and IP, a full 2D Cole-Cole or Constant Phase Angle model parameterization is used


The inversion is done iteratively and it supports:


  • Smooth models: The resistivity model is discretized using several layers (~10-30) with fixed 
    layer boundaries. The regularization penalizes vertical changes in resistivity, resulting in a vertical smooth resistivity model.

  • Sharp models: The resistivity model is discretized using several layers (~10-30) with fixed
    layer boundaries. The regularization penalizes the number of vertical resistivity
    transitions of a certain size, resulting in resistivity models with relatively sharp
    vertical resistivity transitions.

  • Layered models: The resistivity model is characterized by a small number of layers (~4-6). Both layer thickness and resistivity are model parameters.

  • Prior constraints on any model parameter. The prior constraints can be initialized from grids and directly from the GIS map or can be specified at borehole locations with decreasing strength moving away from the borehole locations.

  • Depth of investigation estimates (DOI)


The high performance is obtained by an efficient code parallelization, iterative sparse matrix solvers and sparse matrix data stored for efficient memory handling. AarhusInv is therefore capable of handling vary large airborne EM surveys in a single spatial constraint model setup (SCI) and utilizing multi core CPU’s.

The AarhusInv inversion code is available through the Aarhus Workbench and SPIA program; it cannot be acquired as a standalone program.


© 2019 by Aarhus GeoSoftware

Skolegade 21, 1

8000 Aarhus C, Denmark

VAT: DK37108472

Bank: Jutlander Bank. Reg. no.:8139

Account no.:2072741946 - IBAN no.: DK8281392072741946

SWIFT code: JUTBDK21 - Intermediary banK: KRONDK22

  • LinkedIn Clean