Judging from the amount of time that seismic processors put into phase matching intersecting 2d seismic lines, and overlapping 3d seismic surveys, it appears to be simply a statement of the truth. This article is based on crains seismic petrophysics by e. 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 in both the time domain and with the. Time seriesdata processing and analysis math 587geop 505. Seismic data processing to interpret subsurface features is both computationally and data intensive. Seismic processing often involves several stages of deconvolution, each of a different type and with a different objective.
A process designed to restore a waveshape to the form it had before it underwent a linear filtering action convolution. Pdf seismic data deconvolution using kalman filter based. Predictive deconvolution in seismic data processing in atala prospect of rivers state, nigeria. Young work completed under doe st485d saiidia national laboratories received sponsored by doe aug 1 6 19% osti to support the signal processing and data visualization needs of ctbt related projects at snl, a matlab based gui was developed. Blind seismic deconvolution is the recovery of reflectivity sequence from the seismic records, when the seismic wavelet is unknown. The comprehensive and physically complete linear imaging foundation developed presents new results at the leading edge of seismic processing for target location and identification. Their directivity is due to the use of arrays of source elements and to the free. Compressive sensing in seismic signal processing is a construction of the unknown reflectivity sequence from the incoherent measurements of the seismic records. Digital filtering theory applies to virtually any sampled information in time e.
Seismic data interpretation using digital image processing. Stuart kozola is a product manager at mathworks and focuses on matlab and addon products for data analysis, mathematical modeling, and computational finance. This book serves as a fundamental guide to seismic imaging principles and algorithms and their foundation in inverse scattering theory and is a valuable resource for. In this paper, a seismic blind deconvolution algorithm based on bayesian. Feb 23, 2011 introduce gpu computing for speeding up solutions of the wave equation for seismic analysis. Convolution gives the overlapp of one function with another function reversed in time. The demo also shows how to speed up the solution of the wave equation finite difference pde using a custom cuda kernel. Common procedures to streamline seismic data processing include. 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. Read the readme file to locate the public data sources on the internet. In seismic processing the purpose of deconvolution is to remove or collapse the seismic wavelet. An investigation of the tools of seismic data processing 2 in the frequency domain, the convolution is equivalent to multiplication i. This short book is for students, professors and professionals interested in signal processing of seismic data using matlab tm. An investigation of the tools of seismic data processing.
Working with data files, such as segy, that are too large to fit in system memory. Usually an averaging of 8 makes the best images, if the dye is strong enough. He also received segs awards for best paper in geophysics in 1992 imaging salt with seismic turning waves and best paper presented at the annual meeting in 2002 atomic meshing of seismic images. Deconvolution deconvolution is the reverse of convolution. In seismology we create an impulse seismic energysource and this signal after interacting with the. Durham etheses directional seismic source signature. 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 in both the time domain and with the help of the fourier transform. Deconvolution is basically the opposite of convolution. An ideal seismic source would be a spike maximum amplitude at every frequency, but the best practical one would be minimum phase. Bancroft, thais guirigay, and helen isaac, crewes university of calgary, alberta, canada summary the inversion process to recover rock properties is typically approximated with seismic migration that is a transpose process. The commonest way to perform deconvolution is to design a wiener filter to transform one wavelet into another wavelet in a leastsquares sense. Includes readytorun matlab algorithms for most of the techniques presented. Seismic deconvolution using sparse spike inversion vs. Since wiener deconvolution tries to estimate stationary wavelets from noisy, nonstationary seismic data, we should not be too surprised to see some.
V31v42 january 2016 with 545 reads how we measure reads. Recent progress in numerical methods and computer science allows us today to simulate the propagation of seismic waves through realistically heterogeneous earth models with unprecedented accuracy. Seismic inversion is now commonly used on poststack and prestack seismic data for estimating rock properties used in reservoir characterization. Note that if we need to multiply two polynomials, with coefficients uk and wk, we would use exactly the formula above. In the lecture portion of the course gpgn452561 now gpgn4615 61 advanced seismic methodsseismic processing the student is given a word, picture, and chalkboard introduction of the process of seismic data acquisition and the application of a myriad of processing steps for converting raw seismic data into a scientifically useful picture of. It is quite common to convert a given wavelet source wavelet into its minimum phase equivalent since several processing stages e. While seismic migration is aimed at imaging the reflectors or the interfaces at their correct subsurface locations, seismic inversion attempts to estimate elastic and flow properties of the layers. In this way, the deconvolved seismic record provides an estimate of the reflectivity. Convolution and filtering seismology and geophysics. 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. An illustration of the learning curve for effective use of the dap architecture. Mar 18, 2016 download deconvolution an intuitive image deconvolution and deblur application that can apply multiple different algorithms in order to improve the quality of your pictures. Deconvolution techniques used in conventional processing are based on optimum wiener filtering.
Enhancing a seismic image after migration using deconvolution. Prior to joining mathworks in 2006, stuart worked at pratt. Deghosting benefits traditional seismic processing and is a prerequisite to all inversescatteringseries based processing. Optics cleaning and objective correction collar alignment tip. 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. The seismic convolutional model states that a seismic record is the convolution of the earths reflectivity with the seismic wavelet.
Seismic processing and imaging gain deeper insight with the industrys most accurate subsurface imaging tools. For example, a convolution can be used to model the filtering of seismic energy by the various rock layers in the earth. Deconvolution is a filtering process which removes a wavelet from the recorded seismic trace by reversing the process of convolution. A seismic blind deconvolution algorithm based on bayesian. Introduce gpu computing for speeding up solutions of the wave equation for seismic analysis. The objective of the deconvolution procedure is to remove the obscuring effect of the wavelets replica making up the seismic trace and therefore obtain an esti. Deconvolution compresses the basic wavelet in the recorded seismogram and attenuates reverberations and shortperiod multiples. Doherty society of exploration geophysicists, 2001 seismic prospecting. It removes the basic seismic wavelet the source time function modified by various effects of the earth and recording system from the recorded seismic trace and thereby increases temporal resolution. Convolution general formulae the resulting lake levels can be calculated by the following procedure, called convolution.
Basis pursuit inversion figures 1b and 2b show the synthetic traces, which are a result of convolution between the wavelet and the re. Full modeling software 2d, 3d, and 4d seismic survey design, processing, and modeling are created using industry and proprietary software packages. This webpage version is the ed intellectual property of the author. This is the only way to have usable images for processing never saturate an image averaging noise on image, statistics tip. Seismic exploration of the three major processing part of the deconvolution, is the seismic exploration of the students undergraduate, graduate students must master a program, very classic sponsored links. Convolution of two series, ui, and wi, denoted uw, is. Earthquakes and explosions generate depth phases such as pp and sp as reflections of the underlying p signal generated by the event. Purchase optimal seismic deconvolution 1st edition.
Deconvolution is a process universally applied to seismic data, but is one that is mysterious to many geoscientists. The quality of modern seismic data owes a great deal to the success of deconvolution. Jan 01, 2016 deconvolution is basically the opposite of convolution. Convolution of a wavelet with the reflection coefficients. The mathematical form of the convolution of two functions, a filter ft and a timeseries xt, is. Dave received the virgil kauffman gold medal from the society of exploration geophysics for his work on dipmoveout processing of seismic data. Figures 1cf and 2cf show the results of each of the. Conduct geospatial and seismic analysis mathworks consulting. The objective of deconvolution is to nullify objectionable effects of an earlier filter action and thus improve the recognizability and resolution of.
James brown, and zhihong nancy cao abstract three basic methods for suppressing multiples exist in published literature. Processing, inversion, and interpretation of seismic data, issue 10, volume 2 ozdogan yilmaz, stephen m. Enhancing a seismic image after migration using deconvolution john c. Processing of seismic reflection data using matlab. Download deconvolution an intuitive image deconvolution and deblur application that can apply multiple different algorithms in order to improve the quality of your pictures. Nonparametric deconvolution of seismic depth phases. Seismic data deconvolution using kalman filter based on a new system model article pdf available in geophysics 811. Large data in matlab a seismic data processing case study. The complete software package is available electronically 6. The deconvolution of seismic data as a fluctuation analysis. The book includes essential research and teaching material for digital signal and image processing individuals interested in learning seismic data interpretation from the point of.
Time seriesdata processing and analysis math 587geop. 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. Bioxtas raw is a program for analysis of smallangle xray scattering saxs data. Seismic exploration deconvolution free open source codes. Accurate determination of the source depth of a seismic event is a potentially important goal for better discrimination between deeper earthquakes and more shallow nuclear tests.
Full waveform tomography is a tomographic technique that takes advantage of numerical solutions of the elastic wave equation. Full seismic waveform modelling and inversion advances in. Pdf seismic data deconvolution using kalman filter based on. Large data in matlab a seismic data processing case study search form the following matlab project contains the source code and matlab examples used for large data in matlab a seismic data processing case study. Presents image processing concepts that may be readily applied directly to seismic data. Do not copy or distribute in any form without explicit permission. Hence, it increases resolution and yields a more interpretable seismic section. The stepbystep demo of the full reflection seismic data processing workflow using a complete real seismic data set places itself as a very useful feature of the book.
1168 1186 1293 1196 1331 455 868 825 928 1433 431 155 333 915 1036 1034 714 46 302 403 1486 1555 979 917 177 1487 1370 578 1306 1571 1584 1441 859 1348 835 1053 814 1231 1326 1384 908 779 18 733