improving seismic facies analysis using wtmmla attributes, self-organizing maps and k-mean clustering

نویسندگان

سعید هادیلو

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

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

مشاور، شرکت مشاوران تهران

چکیده

reservoir models are initially generated from estimates of specific rock properties and maps of reservoir heterogeneity. many types of information are used in reservoir model construction. one of the most important sources of information comes from wells, including well logs and core samples. unfortunately well log and core data are local measurements that may not reflect the reservoir behavior as a whole. in addition, well data are not available at the initial phases of exploration. in contrast to sparse well data, 3d seismic data cover large areas. seismic attributes extracted from 3d seismicdata can provide information for the construction of reservoir models. seismic facies analysis can be accomplished through the use of pattern recognition techniques. when the geological information is incomplete or nonexistent, seismic facies analysis can be done using unsupervised learning techniques. one of the most promising mathematical techniques of unsupervised learning is the kohonen's self organizing map (som) (kohonen, 2001). in this paper we use the som and time-frequency analysis to characterize reservoirs. since variations in frequency content are sensitive to subtle changes in reflective information. in this context, we show that the wavelet transform modulus maxima line amplitudes (wtmmla) that extracted from continuous wavelet transforms (cwt) can be applied to detect singularities. these singularities are analyzed and clustered by som. the som networks map points of input space to points in an output space while preserving the topology. topology preservation means that points which are close in the input space should also be close in the output space (map). normally, the input space is of high dimension while the output is two-dimensional. the seismic attributes, can be represented by vectors in the space rn, x = [x1,x2...xn].we assume that the map has pelements; therefore, there will exist p prototype vectors mi, mi= [mi1. . . min], i = 1, 2, . . . ,p, where n is the dimension of the input vector. after the som training, prototype vectors represent the input data set of seismic attributes, the distances between x and all the prototype vectors are computed. the mapunit with the smallest distance mbto the input vector x is called the best matching unit (bmu) and is computed by, . the prototype vector corresponding to the bmu and their neighbors are moved towards the input winner vector in the input space. since one of the main objectives of this work was the identification of data clusters, we displayed the distances between the neighbor prototype vectors to identify similarities among the vector prototypes. we used the u-matrix (ultsch, 1993), to represent these distances. after the som learning, the u-matrix was generated by computing, for each som prototype vector, the distance between the neighbor prototype vectors and their average. for estimation of the number of existing seismic facies in the data, we used a k-means partitive clustering algorithm. we clustered the prototype vectors instead of the original data. in this manner, large data sets formed by the som prototype vectors can be indirectly grouped. results showed that the proposed method not only provides a better understanding about the group formations, but it is also computationally efficient. another benefit of this methodology is noise reduction because the prototype vectors represent local averages of the original data without any loss of resolution. to automate the classification process, we used the davies and bouldin (1979) index (dbi) as means of evaluating the results of the k-means partitioning. transitions, or irregular structures, present in any kind of signals carry information related to its physical phenomena (mallat, 1999). besides the horizon locations, the identified transition characterization in the interpretation is associated with geological processes. in this way, a possible transition classification could be linked to the seismic facies. detection of transitions or singularities in signals is based on simple mathematical concepts. the signal inflection points are associated with the first-derivative extremes which correspond to the second-derivative zero crossings. for the signal inflection-point positions, using the cwt local peak locations, a wavelet should be chosen as the first-derivative of the smoothing function . one of the wavelet functions that fulfill this requirement is the first-derivative of the gaussian function, called the gauss wavelet. we can extract scalogram's local peaks coincide from the signal inflection points. it can be proven that these lines, which are called wtmmla, can be used to characterize the signal irregularity. the signal irregularities can be characterized mathematically through the wtmmla and h?lder exponent (mallat, 1999). the exponent can be obtained from the slope estimation of the curve created by the log2 of the wtmmla coefficients divided by the log2 of the scales. in this study we used wtmmla as a direct seismic attribute. we calculated cwt coefficients and wtmmla for sixteen seismic data samples around the picked reservoir horizon. the extracted wtmmla can show the possible heterogeneity and singularity within the reservoir. we used these attribute as input vector for the som step and obtained the u_matrix. the k-mean and dbi estimate the number of seismic facies.utilizing of cwt to locate events in time through the identification of signal singularities also proved to be useful as an appropriate tool for detection of seismic events. therefor this method proved to be less sensitive to interpretation errors. the performance of the method was tested on kangan formation at one of the iranian oil fields.

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

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

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

منابع مشابه

using game theory techniques in self-organizing maps training

شبکه خود سازمانده پرکاربردترین شبکه عصبی برای انجام خوشه بندی و کوانتیزه نمودن برداری است. از زمان معرفی این شبکه تاکنون، از این روش در مسائل مختلف در حوزه های گوناگون استفاده و توسعه ها و بهبودهای متعددی برای آن ارائه شده است. شبکه خودسازمانده از تعدادی سلول برای تخمین تابع توزیع الگوهای ورودی در فضای چندبعدی استفاده می کند. احتمال وجود سلول مرده مشکلی اساسی در الگوریتم شبکه خودسازمانده به حسا...

Self-organizing Maps as Substitutes for K-Means Clustering

One of the most widely used clustering techniques used in GISc problems is the k-means algorithm. One of the most important issues in the correct use of k-means is the initialization procedure that ultimately determines which part of the solution space will be searched. In this paper we briefly review different initialization procedures, and propose Kohonen’s SelfOrganizing Maps as the most con...

متن کامل

Self-Organizing Maps in Seismic Image Segmentation

Unsupervised neural networks such as the Kohonen Self-Organizing Maps (SOM) have been widely used for searching natural clusters in multidimensional and massive data. One example where the data available for analysis can be extremely large is seismic interpretation for hydrocarbon exploration. In order to assist the interpreter in identifying characteristics of interest confined in the seismic ...

متن کامل

K-Dynamical Self Organizing Maps

Neural maps are a very popular class of unsupervised neural networks that project high-dimensional data of the input space onto a neuron position in a low-dimensional output space grid. It is desirable that the projection effectively preserves the structure of the data. In this paper we present a hybrid model called K-Dynamical Self Organizing Maps (KDSOM ) consisting of K Self Organizing Maps ...

متن کامل

Determining an Optimal Seismic Network Configuration Using Self-Organizing Maps

The Seismic Research Centre, University of the West Indies operates a seismic network that performs suboptimally in detecting, locating, and correctly determining the magnitude of earthquakes due to a diverse constituion of seismometers and the utilization of a site selection process that approximates an educated guess. My work seeks to apply Self-Organizing Maps (SOM) to arrive at the optimal ...

متن کامل

Gene Clustering Using Self-Organizing Maps and Particle Swarm Optimization

Gene clustering, the process of grouping related genes in the same cluster, is at the foundation of different genomic studies that aim at analyzing the function of genes. Microarray technologies have made it possible to measure gene expression levels for thousand of genes simultaneously. For knowledge to be extracted from the datasets generated by these technologies, the datasets have to be pre...

متن کامل

منابع من

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


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

جلد ۳۸، شماره ۲، صفحات ۴۵-۵۶

کلمات کلیدی
[ ' r e s e r v o i r m o d e l s a r e i n i t i a l l y g e n e r a t e d f r o m e s t i m a t e s o f s p e c i f i c r o c k p r o p e r t i e s a n d m a p s o f r e s e r v o i r h e t e r o g e n e i t y . m a n y t y p e s o f i n f o r m a t i o n a r e u s e d i n r e s e r v o i r m o d e l c o n s t r u c t i o n . o n e o f t h e m o s t i m p o r t a n t s o u r c e s o f i n f o r m a t i o n c o m e s f r o m w e l l s ' , ' i n c l u d i n g w e l l l o g s a n d c o r e s a m p l e s . u n f o r t u n a t e l y w e l l l o g a n d c o r e d a t a a r e l o c a l m e a s u r e m e n t s t h a t m a y n o t r e f l e c t t h e r e s e r v o i r b e h a v i o r a s a w h o l e . i n a d d i t i o n ' , ' w e l l d a t a a r e n o t a v a i l a b l e a t t h e i n i t i a l p h a s e s o f e x p l o r a t i o n . i n c o n t r a s t t o s p a r s e w e l l d a t a ' , ' 3 d s e i s m i c d a t a c o v e r l a r g e a r e a s . s e i s m i c a t t r i b u t e s e x t r a c t e d f r o m 3 d s e i s m i c d a t a c a n p r o v i d e i n f o r m a t i o n f o r t h e c o n s t r u c t i o n o f r e s e r v o i r m o d e l s . s e i s m i c f a c i e s a n a l y s i s c a n b e a c c o m p l i s h e d t h r o u g h t h e u s e o f p a t t e r n r e c o g n i t i o n t e c h n i q u e s . w h e n t h e g e o l o g i c a l i n f o r m a t i o n i s i n c o m p l e t e o r n o n e x i s t e n t ' , ' s e i s m i c f a c i e s a n a l y s i s c a n b e d o n e u s i n g u n s u p e r v i s e d l e a r n i n g t e c h n i q u e s . o n e o f t h e m o s t p r o m i s i n g m a t h e m a t i c a l t e c h n i q u e s o f u n s u p e r v i s e d l e a r n i n g i s t h e k o h o n e n & a p o s ; s s e l f o r g a n i z i n g m a p ( s o m ) ( k o h o n e n ' , ' 2 0 0 1 ) . r n i n t h i s p a p e r w e u s e t h e s o m a n d t i m e ' , ' f r e q u e n c y a n a l y s i s t o c h a r a c t e r i z e r e s e r v o i r s . s i n c e v a r i a t i o n s i n f r e q u e n c y c o n t e n t a r e s e n s i t i v e t o s u b t l e c h a n g e s i n r e f l e c t i v e i n f o r m a t i o n . i n t h i s c o n t e x t ' , ' w e s h o w t h a t t h e w a v e l e t t r a n s f o r m m o d u l u s m a x i m a l i n e a m p l i t u d e s ( w t m m l a ) t h a t e x t r a c t e d f r o m 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 s ( c w t ) c a n b e a p p l i e d t o d e t e c t s i n g u l a r i t i e s . t h e s e s i n g u l a r i t i e s a r e a n a l y z e d a n d c l u s t e r e d b y s o m . r n t h e s o m n e t w o r k s m a p p o i n t s o f i n p u t s p a c e t o p o i n t s i n a n o u t p u t s p a c e w h i l e p r e s e r v i n g t h e t o p o l o g y . t o p o l o g y p r e s e r v a t i o n m e a n s t h a t p o i n t s w h i c h a r e c l o s e i n t h e i n p u t s p a c e s h o u l d a l s o b e c l o s e i n t h e o u t p u t s p a c e ( m a p ) . n o r m a l l y ' , ' t h e i n p u t s p a c e i s o f h i g h d i m e n s i o n w h i l e t h e o u t p u t i s t w o ' , ' d i m e n s i o n a l . t h e s e i s m i c a t t r i b u t e s ' , ' c a n b e r e p r e s e n t e d b y v e c t o r s i n t h e s p a c e r n ' , ' x = [ x 1 ' , ' x 2 . . . x n ] . w e a s s u m e t h a t t h e m a p h a s p e l e m e n t s ; t h e r e f o r e ' , ' t h e r e w i l l e x i s t p p r o t o t y p e v e c t o r s m i ' , ' m i = [ m i 1 . . . m i n ] ' , ' i = 1 ' , 2 , ' . . . ' , ' p ' , ' w h e r e n i s t h e d i m e n s i o n o f t h e i n p u t v e c t o r . a f t e r t h e s o m t r a i n i n g ' , ' p r o t o t y p e v e c t o r s r e p r e s e n t t h e i n p u t d a t a s e t o f s e i s m i c a t t r i b u t e s ' , ' t h e d i s t a n c e s b e t w e e n x a n d a l l t h e p r o t o t y p e v e c t o r s a r e c o m p u t e d . t h e m a p u n i t w i t h t h e s m a l l e s t d i s t a n c e m b t o t h e i n p u t v e c t o r x i s c a l l e d t h e b e s t m a t c h i n g u n i t ( b m u ) a n d i s c o m p u t e d b y ' , ' . t h e p r o t o t y p e v e c t o r c o r r e s p o n d i n g t o t h e b m u a n d t h e i r n e i g h b o r s a r e m o v e d t o w a r d s t h e i n p u t w i n n e r v e c t o r i n t h e i n p u t s p a c e . s i n c e o n e o f t h e m a i n o b j e c t i v e s o f t h i s w o r k w a s t h e i d e n t i f i c a t i o n o f d a t a c l u s t e r s ' , ' w e d i s p l a y e d t h e d i s t a n c e s b e t w e e n t h e n e i g h b o r p r o t o t y p e v e c t o r s t o i d e n t i f y s i m i l a r i t i e s a m o n g t h e v e c t o r p r o t o t y p e s . w e u s e d t h e u ' , ' m a t r i x ( u l t s c h ' , ' 1 9 9 3 ) ' , ' t o r e p r e s e n t t h e s e d i s t a n c e s . a f t e r t h e s o m l e a r n i n g ' , ' t h e u ' , ' m a t r i x w a s g e n e r a t e d b y c o m p u t i n g ' , ' f o r e a c h s o m p r o t o t y p e v e c t o r ' , ' t h e d i s t a n c e b e t w e e n t h e n e i g h b o r p r o t o t y p e v e c t o r s a n d t h e i r a v e r a g e . r n f o r e s t i m a t i o n o f t h e n u m b e r o f e x i s t i n g s e i s m i c f a c i e s i n t h e d a t a ' , ' w e u s e d a k ' , ' m e a n s p a r t i t i v e c l u s t e r i n g a l g o r i t h m . w e c l u s t e r e d t h e p r o t o t y p e v e c t o r s i n s t e a d o f t h e o r i g i n a l d a t a . i n t h i s m a n n e r ' , ' l a r g e d a t a s e t s f o r m e d b y t h e s o m p r o t o t y p e v e c t o r s c a n b e i n d i r e c t l y g r o u p e d . r e s u l t s s h o w e d t h a t t h e p r o p o s e d m e t h o d n o t o n l y p r o v i d e s a b e t t e r u n d e r s t a n d i n g a b o u t t h e g r o u p f o r m a t i o n s ' , ' b u t i t i s a l s o c o m p u t a t i o n a l l y e f f i c i e n t . a n o t h e r b e n e f i t o f t h i s m e t h o d o l o g y i s n o i s e r e d u c t i o n b e c a u s e t h e p r o t o t y p e v e c t o r s r e p r e s e n t l o c a l a v e r a g e s o f t h e o r i g i n a l d a t a w i t h o u t a n y l o s s o f r e s o l u t i o n . t o a u t o m a t e t h e c l a s s i f i c a t i o n p r o c e s s ' , ' w e u s e d t h e d a v i e s a n d b o u l d i n ( 1 9 7 9 ) i n d e x ( d b i ) a s m e a n s o f e v a l u a t i n g t h e r e s u l t s o f t h e k ' , ' m e a n s p a r t i t i o n i n g . r n t r a n s i t i o n s ' , ' o r i r r e g u l a r s t r u c t u r e s ' , ' p r e s e n t i n a n y k i n d o f s i g n a l s c a r r y i n f o r m a t i o n r e l a t e d t o i t s p h y s i c a l p h e n o m e n a ( m a l l a t ' , ' 1 9 9 9 ) . b e s i d e s t h e h o r i z o n l o c a t i o n s ' , ' t h e i d e n t i f i e d t r a n s i t i o n c h a r a c t e r i z a t i o n i n t h e i n t e r p r e t a t i o n i s a s s o c i a t e d w i t h g e o l o g i c a l p r o c e s s e s . i n t h i s w a y ' , ' a p o s s i b l e t r a n s i t i o n c l a s s i f i c a t i o n c o u l d b e l i n k e d t o t h e s e i s m i c f a c i e s . d e t e c t i o n o f t r a n s i t i o n s o r s i n g u l a r i t i e s i n s i g n a l s i s b a s e d o n s i m p l e m a t h e m a t i c a l c o n c e p t s . t h e s i g n a l i n f l e c t i o n p o i n t s a r e a s s o c i a t e d w i t h t h e f i r s t ' , ' d e r i v a t i v e e x t r e m e s w h i c h c o r r e s p o n d t o t h e s e c o n d ' , ' d e r i v a t i v e z e r o c r o s s i n g s . f o r t h e s i g n a l i n f l e c t i o n ' , ' p o i n t p o s i t i o n s ' , ' u s i n g t h e c w t l o c a l p e a k l o c a t i o n s ' , ' a w a v e l e t s h o u l d b e c h o s e n a s t h e f i r s t ' , ' d e r i v a t i v e o f t h e s m o o t h i n g f u n c t i o n . o n e o f t h e w a v e l e t f u n c t i o n s t h a t f u l f i l l t h i s r e q u i r e m e n t i s t h e f i r s t ' , ' d e r i v a t i v e o f t h e g a u s s i a n f u n c t i o n ' , ' c a l l e d t h e g a u s s w a v e l e t . w e c a n e x t r a c t s c a l o g r a m & a p o s ; s l o c a l p e a k s c o i n c i d e f r o m t h e s i g n a l i n f l e c t i o n p o i n t s . i t c a n b e p r o v e n t h a t t h e s e l i n e s ' , ' w h i c h a r e c a l l e d w t m m l a ' , ' c a n b e u s e d t o c h a r a c t e r i z e t h e s i g n a l i r r e g u l a r i t y . t h e s i g n a l i r r e g u l a r i t i e s c a n b e c h a r a c t e r i z e d m a t h e m a t i c a l l y t h r o u g h t h e w t m m l a a n d h ? l d e r e x p o n e n t ( m a l l a t ' , ' 1 9 9 9 ) . t h e e x p o n e n t c a n b e o b t a i n e d f r o m t h e s l o p e e s t i m a t i o n o f t h e c u r v e c r e a t e d b y t h e l o g 2 o f t h e w t m m l a c o e f f i c i e n t s d i v i d e d b y t h e l o g 2 o f t h e s c a l e s . i n t h i s s t u d y w e u s e d w t m m l a a s a d i r e c t s e i s m i c a t t r i b u t e . w e c a l c u l a t e d c w t c o e f f i c i e n t s a n d w t m m l a f o r s i x t e e n s e i s m i c d a t a s a m p l e s a r o u n d t h e p i c k e d r e s e r v o i r h o r i z o n . t h e e x t r a c t e d w t m m l a c a n s h o w t h e p o s s i b l e h e t e r o g e n e i t y a n d s i n g u l a r i t y w i t h i n t h e r e s e r v o i r . w e u s e d t h e s e a t t r i b u t e a s i n p u t v e c t o r f o r t h e s o m s t e p a n d o b t a i n e d t h e u _ m a t r i x . t h e k ' , ' m e a n a n d d b i e s t i m a t e t h e n u m b e r o f s e i s m i c f a c i e s . u t i l i z i n g o f c w t t o l o c a t e e v e n t s i n t i m e t h r o u g h t h e i d e n t i f i c a t i o n o f s i g n a l s i n g u l a r i t i e s a l s o p r o v e d t o b e u s e f u l a s a n a p p r o p r i a t e t o o l f o r d e t e c t i o n o f s e i s m i c e v e n t s . t h e r e f o r t h i s m e t h o d p r o v e d t o b e l e s s s e n s i t i v e t o i n t e r p r e t a t i o n e r r o r s . t h e p e r f o r m a n c e o f t h e m e t h o d w a s t e s t e d o n k a n g a n f o r m a t i o n a t o n e o f t h e i r a n i a n o i l f i e l d s . ' ]

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

copyright © 2015-2023