Location of Repository

Parallel computation of optimized arrays for 2-D electrical\ud imaging surveys

By M.H. Loke, P.B. Wilkinson and J.E. Chambers


Modern automatic multi-electrode survey instruments have made it possible to use non-traditional arrays to maximize the subsurface resolution from electrical imaging surveys. Previous studies have shown that one of the best methods for generating optimized arrays is to select the set of array configurations that maximizes the model resolution for a homogeneous earth model. The Sherman–Morrison Rank-1 update is used to calculate the change in the model resolution when a new array is added to a selected set of array configurations. This method had the disadvantage that it required several hours of computer time even for short 2-D survey lines. The algorithm was modified to calculate the change in the model resolution rather than the entire resolution matrix. This reduces the computer time and memory required as well as the computational round-off errors. The matrix–vector multiplications for a single add-on array were replaced with matrix–matrix multiplications for 28 add-on arrays to further reduce the computer time. The temporary variables were stored in the double-precision Single Instruction Multiple Data (SIMD) registers within the CPU to minimize computer memory access. A further reduction in the computer time is achieved by using the computer graphics card Graphics Processor Unit (GPU) as a highly parallel mathematical coprocessor. This makes it possible to carry out the calculations for 512 add-on arrays in parallel using the GPU. The changes reduce the computer time by more than two orders of magnitude. The algorithm used to generate an optimized data set adds a specified number of new array configurations after each iteration to the existing set. The resolution of the optimized data set can be increased by adding a smaller number of new array configurations after each iteration. Although this increases the computer time required to generate an optimized data set with the same number of data points, the new fast numerical routines has made this practical on commonly available microcomputers.\ud \u

Topics: Computer Science, Physics, Earth Sciences, Mathematics
Publisher: Wiley
Year: 2010
DOI identifier: 10.1111/j.1365-246X.2010.04796.x
OAI identifier: oai:nora.nerc.ac.uk:12869

Suggested articles



  1. (2005). 2-D resistivity surveying for hydrocarbons – A primer, CSEG Recorder,
  2. (1996). 2D resistivity surveying for environmental and engineering applications,
  3. (2005). 32/64-bit 80x86 assembly language architecture,
  4. (2006). 3D resistivity in the Athabasca basin with the pole-pole array,
  5. (2002). A 3D resistivity investigation of a contaminated site at Lernacken in Sweden,
  6. (2004). A comparison of automatic techniques for estimating the regularization parameter in non-linear inverse problems,
  7. (2002). A comparison of Gauss-Newton and quasi-Newton methods in resistivity imaging inversion,
  8. (2003). A comparison of smooth and blocky inversion methods in 2D electrical imaging surveys,
  9. (2004). A numerical comparison of 2D resistivity imaging with 10 electrode arrays,
  10. (2006). A survey of current trends in near-surface electrical and electromagnetic methods,
  11. (1956). A tri-potential method of resistivity prospecting,
  12. (2009). Accelerating a three-dimensional finite-difference wave propagation code using GPU graphics cards.
  13. (2010). An introduction to this special section: High performance computing.
  14. (1994). Applied geophysical inversion,
  15. (2000). Aquifer characterization in the Blue Ridge physiographic province using resistivity profiling and borehole geophysics : Geologic analysis,
  16. (2007). Capacitive Resistivity Imaging with Towed Arrays,
  17. (2006). Comparison of the spatial resolution of standard and optimised Electrical Resistivity Tomography arrays, doi
  18. (2003). Electrical resistivity imaging for detecting soil cracking at the centimetric Scale.
  19. (2006). Electrical resistivity tomography applied to geologic, hydrogeologic, and engineering investigations at a former waste-disposal site,
  20. (1996). Examples of resistivity imaging using ME-100 resistivity field acquisition system,
  21. (2004). Experimental design: Electrical resistivity data sets that provide optimum subsurface information, doi
  22. (2010). Fast computation of optimized electrode arrays for 2D resistivity surveys, Computers and Geosciences,
  23. (1989). Geophysical data analysis:
  24. (2006). Improved strategies for the automatic selection of optimized sets of electrical resistivity tomography measurement configurations, doi
  25. (1972). Interpretation of inaccurate, insufficient and inconsistent data,
  26. (1995). Least-squares deconvolution of apparent resistivity pseudosections,
  27. (2010). Leveraging graphics processing units (GPUs) for real-time seismic interpretation, The Leading Edge,
  28. (2010). Massively parallel forward modeling of scalar and tensor gravimetry data,
  29. (1998). Numerical linear algebra for high-performance computers, Society for Industrial and Applied Mathematics,
  30. (1992). Numerical Recipes in C, 2nd edn,
  31. (1990). Occam's inversion to generate smooth, twodimensional models form magnetotelluric data,
  32. (2001). Parallel Programming in OpenMP,
  33. (2003). Properties and effects of measurement errors on 2D resistivity imaging surveying, doi
  34. (2007). Resolution analysis of geophysical images: Comparison between point spread function and region of data influence measures,
  35. (2003). Resolution, stability and efficiency of resistivity tomography estimated from a generalized inverse approach,
  36. (2010). Selecting the right hardware for reverse time migration, The Leading Edge,
  37. (2009). The point-spread function measure of resolution for the 3-D electrical resistivity experiment.
  38. (2002). The Software Optimization Cookbook : High-performance recipes for the Intel architecture, Intel Press,
  39. (2010). Trends for high-performance scientific computing,
  40. (2008). Using OpenMP,

To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.