Nonstationary deconvolution methods perform deconvolution of the seismic data that also compensates for waves attenuation. Time seriesdata processing and analysis math 587geop 505 brian borchers and rick aster november 8, 20 notes on deconvolution we have seen how to perform convolution of discrete and continuous signals. The complete software package is available electronically 6. Image processing comparison of deconvolution software. Deconvolution is performed along the time axis to increase vertical resolution by. Convolution of two series, ui, and wi, denoted uw, is.
Seismically, rock layers are defined by the densities and velocities with which seismic waves propagate through them. The deconvolution demo uses simple synthetic spike traces to illustrate both dereverberation and spiking decon using supef and other tools. The general form of the matrix equation for a filter of length n is represented in equation 1, 6. Center for wave phenomena colorado school of mines golden, co 80401, usa. Center for geophysical investigation of the shallow. Seismic unix user manual free ebook download as pdf file. Processing, inversion and reconstruction of seismic data. Deconvolution is carried out in seismic unix with a driver routine modified from benz, 1999, using weiner predictive filtering, e. Seismic data consists of a series of recorded wiggle traces that describe a set of echoes from interfaces between rock layers in the subsurface that have different rock properties. Landmark promax vsp user training manual xv preface about the manual this two volume manual is intended to accompany the instruction given during the standard promax seismic processing and analysis course. Hence, it increases resolution and yields a more interpretable seismic section.
Estimating hydrogeologic parameters from radar data charles t. Threedimensional sparse seismic deconvolution based on. Seismic unix su package is a free software developed and maintained by the center for wave phenomena cwp at colorado school of mines. Modelbased deconvolution for dominant, thinbed seismic. Mathematics let subcommands deconvolution dataplot reference manual march 19, 1997 337 deconvolution purpose compute the discrete deconvolution of two variables. Predictive deconvolution in seismic data processing in. This course covers practical aspects of signal theory and inverse problems with application to seismic data processing. Principles of 3d seismic interpretation and applications. Introduction deconvolution is a fundamental aspect of seismic signal processing. Reflection seismology or seismic reflection is a method of exploration geophysics that uses the principles of seismology to estimate the properties of the earths subsurface from reflected seismic waves. James brown, and zhihong nancy cao abstract three basic methods for suppressing multiples exist in published literature. The use of l1 penalty in similar problems promotes sparsity of the solution rm 1 see chen et al. In particular, the course stresses regularization methods for inverse problems that arise in the inversion of seismic data, noise elimination and reconstruction of seismic surveys. Software package for passive seismic data preprocessing by chunquan yu, yingcai zheng, and xuefeng shang abstract we introduce an opensource matlab software package, named crazyseismic, for passive seismic data preprocessing.
Deconvolution is basically the opposite of convolution. Deconvolution, or inverse filtering, can improve seismic data that were adversely affected by filtering, or convolution that occurs naturally as seismic energy is filtered by the earth. Marine processing part 6 predictive deconvolution this sequence of blog posts will build up into a complete description of a 2d marine processing sequence and how it is derived. Gabor deconvolution easily beats wiener deconvolution on a nonstationary synthetic. Predictive deconvolution remove shortperiod multiples most notably from relatively flat, shallow waterbottom the periodicity of the multiples is exploited to design a filter that removes the predictable part of the wavelet multiples, leaving only its nonpredictable part signal zerooffset gather after prediction decon.
Similar to the conventional kalman filter model for seismic deconvolution, our new kalman model is also based on the common viewpoint that a reflected seismic record can be regarded as a. Deconvolution techniques used in conventional processing are based on optimum wiener filtering. To understand deconvolution, first we need to examine the constituent elements of a recorded seismic trace the convolutional model. An interactive java application demonstrates the method and can be used to explore the impact of different parameter choices.
The method requires a controlled seismic source of energy, such as dynamite or tovex blast, a specialized air gun or a seismic vibrator, commonly known by the trademark. We have presented a nonstationary deconvolution algorithm based on the gabor transform that effectively corrects seismic data for source signature and attenuation. Seismic 2d reflection processing and interpretation of shallow refraction data abstract posiva oy takes care of the final disposal of spent nuclear fuel in finland. The mathematical form of the convolution of two functions, a filter ft and a timeseries xt, is. Seismic unix su package is a free software developed and maintained by the. Deconvolution is a mathematical process that allows the adjustment of the amplitude and frequency contents of a seismic ground motion applied at the base of the foundation in order to get the desired output at the damfoundation interface. The deconvolution procedure here was adopted from 1 but di erent minimizing criteria and a di erent generation. Index terms data fusion, deconvolution, distributed.
By far the most important application is predictive deconvolution in which a repeating signal e. Note that if we need to multiply two polynomials, with coefficients uk and wk, we would use exactly the formula above. Deconvolution can also be performed on other types of data, such as gravity, magnetic or well log data. Deconvolution is process designed to enhance the vertical. Reflections on the deconvolution of land seismic data. Eq 332 in practice, h is typically a data stream while g is a response function. Cryosphere geophysics research group hydraulic tomography. Seismic deconvolution using sparse spike inversion vs. The deconvolution demo uses simple synthetic spike traces to. Appreciation is expressed to the author, to alistar r.
Deconvolution of seismic data using phase change operators. Laboratory manual for seismic data processing courses at kfupm using the seismic unx. Deconvolution is a process universally applied to seismic data, but is one that is mysterious to many geoscientists. Seismic unix is an open source seismic utilities package which was supported by the center. Intermediateadvanced course description this course delivers techniques related to practical analysis and interpretation of 3d seismic. Convolution and deconvolution university of california. Set your browser font size in the range 12 to 14pt for optimal display. Because of the power and flexibility of promax, it is unreasonable to attempt to cover all possible features and applications in this manual. Builtin core functions such as seismic phase traveltime calculation and multichannel cross correlation significantly. It is possible to use it to manipulate and create your own seismograms, and also to convert them between the su standard file and the industry standard seg y. Convolution gives the overlapp of one function with another function reversed in time. Here you can find a list of the programs that the su package has, with a brief description. Predictive deconvolution in seismic data processing in atala prospect of rivers state, nigeria egbai, j.
Multichannel blind seismic deconvolution using dynamic. In the example above seismic unix will create 100 seismograms in 100 different. A step in seismic signal processing to recover high frequencies, attenuate multiples, equalize amplitudes, produce a zerophase wavelet or for other purposes that generally affect the waveshape. Bowen bullard laboratories, department of earth sciences, university of cambridge, uk. The commonest way to perform deconvolution is to design a wiener filter to transform one wavelet into another wavelet in a leastsquares sense. Seismic reflection methods environmental geophysics us epa. Highresolution seismic data deconvolution by a0 algorithm. We have used the varimax norm as the objective function to measure the simplicity of the seismic trace after.
Our example file is in seismic unix format thus the. Unix translation, header load, qc, shotrecord velocity filter, deconvolution, cdp gather. Sparse spikes deconvolution is one of the oldest inverse problems, which is a stylized version of recovery in seismic imaging. The image below is a single slice taken from a stack before and after deconvolution using these plugins. Deconvolution attempts to compress the wavelet, thereby increasing the resolution of the seismic data. Laboratory manual for seismic data processing courses at kfupm. Multidimensional seismic data processing forms a major component of seismic profiling, a technique used in geophysical exploration. Convolution in coordinate space the image is equivalent to multiplication in frequency space. Seismic data processing with seismic unx course notes. Deconvolution tutorial national radio astronomy observatory. It is also to be noted that not all geological information can be estimated. Convolution and filtering seismology and geophysics. The objective of the deconvolution procedure is to remove the obscuring effect of the wavelets replica making up the seismic trace and consequently.
The best known example of such kind is seismic unix, a free reflection processing system. On real data, gabor deconvolution is very competitive with the combination of wiener deconvolution and. Reflections on the deconvolution of land seismic data cseg. Introduction to deconvolution and image preparation. This document uses html 3 tables, and inline images. Find powerpoint presentations and slides using the power of, find free presentations research about deconvolution seismic ppt.
As for convolution, the reflection coefficient rt is convolved with seismic wavelet wt to obtain the seismic trace. The value of using seismic unix is that it is free and it contains many common seismic waveform processing. The boise hydrogeophysical research site, an in site field laboratory for testing geophysical methods to determine the 3dimensional distribution of permeability in heterogeneous alluvial aquifer systems. Seismic unix has programs for performing many of the operations of geophysical signal processing.
Predictive deconvolution is the use of information from earlier part of a seismic trace to predict and deconvolve the. Early in our careers as geophysicists, most of us took at least one course on seismic signal analysis where we were taught that standard wiener deconvolution converts the minimumphase source wavelet in our seismic data to a wavelet with a phase spectrum that is zero and an amplitude spectrum that. Another goal of deconvolution is to produce a wavelet with a simple phase character, ideally a zerophase wavelet, which is the same for every trace in the seismic dataset. To increase the vertical resolution of data by compressing the source wavelet to a zerostrike if possible which is known as spiking deconvolution. Pdf deconvolution of seismic data using phase change operators. Digital filtering theory applies to virtually any sampled information in time e. Both su and the linux os require working in a commanddriven. If you have little or no practical experience of seismic data processing, it is strongly recommended that you work through all of the exercises from start to finish. Program 3ddeconvolve was developed to provide deconvolution analysis of fmri time series data. Given the input stimulus functions, and the measured fmri signal data, program 3ddeconvolve. Pdf open source seismic processing softwares provides a low cost alternative to. Notes for working with color images are included at the end of this section. Time seriesdata processing and analysis math 587geop 505.
Currently construction of the underground research facility, onkalo, is going on at the olkiluoto site. We performed several deconvolution tests on different kinds of datasets. In section 2, we describe the signal model and brie. Fourier deconvolution is used here to remove the distorting influence of an exponential tailing response function from a recorded signal window 1, top left that is the result of an unavoidable rc lowpass filter action in the electronics. After eight years of the first distribution of the users manual of the seismic unix su package, it is yet a challenge to execute some seismic processing su scripts, especially for inexperienced.
Lines, volume editor society of exploration geophysicists the international society of applied geophysics. Pdf seismic data deconvolution using kalman filter based on. In seismology we create an impulse seismic energysource and this signal after interacting with the. For example, a convolution can be used to model the filtering of seismic energy by the various rock layers in the earth. Among these methods are the gabor deconvolution algorithm pro posed by margrave et al. Several different deconvolution techniques have been developed for different purposes. The seismic unix project is partially funded by the society of exploration geophysicists seg, and by the center for wave phenomena cwp, department of geophysical engineering.
Deconvolution compresses the basic wavelet in the recorded seismogram and attenuates reverberations and shortperiod multiples. The response function window 2, top right must be known and is usually either calculated on the basis. Aug, 20 marine processing part 6 predictive deconvolution this sequence of blog posts will build up into a complete description of a 2d marine processing sequence and how it is derived. Innovative petrotech solutions, inc principles of 3d seismic interpretation and applications instructor. Deconvolution is the reverse process of convolution. However, it is subject to change at any time and does not represent a commitment on the part of improvision. Note that this process can take around three to five days, depending on the level of effort put into the manual analysis stages. Convolution general formulae the resulting lake levels can be calculated by the following procedure, called convolution. Seismic 2d reflection processing and interpretation of. Seismic deconvolution and inversion with erratic data. Deconvolution is a filtering process which removes a wavelet from the recorded seismic trace by reversing the process of convolution.
Each wiggle trace is the record, in time, of when sound from each source. Using fourier mathematics, therefore, greatly decreases the computational cpu expense required to perform the filtering operation. A new iterative procedure for deconvolution of seismic ground. Renaut june 27, 2006 abstract this note provides complete details of the numerical algorithm implemented in the software package for signal resortoration of seismic signals 5. The abcs of seismic exploration and processing enders a. Robinson and sven treitcl geophysical references series no. One of the most common deconvolution techniques is spiking deconvolution, which tries to remove source wavelets from the recorded data by converting them into a spike. The society of exploration geophysicists the gas research institute.
The iterative deconvolution 3d plugin uses a psf image zstack to correct the image contrast vs. Click to view article as pdf gc what is deconvolution. No processing sequence is definitive and techniques vary with time and software, however the idea is to provide a practical guide for applying seismic processing. Here is a page from one of my previous courses describing convolution.
Introduction to deconvolution and image preparation version 3. In year 2001 olkiluoto was selected for the site of final disposal. The earth is composed of layers of rocks with different lithology and physical properties. Center for wave phenomena colorado school of mines golden, co 80401, usa past support was received from.
Seismic data processing geos 469569 spring 2006 geos 469569 is a mix of digital filtering theory and practical applications of digital techniques to assemble and enhance images of subsurface geology. Early in our careers as geophysicists, most of us took at least one course on seismic signal analysis where we were taught that standard wiener deconvolution converts the minimumphase source wavelet in our seismic data to a wavelet with a phase spectrum that is zero and an amplitude spectrum that is broad and flat. Basis pursuit inversion figures 1b and 2b show the synthetic traces, which are a result of convolution between the wavelet and the re. Comparison of deconvolution software a user point of view part 2 deconvolution is an image processing technique that restores the effective object representation 3 4, allowing to improve images analysis steps such as segmentation 1 or colocalization study 2. Seismic unix user manual command line interface parameter. The seismic unix project is supported by these organizations. Multidimensional seismic data processing wikipedia. The technique itself has various applications, including mapping ocean floors, determining the structure of sediments, mapping subsurface currents and hydrocarbon exploration. Seismic processing objective processing systems general concept of cmp processing processing tools segy and similar file formats general features of promax reading. Printing history edition date edition 1 february 1998 edition 2 august 2000 improvision believes the information in this manual to be correct at the date of issue. Time seriesdata processing and analysis math 587geop. Pdf processing seismic data with open source software. Laboratory manual for seismic data processing courses at. Deconvolution may be applied to greyscale images or true color images.