The inversion of data from very large three-dimensional electrical resistivity tomography mobile surveys |
| |
| Authors: | MH Loke N Papadopoulos PB Wilkinson D Oikonomou K Simyrdanis DF Rucker |
| |
| Institution: | 1. Geotomo Software Sdn Bhd, Seri, Iskandar, Penang, 11700 Malaysia;2. Laboratory of Geophysical-Remote Sensing & Archaeoenvironment (GeoSat ReSeArch), Institute for Mediterranean Studies, Foundation for Research and Technology Hellas, Rethymno, Crete, 74 100, Greece;3. British Geological Survey, Keyworth, Nottingham, NG12 5GG UK;4. hydroGEOPHYSICS Inc., Tucson, AZ 85745 USA |
| |
| Abstract: | New developments in mobile resistivity meter instrumentation have made it possible to survey large areas with dense data coverage. The mobile system usually has a limited number of electrodes attached to a cable that is pulled along behind an operator so that a large area can be covered within a short time. Such surveys can produce three-dimensional datasets with hundreds of thousands of electrodes positions and data points. Similarly, the inverse model used to interpret the data can have several hundred thousand cells. It is impractical to model such large datasets within a reasonable time on microcomputers used by many small companies employing standard inversion techniques. We describe a model segmentation technique that subdivides the finite-element mesh used to calculate the apparent resistivity and Jacobian matrix values into a number of smaller meshes. A fast technique that optimizes the calculation of the Jacobian matrix values for multi-channel systems was also developed. A one-dimensional wavelet transform method was then used to compress the storage of the Jacobian matrix, in turn reducing the computer time and memory required to solve the least-squares optimization equation to determine the inverse model resistivity values. The new techniques reduce the calculation time and memory required by more than 80% while producing models that differ by less than 1% from that obtained using the standard inversion technique with a single mesh. We present results using a synthetic model and a field dataset that illustrates the effectiveness of the proposed techniques. |
| |
| Keywords: | Computing aspects 3D Inversion Resistivity Tomography |
|
|
