A comparative study of texture attributes for characterizing subsurface structures in seismic volumes

12/19/2018 ∙ by Zhiling Long, et al. ∙ Georgia Institute of Technology 4

In this paper, we explore how to computationally characterize subsurface geological structures presented in seismic volumes using texture attributes. For this purpose, we conduct a comparative study of typical texture attributes presented in the image processing literature. We focus on spatial attributes in this study and examine them in a new application for seismic interpretation, i.e., seismic volume labeling. For this application, a data volume is automatically segmented into various structures, each assigned with its corresponding label. If the labels are assigned with reasonable accuracy, such volume labeling will help initiate an interpretation process in a more effective manner. Our investigation proves the feasibility of accomplishing this task using texture attributes. Through the study, we also identify advantages and disadvantages associated with each attribute.

READ FULL TEXT VIEW PDF
POST COMMENT

Comments

There are no comments yet.

Authors

page 4

page 8

page 10

page 11

page 15

page 16

page 17

page 18

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

Texture patterns are commonly observed in images of natural scenes Gonzalez and Woods (2006). They have been studied extensively in the image processing literature. In recent years, texture characterization has become an active research area in texture image analysis. In this area, several useful techniques have been developed to extract attributes from a texture pattern that capture the unique spatial distribution of the pixel intensities Haralick et al. (1973); Ojala et al. (2002); Guo et al. (2010); Liu et al. (2012); Zhai et al. (2013); Zhai and Neuhoff (2014); Hu et al. (2016)

. These attributes have been proved effective for various texture image analysis tasks such as classification (categorizing textures into groups of different visual appearances), segmentation (identifying boundaries between different textures in a given image), retrieval (finding images with a texture pattern matching a given image), tracking (following changes of spatial location of a certain texture pattern across different images), etc.

Given that seismic data resembles natural texture images in appearance, we believe that we can analyze them from an image processing perspective. From this perspective, a seismic section can be viewed as a natural image, and the structures contained in the image are visible as areas of certain texture patterns. Then, these patterns can be described and analyzed by the image-based attributes. There were indeed some studies that explored to some extent image-based texture attributes for seismic interpretation. For example, attributes based on the gray level co-occurrence matrix (GLCM) were applied to salt dome detection Gao (2003); Berthelot et al. (2013) and deep-marine facies discrimination Gao (2007); attributes derived from the Hilbert transform were utilized for seismic image segmentation Pitas and Kotropoulos (1992); and Gabor filters were adopted for seismic image segmentation as well Röster and Spann (1998).

However, as successful as these applications were, we believe the exploration has been inadequate considering the tremendous success texture analysis has witnessed in the image processing community. Many advanced techniques developed in recent years have not been examined in the seismic context. More useful applications based on these attributes could be developed to fully exploit their potential. In particular, we believe texture attributes are suitable for serving as descriptors that characterize seismic volumes in terms of the various structures contained therein. Such descriptors are generic, capable of identifying the visible structures all at once, rather than being specific, targeted toward only a certain structure (e.g., salt dome or faults). Such generic descriptors are essential for computationally providing a comprehensive description of the subsurface environment, which can render an initial big picture with spots of possible interest highlighted to expedite the interpretation process.

Therefore, in this paper, we conduct a comparative study examining image-based texture attributes within the context of structure-based characterization of a seismic volume. To fit in the context of the application, we specifically focus on the spatial attributes that belong to the group of local descriptors. The local descriptors capture patterns of variations in visual elements such as intensities and edges in localized scales. The captured patterns are encoded into binary strings, accomplishing robust and computationally efficient texture representations. Such descriptors include the local binary pattern (LBP) Ojala et al. (2002), the completed LBP (CLBP) Guo et al. (2010), the multi-scale CLBP (M-CLBP) Guo et al. (2010), the extended LBP (ELBP) Liu et al. (2012), the completed local derivative pattern (CLDP) Hu et al. (2016), and the local radius index (LRI) Zhai et al. (2013). In addition to the above, we also include in our study two spatial attributes that are more familiar to the seismic interpretation community. One is the classic GLCM, and the other is the semblance, both of which have been used as texture attributes for seismic applications such as salt dome detection Berthelot et al. (2013).

To evaluate the performance of each attribute at characterizing subsurface structures, we adopt a newly-developed framework for interpretation, i.e., seismic volume labeling Alaudah and AlRegib (2016). With labeling, a data volume is automatically segmented into different structures and regions, each assigned with its corresponding label. The segmentation and the label assignment (or classification) are performed based on the extracted texture attributes. Thus, the labeling performance reflects the characterizing capability of the attributes. Figure 1 presents an illustrative example of a manually labeled seismic section. The labels highlight significant structures, which can help pinpoint spots of interest to interpreters to make their exploration more focused and effective. It is worth noting that, the purpose of the labeling is not to accurately delineate the subsurface structures in a given seismic volume. Rather, the labeling is intended to generate an initial map that highlights approximate locations of the structures so that an interpreter can quickly determine where to look further. The labeling simultaneously identifies all structures visible in the data. Techniques customized for a specific structure can then be applied at the will of the interpreter for more focused and refined examination.

The rest of the paper is organized as follows. First, we discuss the various attributes of interest. Then, we introduce the computational framework we use for seismic volume labeling. The attributes will be examined in the labeling experiments afterwards. Finally, we draw conclusions according to our observations in the experiments. The related source code will be available on our website (https://ghassanalregib.com/).

Figure 1: An illustrative example of a manually labeled seismic section

2 Texture Attributes

In this section, we discuss the attributes of interest in this study, divided into two groups. The first group is the traditional attributes including the GLCM-based attributes and the semblance, which are widely used in seismic interpretation. The second group is the local descriptors, which are the newer techniques proposed in the image processing literature in recent years. They are introduced in the order of LBP, LBP variants, and LRI.

2.1 Traditional Attributes

2.1.1 GLCM-based Attributes

The GLCM-based attributes have been widely accepted as useful tools for texture analysis since they were proposed four decades ago Haralick et al. (1973)

. The GLCM is a matrix that describes the co-occurrence pattern between gray levels of two neighboring pixels along a certain direction in an image. In essence, it represents a two-dimensional histogram that approximates the joint probability distribution of the neighboring gray values. It can capture textural patterns for the selected neighborhood along the prescribed direction. For example, high values away from the diagonal in a GLCM reveal sharp changes in gray level, whereas high values close to the diagonal indicate small variations.

Given an 8-bit gray scale image of dimension , where represents the gray value at location . The GLCM (i.e., the co-occurrence matrix of gray values), , is calculated as below:

(1)

where

(2)

, , and represents the number of occurrences of gray level adjacent to gray-level separated by a distance in direction . Based on the GLCM, assuming the total number of gray levels is , the corresponding probability mass function is computed as follows:

(3)

Once is available, various GLCM-based attributes can be generated. For simplicity, we use to represent , and list below the typical attributes:

(4)
(5)
(6)
(7)
(8)

Among these attributes, the GLCM contrast is a measure of the local gray-level variations; the GLCM entropy describes the spatial disorder or complexity in textures; and the GLCM energy measures the pixel pair repetitions, also called texture uniformity or angular second moment. Their values are low for smooth regions, and high for areas with rich texture. On the contrary, both the GLCM homogeneity and the GLCM mutual information

Beghdadi et al. (2015) will show high values for smooth areas and low values for complex textures, with the former being inversely correlated to the GLCM contrast and the latter presenting the dependency among the neighbors.

2.1.2 Semblance

The semblance attribute describes the similarity in a certain spatial neighborhood. Although it is not commonly used for texture image analysis, it has been used for seismic interpretation applications such as detection of faults Wang and AlRegib (2017) and salt domes Berthelot et al. (2013), both of which are structures of interest in this study. In this paper, we adopt the well-known dip-guided semblance developed by Marfurt et al. (1998).

2.2 Local Descriptors

2.2.1 Local Binary Pattern

The local binary pattern (LBP) is a simple and efficient texture attribute, which has become a standard local texture descriptor in the spatial domain Ojala et al. (2002). It describes the intensity difference between a pixel and its local circular neighborhood, denoted by , where defines the number of pixels evenly distributed on the circular neighborhood with radius . To ensure robustness against intensity changes, LBP employs the signs of the differences instead of the exact values to form unique binary codes for the description of local texture patterns.

LBP is calculated as follows:

(9)

where and , represent the intensity of the center pixel and its corresponding neighboring pixels, respectively. Function extracts the sign information of the differences, with the value being for non-negative ones and for negative ones. As the equation shows, LBP encodes the local intensity variation (i.e., the intensity difference between a pixel and its neighbors) into a binary code, resulting in a computationally efficient representation.

The calculation will produce binary patterns. However, a rotation of a texture pattern may lead to different coding results, because is always assigned to an element in a fixed location (e.g., the one to the right of ). To account for the rotation effect, binary patterns with the same circularly shifted code are grouped into one rotation invariant pattern, denoted by . This grouping reduces the total number of possible patterns. As an example, for , the scheme can reduce the total number of patterns from to .

In addition, based on the fact that some binary patterns have higher frequencies of occurrence in texture images than others, Ojala et al. (2002) defined uniform patterns that contain at most two bitwise transitions (i.e., to or to ) when traversed circularly. Such uniform patterns are determined as follows:

(10)

where function counts the bitwise transitions for each rotation invariant pattern, and superscript “” indicates “uniform patterns with rotation invariance.” Introducing the uniform patterns further reduces the number of patterns. When , the number of patterns changes from (for ) to .

2.2.2 LBP Variants

Although LBP is simple and efficient for texture analysis, its performance can be further enhanced by including more local information in addition to the neighbor-center sign information. Following this strategy, a series of LBP variants have been developed. In this study, we select several typical examples among them, including CLBP Guo et al. (2010), M-CLBP Guo et al. (2010), ELBP Liu et al. (2012), and CLDP Hu et al. (2016). As summarized in Table 1, these techniques vary in two aspects: 1) what local information to include in the descriptive components; and 2) how to incorporate such information. To obtain uniform patterns with rotation invariance, the same “” mapping as adopted for LBP is used.

For texture analysis, usually the local descriptors are not examined directly. Instead, histograms are generated from the attributes and used by the algorithms for analysis. In other words, it is the probability distribution of the attributes that helps provide a robust representation of texture patterns. When there are a few descriptive components involved, they are combined to yield either a joint histogram or a concatenated histogram. In Figure 2, histograms obtained for some example subsurface structures using LBP and its variants are given for an illustration.

Techniques Components Notes
LBP Neighbor-center Difference The original form
CLBP
Neighbor-center Difference ()
Neighbor Intensity ()
Center Intensity ()
is identical to the orignal LBP;
intensity of neighbors and center also added
M-CLBP
Combination of CLBP at different scales
(i.e., different values for the radius )
Multiscale implementation of CLBP
for a more comprehensive characterization
ELBP
Neighbor Intensity ()
Center Intensity ()
Radial Difference ()
Ignores the neighbor-center difference;
incorporates cross-scale correlation,
while M-CLBP considers each scale separately
CLDP
Neighbor-center Difference ()
Neighbor Intensity ()
Center Intensity ()
Radial Sign Difference ()
Adds to CLBP cross-scale correlation
(), which is different from
Table 1: Comparison of LBP and Variants
Figure 2: Example histograms calculated using LBP and some variants for typical structures. The left column is for faults, and the right column is for salt dome. The histograms are ordered as (top to bottom): LBP, CLBP, and ELBP.

2.2.3 Local Radius Index

LBP and its variants, as discussed above, all examine the local variation of the pixel intensities. For texture images, edges are also commonly observed. Although LBP-like attributes implicitly capture the edge information, a direct description of the spatial distribution of edges can be more effective. For this purpose, the local radius index (LRI) was proposed by Zhai et al. (2013), which characterizes texture patterns using the local distribution of distances between adjacent edges along a particular angle.

Based on how the local index is computed, there are two variants of LRI: LRI-A and LRI-D. For LRI-A, the inter-edge distance in a given direction is calculated, which represents the width of adjacent smooth regions. In contrast, for LRI-D, the distance is measured from a pixel to its nearest edge, i.e., the boundary of the next smooth region. In this paper, we will only discuss LRI-A, because we did not observe any significant difference in performance between the two in our experiments. An example is shown in Figure 3 illustrating how to compute LRI-A for a small group of pixels. The corresponding procedure is given below.

For pixel and direction , let denote neighboring pixels in direction , where , then

(11)

where we use threshold to define an edge, and to prescribe the size of the texture elements. Consequently, controls the noise sensitivity, while determines the computational complexity. In the volume labeling experiments to be discussed in the next two sections, we used and , where

stands for the standard deviation of the local intensities.

Figure 3: An example illustrating how to compute LRI-A, where black dots refer to edge pixels and white dots are non-edge pixels.

3 Computational Seismic Volume Labeling

To evaluate the capability of the texture attributes for characterizing subsurface structures, we adopt our newly-developed framework for interpretation (i.e., computational seismic volume labeling) Alaudah and AlRegib (2016). The objective is to automatically divide a data volume into segments consisting of various structures and assign each structure with its corresponding label. As demonstrated in Figure 4, this task is accomplished in a setting that combines segmentation, retrieval, and supervised classification. It consists of a training process and a labeling (or testing) process.

Figure 4: An illustration of the framework for computational seismic volume labeling. Alaudah and AlRegib (2016)

Our training process involves three steps. The first step is data extraction, in which, given a small set of manually labeled exemplar patches, a large number of image patches with the same geological structures are automatically extracted to form the training samples. To do this, we use seismic sections extracted along the crossline direction of the widely used Netherlands North Sea Offshore F3 Block dGB Earth Sciences (1987). All sections are normalized to remove contrast and mean variations between sections. Patches containing various subsurface structures are extracted from seismic sections. We define three structures that are of interest, namely, Chaotic layers, Faults, and Salt dome. We also define the Other class for patches that does not contain any of the previous three structures. Examplar patches from each class are shown in Figure 5. Here, the patches are set to a fixed size to make sure the representative texture patterns are captured for each structure of interest, because the labeling is designed to simultaneously identify all structures, not a specific structure. For the same reason, we can only use a square patch, even though rectangular ones may work better for directional structures such as faults.

We extract six exemplar patches (two from Salt dome, two from Other, and one from each of the rest) and manually label them. We then automatically extract 500 image patches from the seismic sections for each type of the defined structures based on their similarity to the exemplar patches, measured according to a recently proposed texture similarity metric Alfarraj et al. (2016). These extracted patches form training samples for a succeeding classification step. The automatic similarity-based patch extraction assists in avoiding the time-consuming manual labeling of the training dataset. Figure 6 presents some examples of the automatically extracted patches, which show excellent consistency comparing to their respective exemplar patches.

(a) (b) (c) (d)
Figure 5: Sample image patches (patch size: ) from the four defined classes of structures Alaudah and AlRegib (2016): (a) Chaotic; (b) Faults; (c) Salt dome; (d) Other.
Figure 6: Illustrative examples of automatically extracted image patches. Here the first column shows four exemplar patches; the second column shows automatically extracted patches, where the similarity scores are among the top 10 of all 500 similarity scores for each exemplar patch; the third column represents extracted patches whose similarity scores are ranked from 11 to 100; and for the fourth column, the similarity scores are ranked from 101 to 500.

The second step of the training process is to extract texture attributes from each patch in the training set. Before applying any of the techniques discussed in this paper to do the attribute extraction from an image patch , we first calculate the Hadamard product of this patch with a two-dimensional Gaussian kernel of the same size, . This kernel gives more weights to the structures at the center of the patch and less to those on the periphery, thus emphasizing local spatial correlations in seismic data. The procedure can be expressed as follow:

(12)

where is the Hadamard product, and the Gaussian kernel is defined as

(13)

where and are the - and - coordinates of the center of , respectively. The value of was set to 25 in our experiments so that pixels in the corners of the patch have weights of less than .

After this pre-processing, one of the eight techniques is applied to the patch to generate texture attributes. For LBP, the direct output generated for a given patch is a map of LBP codes of the same size as the original patch. For texture classification, a histogram is typically calculated from the LBP map and fed into the classifier as the final attributes. For other local descriptors (e.g., CLBP, M-CLBP, ELBP, CLDP, and LRI), multiple components are calculated in the form of multiple attribute maps (e.g., maps for CLBP-S, CLBP-C, and CLBP-M, respectively). In such cases, either a joint histogram or a concatenation of the individual histograms are formed for the classifier to use. To retain consistency, the semblance attribute is converted to a histogram when examined for the labeling task. The only exception is for the GLCM attributes, each of which is already computed from a two-dimensional histogram (i.e., the co-occurrence matrix) based on the patch. In this case, we combine different GLCM attributes into one vector as an input to the classifier.

In the third step of the training process, histograms generated from the texture attributes are used to train the classifier to establish statistical models for each category of subsurface structures of interest. In this research, we choose the support vector machine (SVM)

Vapnik (1999)

as the classifier, which is a powerful binary classification algorithm. It seeks to find the optimal separating hyperplane between two classes by identifying the one with the maximum margin. Since we have a multi-class classification problem, we train four hard-margin SVMs with linear kernels using the one-versus-all (OVA) approach.

Once the training process is complete, the label assignment is fulfilled in the labeling process, which is divided into three steps. First, a segmentation is performed in which each seismic section to be labeled is automatically divided into segments according to structures. To accomplish this purpose, we employ a superpixel-based segmentation approach that groups neighboring image pixels of similar appearance into a cluster or a superpixel. Each superpixel is treated as a single unit in the following processing. In this work, we adopt the superpixel segmentation algorithm based on the simple linear iterative clustering (SLIC) Achanta et al. (2012). In the original SLIC, vectors in the form of are generated for each pixel in an image to be segmented, where , , and are the three components of the color model, and , are the coordinates for each pixel. Then clustering is performed in a space formed by these vectors to obtain the superpixels. Because seismic images are in grayscale, we compute vectors for the pixels in a modified form, i.e., , in which and refer to the gradient along the x- and y- directions, respectively. Afterwards, we generate superpixels by clustering the vectors.

In the second step of labeling, similar to the training process, texture attributes are extracted for each segment or superpixel. Typically, the size of a superpixel is smaller than that of the image patches in the training dataset. To make sure that the attribute extraction is consistent between the training and the labeling processes, in this work, we select a neighborhood centered around the centroid of the superpixel, which is of the same size as that of the training patches. Attributes are extracted from this selected neighborhood (rather than the smaller area covered by the superpixel) to represent the superpixel. Then, in the last step of labeling, histograms of attributes generated for each superpixel are fed into the SVM classifier. They are compared against the trained SVM models to determine which of the four types of structures each superpixel belongs to. Thus, each one of them is assigned with an appropriate label.

The training and labeling processes consist of four components: data extraction (for training), attributes extraction (for both training and labeling), segmentation (for labeling), and classification (for both training and labeling). Training the classifier is considered as part of the classification component.

4 Results

We examined texture attributes within the framework of the computational labeling. As mentioned above, we used seismic sections extracted from the Netherlands North Sea Offshore F3 Block dGB Earth Sciences (1987). We specified four structures to be labeled: Chaotic layers, Faults, Salt dome, and Other. After training the SVM classifier as described in the previous section, four selected seismic sections, such as crossline 61, 211, 231, and 281 were labeled using each attribute. In Table 2, we list parameter settings we used to generate each attribute during the experiments.

Attribute Settings
GLCM Combined: Contrast, Entropy, Energy, Homogeneity, Correlation, and Mutual Information
Semblance Neighborhood:
LBP Radius (): ; Samples (): ; Mapping:
CLBP Radius (): ; Samples (): ; Mapping:
M-CLBP Radius (): ; Samples (): ; Mapping:
ELBP Radius (): ; Samples (): ; Mapping:
CLDP Radius (): ; Samples (): ; Mapping:
LRI LRI-A; Threshold for edge (): ; Threshold for size (): ; Directions:
Table 2: Parameter settings used for each attribute

The labeling performance was evaluated both subjectively by visual inspection and objectively in terms of four metrics commonly accepted in the semantic segmentation literature Long et al. (2015): pixel accuracy (PA), mean class accuracy (MCA), mean intersection over union (MIU), and frequency-weighted intersection over union (FWIU). The four metrics all range from to , with a greater value indicating a better performance. PA counts the overall rate of correctly labeled pixels, and MCA calculates the average rate of such labels among all classes. However, both of them neglect wrongly labeled pixels that have a negative impact on labeling performance. To account for such false labels, MIU and FWIU can be employed. Assuming that manually labeled data from an expert interpreter is available as the ground truth, the metrics are defined as below:

(14)
(15)
(16)
(17)

where refers to the number of pixels that belong to class but misclassified as , indicates the number of classes, and denotes the total number of pixels in class .

Figure 7 shows the crossline 281 labeled using each attribute for illustration purpose. The observations are very similar for all four seismic sections being tested. In general, labeling using most of the attributes is able to locate the main structures in presence. The major errors are in the Faults class. This phenomenon is especially true when only one fault exists in the middle of the patch, which is the case for the large fault in the middle of crosslines 211 and 231 (see Figure 8). For GLCM-based labeling, in addition to this single fault, it also misses part of the multiple faults. It is common to label some of the horizons as faults, as well. We believe the reason for the errors in identifying faults lies in the fact that, compared to other structures of interest, Faults are subtle and thin. Patches consisting of faults are typically mixtures of faults and other structures in the neighborhood. Thus, the labeling algorithm can be confused and assign wrong labels. Another factor is that we are also limited by the single Faults exemplar image that was used. To overcome this problem, a solution is to map the patch-level labels into pixel-level labels so that the labeling can be performed at the pixel level. Details of such an approach can be found in a recent work Alaudah and AlRegib (2017).

Generally, labeling errors in the results displayed include three types. In some cases, falsely labeled patches are actually connected to correctly labeled structures, appearing as an expanded highlighted area. This kind of labeling errors are usually acceptable if the extended area is within a reasonable range so that an interpreter can still easily determine what structure is being highlighted. Some of the falsely labeled patches are not connected to a correctly labeled structure, but they are grouped together as a few large highlighted areas. Such labeling errors are manageable as they are commonly in small quantities. Thus, an interpreter can probably rule them out quickly. The worst type of labeling errors are those scattered, isolated small patches with wrong labels. Not only are they visually annoying, but they also cause significant inconvenience for an interpreter to examine each one of them.

Specifically, the GLCM-based labeling locates the chaotic structure well. It also identifies the salt dome, but shows some scattered errors at the same time. The semblance-based approach yields a more noisy appearance in the labeled section, especially for the salt dome, which is not desirable. For LBP and its variants, labeling errors for the salt dome are mostly of the first type (i.e., connected errors). Therefore, the labeled salt dome areas generally appear reliable. Regarding the limited scattered errors for the salt dome, the labeling quality is the highest with CLDP, followed by M-CLBP, CLBP, ELBP, and LBP. As to the scattered errors for faults, the LBP-based labeling shows the least of them, while M-CLBP goes to the next and CLDP is associated with the highest. For the chaotic structure, in general, ELBP gives the best labeling quality. Finally, the LRI-based labeling locates the true salt dome better than the LBP-like techniques, capturing more of the salt body. However, at the same time, it shows a little more labeling errors for the salt dome. A major drawback with LRI is that faults in the labeled results exhibit many scattered errors.

(a)
(b) (c)
(d) (e)
(f) (g)
(h) (i)
Figure 7: Labeling results for the seismic section at crossline 281, in the Netherlands North Sea F3 block database dGB Earth Sciences (1987), using different texture attributes. The Chaotic class is in blue, Faults is in green, Salt Dome is red, and Other is grey. (a) Manual (manually labeled result overlaid upon the original data); (b) GLCM; (c) Semblance; (d) LBP; (e) CLBP; (f) M-CLBP; (g) ELBP; (h) CLDP; (i) LRI.
(a)
(b) (c)
(d) (e)
(f) (g)
(h) (i)
Figure 8: Labeling results for the seismic section at crossline 231, in the Netherlands North Sea F3 block database dGB Earth Sciences (1987), using different texture attributes. The Chaotic class is in blue, Faults is in green, Salt Dome is red, and Other is grey. (a) Manual (manually labeled result overlaid upon the original data); (b) GLCM; (c) Semblance; (d) LBP; (e) CLBP; (f) M-CLBP; (g) ELBP; (h) CLDP; (i) LRI.

We show objective evaluation results in Table  3 and Figure 9, which match our observation with the actual labeled sections. On average, the two traditional techniques perform at the two extremes, with the GLCM attributes yielding the best results and the semblance being the worst. Among the remaining six local descriptors, ELBP is associated with the highest PA and FWIU, mainly because its background is much less noisy. LRI is the one giving the highest MCA and MIU, indicating that it performs more uniformly across different types of structures. M-CLBP is the only local descriptor that yielded top three performance for all four metrics. As demonstrated in Figure 9, the performance of each attribute is relatively consistent across all four seismic sections being labeled. Considering the very limited manually labeled exemplars that were used, and the challenging nature of the task, the overall performance of the labeling with most texture attributes studied here is very promising.

Attributes GLCM SEMB. LBP CLBP M-CLBP ELBP CLDP LRI
PA 0.7625 0.5066 0.6491 0.7117 0.7171 0.7427 0.6698 0.6814
MCA 0.7465 0.6113 0.6503 0.7008 0.7102 0.6134 0.6975 0.7453
MIU 0.4725 0.2729 0.3601 0.4322 0.4346 0.4112 0.3989 0.4390
FWIU 0.6616 0.3910 0.5408 0.6095 0.6107 0.6322 0.5586 0.5802
Table 3: Evaluation of the labeling performance, averaged over four different seismic sections, with the top three highlighted for each metric
Figure 9: Evaluation of the labeling performance on each seismic section.

In addition to subsurface structure characterization, texture attributes are also important for facies analysis. Naturally, it is of interest to examine if the attributes (i.e., the local descriptors) discussed in this paper can be useful for analyzing facies. Therefore, we applied the framework introduced in this paper to labeling facies instead of structures. We performed a simple labeling experiment to demonstrate the feasibility. In this experiment, we selected three facies in the F3 Block dataset, namely, highstand system tract (HST), lowstand system tract (LST), and transgressive system tract (TST), as described in Illidge et al. (2016). As shown in Figure 10, our labeling tool is capable of producing reasonable results. In particular, the M-CLBP attribute yielded the best labeling output, which is clean and identifies the main locations of each facies very well. Again, this fulfills the purpose of the labeling, which is not to provide accurate delineation but to highlight approximate locations of targets of interest. We also note that the labeling framework is designed towards subsurface structure characterization. Thus, it may need to be adjusted if facies analysis becomes the main application of interest. For example, the segmentation process may be better implemented using the GoT algorithm Shafiq et al. (2017). However, further discussion along this direction is beyond the scope of this paper.

(a) (b)
(c) (d)
Figure 10: Example results labeling facies using inline 370 of F3 Block. Here, for demonstration purpose, only areas consisting of typical HST, LST, and TST are considered. Training samples were extracted from inline 530. (a) Original inline section with the three facies manually labeled (HST: red, LST: green, and TST: blue); (b) Results using GLCM; (c) Results using M-CLBP; (d) Results using LRI. For the experiment, the patch size adopted was .

5 Conclusion

In this paper, we conducted a comparative study of a group of spatial texture attributes, including both traditional attributes commonly used in seismic interpretation and local texture descriptors that gained popularity in recent texture image analysis literature. We examined these attributes in a new framework for seismic volume labeling and demonstrated that most of them can be utilized as a generic attribute to characterize different subsurface structures all at the same time. Thus, combined with a suitable interpretation tool such as the labeling tool, they can help provide an initial interpretation with structures of interest being highlighted so that the following interpretation can be expedited.

With respect to the specific attributes discussed in this paper, we found that in the current workflow, none of these attributes can provide satisfactory labeling results for single faults. The subtlety and thinness of such structures are not well captured by the patch-based training. However, we believe the same local descriptors can perform much better for the faults, given that they are able to keep track of local variations at a small scale. The key will be adjusting the workflow to incorporate a pixel-based training so that the classifier is trained with better-suited data such as delineated fault pixels.

The GLCM-based approach yielded the best overall performance in the labeling experiments, which can be attributed to two factors. First, the GLCM attributes are values derived from a histogram, not the histogram itself as with other attributes. Using the derived values can be advantageous in providing a higher-level description of the texture pattern. Second, the GLCM approach combines different types of attributes, creating a more comprehensive representation of various characteristics of the texture pattern. However, creating a GLCM requires a patch of an enough size, which sets such a limit that GLCM attributes are not suitable for small scale structures.

For the purpose of labeling, different attributes can be selected according to the significance of each structure to be located. If salt dome is of the most importance, then CLDP is the best attribute to use for a reliable representation of salt domes. If faults are more important than the others, then LBP and M-CLBP can be considered. If there is no preference for a specific structure, then GLCM, M-CLBP, ELBP and LRI are all good candidates for the labeling system.

There are two purposes of this paper. The first is to introduce and explore the image-based texture attributes, i.e., the local descriptors. The second is to introduce the seismic volume labeling workflow as a useful interpretation tool. However, the labeling tool is not limited to work with the local descriptors only. We believe it is worthwhile to explore its combination with any other kinds of texture attributes. It can also be customized to label specific structures of interest. For example, the labeling tool may be combined with attributes based on textural orientation variations for unconformity detection Ringdal (2012); Wu and Hale (2015)

; it can also be studied in terms of structure tensors

Bakker (2002); Fehmers and Höcker (2003); Wu and Janson (2017) for better identification of faults and channels Wu (2017). We believe these interesting future research will further prove the value of the labeling tool.

6 Acknowledgments

This work is supported by the Center for Energy and Geo Processing (CeGP) at Georgia Tech and King Fahd University of Petroleum and Minerals (KFUPM).

References

  • Achanta et al. (2012) Achanta, R., A. Shaji, K. Smith, A. Lucchi, P. Fua, and S. Süsstrunk, 2012, SLIC superpixels compared to state-of-the-art superpixel methods: Pattern Analysis and Machine Intelligence, IEEE Transactions on, 34, 2274–2282.
  • Alaudah and AlRegib (2016) Alaudah, Y., and G. AlRegib, 2016, Weakly-supervised labeling of seismic volumes using reference exemplars: Image Processing (ICIP), 2016 IEEE International Conference on, IEEE, 4373–4377.
  • Alaudah and AlRegib (2017) ——–, 2017, A weakly supervised approach to seismic structure labeling, in SEG Technical Program Expanded Abstracts 2017: Society of Exploration Geophysicists, 2158–2163.
  • Alfarraj et al. (2016)

    Alfarraj, M., Y. Alaudah, and G. AlRegib, 2016, Content-adaptive non-parametric texture similarity measure for image retrieval: Multimedia Signal Processing (MMSP), 2016 IEEE Workshop on, IEEE, 1–6.

  • Bakker (2002) Bakker, P., 2002, Image structure analysis for seismic interpretation: Delft University of Technology.
  • Beghdadi et al. (2015) Beghdadi, A., M. A. Qureshi, and M. Deriche, 2015, A critical look to some contrast enhancement evaluation measures: Colour and Visual Computing Symposium (CVCS), 2015, IEEE, 1–6.
  • Berthelot et al. (2013) Berthelot, A., A. H. Solberg, and L.-J. Gelius, 2013, Texture attributes for detection of salt: Journal of Applied Geophysics, 88, 52–69.
  • dGB Earth Sciences (1987) dGB Earth Sciences, B., 1987, The Netherlands Offshore, The North Sea, F3 Block - Complete: https://opendtect.org/osr/pmwiki.php/Main/Netherlands/OffshoreF3BlockComplete4GB,.
  • Fehmers and Höcker (2003) Fehmers, G. C., and C. F. Höcker, 2003, Fast structural interpretation with structure-oriented filtering: Geophysics, 68, 1286–1293.
  • Gao (2003) Gao, D., 2003, Volume texture extraction for 3d seismic visualization and interpretation: GEOPHYSICS, 68, P1294–P1302.
  • Gao (2007) ——–, 2007, Application of three-dimensional seismic texture analysis with special reference to deep-marine facies discrimination and interpretation: Offshore angola, west africa: AAPG bulletin, 91, 1665–1683.
  • Gonzalez and Woods (2006) Gonzalez, R., and R. Woods, 2006, Digital image processing (3rd edition): Prentice-Hall, Inc.
  • Guo et al. (2010) Guo, Z., L. Zhang, and D. Zhang, 2010, A completed modeling of local binary pattern operator for texture classification: Image Processing, IEEE Transactions on, 19, 1657–1663.
  • Haralick et al. (1973) Haralick, R. M., K. S. Shanmugam, and I. Dinstein, 1973, Textural features for image classification.: IEEE Transactions on Systems, Man, and Cybernetics, 3, 610–621.
  • Hu et al. (2016) Hu, Y., Z. Long, and G. AlRegib, 2016, Completed local derivative pattern for rotation invariant texture classification: Proceedings of the International Conference on Image Processing (ICIP), IEEE, 3548–3552.
  • Illidge et al. (2016) Illidge, E., J. Camargo, and J. Pinto, 2016, Turbidites characterization from seismic stratigraphy analysis: Application to the netherlands offshore f3 block: Presented at the AAPG/SEG International Conference & Exhibition.
  • Liu et al. (2012) Liu, L., L. Zhao, Y. Long, G. Kuang, and P. Fieguth, 2012, Extended local binary patterns for texture classification: Image and Vision Computing, 30, 86–99.
  • Long et al. (2015)

    Long, J., E. Shelhamer, and T. Darrell, 2015, Fully convolutional networks for semantic segmentation: Proceedings of the IEEE conference on computer vision and pattern recognition, 3431–3440.

  • Marfurt et al. (1998) Marfurt, K. J., R. L. Kirlin, S. L. Farmer, and M. S. Bahorich, 1998, 3-D seismic attributes using a semblance-based coherency algorithm: Geophysics, 63, 1150–1165.
  • Ojala et al. (2002) Ojala, T., M. Pietikäinen, and T. Mäenpää, 2002, Multiresolution gray-scale and rotation invariant texture classification with local binary patterns: Pattern Analysis and Machine Intelligence, IEEE Transactions on, 24, 971–987.
  • Pitas and Kotropoulos (1992) Pitas, I., and C. Kotropoulos, 1992, A texture-based approach to the segmentation of seismic images: Pattern Recognition, 25, 929–945.
  • Ringdal (2012) Ringdal, K., 2012, Flow-based segmentation of seismic data: PhD thesis, Citeseer.
  • Röster and Spann (1998) Röster, K., and M. Spann, 1998, A system for seismic data processing: Presented at the EUSIPCO 1998.
  • Shafiq et al. (2017) Shafiq, M. A., Z. Wang, G. AlRegib, A. Amin, and M. Deriche, 2017, A texture-based interpretation workflow with application to delineating salt domes: Interpretation, 5, SJ1–SJ19.
  • Vapnik (1999)

    Vapnik, V., 1999, An overview of statistical learning theory: IEEE Transactions on Neural Networks,

    10, 988–999.
  • Wang and AlRegib (2017) Wang, Z., and G. AlRegib, 2017, Interactive fault extraction in 3-d seismic data using the hough transform and tracking vectors: IEEE Transactions on Computational Imaging, 3, 99–109.
  • Wu (2017) Wu, X., 2017, Directional structure-tensor-based coherence to detect seismic faults and channels: Geophysics, 82, A13–A17.
  • Wu and Hale (2015) Wu, X., and D. Hale, 2015, 3d seismic image processing for unconformities: Geophysics, 80, IM35–IM44.
  • Wu and Janson (2017)

    Wu, X., and X. Janson, 2017, Directional structure tensors in estimating seismic structural and stratigraphic orientations: Geophysical Journal International,

    210, 534–548.
  • Zhai and Neuhoff (2014) Zhai, Y., and D. L. Neuhoff, 2014, Rotation-invariant local radius index: A compact texture similarity feature for classification: International Conference on Image Processing (ICIP), IEEE, 5711–5715.
  • Zhai et al. (2013) Zhai, Y., D. L. Neuhoff, and T. N. Pappas, 2013, Local radius index-a new texture similarity feature: International Conference on Acoustics, Speech and Signal Processing (ICASSP), IEEE, 1434–1438.