Optimized arrays for 2-D resistivity survey lines with a large number of electrodes

• Use of symmetry in model resolution values to reduce calculate time by half.
• Improved single-precision computational algorithm.
• Matrix decomposition technique that is 30 times faster than direct updates.
• Symmetrical arrays subsets that reduce calculation time by more than 10 times.
• Advantages of optimized arrays shown by synthetic model and field survey.

Previous studies show that optimized arrays generated using the ‘Compare R’ method have significantly better resolution than conventional arrays. This method determines the optimum set of arrays by selecting those that give the maximum model resolution. The number of possible arrays (the comprehensive data set) increases with the fourth power of the number of electrodes. The optimization method faces practical limitations for 2-D survey lines with more than 60 electrodes where the number of possible arrays exceeds a million. Several techniques are proposed to reduce the calculation time for such survey lines. A single-precision version of the ‘Compare R’ algorithm using a new ranking function reduces the calculation time by two to eight times while providing results similar to the double-precision version. Recent improvements in computer GPU technology can reduce the calculation time by about seven times. The calculation time is reduced by half by using the fact that arrays that are symmetrical about the center of the line produce identical changes in the model resolution values. It is further reduced by more than thirty times by calculating the Sherman–Morrison update for all the possible two-electrode combinations, which are then used to calculate the model resolution values for the four-electrode arrays. The calculation time is reduced by more then ten times by using a subset of the comprehensive data set consisting of only symmetrical arrays. Tests with a synthetic model and field data set show that optimized arrays derived from this subset produce inversion models with differences of less than 10% from those derived using the full comprehensive data set. The optimized data sets produced models that are more accurate than the Wenner–Schlumberger array data sets in all the tests.