Iterative Fault Enhancement Workflow

Similar documents
Enhancing faults and axial planes Program fault_enhancement

PRESTACK STRUCTURE-ORIENTED FILTERING PROGRAM sof_prestack

Improving Dip Estimates Program filter_dip_components

Prestack Workflows: Prestack Structure-Oriented Filtering. The Prestack Structure-Oriented Filtering Workflow. Introduction

Constructing a Geometric Attribute Workflow

(1B) Click Browse and select the AASPI-format file (*.H) that you want to display. After you browse the file, the Colorbar file name (1C), plot title

Seismic imaging PROGRAM aaspi_azim_offset_mig

STRUCTURE-ORIENTED FILTERING OF POSTSTACK DATA Program sof3d

Program flatten reads in an input seismic or attribute volume as well as a picked horizon and generates a flattened output volume:

Seismic Attributes on Frequency-enhanced Seismic Data

MODULATING A POLYCHROMATIC IMAGE BY ONE PLOTTED AGAINST LIGHTNESS: PROGRAM hlplot

Spectral Attributes Program spec_cmp. COMPUTING SPECTRAL COMPONENTS USING A COMPLEX MATCHING PURSUITE ALGORITHM PROGRAM spec_cmp

MULTIATTRIBUTE DISPLAY: PROGRAMS hlplot, hsplot, hlsplot, rgbplot, and crossplot

Well Log Utilities: Program LAS file importer

Which curvature is right for you? Part - 1

MODULATING A POLYCHROMATIC IMAGE BY A 2 ND IMAGE PLOTTED AGAINST SATURATION AND A 3 RD IMAGE PLOTTED AGAINST LIGHTNESS: PROGRAM hlsplot

COMPUTING APPARENT CURVATURE PROGRAM euler_curvature. Euler curvature computation flow chart. Computing Euler curvature

MAXIMUM ENTROPY SPECTRAL ANALYSIS PROGRAM spec_max_entropy

COMPUTING APPARENT DIP AND AMPLITUDE GRADIENTS PROGRAM apparent_cmpt. Apparent curvature computation flow charts

HRS-9 Data Slice Editing and Math. Jon Brown 2013

AUTOMATIC DETECTION AND EXTRACTION OF FAULTS FROM THREE- DIMENSIONAL SEISMIC DATA. Trygve Randen, Stein Inge Pedersen, and Lars Sønneland

Coherence and curvature attributes on preconditioned seismic data

FLATTENING AND GENERATING STRATAL SLICES PROGRAMS flatten, vector_flatten, stratal_slice, and vector_stratal_slice

Observing fracture lineaments with Euler curvature

Attribute Expression of the Mississippi Lime. Attribute expression of the Mississippi Lime

2D Seismic Utilities Programs: aaspi_segy_read_batch aaspi_2d_combine aaspi_pseudo3d_split aaspi_mega2d_split aaspi_segy_write_batch

Volumetric Attributes: Constrained least-squares spectral analysis Program spec_clssa

Computation flow chart

Downloaded 01/23/14 to Redistribution subject to SEG license or copyright; see Terms of Use at

Volumetric curvature attributes adding value to 3D seismic data interpretation

P. Bilsby (WesternGeco), D.F. Halliday* (Schlumberger Cambridge Research) & L.R. West (WesternGeco)

Volumetric Classification: Program pca3d

Structurally oriented coherent noise filtering

Interactive Velocity Analysis (aaspi_iva): Program: Precompute Velocity Semblance, Interactive Velocity Analysis, Velocity Interpolation

Structure-preserving Smoothing for Seismic Amplitude Data by Anisotropic Diffusion using GPGPU

AASPI SOFTWARE PARALLELIZATION

Volumetric Attributes- program curvature3d

Downloaded 10/29/15 to Redistribution subject to SEG license or copyright; see Terms of Use at

We LHR5 06 Multi-dimensional Seismic Data Decomposition for Improved Diffraction Imaging and High Resolution Interpretation

Sharpening through spatial filtering

Using Similarity Attribute as a Quality Control Tool in 5D Interpolation

Abstractacceptedforpresentationatthe2018SEGConventionatAnaheim,California.Presentationtobemadeinsesion

JNTUWORLD. 4. Prove that the average value of laplacian of the equation 2 h = ((r2 σ 2 )/σ 4 ))exp( r 2 /2σ 2 ) is zero. [16]

Chapter 11 Representation & Description

Techniques of Volume Blending for Aiding Seismic Data Interpretation

Coherent partial stacking by offset continuation of 2-D prestack data

Volumetric Curvature-Attribute Applications for Detection of Fracture Lineaments and their Calibration

Image Processing and Analysis

Linear Operations Using Masks

Filtering and Edge Detection. Computer Vision I. CSE252A Lecture 10. Announcement

Background on Kingdom Suite for the Imperial Barrel Competition 3D Horizon/Fault Interpretation Parts 1 & 2 - Fault Interpretation and Correlation

Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines

Seismic Time Processing. The Basis for Modern Seismic Exploration

Downloaded 09/01/14 to Redistribution subject to SEG license or copyright; see Terms of Use at

Lecture 4: Image Processing

Petrel TIPS&TRICKS from SCM

We Fast Beam Migration Using Plane Wave Destructor (PWD) Beam Forming SUMMARY

Curvature. Definition of curvature in 2D and 3D

Correlating Production or Injection Volumes to Volumetric Attributes Program cigar_probe. Computation flow chart

Geometric Seismic Attribute Estimation using Data-adaptive Windows

Kinematic Interactivity

P071 Land Data Regularization and Interpolation Using Azimuth Moveout (AMO)

PSPRO AZIMUTHAL. Big Data Software for Interactive Analysis and Interpretation of Full-Azimuth Seismic Gathers. Power to Decide

OVERVIEW. Photopia for SOLIDWORKS Tutorial. 3 PODT Lens Design

Because they are second-order derivatives, seismic curvature

Fault Detection Using Color Blending and Color Transformations

CoE4TN4 Image Processing

Residual move-out analysis with 3-D angle-domain common-image gathers

Preconditioning seismic data with 5D interpolation for computing geometric attributes

Geology 554 Environmental and Exploration Geophysics II Final Exam

7. Vertical Layering

Common-angle processing using reflection angle computed by kinematic pre-stack time demigration

DMT Manual Version 2.0 December 4, 2006

Algebra Based Physics

Digital Image Processing Fundamentals

Announcements. Edge Detection. An Isotropic Gaussian. Filters are templates. Assignment 2 on tracking due this Friday Midterm: Tuesday, May 3.

Software Installation - Accessing Linux and Checking your Environmental Variables

Volumetric Attributes: Computing Texture Attributes Program glcm3d

Petrel TIPS&TRICKS from SCM

Tu A11 01 Seismic Re-processing and Q-FWI Model Building to Optimise AVO and Resolution in the Shallow Wisting Discovery

FIFI-LS: Basic Cube Analysis using SOSPEX

G009 Multi-dimensional Coherency Driven Denoising of Irregular Data

3D shot-profile migration in ellipsoidal coordinates

Using Kingdom Suite for 3D Horizon/Fault Interpretation Parts 3 & 4 - Horizon Interpretation and Correlation

Th A4 11 Interpretational Aspects of Multispectral Coherence

Petrel TIPS&TRICKS from SCM

Seismic Data Enhancement While Interpretation A Challenge For Better Imagination of Structure and Reservoir Characterization With Confidence

Lecture: Edge Detection

Welcome to Solid Edge University 2015

Azimuthal binning for improved fracture delineation Gabriel Perez*, Kurt J. Marfurt and Susan Nissen

Chapter 26 Geometrical Optics

Examples of GLOBE Claritas Processing

Geology 554 Interpretation Project Big Injun Sand, Trenton/Black River Plays, Central Appalachian Basin, WV

Tu B1 02 Diffraction Modelling and Imaging of Fluid Escape Features, South China Sea

Th Volume Based Modeling - Automated Construction of Complex Structural Models

Petrel TIPS&TRICKS from SCM

Summary. Figure 1: Simplified CRS-based imaging workflow. This paper deals with the boxes highlighted in green.

Waves & Oscillations

Morphological Image Processing

VISTA. Desktop seismic data processing software

Transcription:

Iterative Fault Enhancement Workflow. Introduction A greater degree of filtering can be obtained using three methods: (1) applying more aggressive filter parameters using the original window, (2) applying the original less aggressive filters but using a larger filter window, or (3) cascading the original less aggressive filters using the original filter window in an iterative manner. Iterative application has three advantages over using a larger window size: (1) cascaded filters will implicitly provide a tapered response, (2) cascaded filters adapt to structure, and (3) cascaded filters are more economic to apply than a single filter of the same effective length. Iterative Fault Enhancement Workflow The AASPI iterative fault enhancement workflow GUI can be invoked from the main aaspi_util by first selecting the AASPI Workflows tab: The following workflow GUI will then pop up. Attribute-Assisted Seismic Processing and Interpretation - 24 December 2018 Page 1

1 2 3 4 5 6 7 8 Step 1. Defining input data and data parameters used by all programs in the workflow The input data consists of (1) an attribute that is sensitive to edges, typically one of the similarity (coherence) volumes, or total aberrancy if we wish to enhance faults that appear as flexures. If we wish to enhance axial planes, one would choose the most-positive principal curvature, k1.the confidence in an edge attribute can be sensitive to the amplitude of the seismic reflectivity or to some other measure of signal-to-noise ratio. In this example, we have used energy-ratio coherence as the input edge attribute and (2) total energy as the optional weight. If no weight file is chosen, all edge attribute values are treated equally (i.e. the weight is set to 1.0). The algorithms used in this workflow are sensitive to structure, requiring the use of (3) inline and (4) crossline components of reflector dip. Attribute-Assisted Seismic Processing and Interpretation - 24 December 2018 Page 2

The unique project name, suffix, and verbose options are the same as in most AASPI applications. The use mpi and number of processors options are common to all AASPI programs that run in parallel (which is the case for filter_single_attribute, fault_enhancement, and skeletonize3d used in this workflow). In this example, we are running across (5) three different workstations, each with 16 processors, for a total of 48 processors in all. Improved smoothing and sharpening is achieved by defining (6) the number of iterations. We ve found that for the default, smaller windows, only minimal changes occur beyond three iterations. If this is your first time to run the iterative fault enhancement workflow, you may wish to check (7) keep all generated iteration-generated files. In this case with three iterations, there will be file for fault dip azimuth, fault dip magnitude, and fault probability with the strings 1, 2, and 3 appended to their names. The default is to not keep the intermediate iteration output in order to save disk space. Finally, (8) click Save workflow environment parameters. These parameters will be shared with the successive programs. Once saved, Step 2 becomes activated, application program by program. Step 2. Defining application program parameters Once the parameters are saved, the first application program in the workflow, filter_single_attribute, becomes active: Clicking the filter_single_attribute button invokes the corresponding GUI: Attribute-Assisted Seismic Processing and Interpretation - 24 December 2018 Page 3

Note that the input data at the top of the GUI are grayed out and cannot be changed. They were defined previously by the workflow GUI. The default parameters for this filter are to apply a 50% alpha trimmed-mean filter to the similarity volume using an elliptical window of radius 25 m by 50 m by 0 ms, resulting in a zero-sample thick median filter being applied to the similarity volume along structural dip. This median filter will preserve similarity anomalies parallel to structure and suppress thin anomalies (such as faults) cutting structure at steep angles. The checkmark indicates that this filtered trend will be subtracted from the input similarity volume, resulting in a volume that favor faults, but suppresses anomalies subparallel to structure. Detailed description of this program can be obtained by bringing up the documentation using the Help button in the upper right. After inspecting and perhaps modifying these parameters, simply Save parameters and return to Workflow GUI, to return. At this point, one can select the next application, fault_enhancement, in the iterative fault enhancement workflow GUI: Attribute-Assisted Seismic Processing and Interpretation - 24 December 2018 Page 4

Clicking the fault enhancement button invokes the corresponding GUI: Attribute-Assisted Seismic Processing and Interpretation - 24 December 2018 Page 5

As in the previous application, the input data at the top of the GUI are grayed out and cannot be changed and were defined previously by the workflow GUI. The default parameters for this filter are to apply a directional derivative of a Gaussian that smooths along the hypothesized fault plane using a Gaussian smoother with σ1=75 m, and to sharpen perpendicular to this same hypothesized fault plane using a dlog sharpening/smoothing parameter of σ3=25 m. Detailed description of this program can be obtained by bringing up the documentation using the Help button in the upper right. After inspecting and perhaps modifying these parameters, simply Save parameters and return to Workflow GUI, to return. Attribute-Assisted Seismic Processing and Interpretation - 24 December 2018 Page 6

At this point, one can select the next application, skeletonize3d, in the iterative fault enhancement workflow GUI: Note that there are no parameters to change in this GUI. Simply click Save parameters and return to Workflow GUI and find that the Execute iterative fault enhancement workflow button is now activated. Attribute-Assisted Seismic Processing and Interpretation - 24 December 2018 Page 7

Step 3: Execute the workflow At this point simply click Execute Iterative Fault Enhancement Workflow button to submit the suite of processes. Example Detailed examples on the generation of each of these steps can be found in the documentation for filter_single_attribute, fault_enhancement, and skeletonize3d. In this documentation, we simply display the input similarity volume, the results of subtracting the median filtered similarity from the input, the results of three iterations of fault enhancement, and the final skeletonized image. Vertical slices through these four volumes appear as follows: Attribute-Assisted Seismic Processing and Interpretation - 24 December 2018 Page 8

Attribute-Assisted Seismic Processing and Interpretation - 24 December 2018 Page 9

Time slices at t=1.28 through these four volumes appear as follows: Attribute-Assisted Seismic Processing and Interpretation - 24 December 2018 Page 10

Attribute-Assisted Seismic Processing and Interpretation - 24 December 2018 Page 11

Attribute-Assisted Seismic Processing and Interpretation - 24 December 2018 Page 12