Rapid inversion of data from 2D resistivity surveys with electrode displacements
Tóm tắt
Resistivity monitoring surveys are used to detect temporal changes in the subsurface using repeated measurements over the same site. The positions of the electrodes are typically measured at the start of the survey program and possibly at occasional later times. In areas with unstable ground, such as landslide‐prone slopes, the positions of the electrodes can be displaced by ground movements. If this occurs at times when the positions of the electrodes are not directly measured, they have to be estimated. This can be done by interpolation or, as in recent developments, from the resistivity data using new inverse methods. The smoothness‐constrained least squares optimisation method can be modified to include the electrode positions as additional unknown parameters. The Jacobian matrices with the sensitivity of the apparent resistivity measurements to changes in the electrode positions are then required by the optimisation method. In this paper, a fast adjoint‐equation method is used to calculate the Jacobian matrices required by the least squares method to reduce the calculation time. In areas with large near‐surface resistivity contrasts, the inversion routine sometimes cannot accurately distinguish between electrode displacements and subsurface resistivity variations. To overcome this problem, the model for the initial time‐lapse dataset (with accurately known electrode positions) is used as the starting model for the inversion of the later‐time dataset. This greatly improves the accuracy of the estimated electrode positions compared to the use of a homogeneous half‐space starting model. In areas where the movement of the electrodes is expected to occur in a fixed direction, the method of transformations can be used to include this information as an additional constraint in the optimisation routine.
Từ khóa
Tài liệu tham khảo
Daniels R.W., 1978, An Introduction to Numerical Methods and Optimization Techniques
Golub G., 1996, Matrix Computations
Jennings A., 1992, Matrix Computation
LokeM.H.2000.Topographic modelling in resistivity imaging inversion.62nd EAGE Conference & Technical Exhibition Glasgow Scotland Extended Abstracts D‐2.
LokeM.H.andDahlinT.2010.Methods to reduce banding effects in 3‐D resistivity inversion.16th European Meeting of Environmental and Engineering Geophysics Zurich Switzerland Expanded Abstracts A16.
LokeM.H. WilkinsonP.B.andChambersJ.E.2016.3‐D resistivity inversion with electrodes displacements.ASEG‐PESA‐AIG 25th Geophysical Conference and Exhibition August 21–24 2016 Adelaide Australia Extended Abstracts 809–813.
Press W.H., 1992, Numerical Recipes in C
Schwarz H.R., 1973, Numerical Analysis of Symmetric Matrices
Silvester P.P., 1990, Finite Elements for Electrical Engineers