It is independent of the seismic source signal and there. For instance, unsuccessful interpolation leads to erroneous multiple predictions that adversely affect the performance of multiple elimination, and to imaging artifacts. For reference, the transform used to generate these images is the fdct via wrapping, complexvalued and with curvelets at the finest scale. This paper introduces ground roll suppression of a seismic signal based on a fast discrete curvelet transform fdct and the direct wave removal in ground penetrating radar gpr based on a ct. Gray and color image contrast enhancement by the curvelet transform, ieee transaction on image processing, in press. In particular, finescale basis functions are long ridges. Wavefield reconstruction is a crucial step in the seismic processing flow. After the curvelet representation of the data is computed 350, selected geophysical data attributes and their interdependencies are extracted 355, from which geological features may be identified 360, either from attribute data.
The curvelet transform candes and donoho, 2004 is designed to represent curvelike. Ground roll attenuation based on an empirical curvelet transform. The curvelet transform is a new multiscale transform with strong directional character that provides an optimal representation of objects that have discontinuities along edges. It has many applications in several fields, including denoising, wave propagation in disordered media and pattern recognition. The basic transforms are needed for subsequent seismic migration. The curvelets are localized not only in the spatial domain location and the frequency domain scale, but also in angular orientation, which is a step ahead compared to wavelet transform. Software engineering, king fahad university of petroleum and minerals, dhahran, 2004 a thesis submitted in partial fulfilment of the requirements for the degree of master of science in c fadhel alhashim 2009 october.
A new seismic interpolation and denoising method with a curvelet transform matching filter, employing the fast iterative shrinkage thresholding algorithm fista, is proposed. Application of curvelet denoising to 2d and 3d seismic data. Seismic signal enhancement using 2d fast discrete curvelet. Implementation of the uniform discrete curvelet transform. Transform can be made arbitrarily tight at the cost of oversampling. The secondgeneration curvelet transform 11, 12, has been shown to be a very ecient tool for many di. Formerly, in practice, we usually split 3d seismic data into 2d data in inline, crossline and horizontal directions, and applied the process in sequence with the 2d curvelet transform. For the application of the curvelet transform in seismic denoising. Multiresolution analysis mra using wavelet, ridgelet, and curvelet transforms has been used in the proposed segmentation system. Seismic denoising using curvelet analysis sciencedirect. Curvelet provides a sparser representation of the seismic signal, and it has been successfully used in seismic data recovery, denoising, attenuating random and coherent noise 21, 22, 23, 24.
A sparsitypromoting approach article pdf available in geophysics 756 november 2010 with 187 reads how we measure reads. Curvelet transforms and filtering of seismic attributes for reservoir modeling. Software engineering, king fahad university of petroleum and minerals, dhahran, 2004 a thesis submitted in partial fulfilment of the requirements for the degree of master of science in the faculty of graduate studies geophysics the university of. The specialist will recognize the connection with ideas from microlocal analysis.
Parts of this research were supported by the national science foundation, the alfred p. Seismic or seismoelectric wave fronts can be optimally. We present a new nonparametric transformbased reconstruction method that exploits the compression of seismic data b the recently developed curvelet. Curvelet domain separation for removing multiples from noisy seismic data involves transforming seismic data into the curvelet domain and a process for simultaneously separating multiples and primaries in the curvelet domain. Traditional curvelet transform cannot provide an adaptive basis function to. No accessseg technical program expanded abstracts 2006. Fourier analysis, wavelet transform, inverse theory, spatial analysis. To store a complex valued cartesian grid of size 5123with doubleprecision, we need 2gb of memory.
After the curvelet representation of the data is computed 350, selected geophysical data attributes and their interdependencies are extracted 355, from which geological features. Multiresolution analysis using wavelet, ridgelet, and. The curvelet transform cand e s and donoho, 2004 is designed to represent curvelike singularities optimally by decomposing seismic data into a superposition of localized plane waves, called curvelets. Incoherent noise suppression and deconvolution using curvelet. The aliasfree curvelet coefficients are upscaled to estimate a mask function that is used to constrain the. The mother curvelet define a position vector with as translation parameters. Fast discrete curvelet transforms multiscale modeling. Secondly, we combine the second generation nfdct with the standard fast discrete curvelet transform fdct to form a new curveletbased method, coined nonequispaced curvelet reconstruction with sparsitypromoting inversion ncrsi for the. The first digital transformation is based on unequally spaced fast fourier transforms, while the second is based on the wrapping of specially selected fourier samples.
Seismic noise attenuation using curvelet transform and dip map. This package is an implementation of the uniform discrete curvelet transform as described in uniform discrete curvelet transform by truong t. Application of curvelet denoising to 2d and 3d seismic. A curvelet is a new and effective spectral transform, that allows sparse representations of complex data. Curveletbased seismoelectric data processing archive ouverte hal. This procedure is accomplished applying the fast discrete curvelet transform, as described in. It is a highdimensional generalization of the wavelet transform designed to represent images at different scales and different orientations angles. The curvelet transform is a higher dimensional generalization of the wavelet transform designed to represent images at different scales and different angles. This paper describes two digital implementations of a new mathematical transform, namely, the second generation curvelet transform in two and three dimensions. Incoherent noise present in seismic reflection data corrupts the quality of the signal and can often lead to misinterpretation.
Accelerating seismic interpolation with a gradient. A seismic interpolation and denoising method with curvelet. The aforementioned coefficients are divided into two groups of scales. Geometric seismic attributes capture and enhance spatial features such as faults and channels. Mar 11, 2020 this software was written by darren thomson and gilles hennenfent from slim group at university of british columbia. The inverse curvelet transform is then applied and, after taking the real part, gives the image below. Noise suppression based on a fast discrete curvelet transform. Nonequispaced curvelet transform for seismic data reconstruction. After that, i have used hard thresholding to remove low frequency components and then applied inverse curvelet transform. In the field of seismic exploration, ground roll seriously affects the deep. The secondgeneration curvelet transform 1012 has been shown to be a very efficient tool for many different applications in image processing, seismic data exploration, fluid mechanics, and. Accelerating seismic interpolation with a gradient projection method based on tight frame property of curvelet jingjie cao1,4 yanfei wang2 benfeng wang3 1shijiazhuang university of economics, shijiazhuang, hebei 050031, china. Sometimes these attributes may be hidden within specific spectral bands, and may require extraction by spectral decomposition.
The two implementations essentially differ by the choice of spatial. The curvelet transform is a higher dimensional generalization of the wavelet transform. Mar 21, 2012 hi i want to decomposition multi focus images to cartoon and texture content. Seismic data processing with the parallel windowed.
The curvelet domain lends itself particularly well for denoising because coherent seismic energy maps to a relatively small number of significant curvelet coefficents. The most negative values are in white and the most positive values in black. This process consists, roughly speaking, of obtaining inner products in the fourier domain. Processing seismic laboratory for imaging and modeling. Curvelet transforms and filtering of seismic attributes for reservoir modeling matz haugen and tapan mukerji. Seismic data processing with the parallel windowed curvelet transform by fadhel alhashim b. A curvelet transform differs from other directional wavelet transforms in that the degree of localisation in orientation varies with scale.
Incoherent noise suppression and deconvolution using. Verschuur in their 2008 article entitled adaptive curveletdomain primarymultiple separation, published in geophysics. The university of british columbia vancouver geophysics seismic data processing with the parallel windowed curvelet transform by fadhel alhashim b. The final solution to the seismic data regularization problem, is obtained by applying the weightedinverse curvelet transform to, that is. Curvelet transforms and filtering of seismic attributes for. As mentioned, curvelet transform has its origins in image processing and thus the input data is in a form of an image, i. Nonparametric seismic data recovery with curvelet frames. The latter form, known as the fdct method, has become a useful tool in image processing and signal analysis.
Moreover, the curvelet transform is redundant in general, and the 3d transform described in section 4 has a redundancy factor of about 5. Curvelet transform are widely used in seismic exploration due to its. The curvelet transform is a fairly recent image processing technique that is able to easily approximate curves present in images. We explore curvelet transforms and filtering of coherence attributes based. It is impossible for most computers to store data of this size in the memory. Continuous curvelet transform the continuous curvelet transform has gone through two major revisions. The curvelet transform is a multiscale and multidimensional transform candes et al. Since the wrapping version is faster and invertible up to numerical precision, while the usfft version is only approximately invertible, we use only the wrapping version throughout this paper.
A ground roll is an interferential wave existing widely in seismic data, and it is an important step in seismic data processing to suppress ground roll. This spectral technique is based on directional basis functions that represent objects having discontinuities along smooth curves. Performance evaluation of wavelet, ridgelet, curvelet and. The total number of matrixvector multiplications required by this method is similar to those required by iterativereweighted leastsquares gersztenkorn et al. The quality of seismic data restricts the effective application of. Method for identifying geologic features, such as hydrocarbon indicators, from geophysical data, such as seismic data, by taking a curvelet transform of the data. The curvelet transform for image denoising, ieee transaction on image processing, 11, 6, 2002. Seismic data processing with the parallel windowed curvelet. The irregularities in the acquisition have to be dealt with beforehand resampling of the data, interpolation or filling the empty traces with zeros. Scale and directionguided interpolation of aliased. The curvelet transform naturally exploits the highdimensional and strong geometrical structure of seismic data.
The paper fast discrete curvelet transforms explains the curvelet transforms in detail. Highfidelity adaptive curvelet domain primarymultiple. The sister webpage for the wave atom transform is at. Then go to and register to sign a licensing agreement and download the curvelab 2. Scale and directionguided interpolation of aliased seismic. The discrete curvelet transform was introduced by in two forms, the wrapping version and the unequally spaced fft usfft version. We present a new nonparametric transform based reconstruction method that exploits the compression of seismic data b the recently developed curvelet. Sampling and reconstruction of seismic wavefields in the. The groups software platform for seismic inversion, in python. Software engineering, king fahad university of petroleum and minerals, dhahran, 2004 a thesis submitted in partial fulfilment of the requirements for the degree of master of science in c fadhel alhashim 2009 october, 2009 the faculty of.
Pdf nonequispaced curvelet transform for seismic data. Seismic denoising using curvelet analysis request pdf. The curvelet transform the curvelet transform can be seen as a combination of reversible transformations. Curveletbased separation methods have the advantage of minimizing the damage to the primary events based on the compatible nature of the curvelet transform to seismic data, as described by f. Seismic data recovery from data with missing traces on otherwise regular. Our extension of the curvelet transform for nonuniformly sampled data dashed line retains the. It is particularly a challenging task to classify cancers in human organs in scanners output using shape or graylevel information. Curvelet transforms and filtering of seismic attributes. Ridgelet and curvelet first generation toolbox file.
Curvelet transforms have been widely used for seismic process. Seismic imaging in the curvelet domain and its implications for the. The approach treats the matching filter, seismic interpolation, and denoising all as the same inverse problem using an inversion iteration algorithm. Seismic imaging with the generalized radon transform.
1378 1665 1443 327 765 25 221 202 1535 715 1604 1272 139 696 812 1209 245 1286 1131 1198 503 1493 1652 104 605 1326 672 1113 31 70 249 132 983 965 517 1153 1370 1392 111