Piecewise 1D laterally constrained inversion of resistivity data

In a sedimentary environment, layered models are often capable of representing the actual geology more accurately than smooth minimum structure models. Furthermore, interval thicknesses and resistivities are often the parameters to which non-geophysicist experts can relate and base decisions on when using them in waste site remediation, groundwater modelling and physical planning. We present a laterally constrained inversion scheme for continuous resistivity data based on a layered earth model (1D). All 1D data sets and models are inverted as one system, producing layered sections with lateral smooth transitions. The models are regularized through laterally equal constraints that tie interface depths and resistivities of adjacent layers. Prior information, e.g. originating from electric logs, migrates through the lateral constraints to the adjacent models, making resolution of equivalences possible to some extent. Information from areas with well-resolved parameters will migrate through the constraints in a similar way to help resolve the poorly constrained parameters. The estimated model is complemented by a full sensitivity analysis of the model parameters, supporting quantitative evaluation of the inversion result. Examples from synthetic 2D models show that the model recognition of a sublayered 2D wedge model is improved using the laterally constrained inversion approach when compared with a section of combined 1D models and when compared with a 2D minimum structure inversion. Case histories with data from two different continuous DC systems support the conclusions drawn from the synthetic example.