3D seismic image processing for unconformities

Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. IM35-IM44 ◽  
Author(s):  
Xinming Wu ◽  
Dave Hale
Keyword(s):  
Image Processing ◽  
3D Seismic ◽  
Amplitude Image ◽  
Tensor Method ◽  
Seismic Amplitude ◽  
Seismic Image ◽  
Automatic Methods ◽  
Normal Vectors ◽  
Seismic Reflectors ◽  
Seismic Images

In seismic images, an unconformity can be first identified by reflector terminations (i.e., truncation, toplap, onlap, or downlap), and then it can be traced downdip to its corresponding correlative conformity, or updip to a parallel unconformity; for example, in topsets. Unconformity detection is a significant aspect of seismic stratigraphic interpretation, but most automatic methods work only in 2D and can only detect angular unconformities with reflector terminations. Moreover, unconformities pose challenges for automatic techniques used in seismic interpretation. First, it is difficult to accurately estimate normal vectors or slopes of seismic reflectors at an unconformity with multioriented structures due to reflector terminations. Second, seismic flattening methods cannot correctly flatten reflectors at unconformities that represent hiatuses or geologic age gaps. We have developed a 3D unconformity attribute computed from a seismic amplitude image to detect unconformities by highlighting the angular unconformities and corresponding parallel unconformities or correlative conformities. These detected unconformity surfaces were further used as constraints for a structure-tensor method to more accurately estimate seismic normal vectors at unconformities. Finally, using detected unconformities as constraints and more accurate normal vectors, we could better flatten seismic images with unconformities.

Methods to compute fault images, extract fault surfaces, and estimate fault throws from 3D seismic images

Geophysics ◽  
2013 ◽  
Vol 78 (2) ◽  
pp. O33-O43 ◽  
Author(s):  
Dave Hale
Keyword(s):  
Image Processing ◽  
Data Structure ◽  
Efficient Algorithm ◽  
3D Seismic ◽  
Quadrilateral Mesh ◽  
Fault Surface ◽  
Seismic Image ◽  
Automatic Methods ◽  
Fault Interpretation ◽  
Seismic Images

Fault interpretation enhances our understanding of complex geologic structures and stratigraphy apparent in 3D seismic images. Common steps in this interpretation include image processing to highlight faults, the construction of fault surfaces, and estimation of fault throws. Although all three of these steps have been automated to some extent by others, fault interpretation today typically requires significant manual effort, suggesting that further improvements in automatic methods are feasible and worthwhile. I first used an efficient algorithm to compute images of fault likelihoods, strikes, and dips from a 3D seismic image. From these three fault images, I then automatically extracted fault surfaces as meshes of quadrilaterals that coincide with ridges of fault likelihood. A quadrilateral mesh is a simple data structure alongside which one can easily gather samples of the 3D seismic image. I automatically estimated fault throws by minimizing differences in values of samples gathered from opposite sides of a fault, while constraining the variation of throw within a fault surface. I tested the fidelity of estimated fault throws by using them to undo faulting. After unfaulting, reflectors in 3D seismic images were more continuous than those in the original 3D seismic image. In one example, this unfaulting test supported the observation that some extracted fault surfaces have unusual conical shapes.

Horizon volumes with interpreted constraints

Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. IM21-IM33 ◽  
Author(s):  
Xinming Wu ◽  
Dave Hale
Keyword(s):  
Structural Features ◽  
Control Points ◽  
3D Seismic ◽  
Entire Volume ◽  
Sedimentary Environments ◽  
Seismic Image ◽  
Normal Vectors ◽  
Seismic Reflectors ◽  
Seed Points ◽  
Seismic Images

Horizons are geologically significant surfaces that can be extracted from seismic images. Color coding of horizons based on amplitude or other attributes can help reveal ancient sedimentary environments and structural features. Extracted horizons are also used for building structure models and stratigraphic interpretations. We propose two methods for constructing seismic horizons aligned with reflectors in a 3D seismic image. The first method generates horizons one at a time; the second method generates an entire volume of horizons at once by first computing a relative geologic time volume from seismic normal vectors. Rather than gradually building a horizon by extending one or more seed points to a surface along seismic reflectors, both of our methods generate whole horizons at once by solving partial differential equations derived from seismic normal vectors. The most significant new aspect of both methods is the ability to specify, perhaps interactively during interpretation, a small number of control points that may be scattered throughout a 3D seismic image. Experiments revealed that with our method, control points enable the extraction of more accurate horizons from seismic images in which noise, unconformities, and faults are apparent. These points represent constraints that we implemented as preconditioners in the conjugate gradient method used to construct horizons.

Automatically interpreting all faults, unconformities, and horizons from 3D seismic images

2016 ◽  
Vol 4 (2) ◽  
pp. T227-T237 ◽  
Author(s):  
Xinming Wu ◽  
Dave Hale
Keyword(s):  
Image Processing ◽  
Partial Differential Equations ◽  
Differential Equations ◽  
Array Processing ◽  
3D Seismic ◽  
Partial Differential ◽  
Seismic Image ◽  
Processing Procedure ◽  
Seismic Images

Extracting fault, unconformity, and horizon surfaces from a seismic image is useful for interpretation of geologic structures and stratigraphic features. Although others automate the extraction of each type of these surfaces to some extent, it is difficult to automatically interpret a seismic image with all three types of surfaces because they could intersect with each other. For example, horizons can be especially difficult to extract from a seismic image complicated by faults and unconformities because a horizon surface can be dislocated at faults and terminated at unconformities. We have proposed a processing procedure to automatically extract all the faults, unconformities, and horizon surfaces from a 3D seismic image. In our processing, we first extracted fault surfaces, estimated fault slips, and undid the faulting in the seismic image. Then, we extracted unconformities from the unfaulted image with continuous reflectors across faults. Finally, we used the unconformities as constraints for image flattening and horizon extraction. Most of the processing was image processing or array processing and was achieved by efficiently solving partial differential equations. We used a 3D real example with faults and unconformities to demonstrate the entire image processing.

Moving faults while unfaulting 3D seismic images

Geophysics ◽  
2016 ◽  
Vol 81 (2) ◽  
pp. IM25-IM33 ◽  
Author(s):  
Xinming Wu ◽  
Simon Luo ◽  
Dave Hale
Keyword(s):  
Partial Differential Equations ◽  
Differential Equations ◽  
Extraction Methods ◽  
Seismic Interpretation ◽  
3D Seismic ◽  
Partial Differential ◽  
Multiple Faults ◽  
Seismic Image ◽  
Seismic Reflectors ◽  
Seismic Images

Unfaulting seismic images to correlate seismic reflectors across faults is helpful in seismic interpretation and is useful for seismic horizon extraction. Methods for unfaulting typically assume that fault geometries need not change during unfaulting. However, for seismic images containing multiple faults and, especially, intersecting faults, this assumption often results in unnecessary distortions in unfaulted images. We have developed two methods to compute vector shifts that simultaneously move fault blocks and the faults themselves to obtain an unfaulted image with minimal distortions. For both methods, we have used estimated fault positions and slip vectors to construct unfaulting equations for image samples alongside faults, and we have constructed simple partial differential equations for samples away from faults. We have solved these two different kinds of equations simultaneously to compute unfaulting vector shifts that are continuous everywhere except at faults. We have tested both methods on a synthetic seismic image containing normal, reverse, and intersecting faults. We also have applied one of the methods to a real 3D seismic image complicated by numerous intersecting faults.

3D seismic image processing for faults

Geophysics ◽  
2016 ◽  
Vol 81 (2) ◽  
pp. IM1-IM11 ◽  
Author(s):  
Xinming Wu ◽  
Dave Hale
Keyword(s):  
Image Processing ◽  
Data Structure ◽  
Linked Data ◽  
Difficult Problem ◽  
Fault Slip ◽  
3D Seismic ◽  
Fault Surface ◽  
Subsequent Processing ◽  
Seismic Image ◽  
Seismic Images

Numerous methods have been proposed to automatically extract fault surfaces from 3D seismic images, and those surfaces are often represented by meshes of triangles or quadrilaterals. However, extraction of intersecting faults is still a difficult problem that is not well addressed. Moreover, mesh data structures are more complex than the arrays used to represent seismic images, and they are more complex than necessary for subsequent processing tasks, such as that of automatically estimating fault slip vectors. We have represented a fault surface using a simpler linked data structure, in which each sample of a fault corresponded to exactly one seismic image sample, and the fault samples were linked above and below in the fault dip directions, and left and right in the fault strike directions. This linked data structure was easy to exchange between computers and facilitated subsequent image processing for faults. We then developed a method to construct complete fault surfaces without holes using this simple data structure and to extract multiple intersecting fault surfaces from 3D seismic images. Finally, we used the same structure in subsequent processing to estimate fault slip vectors and to assess the accuracy of estimated slips by unfaulting the seismic images.

Multitask learning for local seismic image processing: fault detection, structure-oriented smoothing with edge-preserving, and seismic normal estimation by using a single convolutional neural network

2019 ◽  
Vol 219 (3) ◽  
pp. 2097-2109 ◽  
Author(s):  
Xinming Wu ◽  
Luming Liang ◽  
Yunzhi Shi ◽  
Zhicheng Geng ◽  
Sergey Fomel
Keyword(s):  
Neural Network ◽  
Image Processing ◽  
Fault Detection ◽  
Convolutional Neural Network ◽  
Basic Step ◽  
Edge Preserving ◽  
Structural Interpretation ◽  
Seismic Image ◽  
Normal Vectors ◽  
Seismic Images

Summary Fault detection in a seismic image is a key step of structural interpretation. Structure-oriented smoothing with edge-preserving removes noise while enhancing seismic structures and sharpening structural edges in a seismic image, which, therefore, facilitates and accelerates the seismic structural interpretation. Estimating seismic normal vectors or reflection slopes is a basic step for many other seismic data processing tasks. All the three seismic image processing tasks are related to each other as they all involve the analysis of seismic structural features. In conventional seismic image processing schemes, however, these three tasks are often independently performed by different algorithms and challenges remain in each of them. We propose to simultaneously perform all the three tasks by using a single convolutional neural network (CNN). To train the network, we automatically create thousands of 3-D noisy synthetic seismic images and corresponding ground truth of fault images, clean seismic images and seismic normal vectors. Although trained with only the synthetic data sets, the network automatically learns to accurately perform all the three image processing tasks in a general seismic image. Multiple field examples show that the network is significantly superior to the conventional methods in all the three tasks of computing a more accurate and sharper fault detection, a smoothed seismic volume with better enhanced structures and structural edges, and more accurate seismic normal vectors or reflection slopes. Using a Titan Xp GPU, the training processing takes about 8 hr and the trained model takes only half a second to process a seismic volume with $128\, \times \, 128\, \times \, 128$ image samples.

Detecting faults and channels while enhancing seismic structural and stratigraphic features

2019 ◽  
Vol 7 (1) ◽  
pp. T155-T166 ◽  
Author(s):  
Xinming Wu ◽  
Zhenwei Guo
Keyword(s):  
Iterative Scheme ◽  
3D Seismic ◽  
Linear Features ◽  
Smoothing Methods ◽  
Seismic Image ◽  
Diffusion Scheme ◽  
Seismic Reflections ◽  
2D And 3D ◽  
Seismic Images

A 3D seismic image contains structural and stratigraphic features such as reflections, faults, and channels. When smoothing such an image, we want to enhance all of these features so that they are easier to interpret. Most smoothing methods aim to enhance reflections but may blur faults and channels in the image. A few methods smooth seismic reflections while preserving faults and channel boundaries. However, it has not well-discussed to smooth simultaneously along the seismic reflections and channels, which are linear features apparent within dipping reflections. In addition, to interpret faults and channels, extra steps are required to compute attributes or mappings of faults and channels from a seismic image. Such fault and channel attributes are often sensitive to noise because they are typically computed as discontinuities of seismic reflections. In this paper, we have developed methods to simultaneously enhance seismic reflections, faults, and channels while obtaining mappings of the faults and channels. In these methods, we first estimate the orientations of the reflections, faults, and channels directly in a seismic image. We then use the estimated orientations to control the smoothing directions in an efficient iterative diffusion scheme to smooth a seismic image along the reflections and channels. In this iterative scheme, we also efficiently compute mappings of faults and channels, which are used to control smoothing extents in the diffusion to stop smoothing across them. This diffusion scheme iteratively smooths a seismic image along reflections and channels while updating the mappings of faults and channels. By doing this, we will finally obtain an enhanced seismic image (with enhanced reflections and channels and sharpened faults) and cleaned mappings of faults and channels (discontinuities related to noise are cleaned up). We have examined the methods using 2D and 3D real seismic images.

Edge-aware filtering with Siamese neural networks

2020 ◽  
Vol 39 (10) ◽  
pp. 711-717
Author(s):  
Mehdi Aharchaou ◽  
Michael Matheney ◽  
Joe Molyneux ◽  
Erik Neumann
Keyword(s):  
Neural Networks ◽  
Seismic Data ◽  
Learning Ability ◽  
Seismic Exploration ◽  
Data Sets ◽  
3D Seismic ◽  
New Class ◽  
Seismic Image ◽  
Siamese Networks ◽  
Seismic Images

Recent demands to reduce turnaround times and expedite investment decisions in seismic exploration have invited new ways to process and interpret seismic data. Among these ways is a more integrated collaboration between seismic processors and geologist interpreters aiming to build preliminary geologic models for early business impact. A key aspect has been quick and streamlined delivery of clean high-fidelity 3D seismic images via postmigration filtering capabilities. We present a machine learning-based example of such a capability built on recent advances in deep learning systems. In particular, we leverage the power of Siamese neural networks, a new class of neural networks that is powerful at learning discriminative features. Our novel adaptation, edge-aware filtering, employs a deep Siamese network that ranks similarity between seismic image patches. Once the network is trained, we capitalize on the learned features and self-similarity property of seismic images to achieve within-image stacking power endowed with edge awareness. The method generalizes well to new data sets due to the few-shot learning ability of Siamese networks. Furthermore, the learning-based framework can be extended to a variety of noise types in 3D seismic data. Using a convolutional architecture, we demonstrate on three field data sets that the learned representations lead to superior filtering performance compared to structure-oriented filtering. We examine both filtering quality and ease of application in our analysis. Then, we discuss the potential of edge-aware filtering as a data conditioning tool for rapid structural interpretation.

Seismic fault detection with convolutional neural network

Geophysics ◽  
2018 ◽  
Vol 83 (5) ◽  
pp. O97-O103 ◽  
Author(s):  
Wei Xiong ◽  
Xu Ji ◽  
Yue Ma ◽  
Yuxiang Wang ◽  
Nasher M. AlBinHassan ◽  
...  
Keyword(s):  
Neural Network ◽  
Convolutional Neural Network ◽  
Field Data ◽  
Data Sets ◽  
3D Seismic ◽  
Learning Tools ◽  
Seismic Fault ◽  
Seismic Image ◽  
Coherence Estimation ◽  
Seismic Images

Mapping fault planes using seismic images is a crucial and time-consuming step in hydrocarbon prospecting. Conventionally, this requires significant manual efforts that normally go through several iterations to optimize how the different fault planes connect with each other. Many techniques have been developed to automate this process, such as seismic coherence estimation, edge detection, and ant-tracking, to name a few. However, these techniques do not take advantage of the valuable experience accumulated by the interpreters. We have developed a method that uses the convolutional neural network (CNN) to automatically detect and map fault zones using 3D seismic images in a similar fashion to the way done by interpreters. This new technique is implemented in two steps: training and prediction. In the training step, a CNN model is trained with annotated seismic image cubes of field data, where every point in the seismic image is labeled as fault or nonfault. In the prediction step, the trained model is applied to compute fault probabilities at every location in other seismic image cubes. Unlike reported methods in the literature, our technique does not require precomputed attributes to predict the faults. We verified our approach on the synthetic and field data sets. We clearly determined that the CNN-computed fault probability outperformed that obtained using the coherence technique in terms of exhibiting clearer discontinuities. With the capability of emulating human experience and evolving through training using new field data sets, deep-learning tools manifest huge potential in automating and advancing seismic fault mapping.

Sign in / Sign up

Export Citation Format

Share Document