seismic-interpretation

name	seismic-interpretation
skill_type	workflow
description	End-to-end seismic interpretation workflow from SEG-Y loading through signal processing, rock physics, and visualization. Use when working with seismic data analysis pipelines.
version	1.0.0
author	Geoscience Skills
license	MIT
tags	["Seismic","Interpretation","SEG-Y","Rock Physics","Workflow","AVO","Processing"]
dependencies	["segyio","obspy","bruges","pyvista"]
complements	["segyio","obspy","bruges","disba","pyvista"]
workflow_role	processing

Seismic Interpretation Workflow

End-to-end pipeline for seismic data analysis, from loading SEG-Y files through signal processing, rock physics modelling, and 3D visualization.

Skill Chain

segyio          obspy           bruges          disba           pyvista
[SEG-Y I/O] --> [Signal Proc] --> [Rock Physics] --> [Dispersion] --> [3D Viz]
  |               |                |                 |                |
  Load traces     Filter/FFT       AVO modelling     Surface waves    Volume render
  Read headers    Instrument resp   Fluid sub         Phase velocity   Slice display
  3D geometry     Spectral anal    Synthetics        Inversion        Horizon pick

Decision Points

Question	If Yes	If No
Working with SEG-Y files?	Start with segyio	Use obspy for miniSEED/SAC
Need frequency filtering or spectral analysis?	Use obspy signal tools	Skip to rock physics
Performing AVO or fluid substitution?	Use bruges	Skip to visualization
Analysing surface waves (MASW/SASW)?	Use disba for dispersion	Skip disba
Need 3D volume rendering?	Use pyvista	Use matplotlib for 2D

Step-by-Step Orchestration

Stage 1: Data Loading (segyio)

import segyio
import numpy as np

# Load and inspect SEG-Y
with segyio.open('survey.sgy', 'r', iline=189, xline=193) as f:
    n_traces = f.tracecount
    n_samples = len(f.samples)
    dt = f.samples[1] - f.samples[0]  # ms
    cube = segyio.tools.cube(f)        # (ilines, xlines, samples)

    
    cdp_x = f.attributes(segyio.TraceField.CDP_X)[:]
    cdp_y = f.attributes(segyio.TraceField.CDP_Y)[:]