Pdf processing seismic data with open source software. Deconvolution is performed along the time axis to increase vertical resolution by. Nonstationary deconvolution methods perform deconvolution of the seismic data that also compensates for waves attenuation. Image processing comparison of deconvolution software. In seismology we create an impulse seismic energysource and this signal after interacting with the. Seismic unix su package is a free software developed and maintained by the. Although deconvolution usually is applied to prestack data trace by trace, it is not uncommon to design a single deconvolution operator and apply it to all the traces on a shot record. Processing, inversion and reconstruction of seismic data. 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. 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. Here is a page from one of my previous courses describing convolution. Predictive deconvolution in seismic data processing in. Cryosphere geophysics research group hydraulic tomography.
Gabor deconvolution easily beats wiener deconvolution on a nonstationary synthetic. The image below is a single slice taken from a stack before and after deconvolution using these plugins. Seismically, rock layers are defined by the densities and velocities with which seismic waves propagate through them. Laboratory manual for seismic data processing courses at. In section 2, we describe the signal model and brie. Reflections on the deconvolution of land seismic data. 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. This course covers practical aspects of signal theory and inverse problems with application to seismic data processing. Seismic 2d reflection processing and interpretation of. 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. 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. Convolution of two series, ui, and wi, denoted uw, is. Set your browser font size in the range 12 to 14pt for optimal display.
An interactive java application demonstrates the method and can be used to explore the impact of different parameter choices. Mathematics let subcommands deconvolution dataplot reference manual march 19, 1997 337 deconvolution purpose compute the discrete deconvolution of two variables. Estimating hydrogeologic parameters from radar data charles t. The seismic unix project is supported by these organizations. Using fourier mathematics, therefore, greatly decreases the computational cpu expense required to perform the filtering operation. Deconvolution attempts to compress the wavelet, thereby increasing the resolution of the seismic data. Multidimensional seismic data processing wikipedia. Seismic unix has programs for performing many of the operations of geophysical signal processing. Appreciation is expressed to the author, to alistar r. Each wiggle trace is the record, in time, of when sound from each source. The deconvolution demo uses simple synthetic spike traces to illustrate both dereverberation and spiking decon using supef and other tools. 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. Lines, volume editor society of exploration geophysicists the international society of applied geophysics. To increase the vertical resolution of data by compressing the source wavelet to a zerostrike if possible which is known as spiking deconvolution.
This website beta version contains information on geophysical methods, references to geophysical citations, and a glossary of geophysical terms related to environmental applications. Deconvolution is process designed to enhance the vertical. Deconvolution is basically the opposite of convolution. Center for wave phenomena colorado school of mines golden, co 80401, usa. 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. Deconvolution is carried out in seismic unix with a driver routine modified from benz, 1999, using weiner predictive filtering, e.
Introduction a software developed by wolfgang stefan in cooperation with rosie renaut and ed garnero to deblur seismograms using total variation regularized deconvolution. Center for geophysical investigation of the shallow. 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. To understand deconvolution, first we need to examine the constituent elements of a recorded seismic trace the convolutional model. Reflections on the deconvolution of land seismic data cseg. 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. Deconvolution is a filtering process which removes a wavelet from the recorded seismic trace by reversing the process of convolution. 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. Deconvolution tutorial national radio astronomy observatory. 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. Deconvolution compresses the basic wavelet in the recorded seismogram and attenuates reverberations and shortperiod multiples. Convolution and deconvolution university of california. The iterative deconvolution 3d plugin uses a psf image zstack to correct the image contrast vs.
Cryosphere geophysics research group hydraulic tomography self potential electrical. 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. Laboratory manual for seismic data processing courses at kfupm using the seismic unx. Seismic unix su package is a free software developed and maintained by the center for wave phenomena cwp at colorado school of mines. Sparse spikes deconvolution is one of the oldest inverse problems, which is a stylized version of recovery in seismic imaging. Pdf deconvolution of seismic data using phase change operators. Deconvolution may be applied to greyscale images or true color images. Introduction to deconvolution and image preparation.
Basis pursuit inversion figures 1b and 2b show the synthetic traces, which are a result of convolution between the wavelet and the re. Seismic unix is an open source seismic utilities package which was supported by the center. We have presented a nonstationary deconvolution algorithm based on the gabor transform that effectively corrects seismic data for source signature and attenuation. The complete software package is available electronically 6. 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. Principles of 3d seismic interpretation and applications. Convolution general formulae the resulting lake levels can be calculated by the following procedure, called convolution. Innovative petrotech solutions, inc principles of 3d seismic interpretation and applications instructor. Given the input stimulus functions, and the measured fmri signal data, program 3ddeconvolve. Our example file is in seismic unix format thus the.
Modelbased deconvolution for dominant, thinbed seismic. The use of l1 penalty in similar problems promotes sparsity of the solution rm 1 see chen et al. In the example above seismic unix will create 100 seismograms in 100 different. Notes for working with color images are included at the end of this section. 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.
The general form of the matrix equation for a filter of length n is represented in equation 1, 6. For example, a convolution can be used to model the filtering of seismic energy by the various rock layers in the earth. Note that if we need to multiply two polynomials, with coefficients uk and wk, we would use exactly the formula above. Convolution and filtering seismology and geophysics. 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. Intermediateadvanced course description this course delivers techniques related to practical analysis and interpretation of 3d seismic. Pdf open source seismic processing softwares provides a low cost alternative to. 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. Since geophysical data obtained in such techniques is a function of. Deconvolution of seismic data using phase change operators. The best known example of such kind is seismic unix, a free reflection processing system. Seismic reflection methods environmental geophysics us epa.
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. 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. View and download powerpoint presentations on deconvolution seismic ppt. Currently construction of the underground research facility, onkalo, is going on at the olkiluoto site. Seismic deconvolution and inversion with erratic data. Among these methods are the gabor deconvolution algorithm pro posed by margrave et al. Hence, it increases resolution and yields a more interpretable seismic section. The spiking deconvolution in seismic data processing is routinely applied to compress the source wavelet included in the seismic traces to improve temporal resolution. Introduction deconvolution is a fundamental aspect of seismic signal processing. The mathematical form of the convolution of two functions, a filter ft and a timeseries xt, is. The deconvolution demo uses simple synthetic spike traces to. Seismic processing objective processing systems general concept of cmp processing processing tools segy and similar file formats general features of promax reading.
Threedimensional sparse seismic deconvolution based on. Deconvolution is a process universally applied to seismic data, but is one that is mysterious to many geoscientists. Center for wave phenomena colorado school of mines golden, co 80401, usa past support was received from. 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. Program 3ddeconvolve was developed to provide deconvolution analysis of fmri time series data. Deconvolution is the reverse process of convolution. Time seriesdata processing and analysis math 587geop.
Deconvolution techniques used in conventional processing are based on optimum wiener filtering. We performed several deconvolution tests on different kinds of datasets. 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. Predictive deconvolution is the use of information from earlier part of a seismic trace to predict and deconvolve the. Seismic data processing with seismic unx course notes.
By far the most important application is predictive deconvolution in which a repeating signal e. 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. Bowen bullard laboratories, department of earth sciences, university of cambridge, uk. Multidimensional seismic data processing forms a major component of seismic profiling, a technique used in geophysical exploration. 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. It is also to be noted that not all geological information can be estimated. Here you can find a list of the programs that the su package has, with a brief description. Highresolution seismic data deconvolution by a0 algorithm.
The deconvolution procedure here was adopted from 1 but di erent minimizing criteria and a di erent generation. 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. Multichannel blind seismic deconvolution using dynamic. Note that this process can take around three to five days, depending on the level of effort put into the manual analysis stages. Builtin core functions such as seismic phase traveltime calculation and multichannel cross correlation significantly. We have used the varimax norm as the objective function to measure the simplicity of the seismic trace after.
Both su and the linux os require working in a commanddriven. A new iterative procedure for deconvolution of seismic ground. As for convolution, the reflection coefficient rt is convolved with seismic wavelet wt to obtain the seismic trace. This document uses html 3 tables, and inline images. 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. 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. James brown, and zhihong nancy cao abstract three basic methods for suppressing multiples exist in published literature.
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. 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. 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. 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. In year 2001 olkiluoto was selected for the site of final disposal. Eq 332 in practice, h is typically a data stream while g is a response function. However, it is subject to change at any time and does not represent a commitment on the part of improvision.
The commonest way to perform deconvolution is to design a wiener filter to transform one wavelet into another wavelet in a leastsquares sense. Deconvolution can also be performed on other types of data, such as gravity, magnetic or well log data. Introduction to deconvolution and image preparation version 3. Convolution gives the overlapp of one function with another function reversed in time.
Index terms data fusion, deconvolution, distributed. On real data, gabor deconvolution is very competitive with the combination of wiener deconvolution and. The abcs of seismic exploration and processing enders a. Because of the power and flexibility of promax, it is unreasonable to attempt to cover all possible features and applications in this manual. Click to view article as pdf gc what is deconvolution. Laboratory manual for seismic data processing courses at kfupm. Robinson and sven treitcl geophysical references series no.
Seismic deconvolution using sparse spike inversion vs. Unix translation, header load, qc, shotrecord velocity filter, deconvolution, cdp gather. Convolution in coordinate space the image is equivalent to multiplication in frequency space. Seismic unix user manual free ebook download as pdf file. The response function window 2, top right must be known and is usually either calculated on the basis. The objective of the deconvolution procedure is to remove the obscuring effect of the wavelets replica making up the seismic trace and consequently. The technique itself has various applications, including mapping ocean floors, determining the structure of sediments, mapping subsurface currents and hydrocarbon exploration. The society of exploration geophysicists the gas research institute. Pdf seismic data deconvolution using kalman filter based on. Seismic unix user manual command line interface parameter. The value of using seismic unix is that it is free and it contains many common seismic waveform processing.
969 505 1292 1193 1256 60 159 513 1513 1115 1295 531 237 1075 729 694 1356 468 331 428 908 1247 190 418 103 296 319 992 929 1047 959 1362 832 505 1141 1063 1410 617 1105 912 1317 885 1117 1333 19 836 795 647 159