The point-spread function measure of resolution for the 3-D electrical resistivity experiment

被引:42
作者
Oldenborger, Greg A. [1 ]
Routh, Partha S. [2 ]
机构
[1] Univ British Columbia, Dept Earth & Ocean Sci, Vancouver, BC V6T 1Z4, Canada
[2] Boise State Univ, Dept Geosci, Boise, ID 83725 USA
基金
加拿大自然科学与工程研究理事会;
关键词
Inverse theory; Tomography; Electrical properties; Hydrogeophysics; Hydrology; GROSS EARTH DATA; INVERSE PROBLEMS; TOMOGRAPHY; APPRAISAL;
D O I
10.1111/j.1365-246X.2008.04003.x
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
The solution appraisal component of the inverse problem involves investigation of the relationship between our estimated model and the actual model. However, full appraisal is difficult for large 3-D problems such as electrical resistivity tomography (ERT). We tackle the appraisal problem for 3-D ERT via the point-spread functions (PSFs) of the linearized resolution matrix. The PSFs represent the impulse response of the inverse solution and quantify our parameter-specific resolving capability. We implement an iterative least-squares solution of the PSF for the ERT experiment, using on-the-fly calculation of the sensitivity via an adjoint integral equation with stored Green's functions and subgrid reduction. For a synthetic example, analysis of individual PSFs demonstrates the truly 3-D character of the resolution. The PSFs for the ERT experiment are Gaussian-like in shape, with directional asymmetry and significant off-diagonal features. Computation of attributes representative of the blurring and localization of the PSF reveal significant spatial dependence of the resolution with some correlation to the electrode infrastructure. Application to a time-lapse ground-water monitoring experiment demonstrates the utility of the PSF for assessing feature discrimination, predicting artefacts and identifying model dependence of resolution. For a judicious selection of model parameters, we analyse the PSFs and their attributes to quantify the case-specific localized resolving capability and its variability over regions of interest. We observe approximate interborehole resolving capability of less than 1-1.5 m in the vertical direction and less than 1-2.5 m in the horizontal direction. Resolving capability deteriorates significantly outside the electrode infrastructure.
引用
收藏
页码:405 / 414
页数:10
相关论文
共 28 条
[1]   Image appraisal for 2-D and 3-D electromagnetic inversion [J].
Alumbaugh, DL ;
Newman, GA .
GEOPHYSICS, 2000, 65 (05) :1455-1467
[2]   UNIQUENESS IN INVERSION OF INACCURATE GROSS EARTH DATA [J].
BACKUS, G ;
GILBERT, F .
PHILOSOPHICAL TRANSACTIONS OF THE ROYAL SOCIETY OF LONDON SERIES A-MATHEMATICAL AND PHYSICAL SCIENCES, 1970, 266 (1173) :123-&
[3]   RESOLVING POWER OF GROSS EARTH DATA [J].
BACKUS, G ;
GILBERT, F .
GEOPHYSICAL JOURNAL OF THE ROYAL ASTRONOMICAL SOCIETY, 1968, 16 (02) :169-&
[4]   Hierarchical geostatistics and multifacies systems: Boise Hydrogeophysical Research Site, Boise, Idaho [J].
Barrash, W ;
Clemo, T .
WATER RESOURCES RESEARCH, 2002, 38 (10)
[5]   Finite element three-dimensional direct current resistivity modelling: accuracy and efficiency considerations [J].
Bing, Z ;
Greenhalgh, SA .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2001, 145 (03) :679-688
[6]   RESISTIVITY MODELING FOR ARBITRARILY SHAPED 3-DIMENSIONAL STRUCTURES [J].
DEY, A ;
MORRISON, HF .
GEOPHYSICS, 1979, 44 (04) :753-780
[7]   Resolution, stability and efficiency of resistivity tomography estimated from a generalized inverse approach [J].
Friedel, S .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2003, 153 (02) :305-316
[8]  
Golub G. H., 2012, Matrix computations, V4th
[9]   3-D inversion of induced polarization data [J].
Li, YG ;
Oldenburg, DW .
GEOPHYSICS, 2000, 65 (06) :1931-1945
[10]  
Menke W., 1989, GEOPHYS DATA ANAL DI