Tuesday, December 1, 2020

Table of Contents

Seismic Fundamental
AAA
Aliasing
Amplitude
Angle mute
Angle stack
Anisotropic NMO
Anisotropic parameters
Anisotropy
Anisotropy classification
Artificial Neural Network
Attenuation
Auto correlation
Bowtie
Brute Stack
CDP
Converted Waves
Crooked Line
Curvelet Transform
Dip Moveout
Dephasing
Detectability and Fidelity
Diffraction
Dogleg
Fermat principle
First Break
Frequency
Fresnel Zone
Gain
Groundroll
Higher order NMO
Hockey Stick
Huygens principle
Interval velocity
Prism wave
P wave velocity
Love wave
Migration
Migration aperture
Mistie
Multiple
NMO
Near offset far offset
Noise and data
Non zero apex
P wave
Polarity
PSDM
PSTM
Preserve Amplitude
Prewhitening
Q Migration
RMS Velocity
Ray Tracing
Reflectivity
Seismic Resolution
SELMA
SRME
Seismic record
Seismic reference datum
Seismic refraction
Seismic Refraction and ray tracing
Seismic tomography
Seismic trace
Seismic multicomponent
Semblance
Shear wave
Shear wave splitting
Slant Stack
Source signature
Spectral balancing
Spectral broadening
Spectral decomposition
Spectral ratio method
Spike
Stacking
Static correction
Super gather
Sweep vibroseis
TWT
Tau-p
Time to depth conversion
True amplitude recovery
Tuning thickness
Tuning Cube
Velocity Analysis
Velocity Dispersion
Vertical seismic resolution
Vp/Vs
Water column static
Wavelet
Wavelet polarity onshore
Zero crossing


Seismic Advanced
Anisotropic PSDM
Broadband Seismic 
Controlled Beam Migration
Full wave inversion
Reverse time migration


Signal Processing
Bandwidth
Convolution
Cross Correlation
dB/Octave
DC Component
Deconvolution
Dynamic range and BIT
Fourier Transform
Filtering
Gabor Transform
Hilbert Transform
Main lobe side lobe
Minimum phase zero phase
Multi Taper Method
Nyquist frequency


Seismic Acquisition, Well Logging, Drilling and Non Seismic
2D Marine Acquisition
3D Land survey
3D Seismic acquisition QC
4D Seismic
Checkshot survey
Gamma ray log
Horizontal Seismic Profiling
Pore Pressure Prediction
Remote Resistivity Mapping
Resistivity Log
Rich Azimuth Seismic 
Seismoelectric 
Kelly Bushing
Marine Analog Survey
Marine Geohazard survey
Neutron porosity
VSP
Wide Azimuth Seismic


Seismic Interpretation and Petroleum Geology
Bottom Simulating Reflector
Coal Bed Methane Seismic
Curvature
Envelope
Instantaneous Frequency 1
Instantaneous Frequency 2
Instantaneous phase 1
Instantaneous phase 2
Optical Stacking
Perigram
Perigram x Cosines phase
Petroleum System
Reflection Strength
Seismic Attribute
Seismic Facies 
Texture Attribute 
Unconformity
Velocity sag
Volume Assessment
Well seismic tie


Seismic Processing with Seismic UNIX (Basic)
1. Seismic Unix Installation
2. Data Download, Segy2SU, QC
3. Kill traces and gain
4. Groundroll Removal
5. Bandpass Filter and Decon
6. Geometry and CMP Sorting
7. NMO and Brute Stack
8. Velocity Analysis and Stack
9. Statics
10. Migration (STOLT)


Seismic Processing with Seismic UNIX (Advanced)
Attribute Seismic
Continuous Wavelet Transform
FX Decon
Interactive velocity analysis 
Mute Design
PreStack Time Migration
PreStack Depth Migration
PostStack Depth Migration
Radon Demultiple
Ray Tracing
Super gather 
Synthetic Gathers Generator
Velocity Medeling
VSP Processing
Waveform modeling


Reservoir Characterization
AVA
AVFVO
AVO Attributes
AVO Azimuth
AVO Classification
Acoustic Impedance
Anisotropic AVO
Aki and Richards Approximation
Coherence
Continuous Wavelet Transform
Elastic Impedance
Extended Elastic Impedance
Gassmann Fluid Substitution
Faust Equation
Intercept and Gradient AVO
Inversion Maximum Likelihood
Inversion Sparse Spike
Lambda Mu Rho
Mobil Approximation
Poisson Impedance
Recursive Inversion 
Reflection Impedance
Rock density
Rock Physic Analysis
Seismic Inversion
Seismic Multiattribute
Shuey Approximation
Smith and Gidlow 
Thinbed reflectivity inversion
Tuning Cube
Zoeppritz Elastic Impedance
Zoeppritz Equation


Computational Geophysics
Bresenham Algorithm
Madagascar
Marquardt Inversion
Matlab for horizon snap
SVD 1
SVD 2
SVDLib C
Seislab
SEGY2ASCII Matlab


Perl
Introduction to Perl
Hello World
Bit Reverse
Call GNUPLOT from Perl
Compare files 
Left of right justify
Matrix Multiplication (subroutine) 
Matrix Transpose
Print n-times according to integer in a specific column
Read text file and print on console
Read and write matrix file
Select Columns
Select rows 
Sort by key field
Sort by key fields
DEC2IEEE


Perl for Geoscientists
Ascii to SEGY
Automatic First Break Picks
Boundary Detection 
Compare two files and print closest values
Compare two files and search match columns
Copy header from another omega velocity file
Convert CGG to Omega velocity format
Discrete Fourier Transform 
Edit omega velocity file
Merge two omega velocity files
SEGY Header Dumper trace by trace
SEGY Header Dumper min-max 
SEGY Header Dumper user input 
SEGY to Ascii 
SEGD (Revision 1) to Ascii 
Select data inside polygon
Smoothing with convolution
SPS File Generator
Text to SEGY header
XPS File Generator


Perl Graphical User Interface (GUI)
Data Selector
Perl/Tk-GNUPLOT
Spectrum and RAAC Editor

C Language for Geoscientists
Hello World
Fahrenheit to Celsius Conversion
Progress indicator
Read input from keyboard entry and print on console
Read input from file perform operation and print on console
Read file perform operation and output file
Addition function
Temperature conversion function
Matrix addition function with keyboard entry
Matrix addition function load file
Matrix multiplication
Matrix transpose
Matrix identity
Matrix Inverse
Inverse Matrix with Marquardt method
Matrix operation with malloc
Matrix transpose with malloc
Fast Fourier Transform (FFT)
Format
SEGY Header Dumper


Python for Geophysicist
Python Installation
Plot XY Python
Plot Surface 2D Python
Plot XYZ 3D View Python
Plot XYZ as Points Python
Plot with Contour and Contourf Python
Map Extrapolation Python
1D Smoothing Savitzky Golay 
2D Smoothing Savitzky Golay
2D Smoothing Gaussian Python
2D Smoothing Moving Average and Save Result to file Python
Plot Map Customize Colorscale 
Plot 3D 
Delaunay Python 
Delaunay Triangulaton Python
Triangulation Grid Plot Python
High Resolution Tricontouring
Grid Rotation Smoothing and Plot Python
Extract Surface (Z) at any XY input Python 
Fast Fourier Transform

Miscellaneous
Animated GIF
AWK 1
AWK 2
AWK 3
BP Velocity Benchmark
Data Marine download
GNU Plot  
GNU Plot transparent filled polygon 
GNUPlot for contour and surface
GNUPlot for smoothing
GNUPlot 4D
GNUPlot and Perl to select based on polygon 
IBM and IEEE 
Perl One Liner
Promax 1
Promax 2
Promax 3
Promax 4
SEGY
Unix for geoscientists
VI

Friday, January 29, 2016

Plot XYZ as Points Python

#!/usr/bin/python
import numpy as np
import matplotlib.pyplot as plt
data = np.genfromtxt('input2d.txt')
x = data[:,0]
y = data[:,1]
z = data[:,2]

plt.scatter(x, y, c = z, lw = 0)
plt.clim(vmin=-1200,vmax=-200) #color limit
plt.colorbar()
plt.show()

Thursday, January 28, 2016

Broadband Seismic

Terminologi broadband seismic mengacu pada band-width frekuensi gelombang seismik yang lebih lebar dibandingkan dengan seismic konvensional. 

Keuntungan eksplorasi dengan menggunakan broadband seismic diantaranya meningkatkan resolusi event seismic dan berperan dalam menekan side lobes.

 Hubungan frekuensi dan wavelet gelombang seismic (courtesy GEO ExPro).

Sebagaimana yang terlihat dari gambar di atas (kiri) jika frekuensi semakin melebar ke arah frekuensi tinggi, resolusi temporal wavelet akan meningkat namun tidk berhubungan langsung dengan kemampuan menekan sidelobes.
Sedangkan (kanan), jika frekuensi semakin melebar ke arah frekuensi rendah, sidelobes akan semakin tertekan.


 Courtesy GEO ExPro

Gambar di atas menunjukkan model wedge serta seismogram sintetik dengan rentang frekuensi wavelet 2.5-50Hz dengan 10-50Hz. Terlihat bahwa jika frekuensi melebar ke arah frekuensi rendah 2.5-50Hz, definisi event wedge model jauh lebih sempurna.

Gambar di bawah ini menunjukkan stack data seismic konvensional dibandingkan dengan broadband seismic. Perhatikan bahwa event pada sedimen infill lebih terdefinisikan. disamping itu kandungan frekuensi rendah sebagai 'real' event masih bisa dipertahankan.


Courtesy CGG

Gambar di bawah ini menunjukkan perbandingan spectrum seismic konvensional dibandingkan dengan broadband seismic.
 
 Courtesy CGG

Terdapat beberapa faktor yang menyebabkan rekaman seismic memiliki band width yang lebih lebar diantaranya: sumber gelombang: type vibroseis, gag, teknik akuisisi, sampling interval ataupun spatial interval dalam akuisisi data seismik.