bruges

name	bruges
description	Geophysical equations and rock physics calculations for seismic analysis. Use when Claude needs to: (1) Calculate AVO responses (Zoeppritz, Shuey, Aki-Richards), (2) Perform Gassmann fluid substitution, (3) Generate seismic wavelets (Ricker, Ormsby), (4) Compute reflectivity and synthetic seismograms, (5) Calculate elastic moduli from velocities, (6) Apply Gardner/Castagna empirical relations, (7) Model rock physics effects.
version	1.0.0
author	Geoscience Skills
license	MIT
tags	["Rock Physics","AVO","Gassmann","Wavelets","Seismic Modelling"]
dependencies	["bruges>=0.5.0","numpy","scipy"]
complements	["segyio","obspy","lasio","welly"]
workflow_role	analysis

bruges - Geophysics Equations

Quick Reference

import bruges
import numpy as np

# AVO - Zoeppritz reflectivity
from bruges.reflection import zoeppritz
theta = np.arange(0, 45, 1)
Rpp = zoeppritz(vp1, vs1, rho1, vp2, vs2, rho2, theta)

# Wavelet - Ricker
from bruges.filters import ricker
t, w = ricker(duration=0.128, dt=0.001, f=25)

# Fluid substitution
from bruges.rockphysics import gassmann
vp_new, vs_new, rho_new = gassmann(vp, vs, rho, k_min, rho_min,
                                    k_fl1, rho_fl1, k_fl2, rho_fl2, phi)

Module Overview

Module	Purpose
bruges.reflection	Zoeppritz, Shuey, Aki-Richards AVO equations
bruges.rockphysics	Gassmann, Gardner, Castagna, elastic moduli
bruges.filters	Ricker, Ormsby wavelets, convolution

AVO Analysis

from bruges.reflection import zoeppritz, shuey, akirichards
theta = np.arange(0, 45, 1)  # degrees

Rpp = zoeppritz(vp1, vs1, rho1, vp2, vs2, rho2, theta)  # Exact
Rpp = shuey(vp1, vs1, rho1, vp2, vs2, rho2, theta)      # Fast (~30 deg)
Rpp, G = akirichards(vp1, vs1, rho1, vp2, vs2, rho2, theta, return_gradient=True)

Fluid Substitution

from bruges.rockphysics import gassmann

vp_new, vs_new, rho_new = gassmann(
    vp_sat, vs_sat, rho_sat,   
    k_min, rho_min,             
    k_fl1, rho_fl1,             
    k_fl2, rho_fl2,             
    phi                         
)