Temporal Orthogonal Projection Inversion for EMI Sensing of UXO

被引:14
作者
Song, Lin-Ping [1 ]
Oldenburg, Douglas W. [1 ]
Pasion, L. R. [2 ]
Billings, S. D. [2 ]
Beran, L. [2 ]
机构
[1] Univ British Columbia, Dept Earth Ocean & Atmospher Sci, Vancouver, BC V6T 1Z4, Canada
[2] Black Tusk Geophys Inc, Vancouver, BC V6J 4S5, Canada
来源
IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING | 2015年 / 53卷 / 02期
关键词
Electromagnetic induction (EMI); magnetic dipole polarization; nonlinear inversion; orthogonal projection; subspace; unexploded ordnance (UXO); UNEXPLODED ORDNANCE; DISCRIMINATION; CLASSIFICATION; IDENTIFICATION; POLARIZABILITIES; OBJECTS; MODEL;
D O I
10.1109/TGRS.2014.2332992
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
We present a new approach for inverting time-domain electromagnetic data to recover the location and magnetic dipole polarizations of a limited number of buried objects. We form the multichannel electromagnetic induction (EMI) sensor data as a spatial-temporal response matrix (STRM). The rows of the STRM correspond to measurements sampled at different time channels from one sensor and the columns correspond to measurements sampled at the same time channel from different sensors. The singular value decomposition of the STRM produces the left and right singular vectors that are related to the sensor and the temporal spaces, respectively. If the effective rank of the STRM is r, then the first r singular vectors span signal subspaces (SS), and the remaining singular vectors span the noise subspaces. The original data are projected onto the SS, and the temporal orthogonal projection inversion (TOPI) uses these data in a nonlinear inverse problem to solve for source locations of the objects. The polarizations of the targets are then obtained by solving a linear optimization problem in the original data domain. We present theoretical and numerical analyses to investigate the singular value system of the STRM and the sensitivity of the TOPI to the size of an SS. Only a few subspace vectors are required to generate locations of the objects. The results are insensitive to the exact choice of rank, and this differs from usual methods that involve selecting the number of time channels to be used in the inversion and carefully estimating associated uncertainties. The proposed approach is evaluated using the synthetic and real multistatic EMI data.
引用
收藏
页码:1061 / 1072
页数:12
相关论文
共 50 条