Search

Abstract: Four different methods to automatically select an optimal set of array configurations that gives the maximum subsurface resolution with a limited number of measurements for 2D electrical imaging surveys were tested. The first (CR) method directly calculates the change in the model resolution for each new array added to the base data set, and uses this to select array configurations that gave the maximum model resolution. However this method is the slowest. The algorithm used by the CR method for calculating rank-one updates was optimized to reduce computational time by a factor of eighty. The sequence of calculations was modified to reduce the traffic between the computer main memory and the CPU registers. Further code optimizations were made to take advantage of the parallel processing capabilities of modern CPUs. The second (ETH) and third (BGS) methods use approximations based on the sensitivity values to estimate the change in the model resolution matrix. The ETH and BGS methods, respectively, use the first and second power of the sensitivity values to calculate approximations of the model resolution. Both methods are about an order of magnitude faster than the CR method. The results obtained by the BGS method are significantly better than the ETH method, and it approaches that of the CR method. The fourth method (BGS–CR) uses a combination of the BGS and CR methods. It produces results that are almost identical to the CR method but is several times faster. The different methods were tested using data from synthetic models and field surveys. The models obtained from the inversion of the data sets generated by the four different methods confirm that the models generated by the CR method have the best resolution, followed by the BGS–CR, BGS and ETH methods. [ABSTRACT FROM AUTHOR]

The smoothness-constrained least-squares method is widely used for two-dimensional (2D) and three-dimensional (3D) inversion of apparent resistivity data sets. The Gauss–Newton method that recalculates the Jacobian matrix of partial derivatives for all iterations is commonly used to solve the least-squares equation. The quasi-Newton method has also been used to reduce the computer time. In this method, the Jacobian matrix for a homogeneous earth model is used for the first iteration, and the Jacobian matrices for subsequent iterations are estimated by an updating technique. Since the Gauss–Newton method uses the exact partial derivatives, it should require fewer iterations to converge. However, for many data sets, the quasi-Newton method can be significantly faster than the Gauss–Newton method. The effectiveness of a third method that is a combination of the Gauss–Newton and quasi-Newton methods is also examined. In this combined inversion method, the partial derivatives are directly recalculated for the first two or three iterations, and then estimated by a quasi-Newton updating technique for the later iterations.The three different inversion methods are tested with a number of synthetic and field data sets. In areas with moderate (less than 10:1) subsurface resistivity contrasts, the inversion models obtained by the three methods are similar. In areas with large resistivity contrasts, the Gauss–Newton method gives significantly more accurate results than the quasi-Newton method. However, even for large resistivity contrasts, the differences in the models obtained by the Gauss–Newton method and the combined inversion method are small. As the combined inversion method is faster than the Gauss–Newton method, it represents a satisfactory compromise between speed and accuracy for many data sets. [Copyright &y& Elsevier]

The optimized ‘Compare R’ (CR) method is considered one of the best algorithms to automatically generate a set of electrode arrays that maximizes the model resolution for electrical resistivity imaging. In its original version, the calculation time to generate a set of arrays, even for a short resistivity line, required several hours. Recently, researchers spent considerable effort to reduce execution time. In addition, other modifications have been made to improve the practical aspects, such as minimizing the electrode polarization effects and incorporating errors in the data. In line with this approach, the effectiveness of the optimized arrays that are generated based on a modified version of ‘Compare-R' optimization algorithm is demonstrated through a comprehensive and quantitative comparison with the conventional Dipole-Dipole (DD) and Wenner-Schlumberger (WS) arrays. Four different synthetic models were generated and noisy data cases were considered. The results show that the optimized array has significantly better resolution and target detail than conventional arrays. This is clear in enhancing the detectability accuracy in the deepest part of the inversion models. Of the conventional arrays, the DD array provides inversion images mostly closely related to the performance of the optimized array models. [ABSTRACT FROM AUTHOR]