study of efficiency of seismic time-frequency spectral decomposition by matching pursuit for detecting thin layers

نویسندگان

محمد رداد

دانشجوی دکتری ژئوفیزیک، گروه فیزیک زمین، مؤسسة ژئوفیزیک دانشگاه تهران، ایران عبدالرحیم جواهریان

استاد بازنشسته، گروه فیزیک زمین، مؤسسه ژئوفیزیک، دانشگاه تهران و استاد دانشکده مهندسی نفت دانشگاه صنعتی امیرکبیر، تهران، ایران

چکیده

most geologic changes have a seismic response but sometimes this is expressed only in certain spectral ranges hidden within the broadband data. spectral decomposition is one of the methods which can be utilized to help interpreting such cases. there are several time-frequency methods including: short-time fourier transform (stft), continuous wavelet transform (cwt), wigner-ville distribution (wvd), s-transform, and matching pursuit decomposition (mpd). in this paper, we use the mpd method. this method is newer than the other time-frequency methods used in exploration seismology. mallat and zhang (1993) have improved time and frequency resolution simultaneously by using mpd method. in this method, a signal is decomposed into constructive wavelets. time and frequency properties of wavelets are used locally for spectral decomposition. pursuits are the algorithms which search the best time-frequency matching between the signal and a linear combination of selected wavelets from wavelet dictionary. matching pursuit which is an iterative procedure optimizes signal estimation by each new wavelet chosen from a dictionary. these wavelets combined linearly to obtain the best match with the signal. a signal should expand to waveforms which their time-frequency properties could be matched to local structures. such waveforms are called time-frequency atoms. there are many approaches to match wavelets of dictionary to a seismic signal and to obtain time-frequency spectrum in matching pursuit decomposition. the base of all approaches is the mallat and zhang’s algorithm. however, computing time of the original algorithm is very high due to many iterations and that is why particular conditions have been applied in different researches to limit the matching pursuit algorithm for obtaining the lower performing time. in this work, particular conditions are rather similar to wang’s (2007) method. on seismic data, layer thickness is described on the basis of the seismic travel time. when a layer with different properties has a thickness by one-fourth wavelength, top and base reflections will interfere constructively. for thin layers less than tuning thickness, combined seismic amplitude decreases with thickness (when reflection coefficients of the top and of the base are opposite). generally, spectral decomposition has many applications in interpretation of seismic sections and so there will be extra needs to study and develop them. thin layers and many stratigraphic hydrocarbon reservoirs are beneath the threshold of the vertical seismic resolution (tuning thickness) and because of their low thicknesses, they are not resolvable. for this reason, mapping the small-scale geological structures is one of the important interpretational cases. when the thickness of a thin layer decreases pick frequency slightly increases. in this work, this issue has been used to detect thin layers by time-frequency spectrum and by single-frequency sections obtained from mpd. in this paper, we investigated the performance of the matching pursuit decomposition for time-frequency analysis of seismic sections to delineate and detect thin layers on synthetic data (including simple thin layer model and also wedge model) and real data. it is observed that the interpretation of thin layers is simpler by single-frequency sections. it is shown that for a simple thin layer if considerable frequency in single-frequency sections increases, ability in resolving layer interfaces would be increased. in the wedge model, as the frequency increases resolution threshold of layer interfaces moves to a lower thickness and therefore it would be possible to detect lower thickness layers. the tuning thickness has been decreased from 19 meters in original seismic section to 12 meters in 80 hz frequency section. in the real data, it is shown that when a thin layer is not resolvable in a seismic section it might be detected using the mpd method. in this case, by providing 20, 40, 60, 80 and 100 hz single-frequency sections when high frequency sections are studied, interfaces of thin layer are appeared gradually. it is concluded that time-frequency sections are useful instruments to detect and delineate thin layers.

برای دانلود باید عضویت طلایی داشته باشید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Basis Pursuit for Seismic Spectral decomposition

Spectral decomposition is a powerful analysis tool used to identify the frequency content of seismic data. Many spectral decomposition techniques have been developed, each with their own advantages and disadvantages. The basis pursuit technique produces a high time frequency resolution map through formulating the problem as an inversion scheme. This techniques differs from conventional spectral...

متن کامل

study of cohesive devices in the textbook of english for the students of apsychology by rastegarpour

this study investigates the cohesive devices used in the textbook of english for the students of psychology. the research questions and hypotheses in the present study are based on what frequency and distribution of grammatical and lexical cohesive devices are. then, to answer the questions all grammatical and lexical cohesive devices in reading comprehension passages from 6 units of 21units th...

Matching Pursuit With Time - Frequency

We introduce an algorithm, called matching pursuit, that decomposes any signal into a linear expansion of waveforms that are selected from a redundant dictionary of functions. These waveforms are chosen in order to best match the signal structures. Matching pursuits are general procedures to compute adaptive signal representations. With a dictionary of Gabor functions , a matching pursuit deene...

متن کامل

Sparse Time-Frequency decomposition by adaptive basis pursuit

In our recent paper, we proposed a data-driven time-frequency analysis method. We also presented convincing numerical evidence to demonstrate the effectiveness of our method. Convergence analysis has been also carried out for periodic signals that satisfy certain scale separation property. In this paper, we propose an improved time-frequency analysis method that can be used to decompose signals...

متن کامل

imaging thin beds using attributes achieved in spectral decomposition by short time fourier transforms

one of the most fundamental reservoir characteristics is the thickness. the analysis of thin bed tuning on seismic reflectivity has been studied extensively by wides (1973) and neidel and pogiaggliomi (1977), who discussed the limits of seismic resolution. during the past decade, the industry has developed a plethora of new attributes in studying thin beds by employing spectral decomposition (p...

متن کامل

Empirical mode decomposition for seismic time-frequency analysis

Time-frequency analysis plays a significant role in seismic data processing and interpretation. Complete ensemble empirical mode decomposition decomposes a seismic signal into a sum of oscillatory components, with guaranteed positive and smoothly varying instantaneous frequencies. Analysis on synthetic and real data demonstrates that this method promises higher spectral-spatial resolution than ...

متن کامل

منابع من

با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید


عنوان ژورنال:
فیزیک زمین و فضا

جلد ۳۸، شماره ۱، صفحات ۱۱۳-۱۳۱

کلمات کلیدی
[ ' m o s t g e o l o g i c c h a n g e s h a v e a s e i s m i c r e s p o n s e b u t s o m e t i m e s t h i s i s e x p r e s s e d o n l y i n c e r t a i n s p e c t r a l r a n g e s h i d d e n w i t h i n t h e b r o a d b a n d d a t a . s p e c t r a l d e c o m p o s i t i o n i s o n e o f t h e m e t h o d s w h i c h c a n b e u t i l i z e d t o h e l p i n t e r p r e t i n g s u c h c a s e s . t h e r e a r e s e v e r a l t i m e ' , ' f r e q u e n c y m e t h o d s i n c l u d i n g : s h o r t ' , ' t i m e f o u r i e r t r a n s f o r m ( s t f t ) ' , ' c o n t i n u o u s w a v e l e t t r a n s f o r m ( c w t ) ' , ' w i g n e r ' , ' v i l l e d i s t r i b u t i o n ( w v d ) ' , ' s ' , ' t r a n s f o r m ' , ' a n d m a t c h i n g p u r s u i t d e c o m p o s i t i o n ( m p d ) . i n t h i s p a p e r ' , ' w e u s e t h e m p d m e t h o d . t h i s m e t h o d i s n e w e r t h a n t h e o t h e r t i m e ' , ' f r e q u e n c y m e t h o d s u s e d i n e x p l o r a t i o n s e i s m o l o g y . m a l l a t a n d z h a n g ( 1 9 9 3 ) h a v e i m p r o v e d t i m e a n d f r e q u e n c y r e s o l u t i o n s i m u l t a n e o u s l y b y u s i n g m p d m e t h o d . i n t h i s m e t h o d ' , ' a s i g n a l i s d e c o m p o s e d i n t o c o n s t r u c t i v e w a v e l e t s . t i m e a n d f r e q u e n c y p r o p e r t i e s o f w a v e l e t s a r e u s e d l o c a l l y f o r s p e c t r a l d e c o m p o s i t i o n . p u r s u i t s a r e t h e a l g o r i t h m s w h i c h s e a r c h t h e b e s t t i m e ' , ' f r e q u e n c y m a t c h i n g b e t w e e n t h e s i g n a l a n d a l i n e a r c o m b i n a t i o n o f s e l e c t e d w a v e l e t s f r o m w a v e l e t d i c t i o n a r y . m a t c h i n g p u r s u i t w h i c h i s a n i t e r a t i v e p r o c e d u r e o p t i m i z e s s i g n a l e s t i m a t i o n b y e a c h n e w w a v e l e t c h o s e n f r o m a d i c t i o n a r y . t h e s e w a v e l e t s c o m b i n e d l i n e a r l y t o o b t a i n t h e b e s t m a t c h w i t h t h e s i g n a l . a s i g n a l s h o u l d e x p a n d t o w a v e f o r m s w h i c h t h e i r t i m e ' , ' f r e q u e n c y p r o p e r t i e s c o u l d b e m a t c h e d t o l o c a l s t r u c t u r e s . s u c h w a v e f o r m s a r e c a l l e d t i m e ' , ' f r e q u e n c y a t o m s . t h e r e a r e m a n y a p p r o a c h e s t o m a t c h w a v e l e t s o f d i c t i o n a r y t o a s e i s m i c s i g n a l a n d t o o b t a i n t i m e ' , ' f r e q u e n c y s p e c t r u m i n m a t c h i n g p u r s u i t d e c o m p o s i t i o n . t h e b a s e o f a l l a p p r o a c h e s i s t h e m a l l a t a n d z h a n g ’ s a l g o r i t h m . h o w e v e r ' , ' c o m p u t i n g t i m e o f t h e o r i g i n a l a l g o r i t h m i s v e r y h i g h d u e t o m a n y i t e r a t i o n s a n d t h a t i s w h y p a r t i c u l a r c o n d i t i o n s h a v e b e e n a p p l i e d i n d i f f e r e n t r e s e a r c h e s t o l i m i t t h e m a t c h i n g p u r s u i t a l g o r i t h m f o r o b t a i n i n g t h e l o w e r p e r f o r m i n g t i m e . i n t h i s w o r k ' , ' p a r t i c u l a r c o n d i t i o n s a r e r a t h e r s i m i l a r t o w a n g ’ s ( 2 0 0 7 ) m e t h o d . r n o n s e i s m i c d a t a ' , ' l a y e r t h i c k n e s s i s d e s c r i b e d o n t h e b a s i s o f t h e s e i s m i c t r a v e l t i m e . w h e n a l a y e r w i t h d i f f e r e n t p r o p e r t i e s h a s a t h i c k n e s s b y o n e ' , ' f o u r t h w a v e l e n g t h ' , ' t o p a n d b a s e r e f l e c t i o n s w i l l i n t e r f e r e c o n s t r u c t i v e l y . f o r t h i n l a y e r s l e s s t h a n t u n i n g t h i c k n e s s ' , ' c o m b i n e d s e i s m i c a m p l i t u d e d e c r e a s e s w i t h t h i c k n e s s ( w h e n r e f l e c t i o n c o e f f i c i e n t s o f t h e t o p a n d o f t h e b a s e a r e o p p o s i t e ) . g e n e r a l l y ' , ' s p e c t r a l d e c o m p o s i t i o n h a s m a n y a p p l i c a t i o n s i n i n t e r p r e t a t i o n o f s e i s m i c s e c t i o n s a n d s o t h e r e w i l l b e e x t r a n e e d s t o s t u d y a n d d e v e l o p t h e m . t h i n l a y e r s a n d m a n y s t r a t i g r a p h i c h y d r o c a r b o n r e s e r v o i r s a r e b e n e a t h t h e t h r e s h o l d o f t h e v e r t i c a l s e i s m i c r e s o l u t i o n ( t u n i n g t h i c k n e s s ) a n d b e c a u s e o f t h e i r l o w t h i c k n e s s e s ' , ' t h e y a r e n o t r e s o l v a b l e . f o r t h i s r e a s o n ' , ' m a p p i n g t h e s m a l l ' , ' s c a l e g e o l o g i c a l s t r u c t u r e s i s o n e o f t h e i m p o r t a n t i n t e r p r e t a t i o n a l c a s e s . w h e n t h e t h i c k n e s s o f a t h i n l a y e r d e c r e a s e s p i c k f r e q u e n c y s l i g h t l y i n c r e a s e s . i n t h i s w o r k ' , ' t h i s i s s u e h a s b e e n u s e d t o d e t e c t t h i n l a y e r s b y t i m e ' , ' f r e q u e n c y s p e c t r u m a n d b y s i n g l e ' , ' f r e q u e n c y s e c t i o n s o b t a i n e d f r o m m p d . r n i n t h i s p a p e r ' , ' w e i n v e s t i g a t e d t h e p e r f o r m a n c e o f t h e m a t c h i n g p u r s u i t d e c o m p o s i t i o n f o r t i m e ' , ' f r e q u e n c y a n a l y s i s o f s e i s m i c s e c t i o n s t o d e l i n e a t e a n d d e t e c t t h i n l a y e r s o n s y n t h e t i c d a t a ( i n c l u d i n g s i m p l e t h i n l a y e r m o d e l a n d a l s o w e d g e m o d e l ) a n d r e a l d a t a . i t i s o b s e r v e d t h a t t h e i n t e r p r e t a t i o n o f t h i n l a y e r s i s s i m p l e r b y s i n g l e ' , ' f r e q u e n c y s e c t i o n s . i t i s s h o w n t h a t f o r a s i m p l e t h i n l a y e r i f c o n s i d e r a b l e f r e q u e n c y i n s i n g l e ' , ' f r e q u e n c y s e c t i o n s i n c r e a s e s ' , ' a b i l i t y i n r e s o l v i n g l a y e r i n t e r f a c e s w o u l d b e i n c r e a s e d . i n t h e w e d g e m o d e l ' , ' a s t h e f r e q u e n c y i n c r e a s e s r e s o l u t i o n t h r e s h o l d o f l a y e r i n t e r f a c e s m o v e s t o a l o w e r t h i c k n e s s a n d t h e r e f o r e i t w o u l d b e p o s s i b l e t o d e t e c t l o w e r t h i c k n e s s l a y e r s . t h e t u n i n g t h i c k n e s s h a s b e e n d e c r e a s e d f r o m 1 9 m e t e r s i n o r i g i n a l s e i s m i c s e c t i o n t o 1 2 m e t e r s i n 8 0 h z f r e q u e n c y s e c t i o n . i n t h e r e a l d a t a ' , ' i t i s s h o w n t h a t w h e n a t h i n l a y e r i s n o t r e s o l v a b l e i n a s e i s m i c s e c t i o n i t m i g h t b e d e t e c t e d u s i n g t h e m p d m e t h o d . i n t h i s c a s e ' , ' b y p r o v i d i n g 2 0 ' , 4 0 , 6 0 , ' 8 0 a n d 1 0 0 h z s i n g l e ' , ' f r e q u e n c y s e c t i o n s w h e n h i g h f r e q u e n c y s e c t i o n s a r e s t u d i e d ' , ' i n t e r f a c e s o f t h i n l a y e r a r e a p p e a r e d g r a d u a l l y . i t i s c o n c l u d e d t h a t t i m e ' , ' f r e q u e n c y s e c t i o n s a r e u s e f u l i n s t r u m e n t s t o d e t e c t a n d d e l i n e a t e t h i n l a y e r s . ' ]

میزبانی شده توسط پلتفرم ابری doprax.com

copyright © 2015-2023