- Vaclav Bucha
- Department of Geophysics, Charles University Prague, Ke Karlovu 3, 121 16 Praha 2, Czech Republic
- Petr Bulant
- Department of Geophysics, Charles University Prague, Ke Karlovu 3, 121 16 Praha 2, Czech Republic
- Ludek Klimes
- Department of Geophysics, Charles University Prague, Ke Karlovu 3, 121 16 Praha 2, Czech Republic

This package is property of its authors, and may be used free of charge. If any result obtained with help of any part of this package is published, we appreciate citation of all papers related to the used part of this package. Refer to the references related to this package.

**Model with a lenticular inclusion (LEN)****Preliminary Reference Earth Model (PREM)****Model with homogeneous layers (ELF1, MI)****Salt dome models (SD1 and SD2)****Model with interface slitted by a curved vertical fault (SLIT)****Western Bohemia a priori model (WB2)****Model of a real geological structure (L7)****Gridded 2-D Marmousi model (MAR)****Unconformity 2-D model (U2D)****1-D model (RM)****Weakly anisotropic models (QI, QIH, QI2, QI4, QI8, SC1_I, SC1_II, KISS and ORT)****Sample synthetic 3-D model (98)****1-D anisotropic "twisted crystal" model (TC)****Simple smooth 3-D model (N2)****2-D representation of the Kummer random medium (RAN)****2-D model with a salt body (HESS)****1-D constant velocity gradient model (VGR)****SEG/EAGE Salt Model (SAL)****2-D model Pluto 1.5 (PLU)****Computation of R/T coefficients****Models of the Cotton Valley site (CV)****Model with planar interfaces (EX)****Correlation functions of random media (CORFUN)****2-D model P1I with reflection surfaces (P1, P1I)****Simple models to test the Born approximation**

- 3-D model with a lenticular inclusion is a synthetic demo model for SW3D package MODEL.
- VRML and GOCAD representations and snapshots of the model
*Examples with history files:*- Transforming surfaces and velocity sections into VRML or GOCAD. Computing 3-D P-wave velocity cube (GOCAD Voxet). len-mod.h
- Generating colour PostScript files with the velocity sections across the model. len-grd.h
- Generating 3-D grids of values of velocity and of indices of complex blocks. len-grid.h
- Performing initial-value ray tracing of rays shot in a single vertical plane from the point source. len-crti.h
- Performing 3-D two-point ray tracing for a point source. Adding source, receivers and two-point rays into VRML or GOCAD representation of the model with surfaces. len-crt.h
- Transforming wavefronts and two-point rays into VRML. len-wf.h
- 3-D two-point ray tracing for the plane wave incident at the bottom of the model. lenb-crt.h
- 3-D two-point ray tracing for exploding reflector (zero-offset rays). lenz-crt.h
- Controlled initial-value ray tracing in the model followed by 3-D interpolation of travel times within ray cells. leni-mtt.h
- Initial-value ray tracing of the cone of rays in the model followed by calculation of the derivatives of the travel times with respect to the model B-spline coefficients. Adding source, receivers and rays into VRML or GOCAD representation of the model with surfaces. len-cone.h
- Performing shortest path calculation of the first arrival travel times and the corresponding rays on a rectangular grid in the model (network ray tracing). len-net.h
- Simple demonstration history file used to debug and test the inversion programs. The data to be inverted are created by discretization of model LEN. The same model is then obtained by the overdetermined least-square inversion. len-inv.h

- Brief description of data files for packages FORMS, MODEL, CRT, NET
- Data files for packages FORMS, MODEL, CRT, NET
- Related papers

**Preliminary Reference Earth Model (PREM)** - 3-D B-spline approximation of the model PREM (Dziewonski and Anderson 1981), specified in geographical spherical coordinates.
- Snapshot of the model
*Examples with history files:*- Generating PostScript files with the sections across the model PREM displaying the velocities and indices of complex blocks. prem-grd.h
- Transforming surfaces and a velocity section into VRML or GOCAD. prem-mod.h
- Tracing pPKP+pPKiKP, pPKIKP, PKP+PKiKP and PKIKP rays in spherical coordinates. prem-crt.h

- Guide to the input data for the complete ray tracing program.
- Data files for packages FORMS, MODEL, CRT

**Model with homogeneous layers (ELF1, MI)** - 3-D model with four homogeneous layers separated by dipping plane interfaces.
- VRML and GOCAD representations and snapshots of the model
*Examples with history files and perl scripts:*- Transforming model with homogeneous layers separated by plane interfaces (ELF1) into VRML or GOCAD. elf1-mod.h
- Performing 3-D two-point ray tracing for moving source and receiver profile. elf1a.pl, elf1b.pl, elf1c.pl

- Data files for packages FORMS, MODEL, CRT
- Related papers

**Salt dome model (SD1) and dislocated salt dome model (SD2)** - Sample 3-D synthetic models of a salt dome.
- VRML and GOCAD representations and snapshots of the model
*Examples with history files:*- Data files for packages FORMS, MODEL
- Related papers

**Model with interface slitted by a curved vertical fault (SLIT)** - Image of the model
- Data files for packages FORMS, MODEL
- Related papers

**Western Bohemia a priori model (WB2)** - 3-D Western Bohemia a priori model prepared as a test model resembling the conditions of refraction seismic measurements performed in Western Bohemia.
- VRML and GOCAD representations and snapshots of the model
*Examples with history files:*- Data files for packages FORMS, MODEL
- Related papers

**Model of a real geological structure (L7)** - 3-D model of a real geological structure provided by Institut Francais du Petrole. Present version contains only two of seven structural interfaces.
- Image of the model
- VRML and GOCAD representations and snapshots of the model
*Examples with history files:*- Transforming interfaces in model L7 into VRML or GOCAD. l7-mod.h
- Generating PostScript files with the sections across model L7 displaying the P-wave velocities.
- Two-point ray tracing in model L7, adding rays to the VRML or GOCAD files with interfaces.

- Data files for packages FORMS, MODEL, CRT
- Related papers

**Gridded 2-D Marmousi model and INRIA bench-mark versions of the 2-D Marmousi model (MAR)** - 2-D Marmousi model is based on real geological structure (Versteeg & Grau 1991).
*Examples with history files:*- INRIA bench-mark versions of the Marmousi model
- Generating PostScript files showing the P wave velocity sections in both INRIA bench-mark versions of the Marmousi model. mar-grd.h
- Calculating travel times in the INRIA Marmousi models
- Smoothing the Marmousi model for Gaussian-packet migrations.
- Optimization of the Shape of Gaussian Beams in the smoothed Marmousi model.
- Decomposition of the wave field into optimized Gaussian packets.
- Gaussian packet migration in the smoothed Marmousi model.
- Zero-offset ray tracing in the smoothed Marmousi model.

- Data files for packages FORMS, MODEL, CRT
- Related papers
- Reference:
- Versteeg, R.J. & Grau, G. (1991): The Marmousi experience. Proc. EAGE workshop on Practical Aspects of Seismic Data Inversion (Copenhagen, 1990), Eur. Assoc. Explor. Geophysicists, Zeist.

**Unconformity 2-D model (U2D)** - Unconformity 2-D model with attenution by Cormier and Mellen (1984).
- Image of the model
*Examples with history files:*- Generating PostScript file with the section across the model U2D displaying the P-wave velocity. u2d-grd.h
- Performing 2-D two-point ray tracing and generating synthetic seismograms with no attenuation, noncausal attenuation, and causal attenuation. u2d-ss.h
- Performing 2-D two-point ray tracing and plotting travel times. u2d-tt.h
- Comparing ray-theory and finite-difference synthetic seismograms. u2d-fd.h
- Generating system of equations for travel-time inversion. u2d-inv.h

- Data files for packages FORMS, MODEL, CRT, FD
- Related papers
- Reference:
- Cormier, V.F. & Mellen, M.H. (1984): Application of asymptotic ray theory to vertical seismic profiling. In: Toksoz, M.N. & Stewart, R.R. (eds.), Vertical Seismic Profiling: Advanced Concepts, pp. 28-44, Geophysical Press, London.

**1-D model (RM)** - 1-D model RM consists of a thin homogeneous low-velocity layer covering homogeneous halfspace.
*Examples with history files:*- Comparing ray-matrix and 2-D finite differences seismograms in 1-D model RM. rm-fd.h

- Data files for packages FORMS, MODEL, CRT, RM, FD
- Related papers

**Weakly anisotropic models (QI, QIH, QI2, QI4, QI8, SC1_I, SC1_II, KISS and ORT)** - Weakly anisotropic models QI, QIH, QI2, QI4, QI8, SC1_I, SC1_II, KISS and ORT for tests of coupling ray theory and for comparison of synthetic seismograms with package ANRAY, FORT and with Fourier pseudo-spectral method.
*Examples with history files:*- Computing coupling ray theory and four kinds of quasi-isotropic ray theory synthetic seismograms. Calculating second-order perturbations of travel times. qi-ss.h
- Running history file 'qi-ss.h' for four different models, which differ in the degree of anisotropy. Calculating travel times by program ANRAY, and comparing the results of ANRAY and CRT. qi-ss4m.h
- Test of anisotropic common S wave ray tracing in different models. qi-acr.h
- Computing isotropic, anisotropic, coupling ray theory, and prevailing-frequency approximation of the coupling ray theory synthetic seismograms in the model with a very weak anisotropy. qih-ss.h
- Computing anisotropic, isotropic, coupling-ray theory, and prevailing-frequency approximation of the coupling ray theory S-wave synthetic seismograms in models QI, QI2 and QI4. Comparing the synthetic seismograms with package FORT and with Fourier pseudo-spectral method. qi1-ss.h, qi2-ss.h, qi4-ss.h
- Computing coupling-ray theory and prevailing-frequency approximation of the coupling ray theory S-wave synthetic seismograms in models SC1_I, SC1_II and ORT. Comparing the synthetic seismograms with package FORT and with Fourier pseudo-spectral method. sc-ss.h,
- Computing coupling-ray theory and prevailing-frequency approximation of the coupling ray theory S-wave synthetic seismograms in model KISS. Comparing the synthetic seismograms with package FORT and with Fourier pseudo-spectral method. kis-ss.h
- Computing prevailing-frequency approximation of the coupling-ray theory S-wave Green tensors along common reference anisotropic S-wave rays in models QIH, QI1, QI2, QI4, KISS, SC1_I, SC1_II and ORT using the program CRTPFA. Calculating synthetic seismograms from the obtained Green tensors and comparing them with synthetic seismograms obtained by program GREEN. Interpolation of prevailing-frequency approximation of the coupling-ray theory within ray tubes. qih-pfa.h, qi1-pfa.h, qi2-pfa.h, qi4-pfa.h, kis-pfa.h, sc-pfa.h

- Data files for packages FORMS, MODEL, CRT
- Related papers

**Sample synthetic 3-D model (98)** - Sample synthetic 3-D model created during the technical session of the SW3D consortium meeting in June 1998. The model is composed of three layers separated by two interfaces. A salt body in the middle layer is bounded by a spherical interface.
- VRML and GOCAD representations and snapshots of the model
*Examples with history files:*- Sample computations in the model 98 (generating PostScript figures of velocity sections, transforming surfaces into VRML or GOCAD, calculating multi-valued travel times, generating synthetic seismograms)

- Data files for packages FORMS, MODEL, CRT
- Related papers

**1-D anisotropic "twisted crystal" model (TC)** - 1-D anisotropic "simplified twisted crystal" and "oblique twisted crystal" models.
*Examples with history files:*- Analytically computing propagator matrices in model "simplified twisted crystal", and comparing the results with Green functions computed by CRT. Comparing the analytical solution with finite-difference solution in the frequency domain. Generating synthetic seismograms and plotting a figure of comparison of different ray-tracing methods. tc.h
- Analytically computing exact propagator matrix in model "oblique twisted crystal", and comparing the results with different quasi-isotropic approximations of the coupling ray theory computed by CRT. Generating synthetic seismograms and plotting a figure of the comparison. tco.h

- Data files: one-purpose Fortran codes, data for packages FORMS, MODEL, CRT and for the one-purpose codes
- Description of files related to the finite-difference solution
- Related papers

**Simple smooth 3-D model (N2)** - Simple smooth 3-D model to demonstrate network ray tracing in the paper by Klimes and Kvasnicka (1994).
*Examples with history files:*- Calculating velocity isolines in 2-D sections. Computing rays and travel times by network ray tracing. n2-net.h

- Brief description of data files for packages MODEL, NET
- Data files for packages MODEL, NET
- Related papers
- Reference:
- Klimes, L. and Kvasnicka, M. (1994): 3-D network ray tracing.
Geophys. J. int.,
**116**, 726-738.

**2-D representation of the Kummer random medium (RAN)** - Klimes, L. and Kvasnicka, M. (1994): 3-D network ray tracing.
Geophys. J. int.,
*Examples with history files:*- Generating a 2-D representation of the Kummer random medium and calculating first arrival travel times. ran-net.h

- Data files for packages FORMS, NET
- Related papers

**2-D model with a salt body (HESS)** - 2-D model of a real geological structure with a salt body by Scott Morton of the Amerada Hess Corporation.
*Examples with history files:*- Data files for packages FORMS, MODEL, CRT
- Related papers

**1-D constant velocity gradient model (VGR)** - 1-D constant velocity gradient model by Pascal Podvin (Podvin and Lecomte 1991).
*Examples with history files:*- Brief description of data files for packages FORMS, CRT, NET
- Data files for packages FORMS, CRT, NET
- Related papers
- Reference:
- Podvin, P. & Lecomte, I. (1991): Finite difference computation
of traveltimes in very contrasted velocity models: a massively
parallel approach and its associated tools.
Geophys.J.int.,
**105**, 271-284.

**SEG/EAGE Salt Model (SAL)** - Podvin, P. & Lecomte, I. (1991): Finite difference computation
of traveltimes in very contrasted velocity models: a massively
parallel approach and its associated tools.
Geophys.J.int.,
- 3-D SEG/EAGE Salt Model with features similar to those of the U.S. Gulf of Mexico salt structure (Aminzadeh et al. 1997).
- VRML and GOCAD representations and snapshots of the model smoothed for ray tracing
*Examples with history files:*- Data files for packages FORMS, MODEL, CRT
- Related papers
- Reference:
- Aminzadeh, F., Brac, J., Kunz, T. (1997): 3-D Salt and Overthrust Models. SEG/EAGE 3-D Modeling Series No.1., Soc. Explor. Geophysicists, Tulsa.

**2-D model Pluto 1.5 (PLU)** - 2-D model Pluto 1.5 of a real geological structure with a salt bodies created by Smaart Joint Venture.
*Examples with history files:*- Data files for packages FORMS, MODEL, CRT
- Related papers

**Computation of R/T coefficients** *Examples with history files:*- Data files for package FORMS
- Related papers

**Models of the Cotton Valley site (CV)** - Several models of a real structure in the Cotton Valley site based on a simultaneous inversion of the data by James Rutledge.
*Examples with history files:*- Data files for packages FORMS, MODEL, CRT

**Model with planar interfaces (EX)** - 3-D analytical model with five homogeneous layers separated by dipping plane interfaces.
- Sample synthetic 3-D model created during the technical session of the SW3D consortium meeting in June 2004. The model is composed of five layers separated by three planar interfaces and one interface composed of two horizontal and two tilted planes.
- VRML and GOCAD representations and snapshots of the reflected two-point rays calculated in the model
*Examples with history files:*- Data files for packages FORMS, MODEL, CRT

**Correlation functions of random media (CORFUN)** - History files to calculate different realizations of random media in terms of a white noise filtered by spectral filters enabling the random media used in geophysics to be specified (Gaussian, exponential, von Karman, self-affine, Kummer).
*Examples with history files:*- Generating PostScript figures of random functions with different correlation functions using 2-D program 'grdran2d.for'. corfun.h
- Generating and plotting representations of random functions with different correlation functions. Calculates the same figures as 'corfun.h' but using 3-D programs 'grdran.for', 'grdfft.for', 'grdcor.for', 'grdcal.for' and 'grdstat.for'. corfft.h
- Calculating the von Karman medium correlation function by FFT (programs 'grdfft.for', 'grdcor.for', 'grdcal.for') and comparing the result with analytic solution calculated by 'grdckn.for'. cknfft.h

- Data files for package FORMS
- Related papers

**2-D model P1I with reflection surfaces (P1, P1I)** - 2-D model with several reflection surfaces, data provided by Stream Oil & Gas Ltd.
*Examples with history files:*- Constructing the smooth model - model P1. p1-inv.h
- Constructing the model with interfaces - model P1I. p1i-inv.h
- Calculating the velocity differences in the models. p1-dv.h
- Calculating synthetic seismograms simulating a reflection measurement. p1i-ss.h
- Computing 3D seismograms using 2D ray-based Born approximation
in heterogenous model P1-8-10%. The background (unperturbed) model
is heterogenous smooth model P1 without interfaces.
Main history file to run history files 'p1-grid.h', 'p1-mtt.h',
'p1-iwf.h', 'p1-born.h', 'p1-ss2.h' and 'p1-ss.h'.
p1-borna.h
- Computing the elastic parameters at the gridpoints in both unperturbed and perturbed model and the perturbations of these elastic parameters. p1-grid.h
- Computing the Green functions from the source and from the receiver, slowness vectors, travel times and other quantities at the gridpoints. p1-mtt.h
- Computing the incident wavefield necessary for the Born approximation using the representation theorem. p1-iwf.h
- Computing the 2D Born approximation. p1-born.h
- Generating the ray-theory seismogram. p1-ss2.h
- Generating the Born seismogram from RF format file computed by history file 'p1-born.h' and visualizing it together with the ray-theory seismogram computed by history file 'p1-ss2.h'. p1-ss.h

- Data files for packages FORMS, MODEL, CRT
- Related papers

**Simple models to test the Born approximation** - Three simple models for computation of 3D seismograms using the 2D ray-based Born approximation. Each model contains two homogenous isotropic layers. The models differ only by the shape of the interface between the layers. The background is a homogenous isotropic model without interfaces.
*Examples with history files:*- Running history files 'm2d-mtt.h', 'm2d-iwf.h', 'm2d-born.h',
'm2d-ss2.h' and 'm2d-ss.h'.
m2d-main.h
- Computing the elastic perturbations, the Green functions from the source and from the receiver, slowness vectors, travel times and other quantities at the gridpoints. m2d-mtt.h
- Computing the incident wavefield necessary for the Born approximation using the representation theorem. m2d-iwf.h
- Computing the Born approximation. m2d-born.h
- Generating the ray-theory seismogram. m2d-ss2.h
- Generating the Born seismogram from RF format file computed by history file 'm2d-born.h' and visualizing it together with the ray-theory seismogram computed by history file 'm2d-ss2.h'. m2d-ss.h

- Running history files 'm2d-mtt.h', 'm2d-iwf.h', 'm2d-born.h',
'm2d-ss2.h' and 'm2d-ss.h'.
m2d-main.h
- Data files for packages FORMS, MODEL, CRT
- Related papers