Optimal design of simultaneous source encoding

A broad range of parameter estimation problems involve the collection of an excessively large number of observations N. Typically, each such observation involves excitation of the domain through injection of energy at some predefined sites and recording of the response of the domain at another set of locations. It has been observed that similar results can often be obtained by considering a far smaller number K of multiple linear superpositions of experiments with . This allows the construction of the solution to the inverse problem in time instead of . Given these considerations it should not be necessary to perform all the N experiments but only a much smaller number of K experiments with simultaneous sources in superpositions with certain weights. Devising such procedure would results in a drastic reduction in acquisition time. The question we attempt to rigorously investigate in this work is: what are the optimal weights? We formulate the problem as an optimal experimental design problem and show that by leveraging techniques from this field an answer is readily available. Designing optimal experiments requires some statistical framework and therefore the statistical framework that one chooses to work with plays a major role in the selection of the weights.

[1]  Hansruedi Maurer,et al.  Recent advances in optimized geophysical survey design , 2010 .

[2]  L. Tenorio,et al.  Data analysis tools for uncertainty quantification of inverse problems , 2011 .

[3]  E. Miller,et al.  A shape reconstruction method for electromagnetic tomography using adjoint fields and level sets , 2000 .

[4]  Darrell Coles,et al.  A method of fast, sequential experimental design for linearized geophysical inverse problems , 2009 .

[5]  F. Herrmann,et al.  Compressive simultaneous full-waveform simulation , 2009 .

[6]  E. Haber,et al.  Numerical methods for the design of large-scale nonlinear discrete ill-posed inverse problems , 2010 .

[7]  M. Hutchinson.  A stochastic estimator of the trace of the influence matrix for laplacian smoothing splines , 1989 .

[8]  E. Haber,et al.  Optimal Experimental Design for the Large‐Scale Nonlinear Ill‐Posed Problem of Impedance Imaging , 2010 .

[9]  A. R. Borges,et al.  A quantitative algorithm for parameter estimation in magnetic induction tomography , 2004 .

[10]  S. Arridge.  Optical tomography in medical imaging , 1999 .

[11]  Tom Lahmer,et al.  Optimal experimental design for nonlinear ill-posed problems applied to gravity dams , 2011 .

[12]  Max Deffenbaugh,et al.  Efficient seismic forward modeling using simultaneous random sources and sparsity , 2010 .

[13]  Per Christian Hansen,et al.  Rank-Deficient and Discrete Ill-Posed Problems , 1996 .

[14]  Eldad Haber,et al.  An Effective Method for Parameter Estimation with PDE Constraints with Multiple Right-Hand Sides , 2012, SIAM J. Optim..

[15]  Curtis C. Ober,et al.  Faster shot-record depth migrations using phase encoding , 1998 .

[16]  D K Smith,et al.  Numerical Optimization , 2001, J. Oper. Res. Soc..

[17]  Craig J. Beasley,et al.  A new look at marine simultaneous sources , 2008 .

[18]  Max Deffenbaugh,et al.  Efficient Seismic Forward Modeling and Acquisition using Simultaneous Random Sources and Sparsity , .

[19]  William W. Symes,et al.  Source Synthesis For Waveform Inversion , 2010 .

[20]  Louis A. Romero,et al.  Phase encoding of shot records in prestack migration , 2000 .

[21]  Kurt M. Strack,et al.  Society of Exploration Geophysicists , 2007 .

[22]  E. Haber,et al.  Numerical methods for experimental design of large-scale linear ill-posed inverse problems , 2008 .

[23]  Per Christian Hansen,et al.  Truncated Singular Value Decomposition Solutions to Discrete Ill-Posed Problems with Ill-Determined Numerical Rank , 1990, SIAM J. Sci. Comput..

[24]  Aria Abubakar,et al.  Full‐waveform seismic inversion using the source‐receiver compression approach , 2010 .