Curvelet imaging and processing: adaptive multiple elimination Felix J. Herrmann1 and Eric Verschuur21 Department of Earth and Ocean Sciences, University of British Columbia, Canada2 Faculty of Applied Sciences, Delft University of Technology, The NetherlandsAbstractPredictive multiple suppression methods consist of two main steps: a prediction step, in which multiples are predicted from the seismic data, and a subtraction step, in which the predicted multiples are matched with the true multiples in the data. The last step appears crucial in practice: an incorrect adaptive subtraction method will cause multiples to be sub-optimally subtracted or primaries being distorted, or both. Therefore, we propose a new domain for separation of primaries and multiples via the Curvelet transform. This transform maps the data into almost orthogonal localized events with a directional and spatial-temporal component. The multiples are suppressed by thresholding the input data at those Curvelet components where the predicted multiples have large amplitudes. In this way the more traditional filtering of predicted multiples to fit the input data is avoided. An initial field data example shows a considerable improvement in multiple suppression.IntroductionIn complex areas move-out filtering multiple suppression techniques may fail because underlying assumptions are not met. Several attempts have been made to address this problem by either extending move-out discrimination methods towards 3D complexities (e.g. by introducing apex-shifted hyperbolic transforms [7]) or by coming up with matching techniques in the wave-equation based predictive methods [see e.g. 17,1]. Least-squares matching the predicted multiples in time and space overlapping windows, [16] provides a straightforward subtraction method, where the predicted multiples are matched to the true multiples for 2-D input data. Unfortunately, this matching procedure fails when the underlying 2D assumption are severely violated. There have been several attempts to address this issue and the proposed solutions range from including surrounding shot positions [11] to methods based on model- [14], data-driven [10] time delays and separation of predicted multiples into (in)-coherent parts [12]. Even though these recent advances in adaptive subtraction and other techniques have improved the attenuation of multiples, these methods continue to suffer from (i) a relative strong sensitivity to the accuracy of the predicted multiples; (ii) creation of spurious artifacts or worse (iii) a possible distortions of the primary energy. For these situations, subtraction techniques based on a different concept are needed to complement the processor?s tool box.The method we are proposing here holds the middle between two complementary approaches common in multiple elimination: prediction in combination with subtraction and filtering [15,7]. Whereas the first approach aims to predict the multiples and then subtract, the second approach tries to find a domain in which the primaries and multiples separate, followed by some filtering operation and reconstruction. Our method is not distant from either since it uses the predicted multiples to non-linearly filter data in a domain spanned by almost orthogonal and local basis functions. We use the recently developed Curvelet transform [see e.g. 4], that decomposes data into basis functions that not only obtain optimal sparseness on the coefficients and hence reduce the dimensionality of the problem but which are also local in both location and angle/dip, facilitating the definition of non-linear estimators based on thresholding. Main assumption of this proposal is that multiples and primaries have locally a different temporal, spatial and dip behavior, and therefore map into different areas in the Curvelet domain. Multiples give rise to large Curvelet coefficients in the input and these coefficients can be muted by our estimation procedure when the threshold is set according to the Curvelet transform of the predicted multiples. As such, our suppression technique has at each location in the transformed domain one parameter, namely the threshold yielded by the predicted multiple, beyond which the input data is suppressed. In that sense, our procedure is similar to the ones proposed by [18] and [15], although the latter use the non-localized FK/Radon domains for their separation while we use localized basis functions and non-linear estimation by thresholding. Non-locality and non-optimality in their approximation renders the first filtering techniques less effective because primaries and multiples will still have a considerable overlap. The Curvelet transform is able to make a local discrimination between interfering events with different temporal and spatial characteristics.Adaptive subtraction by non-linear estimation of the Curvelet coefficientsThe denoising problemRemoving predicted multiples can be seen as a particular instance of a generic denoising problem that estimates the model (primaries) m from noisy data (data including multiples) [See e.g., 13] EMBED Equation.3 (1)with n white Gaussian noise. Main difference between this problem and ours is that in our case the noise is coherent and considered to be the multiples to which we have access through prediction.Question now is: how can we solve this denoising problem effectively? In other words, how can we construct a diagonal decision operator that minimizes the energy difference between the estimate and the true model. It appears from the work of [6], [13] and others that, for a certain class of models, one can obtain nearly optimal denoising results, i.e. near optimal SNR for denoised data, by projecting noisy data onto a basis-function representation that is optimal for that particular class of models. In that case, most of the model?s energy resides in only a few coefficients, allowing for the definition of a shrinkage estimator that separates noise from the model. For basis functions that are also local, one can show that soft thresholding on the coefficients suffices to approximately solve the denoising problem, i.e. EMBED Equation.3 (2)In this expression, EMBED Equation.3 refers to the (pseudo)-inverse of EMBED Equation.3 , which is the basis-function expansion. EMBED Equation.3 is a soft/hard thresholding operator with a threshold that for orthonormal basis functions equals [13,6, ] EMBED Equation.3 with EMBED Equation.3 the standard deviation of the noise and N the number of data samples. Question now is: can we extend these results to colored noise and to be specific to noise given by the predicted multiples? Before answering this question let us first be more specific with respect to the choice of the appropriate basis functions for seismic data, primaries as well as multiples.The basis functionsCurvelets as proposed by [3], constitute a relatively new family of non-separable wavelet bases that are designed to effectively represent seismic data with reflectors that generally tend to lie on piece-wise smooth curves. This property makes Curvelets suitable to represent events in seismic whether these are located in shot records or time slices. For these type of signals, Curvelets obtain nearly optimal sparseness, because of (i) the rapid decay for the reconstruction error as a function of the largest coefficients; (ii) the ability to concentrate the signal?s energy in a limited number of coefficients; (iii) the ability to map noise and signal to different areas in the Curvelet domain. So how do Curvelets obtain such a high non-linear approximation rate? Without being all inclusive [see for details, 2,4,3,5], the answer to this question lies in the fact that Curvelets are ? multi-scale, i.e. they live in different dyadic corona (see for more detail [3] or the other contributions of the fist author to the proceedings of this conference) in the FK-domain. ? multi-directional, i.e. they live on wedges within these corona. ? anisotropic, i.e.?they obey the following scaling law width ?length2. ? directional selective with EQ # orientations? \F(1, \R(,scale)). ? local both in (x,t)) and KF. ? almost ortho rectional selective and anisotropic for the higher frequencies. They are localized in both the space (or (x,t)) and spatial KF-domains and have, as consequence of their partitioning, the tendency to align themselves with curves/wavefronts. As such they can be more flexible then a representation yielded by high-resolution Radon [as described by e.g. Trad [15] because they are local and able to follow any piece-wise smooth curve.A tantalizing perspectiveAs the examples in the next section clearly demonstrate, the optimal denoising capabilities for incoherent noise carry over to coherent noise removal provided we have reasonable accurate predictions for the noise, the multiples in this case. By choosing a threshold defined by the predicted multiples, i.e. EMBED Equation.3 (3)we are able to awe are able to adaptively decide whether a certain event belongs to primary or multiple energy. The n contains the predicted multiples and h represents an additional control parameter which sets the confidence interval (e.g. 95 % for EMBED Equation.3 =3) (de)-emphasizing the thresholding.Results of multiple subtraction via thresholding in the Curvelet domainThe above methodology is tested on a field dataset from offshore Scotland. This is a 2D line, which is known to suffer from strong surface-related multiples. Because the geology of the shallow sub-bottom layers is laterally complex, multiples have been observed to exhibit 3D characteristics. Applying surface-related multiple elimination to this data gives only an attenuation of surface multiples, but not a complete suppression. This lack of suppression can be observed for a shot record shown in Fig. 1. When comparing the predicted multiples with the input data a good resemblance is observed in a global sense. However, when a matching filter is calculated, it appears that the predicted multiples do not coincide well enough with the true multiples, mainly due to 3D effects. Locally the temporal and spatial shape of the predicted multiples - especially for the higher order multiples - differs from the true events. Estimating shaping filters in spatially and temporally varying, overlapping windows, as described by [16] can not resolve this mismatch completely. The adaptive subtraction result is displayed in Fig. 1 in the third panel. Figure 1: Comparison between adaptive least-squares (using time and space varying least-squares filters)and Curvelet subtraction of predicted multiples for one shot record. Input data (left), predicted multiples (second), least-squares subtraction (third) and the results by thresholding Curvelet (right).Figure 2: Time slice through the prestack data volume of all shot records. Top left: Time slice through the input data with all multiples. Top right: Slice through the 2D predicted surface-related multiples. Bottom left: Time slice through the shot records after adaptive subtraction per shot record of the predicted multiples. Bottom right: The Curvelet equivalent of the adaptive subtraction computed for a single time horizon.To improve this result, more freedom can be incorporated in the subtraction process. However, note that the more the predicted multiples can adapt to the input data, the higher the chance that primary energy will be distorted as well, as the primary and multiple energy are not orthogonal to each other in the least-squares sense. Multi-gather subtraction using 9 surrounding multiple panels gives better suppression of the multiples but leaves a lot of multiple remnants to be observed, especially for the large offset at large travel times. The output gather of the Curvelet method (see fourth panel in Fig. 1) looks much cleaner. Furthermore, clear events have been restored from interference with the multiples in the lower left area. Also note the good preservation of a primary event around -1000 meter offset and 1.1 seconds. The amplitudes seem to be identical to the ones in the input data, whereas the adaptive subtraction has distorted these amplitudes to some extend. For the Curvelet-domain procedure, the threshold value can be set by the user, thus creating more or less suppression. Too weak thresholding leaves too much multiples in the data and that a too strong thresholding procedure seems to remove too much energy. The filtering via the Curvelet domain can also be applied within each 2D cross-section through the seismic data volume. Time slices at t=1.48 through the prestack shot-offset data and predicted-multiple volumes are displayed in Fig. 2 (top row). Results after application of the adaptive subtraction per shot record and the time-slice only Curvelet subtraction are included in Fig. 2 (bottom row). Again, the Curvelet results look much cleaner than the least-squares subtraction result. Much of the noisy remnants at small offsets in the least-squares subtraction result have been removed. Notice also the improved suppression of the higher order peg-legs in the area around shot 900 and offset 1500 meter.ConclusionsIn this paper a new concept related to multiple subtraction has been described, based on the Curvelet transform. The Curvelet transform is an almost orthogonal transformation into local basis functions parameterized by their relate temporal- and spatial-frequency content. Because of their anisotropic shape, Curvelets are directional selective, i.e. they have local angle-discrimination capabilities. Our method uses the predicted multiples from the surface-related multiple prediction method, as a guide to suppress the Curvelet coefficients related to the multiple events directly in the original data. Therefore, the multiples are not actually subtracted, but areas in the Curvelet domain related to multiple energy are muted. The success of this method depends on the assumption that primaries and multiples map into different areas in the Curvelet domain. Based on a field data example, we can conclude that the Curvelet-based multiple filtering is effective and is able to suppress multiples, while preserving primary energy. Especially in situations with clear 3D effects in the 2D seismic data, it appears to perform better than the more traditional least-squares adaptive subtraction methods. This success does not really come as a surprise given the successful application of these techniques to the removal of noise colored by migration [8,9] and to the computation of 4D difference cubes (see for details in both other contributions of the first to proceedings of this conference). Leaves us to hope for future higher dimensional implementation of the Curvelet transform possibly supplemented by constraint optimization imposing sparseness constraints (see the migration paper in these proceedings).AcknowledgmentsThe authors thank Elf Caledonia Ltd (now part of Total) for providing the field data and Emmanuel Cand?s and David Donoho for making an early version of their Curvelet code (Digital Curvelet Transforms via Unequispaced Fourier Transforms, presented at the ONR Meeting, University of Minnesota, May, 2003) available for evaluation. This work was in part financially supported by a NSERC Discovery Grant.References[1] A.?J. Berkhout and D.?J. Verschuur. Estimation of multiple scattering by iterative inversion, part?I: ?theoretical considerations. Geophysics, 62(5): 1586?1595, 1997.[2] E.?J. Cand?s and D.?L. Donoho. Curvelets ? a surprisingly effective nonadaptive representation for objects with edges. Curves and Surfaces. Vanderbilt University Press, 2000.[3] E.?J. Cand?s and F.?Guo. New multiscale transforms, minimum total variation synthesis: Applications to edge-preserving image reconstruction. Signal Processing, pages 1519?1543, 2002.[4] E.?J. E.?J. Cand?s and D.?Donoho. New tight frames of curvelets and optimal representations of objects with smooth singularities. Technical report, Stanford, 2002. submitted.[5] M.?Do and M.?Vetterli. Beyond wavelets, chapter Contourlets. Academic Press, 2002.[6] D.?L. Donoho and I.?M. Johnstone. Minimax estimation via wavelet shrinkage. Annals of Statistics, 26(3): 879?921, 1998.[7] N.?Hargreaves, B.?verWest, R.?Wombell, and D.?Trad. Multiple attenuation using apex-shifted radon transform. pages 1929?1932, Dallas, 2003. SEG, Soc. Expl. Geophys., Expanded abstracts.[8] F.?J. Herrmann. Multifractional splines: application to seismic imaging. In A.?F. L.?E. Michael A.?Unser, Akram?Aldroubi, editor, Proceedings of SPIE Technical Conference on Wavelets: Applications in Signal and Image Processing X, volume 5207, pages 240?258. SPIE, 2003.[9] F.?J. Herrmann. Optimal seismic imaging with curvelets. In Expanded Abstracts, Tulsa, 2003. Soc. Expl. Geophys.[10] L.?T. Ikelle and S.?Yoo. An analysis of 2D and 3D inverse scattering multiple attenuation. pages 1973?1976, Calgary, 2000. SEG, Soc. Expl. Geophys., Expanded abstracts.[11] H.?Jakubowicz. Extended subtraction of multiples. Private communication, Delft, 1999.[12] M.?M.?N. Kabir. Weighted subtraction for diffracted multiple attenuation. pages 1941?1944, Dallas, 2003. SEG, Soc. Expl. Geophys., Expanded abstracts.[13] S.?G. Mallat. A wavelet tour of signal processing. Academic Press, 1997.[14] W.?S. Ross. Multiple suppression: beyond 2-D. part I: theory. Delft Univ. Tech., pages 1387?1390, Dallas, 1997. Soc. Expl. Geophys., Expanded abstracts.[15] D.?O. Trad . Latest views of the sparse Radon transform. Geophysics, 68(1) : 386 ? 399, 2003[16] D.?J. Verschuur and A.?J. Berkhout. Estimation of multiple scattering by iterative inversion, part?II: ?practical?aspects?and?examples. Geophysics, 62(5): 1596?1611, 1997.[17] D.?J. Verschuur, A.?J. Berkhout, and C.?P.?A. Wapenaar. Adaptive surface-related multiple elimination. Geophysics, 57(9): 1166?1177, 1992.[18] B.?Zhou and S.?A. Greenhalgh. Wave-equation extrapolation-based multiple attenuation: 2-d filtering in the f-k domain. Geophysics, 59(10): 1377?1391, 1994.Great Explorations ? Canada and BeyondPAGE 4 Great Explorations ? Canada and BeyondPAGE 3 Curvelet imaging and processing: adaptive multiple elimination Felix J. Herrmann1 and Eric Verschuur2 1 Department of Earth and Ocean Sciences, University of British Columbia, Canada 2 Faculty of Applied Sciences, Delft University of Technology, The Netherlands Abstrac t Predictive multiple suppression methods consist of two main steps: a prediction step, in which multiples are predicted from the seismic data, and a subtraction step, in which the predicted multiples are matched with the true multiples in the data. The last step appears crucial in practice: an incorrect adaptive subtraction method will cause multiples to be sub-optimally subtracted or primaries being distorted, or both. Therefore, we propose a new domain for separation of primaries and multiples via the Curvelet transform. This transform maps the data into almost orthogonal localized events with a directional and spatial-temporal component. The multiples are suppressed by thresholding the input data at those Curvelet components where the predicted multiples have large amplitudes. In this way the more traditional filtering of predicted multiples to fit the input data is avoided. An initial field data example shows a considerable improvement in multiple suppression. In t r oduc t i on In complex areas move-out filtering multiple suppression techniques may fail because underlying assumptions are not met. Several attempts have been made to address this problem by either extending move-out discrimination methods towards 3D complexities (e.g. by introducing apex-shifted hyperbolic transforms [7]) or by coming up with matching techniques in the wave-equation based predictive methods [see e.g. 17,1]. Least-squares matching the predicted multiples in time and space overlapping windows, [16] provides a straightforward subtraction method, where the predicted multiples are matched to the true multiples for 2-D input data. Unfortunately, this matching procedure fails when the underlying 2D assumption are severely violated. There have been several attempts to address this issue and the proposed solutions range from including surrounding shot positions [11] to methods based on model- [14], data-driven [10] time delays and separation of predicted multiples into (in)-coherent parts [12]. Even though these recent advances in adaptive subtraction and other techniques have improved the attenuation of multiples, these methods continue to suffer from (i) a relative strong sensitivity to the accuracy of the predicted multiples; (ii) creation of spurious artifacts or worse (iii) a possible distortions of the primary energy. For these situations, subtraction techniques based on a different concept are needed to complement the processor’s tool box. The method we are proposing here holds the middle between two complementary approaches common in multiple elimination: prediction in combination with subtraction and filtering [15,7]. Whereas the first approach aims to predict the multiples and then subtract, the second approach tries to find a domain in which the primaries and multiples separate, followed by some filtering operation and reconstruction. Our method is not distant from either since it uses the predicted multiples to non-linearly filter data in a domain spanned by almost orthogonal and local basis functions. We use the recently developed Curvelet transform [see e.g. 4], that decomposes data into basis functions that not only obtain optimal sparseness on the coefficients and hence reduce the dimensionality of the problem but which are also local in both location and angle/dip, facilitating the definition of non-linear estimators based on thresholding. Main assumption of this proposal is that multiples and primaries have locally a different temporal, spatial and dip behavior, and therefore map into different areas in the Curvelet domain. Multiples give rise to large Curvelet coefficients in the input and these coefficients can be muted by our estimation procedure when the threshold is set according to the Curvelet transform of the predicted multiples. As such, our suppression technique has at each location in the transformed domain one parameter, namely the threshold yielded by the predicted multiple, beyond which the input data is suppressed. In that sense, our procedure is similar to the ones proposed by [18] and [15], although the latter use the non-localized FK/Radon domains for their separation while we use localized basis functions and non-linear estimation by thresholding. Non-locality and non-optimality in their approximation renders the first filtering techniques less effective because primaries and multiples will still have a considerable overlap. The Curvelet transform is able to make a local discrimination between interfering events with different temporal and spatial characteristics. Adapt ive subtrac t i on by non- l i nea r est imat i on of the Curvele t coef f i c i e n t s The denoising problem Removing predicted multiples can be seen as a particular instance of a generic denoising problem that estimates the model (primaries) m from noisy data (data including multiples) [See e.g., 13] d = m+ n (1) with n white Gaussian noise. Main difference between this problem and ours is that in our case the noise is coherent and considered to be the multiples to which we have access through prediction. Question now is: how can we solve this denoising problem effectively? In other words, how can we construct a diagonal decision operator that minimizes the energy difference between the estimate and the true model. It appears from the work of [6], [13] and others that, for a certain class of models, one can obtain nearly optimal denoising results, i.e. near optimal SNR for denoised data, by projecting noisy data onto a basis-function representation that is optimal for that particular class of models. In that case, most of the model’s energy resides in only a few Grea t Exp lo r a t i o n s – Canada and Beyond 1 coefficients, allowing for the definition of a shrinkage estimator that separates noise from the model. For basis functions that are also local, one can show that soft thresholding on the coefficients suffices to approximately solve the denoising problem, i.e. ö = B- 1 Q m (Bd) m (2) In this expression, B- 1 refers to the (pseudo)-inverse of B , which is the basis-function expansion. Θ µ is a soft/hard thresholding operator with a threshold that for orthonormal basis functions equals [13,6, ] µ = σ 2log e N with σ the standard deviation of the noise and N the number of data samples. Question now is: can we extend these results to colored noise and to be specific to noise given by the predicted multiples? Before answering this question let us first be more specific with respect to the choice of the appropriate basis functions for seismic data, primaries as well as multiples. The basis functions Curvelets as proposed by [3], constitute a relatively new family of non-separable wavelet bases that are designed to effectively represent seismic data with reflectors that generally tend to lie on piece-wise smooth curves. This property makes Curvelets suitable to represent events in seismic whether these are located in shot records or time slices. For these type of signals, Curvelets obtain nearly optimal sparseness, because of (i) the rapid decay for the reconstruction error as a function of the largest coefficients; (ii) the ability to concentrate the signal’s energy in a limited number of coefficients; (iii) the ability to map noise and signal to different areas in the Curvelet domain. So how do Curvelets obtain such a high non-linear approximation rate? Without being all inclusive [see for details, 2,4,3,5], the answer to this question lies in the fact that Curvelets are • multi-scale, i.e. they live in different dyadic corona (see for more detail [3] or the other contributions of the fist author to the proceedings of this conference) in the FK-domain. • multi-directional, i.e. they live on wedges within these corona. • anisotropic, i.e. they obey the following scaling law width∝length2. • directional selective with . • local both in (x,t)) and KF. • almost orthogonal, they are tight frames with a moderate redundancy. Curvelets live in a wedges of the 2-D Fourier plane and become more directional selective and anisotropic for the higher frequencies. They are localized in both the space (or (x,t)) and spatial KF-domains and have, as consequence of their partitioning, the tendency to align themselves with curves/wavefronts. As such they can be more flexible then a representation yielded by high-resolution Radon [as described by e.g. Trad [15] because they are local and able to follow any piece-wise smooth curve. A tantalizing perspective As the examples in the next section clearly demonstrate, the optimal denoising capabilities for incoherent noise carry over to coherent noise removal provided we have reasonable accurate predictions for the noise, the multiples in this case. By choosing a threshold defined by the predicted multiples, i.e. µ ≈ η Bn (3) we are able to awe are able to adaptively decide whether a certain event belongs to primary or multiple energy. The n contains the predicted multiples and h represents an additional control parameter which sets the confidence interval (e.g. 95 % for η =3) (de)-emphasizing the thresholding. Results of multiple subtraction via thresholding in the Curvelet domain The above methodology is tested on a field dataset from offshore Scotland. This is a 2D line, which is known to suffer from strong surfacerelated multiples. Because the geology of the shallow sub-bottom layers is laterally complex, multiples have been observed to exhibit 3D characteristics. Applying surface-related multiple elimination to this data gives only an attenuation of surface multiples, but not a complete suppression. This lack of suppression can be observed for a shot record shown in Fig. 1. When comparing the predicted multiples with the input data a good resemblance is observed in a global sense. However, when a matching filter is calculated, it appears that the predicted multiples do not coincide well enough with the true multiples, mainly due to 3D effects. Locally the temporal and spatial shape of the predicted multiples - especially for the higher order multiples - differs from the true events. Estimating shaping filters in spatially and temporally varying, overlapping windows, as described by [16] can not resolve this mismatch completely. The adaptive subtraction result is displayed in Fig. 1 in the third panel. Grea t Exp lo r a t i o n s – Canada and Beyond 2 Figure 1: Comparison between adaptive least-squares (using time and space varying least-squares filters)and Curvelet subtraction of predicted multiples for one shot record. Input data (left), predicted multiples (second), least-squares subtraction (third) and the results by thresholding Curvelet (right). Figure 2: Time slice through the prestack data volume of all shot records. Top left: Time slice through the input data with all multiples. Top right: Slice through the 2D predicted surface-related multiples. Bottom left: Time slice through the shot records after adaptive subtraction per shot record of the predicted multiples. Bottom right: The Curvelet equivalent of the adaptive subtraction computed for a single time horizon. To improve this result, more freedom can be incorporated in the subtraction process. However, note that the more the predicted multiples can adapt to the input data, the higher the chance that primary energy will be distorted as well, as the primary and multiple energy are not orthogonal to each other in the least-squares sense. Multi-gather subtraction using 9 surrounding multiple panels gives better suppression of the multiples but leaves a lot of multiple remnants to be observed, especially for the large offset at large travel times. The output gather of the Curvelet method (see fourth panel in Fig. 1) looks much cleaner. Furthermore, clear events have been restored from interference with the multiples in the lower left area. Also note the good preservation of a primary event around -1000 meter offset and 1.1 seconds. The amplitudes seem to be identical to the ones in the input data, whereas the adaptive subtraction has distorted these amplitudes to some extend. For the Curvelet-domain procedure, the threshold value can be set by the user, thus creating more or less suppression. Too weak thresholding leaves too much multiples in the data and that a too strong thresholding procedure seems to remove too much energy. The filtering via the Curvelet domain can also be applied within each 2D cross-section through the seismic data volume. Time slices at t=1.48 through the prestack shot-offset data and predicted-multiple volumes are displayed in Fig. 2 (top row). Results after application of the adaptive subtraction per shot record and the time-slice only Curvelet subtraction are included in Fig. 2 (bottom row). Again, the Curvelet results look much cleaner than the least-squares subtraction result. Much of the noisy remnants at small offsets in the least-squares subtraction result have been removed. Notice also the improved suppression of the higher order peg-legs in the area around shot 900 and offset 1500 meter. Grea t Exp lo r a t i o n s – Canada and Beyond 3 Conclusions In this paper a new concept related to multiple subtraction has been described, based on the Curvelet transform. The Curvelet transform is an almost orthogonal transformation into local basis functions parameterized by their relate temporal- and spatial-frequency content. Because of their anisotropic shape, Curvelets are directional selective, i.e. they have local angle-discrimination capabilities. Our method uses the predicted multiples from the surface-related multiple prediction method, as a guide to suppress the Curvelet coefficients related to the multiple events directly in the original data. Therefore, the multiples are not actually subtracted, but areas in the Curvelet domain related to multiple energy are muted. The success of this method depends on the assumption that primaries and multiples map into different areas in the Curvelet domain. Based on a field data example, we can conclude that the Curvelet-based multiple filtering is effective and is able to suppress multiples, while preserving primary energy. Especially in situations with clear 3D effects in the 2D seismic data, it appears to perform better than the more traditional least-squares adaptive subtraction methods. This success does not really come as a surprise given the successful application of these techniques to the removal of noise colored by migration [8,9] and to the computation of 4D difference cubes (see for details in both other contributions of the first to proceedings of this conference). Leaves us to hope for future higher dimensional implementation of the Curvelet transform possibly supplemented by constraint optimization imposing sparseness constraints (see the migration paper in these proceedings). Acknowledgments The authors thank Elf Caledonia Ltd (now part of Total) for providing the field data and Emmanuel Candés and David Donoho for making an early version of their Curvelet code (Digital Curvelet Transforms via Unequispaced Fourier Transforms, presented at the ONR Meeting, University of Minnesota, May, 2003) available for evaluation. This work was in part financially supported by a NSERC Discovery Grant. Re fe r ences [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] A. J. Berkhout and D. J. Verschuur. Estimation of multiple scattering by iterative inversion, part I: theoretical considerations. Geophysics, 62(5): 1586– 1595, 1997. E. J. Candès and D. L. Donoho. Curvelets – a surprisingly effective nonadaptive representation for objects with edges. Curves and Surfaces. Vanderbilt University Press, 2000. E. J. Candès and F. Guo. New multiscale transforms, minimum total variation synthesis: Applications to edge-preserving image reconstruction. Signal Processing, pages 1519–1543, 2002. E. J. E. J. Candès and D. Donoho. New tight frames of curvelets and optimal representations of objects with smooth singularities. Technical report, Stanford, 2002. submitted. M. Do and M. Vetterli. Beyond wavelets, chapter Contourlets. Academic Press, 2002. D. L. Donoho and I. M. Johnstone. Minimax estimation via wavelet shrinkage. Annals of Statistics, 26(3): 879–921, 1998. N. Hargreaves, B. verWest, R. Wombell, and D. Trad. Multiple attenuation using apex-shifted radon transform. pages 1929–1932, Dallas, 2003. SEG, Soc. Expl. Geophys., Expanded abstracts. F. J. Herrmann. Multifractional splines: application to seismic imaging. In A. F. L. E. Michael A. Unser, Akram Aldroubi, editor, Proceedings of SPIE Technical Conference on Wavelets: Applications in Signal and Image Processing X, volume 5207, pages 240–258. SPIE, 2003. F. J. Herrmann. Optimal seismic imaging with curvelets. In Expanded Abstracts, Tulsa, 2003. Soc. Expl. Geophys. L. T. Ikelle and S. Yoo. An analysis of 2D and 3D inverse scattering multiple attenuation. pages 1973–1976, Calgary, 2000. SEG, Soc. Expl. Geophys., Expanded abstracts. H. Jakubowicz. Extended subtraction of multiples. Private communication, Delft, 1999. M. M. N. Kabir. Weighted subtraction for diffracted multiple attenuation. pages 1941–1944, Dallas, 2003. SEG, Soc. Expl. Geophys., Expanded abstracts. S. G. Mallat. A wavelet tour of signal processing. Academic Press, 1997. W. S. Ross. Multiple suppression: beyond 2-D. part I: theory. Delft Univ. Tech., pages 1387–1390, Dallas, 1997. Soc. Expl. Geophys., Expanded abstracts. D. O. Trad . Latest views of the sparse Radon transform. Geophysics, 68(1) : 386 – 399, 2003 D. J. Verschuur and A. J. Berkhout. Estimation of multiple scattering by iterative inversion, part II: practical aspects and examples. Geophysics, 62(5): 1596–1611, 1997. D. J. Verschuur, A. J. Berkhout, and C. P. A. Wapenaar. Adaptive surface-related multiple elimination. Geophysics, 57(9): 1166–1177, 1992. B. Zhou and S. A. Greenhalgh. Wave-equation extrapolation-based multiple attenuation: 2-d filtering in the f-k domain. Geophysics, 59(10): 1377–1391, 1994. Grea t Exp lo r a t i o n s – Canada and Beyond 4
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Faculty Research and Publications /
- Curvelet imaging and processing : adaptive multiple...
Open Collections
UBC Faculty Research and Publications
Curvelet imaging and processing : adaptive multiple elimination Herrmann, Felix J.; Verschuur, Eric 2004
Page Metadata
Item Metadata
Title | Curvelet imaging and processing : adaptive multiple elimination |
Creator |
Herrmann, Felix J. Verschuur, Eric |
Publisher | Canadian Society of Exploration Geophysicists |
Date Issued | 2004 |
Description | Predictive multiple suppression methods consist of two main steps: a prediction step, in which multiples are predicted from the seismic data, and a subtraction step, in which the predicted multiples are matched with the true multiples in the data. The last step appears crucial in practice: an incorrect adaptive subtraction method will cause multiples to be sub-optimally subtracted or primaries being distorted, or both. Therefore, we propose a new domain for separation of primaries and multiples via the Curvelet transform. This transform maps the data into almost orthogonal localized events with a directional and spatial-temporal component. The multiples are suppressed by thresholding the input data at those Curvelet components where the predicted multiples have large amplitudes. In this way the more traditional filtering of predicted multiples to fit the input data is avoided. An initial field data example shows a considerable improvement in multiple suppression. |
Extent | 1092096 bytes |
Subject |
adaptive subtraction curvelets denoising curvelet transform 3D 4D |
Genre |
Conference Paper |
Type |
Text |
File Format | application/msword |
Language | eng |
Collection |
Seismic Laboratory for Imaging and Modeling |
Date Available | 2008-03-26 |
Provider | Vancouver : University of British Columbia Library |
DOI | 10.14288/1.0107421 |
Affiliation |
Science, Faculty of Earth and Ocean Sciences, Department of |
Citation | Herrmann, Felix J., Verschuur, Eric. 2004. Curvelet imaging and processing: adaptive multiple elimination. CSEG National Convention. |
Peer Review Status | Unreviewed |
Scholarly Level | Faculty |
Copyright Holder | Herrmann Felix J. |
URI | http://hdl.handle.net/2429/612 |
Aggregated Source Repository | DSpace |
Download
- Media
- [if-you-see-this-DO-NOT-CLICK]
- Herrmann_F_Curvelet_multiple.doc [ 1.04MB ]
- [if-you-see-this-DO-NOT-CLICK]
- Metadata
- JSON: 1.0107421.json
- JSON-LD: 1.0107421+ld.json
- RDF/XML (Pretty): 1.0107421.xml
- RDF/JSON: 1.0107421+rdf.json
- Turtle: 1.0107421+rdf-turtle.txt
- N-Triples: 1.0107421+rdf-ntriples.txt
- Original Record: 1.0107421 +original-record.json
- Full Text
- 1.0107421.txt
- Citation
- 1.0107421.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Country | Views | Downloads |
---|---|---|
China | 10 | 11 |
United States | 5 | 0 |
Japan | 3 | 0 |
Germany | 3 | 0 |
Canada | 1 | 0 |
City | Views | Downloads |
---|---|---|
Shenzhen | 8 | 11 |
Ashburn | 5 | 0 |
Unknown | 3 | 0 |
Tokyo | 3 | 0 |
Beijing | 2 | 0 |
Ottawa | 1 | 0 |
{[{ mDataHeader[type] }]} | {[{ month[type] }]} | {[{ tData[type] }]} |
Share
Embed
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"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.340.1-0107421/manifest