Applying Compactness Constraints to Differential Traveltime Tomography

Tomographic imaging problems are typically ill-posed and often require the use of regularization techniques to guarantee a stable solution. Minimization of a weighted norm of model length is one commonly used secondary constraint. Tikhonov methods exploit low-order differential operators to select for solutions that are small, flat, or smooth in one or more dimensions. This class of regularizing functionals may not always be appropriate, particularly in cases where the anomaly being imaged is generated by a nonsmooth spatial process. Time-lapse imaging of flow-induced velocity anomalies is one such case; flow features are often characterized by spatial compactness or connectivity. By performing inversions on differenced arrival time data, the properties of the time-lapse feature can be directly constrained. We develop a differential traveltime tomography algorithm whichselects for compact solutions, i.e., models with a minimum area of support, through application of model-space iteratively reweighted leas...

[1]  K. Kubik,et al.  Compact gravity inversion , 1983 .

[2]  C. Vogel,et al.  Analysis of bounded variation penalty methods for ill-posed problems , 1994 .

[3]  T. Davis,et al.  Multicomponent seismic characterization and monitoring of the CO2 flood at Weyburn Field, Saskatchewan , 2003 .

[4]  P. N. Sen,et al.  A self-similar model for sedimentary rocks with application to the dielectric constant of fused glass beads , 1981 .

[5]  Alain Bourgeat,et al.  TWO-PHASE FLOW IN HETEROGENEOUS POROUS MEDIA , 1991 .

[6]  Gary R. Olhoeft,et al.  Geophysics and solvents: The Borden experiment , 1993 .

[7]  Michael E. Cates,et al.  Seismic monitoring of a CO2 flood in a carbonate reservoir: A rock physics study , 1998 .

[8]  D. Oldenburg,et al.  NON-LINEAR INVERSION USING GENERAL MEASURES OF DATA MISFIT AND MODEL STRUCTURE , 1998 .

[9]  S. Treitel,et al.  Fast l p solution of large, sparse, linear systems: application to seismic travel time tomography , 1988 .

[10]  Amos Nur,et al.  Effects of CO2 Flooding on Wave Velocities in Rocks With Hydrocarbons , 1989 .

[11]  J. Nath,et al.  BINARY SYSTEMS OF TETRACHLOROETHYLENE WITH BENZENE, TOLUENE, P‐XYLENE, CARBON TETRACHLORIDE, AND CYCLOHEXANE. 1. ULTRASONIC VELOCITIES AND ADIABATIC COMPRESSIBILITIES AT 293.15 AND 303.15 K, DIELECTRIC CONSTANTS AT 298.15 AND 308.15 K, AND REFRACTIVE INDEXES AT 298.15 K , 1982 .

[12]  Terrance J. Fulp,et al.  AAPG Memoir 42 and SEG Investigations in Geophysics, No. 9, Chapter 9 (Case Histories of Three-Dimensional Seismic Surveys) -- Case History 4: Three-Dimensional Seismic Monitoring of An Enhanced Oil Recovery Process , 1987 .

[13]  Michael S. Zhdanov,et al.  Minimum support nonlinear parametrization in the solution of a 3D magnetotelluric inverse problem , 2004 .

[14]  Per Christian.  ANALYSIS OF DISCRETE ILL-POSED PROBLEMS BY MEANS OF THE L-CURVE* , 1992 .

[15]  Bruce P. Marion,et al.  Crosswell seismic imaging of reservoir changes caused by CO 2 injection , 1997 .

[16]  D. Oldenburg,et al.  A comparison of automatic techniques for estimating the regularization parameter in non-linear inverse problems , 2004 .

[17]  Per Christian Hansen,et al.  Analysis of Discrete Ill-Posed Problems by Means of the L-Curve , 1992, SIAM Rev..

[18]  Michael S. Zhdanov,et al.  Focusing geophysical inversion images , 1999 .

[19]  R. Parker,et al.  Occam's inversion; a practical algorithm for generating smooth models from electromagnetic sounding data , 1987 .

[20]  David E. Dougherty,et al.  Modified total variation methods for three‐dimensional electrical resistance tomography inverse problems , 2000 .

[21]  Emil O. Frind,et al.  Two‐phase flow in heterogeneous porous media: 1. Model development , 1991 .

[22]  Variability of point source infiltration rates for two-phase flow in heterogeneous porous media , 1996 .

[23]  A. Guillén,et al.  Gravity and magnetic inversion with minimization of a specific functional , 1984 .

[24]  J. Claerbout,et al.  Robust Modeling With Erratic Data , 1973 .

[25]  A. Peter Annan,et al.  Ground‐penetrating radar monitoring of a controlled DNAPL release: 200 MHz radar , 1994 .

[26]  J. Harris,et al.  High-resolution crosswell imaging of a West Texas carbonate reservoir; Part 5, Core analysis , 1995 .

[27]  G. Wahba.  Spline models for observational data , 1990 .

[28]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[29]  Clifford H. Thurber,et al.  Parameter estimation and inverse problems , 2005 .

[30]  Jonathan B. Ajo-Franklin,et al.  The dielectric properties of granular media saturated with DNAPL/water mixtures , 2004 .

[31]  A. N. Tikhonov,et al.  Solutions of ill-posed problems , 1977 .

[32]  M. Sacchi,et al.  Edge Preserving Imaging , 2007 .

[33]  Sally M. Benson,et al.  Measuring permanence of CO2 storage in saline formations: the Frio experiment , 2005 .

[34]  K. Broholm,et al.  Dense chlorinated solvents and other DNAPLs in groundwater , 1996 .

[35]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[36]  Elena Cherkaev,et al.  Non-smooth gravity problem with total variation penalization functional , 2002 .

[37]  Frederick D. Day-Lewis,et al.  Time‐lapse imaging of saline‐tracer transport in fractured rock using difference‐attenuation radar tomography , 2003 .

[38]  T. M. Daley,et al.  Acquisition of time-lapse, 6-component, P- and S-wave, crosswell seismic survey with orbital vibrator and of time-lapse VSP for CO2 injection monitoring , 2004 .

[39]  Harris,et al.  High Resolution Crosswell Imaging of a West Texas Carbonate Reservoir , 1993 .

[40]  S. Benson.  Geologic Sequestration of Carbon Dioxide , 2003 .

[41]  K. Bube,et al.  Hybrid l 1 /l 2 minimization with applications to tomography , 1997 .

[42]  Jason I. Gerhard,et al.  Variability of point source infiltration rates for two‐phase flow in heterogeneous porous media , 1995 .

[43]  John A. Cherry,et al.  Dense Chlorinated Solvents and other DNAPLs in Groundwater , 1996 .

[44]  David Redman,et al.  A Field Experiment to Study the Behavior of Tetrachloroethylene Below the Water Table: Spatial Distribution of Residual and Pooled DNAPL , 1993 .

[45]  G. Wahba.  Practical Approximate Solutions to Linear Operator Equations When the Data are Noisy , 1977 .

[46]  R. J. Greaves,et al.  Three-dimensional seismic monitoring of enhanced oil recovery project , 1983 .

[47]  Brian R. Strazisar,et al.  Preliminary modeling of the long-term fate of CO2 following injection into deep geological formations , 2006 .

[48]  Dianne P. O'Leary,et al.  The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems , 1993, SIAM J. Sci. Comput..

[49]  John E. Peterson,et al.  Cost-effective imaging of CO2 injection with borehole seismic methods , 2006 .