UBC Faculty Research and Publications

Seismic imaging and processing with curvelets Herrmann, Felix J.; Hennenfent, Gilles; Moghaddam, Peyman P. 2007

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Item Metadata


52383-herrmann07b.pdf [ 944.97kB ]
JSON: 52383-1.0107404.json
JSON-LD: 52383-1.0107404-ld.json
RDF/XML (Pretty): 52383-1.0107404-rdf.xml
RDF/JSON: 52383-1.0107404-rdf.json
Turtle: 52383-1.0107404-turtle.txt
N-Triples: 52383-1.0107404-rdf-ntriples.txt
Original Record: 52383-1.0107404-source.json
Full Text

Full Text

ABSTRACT In this paper, we present a nonlinear curvelet-based sparsity-promoting formulation for three problems in seismic processing and imaging namely, seismic data regularization from data with large percentages of traces missing; seismic amplitude recovery for subsalt images obtained by reverse-time migration and primary-multiple separation, given an inaccurate multiple prediction. We argue why these nonlinear formulations are beneficial.  EAGE 69th Conference & Exhibition — London, UK, 11 - 14 June 2007  In this paper, we report recent developments on the application of the curvelet transform (see e.g. Candes et al., 2005a; Hennenfent and Herrmann, 2006b) to seismic processing and imaging. Our approach banks on two fundamental properties of curvelets, namely the • detection of wave-fronts without requiring prior information on the dips or on the velocity model; • invariance under the action of wave propagation. These two properties make this transform suitable for a robust formulation of problems such as seismic data regularization (Hennenfent and Herrmann, 2006a; Herrmann and Hennenfent, 2007), migration-amplitude recovery (Herrmann et al., 2006b) and primary multiple separation (Herrmann et al., 2006a). All these methods exploit sparsity in the curvelet domain that is a direct consequence of the above two properties and corresponds to a rapid decay for the magnitude-sorted curvelet coefficients. Curvelets: As can be observed from Fig. 1, curvelets are localized functions that oscillate in one direction and that are smooth in the other directions. They are multiscale and multi-directional and because of their anisotropic shape (they obey the so-called parabolic scaling relationship, yielding a width ∝ 2j/2 and a length ∝ 2j with j the scale), curvelets are optimal for detecting wavefronts. This explains their high compression rates for seismic data and images as reported in the literature.  (a)  (b)  Figure 1: Example of a 3-D curvelet. Notice the oscillations in one direction and the smoothness in the other two directions.  Sparsity promoting inversion: High compression rates for signal representations are a prerequisite for the robust formulation of stable signal recovery problems and other inverse problems. These compression rates allow for a nonlinear sparsity promoting solution. As such sparsity-promoting norm-one penalty functionals are not new to the geosciences (see for instance the seminal work of Claerbout and Muir (1973), followed by many others), where sparsity is promoted on the model. What is different in the current surge of interest in sparsity-promoting inversion, known as ’compressed sensing’ (Candes et al., 2005b; Donoho et al., 2006), is (i) the existence of sparsity promoting transforms such as the curvelet transform; (ii) the deep theoretical understanding on what the conditions are for a successful solution. This work can be seen as the application of these recent ideas to the seismic situation and involves the solution of the following norm-one nonlinear program, P :  x = arg minx x d = ST x  1  s.t.  Ax − y  2  ≤  EAGE 69th Conference & Exhibition — London, UK, 11 - 14 June 2007  (1)  in which y is the (incomplete) data, A the synthesis matrix and ST the inverse sparsity transform. Both these matrices consist of the inverse curvelet transform matrix, ST (the symbol T denoting the transpose) compounded with other operators depending on the problem. The above constrained optimization problem is solved to an accuracy of that depends on the noise level. Finally, d stands for the recovered vector with the symbol reserved for optimized quantities. CRSI: An important topic in seismic processing is the seismic regularization problem, where attempts are made to recover fully-sampled seismic data volumes from incomplete data, i.e., data with large percentages (> 80 %) of traces missing. By choosing A := RCT , S := C and y = Rd for the incomplete data, one arrives at the formulation for curvelet recovery by sparsity-promoting inversion (CRSI), which has successfully been applied to the recovery of incomplete seismic data (see e.g. Hennenfent and Herrmann, 2006a). In this formulation, R is the restriction matrix, selecting the rows from the curvelet transform matrix that correspond to active traces. As opposed to other recovery methods, such as sparse Fourier recovery and plane wave destruction, curveletbased methods have the advantage of working in situations where there are conflicting dips without stationarity assumptions. The method exploits the high-dimensional continuity of wavefronts and as Fig. 2 demonstrates, recovery results improve when using the 3-D curvelet transform compared to the 2-D transform. For application of this method to real data, refer to other contribution by the authors to the proceedings of this meeting. xr  xr  image interpolation  windowing t  t  xr xs  comparison xr  t xs  data  xr  t  t  volume interpolation  windowing  Figure 2: Illustration of sliced versus volumetric interpolation. Migration amplitude recovery: Because of the ’alleged’ invariance of curvelets under wave propagation, there has been a substantial interest in deriving migration operators in the curvelet domain (Douma and de Hoop, 2006; Chauris, 2006). In these approaches, one comes to benefit when strict sparsity is preserved. Strict sparsity is a significant stronger assumption than the preservation of high decay rates for the sorted coefficients. Curvelets are discrete and hence move around on grids and this makes it a challenge to define fast migration operators. Curvelets, however, prove to be very useful for solving the seismic amplitude recovery problem, during which curvelets are being imaged. On theoretical grounds (Herrmann et al., 2006b), one can expect the following identity to approximately hold AAT r Ψr (2) with r an appropriately chosen discrete reference vector and Ψ the discrete normal operator, formed by compounding the discrete scattering and its transpose, the migration  EAGE 69th Conference & Exhibition — London, UK, 11 - 14 June 2007  operator. The synthesis operator in this case is defined as A := CT Γ with Γ a diagonal weighting matrix. This identity diagonalizes the normal operator and allows for a stable recovery of the migration amplitudes by setting y = KT d, with KT the migration † operator and d the seismic data, and ST := AT with † the pseudo inverse. Results of this procedure on the SEG AA’ dataset with a reverse-time migration operator, are summarized in Fig. 3. The resulting image shows a nice recovery of the amplitudes. For more details on this method and on a method to remove remaining artifacts, refer to another contributions by the authors to the proceedings of this meeting.  (a)  (b)  Figure 3: Reverse-time migration on the SEG AA’ salt model. (a) Conventional migrated image y = KT d. (b) Seismic amplitude recovery through P . Notice the improved amplitudes. Some minor artifacts remain.  Primary-multiple separation: So far, we looked at exploiting the sparsity of curvelets in the data and image domains for the purpose of recovery. The ability of curvelets to detect wavefronts with conflicting dips, allows for a formulation of a coherent signal separation method that uses inaccurate predictions as weightings. By defining the synT thesis matrix as A := CT W1 CT W2 , x = x1 x2 and y = d and by setting the diagonal weighting matrices W1,2 in terms of predictions for the primaries and multiples, the solution of P separates primaries and multiples (Herrmann et al., 2006a) even for inaccurate predictions for which least-squares subtractions fails (see Fig. 4).  Discussion The methodology presented in this paper banks on two favorable properties of curvelets, namely their ability to detect wavefronts and their approximate invariance under wave propagation. These properties allow for a formulation of seismic processing and imaging problems that promote sparsity through the nonlinear optimization problem P . We showed that by compounding the curvelet transform with certain matrices each problem can be cast into one and the same optimization problem. The results show that (i) exploiting the multi-dimensional structure of seismic data with 3-D curvelets leads to a recovery scheme that is able to reconstruct fully sampled data volumes from data with > 80 % traces missing; (ii) the invariance of curvelets under the demigrationmigration operator can be used to recover the seismic amplitudes by inverting a diagonally weighted curvelet matrix. This latter approach is an extension of recent work by Symes (2006). Finally, we also showed that (iii) the sparsity of curvelet can be used to separate coherent signal components given an (inaccurate) prediction. The successful application to this wide spectrum of problems opens an enticing perspective of extending this framework to multiple prediction and the compression of (migration) operators.  EAGE 69th Conference & Exhibition — London, UK, 11 - 14 June 2007  (a)  (b)  (c)  (d)  Figure 4: Example of primary-multiple separation through P for predicted multiples with moveout errors. (a) the total data with primaries and multiples. (b) the true multiples used for the prediction. (c) the result obtained with least-squares adaptive subtraction with localized windows. (d) the result obtained with a single curvelet-domain soft thresholding with λ = 1.4. Notice that least-squares subtractions fails.  Acknowledgments: The authors would like to thank the authors of CurveLab for making their codes available. We also would like to thank Dr. William Symes for making his reverse-time migration code available to us. This work was in part financially supported by the Natural Sciences and Engineering Research Council of Canada Discovery Grant (22R81254) and Collaborative Research and Development Grant DNOISE (334810-05) of Felix J. Herrmann and was carried out as part of the SINBAD project with support, secured through ITF (the Industry Technology Facilitator), from the following organizations: BG Group, BP, Chevron, ExxonMobil and Shell.  REFERENCES Candes, E., L. Demanet, D. Donoho, and L. Ying, 2005a, Fast discrete curvelet transforms: SIAM Multiscale Model. Simul. Candes, E., J. Romberg, and T. Tao, 2005b, Stable signal recovery from incomplete and inaccurate measurements. to appear in Comm. Pure Appl. Math. Chauris, H., 2006, Seismic imaging in the curvelet domain and its implications for the curvelet design: Presented at the 76th Ann. Internat. Mtg., SEG, Soc. Expl. Geophys., Expanded abstracts. Claerbout, J. and F. Muir, 1973, Robust modeling with erratic data: Geophysics, 38, 826–844. Donoho, D., M. Elad, and V. Temlyakov, 2006, Stable recovery of sparse overcomplete representations in the presence of noise: IEEE Trans. Inform. Theory. Douma, H. and M. de Hoop, 2006, Leading-order seismic imaging using curvelets: Presented at the 76th Ann. Internat. Mtg., SEG, Soc. Expl. Geophys., Expanded abstracts. Hennenfent, G. and F. Herrmann, 2006a, Application of stable signal recovery to seismic interpolation: Presented at the SEG International Exposition and 76th Annual Meeting. Hennenfent, G. and F. J. Herrmann, 2006b, Seismic denoising with non-uniformly sampled curvelets: IEEE Comp. in Sci. and Eng., 8, 16–25. Herrmann, F. J., U. Boeniger, and D.-J. E. Verschuur, 2006a, Nonlinear primary-multiple separation with directional curvelet frames: Geoph. J. Int. To appear. Herrmann, F. J. and G. Hennenfent, 2007, Non-parametric seismic data recovery with curvelet frames. Submitted for publication. Herrmann, F. J., P. P. Moghaddam, and C. Stolk, 2006b, Sparsety- and continuity-promoting seismic imaging with curvelet frames. In revision. Symes, W. W., 2006, Optimal scaling for reverse time migration: Technical Report TR 06-19, Department of Computational and Applied Mathematics, Rice University, Houston, Texas, USA.  EAGE 69th Conference & Exhibition — London, UK, 11 - 14 June 2007  


Citation Scheme:


Citations by CSL (citeproc-js)

Usage Statistics



Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            async >
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:


Related Items