PyMERRY: A Python solution for an improved interpretation of electrical resistivity tomography images

Abstract

Electrical resistivity tomography (ERT) is a widely used geophysical method for studying geologic hazards, civil engineering, and environmental remediation. It provides information about the subsurface’s resistivity distribution by analyzing electrical data collected at the surface or in boreholes. However, interpreting ERT images can be complex due to ambiguities in their resolution. To address this issue, we develop a postprocessing method called Python iMprovement of Electrical Resistivity tomography ReliabilitY (PyMERRY) to improve the reliability of ERT images. The PyMERRY code can be applied to any 2D resistivity model obtained from ERT inversion software. It computes a coverage mask that defines the domain well constrained by the data and the inversion process. It also evaluates the resistivity uncertainties in the ERT models. In addition to the sensitivity approaches, PyMERRY provides low- and high-resistivity values for all covered cells. Synthetic tests indicate that the approach is efficient and highlight the importance of resistivity contrasts, mesh selection, electrode spacing, and profile length in the reliability of ERT images. Compared with previous studies, using PyMERRY in south-central Bhutan allows a more accurate interpretation of ERT images. It confirms a high-resistivity contrast across the topographic frontal thrust and reveals the existence of small-scale variations.