RESEARCH PROGRAMME OF SEISMIC WAVES IN COMPLEX 3-D STRUCTURES
The research is focused primarily on the fundamental issues of
high-frequency seismic wave propagation in complex 3-D isotropic
and anisotropic structures, which go beyond the traditional approaches.
The ray method and its extensions, as well as its combination
with other methods are mainly applied and investigated.
The emphasis is put on new, stable, more efficient and flexible
algorithms for both forward numerical modelling and inversion
of seismic wave fields in 3-D inhomogeneous, isotropic
and anisotropic structures.
Considerable attention is also devoted to applications involving
S waves, converted waves, S-wave splitting and coupling in
anisotropic media, particle ground motions, etc.
Much more detailed information can be obtained at
"http://sw3d.mff.cuni.cz".
The research programme was begun on October 1, 1993.
RESEARCH PROGRAMME FOR THE NINTH YEAR
October 1, 2001 - September 30, 2002
1. Program packages
Package MODEL:
Model: General 3-D layered and block isotropic or anisotropic
structures, containing isolated bodies, pinchouts, etc.
Inside the layers and blocks, the elastic parameters may vary
in all three dimensions.
Dissipation and non-planar topography can be considered.
Possibility of model smoothing, data fitting by inversion
including fitting and smoothing GOCAD models,
conversion of model parametrization, structural interface triangulation,
VRML and GOCAD visualization.
Package CRT:
Model: Using package MODEL.
Type of waves: Arbitrary type of elementary seismic body wave
corresponding to the zero-order ray theory (P,S, converted).
Computations: Arbitrary position and shape of the source,
initial-value ray tracing by numerical integration of ray equations,
two-point ray tracing by the shooting method, travel-time computation,
dynamic ray tracing, paraxial-ray propagator matrix,
geometrical spreading, vectorial amplitudes, polarization vectors.
The package may be applied to the
evaluation of the elastodynamic ray-theory Green function, and to
the computation of synthetic seismograms,
including the coupling ray theory for weak anisotropy and the
response of fine layers at receiver sites
(program package RMATRIX by C.J. Thomson, linked to the CRT package).
Least-square travel-time tomography with smoothing using Sobolev
scalar products.
Aquisition schemes: Surface seismics (land and marine),
VSP, cross-hole, OBS, OBC.
Planned innovations:
Anisotropic common-ray approximation of the coupling ray theory for S waves.
The package will be extended to solve various inverse problems,
travel-time tomography in particular.
Package ANRAY:
Model: 3-D laterally varying structures containing isotropic and
anisotropic non-vanishing layers. Specification of elastic parameters
inside individual layers either by linear interpolation between
isosurfaces of elastic parameters, or by B-spline interpolation
within a 3-D rectangular grid of elastic parameters.
Possibility of VRML visualization.
Types of waves: Arbitrary type of elementary seismic body wave
(P, S, qP, qS1, qS2, any converted wave, coupled qS waves).
Computations: Numerical integration of ray tracing and dynamic
ray tracing equations, calculation of ray vectorial amplitudes,
ray-theory Green function including the Green function in the quasi-isotropic
approximation for qS waves, ray synthetic seismograms, particle
ground motions.
Aquisition schemes: Surface seismics (land and marine),
VSP, cross-hole, OBS, OBC.
Planned innovations:
(a) Study of possibilities to increase accuracy of the
quasi-isotropic (QI) approach by increasing accuracy of the travel
time approximation; testing the QI approach by comparing the QI
results with results of more precise methods in specific situations,
like vicinity of singularities, transition from anisotropy to
isotropy, reflection/transmission in weakly anisotropic media. (b)
Modifications enabling calculation of Gaussian beams. (c) Alternative
ray tracing with RHS specified by the polarization vectors.
(d) GOCAD visualization.
Package NET:
Model: Using package MODEL or using gridded velocities.
Types of waves: First arrivals, constrained first arrivals.
Computations: Arbitrary position and shape of the source.
First-arrival travel times in the whole model are computed.
The algorithm of computation is independent of the model's complexity.
Aquisition schemes: Surface seismics (land and marine),
VSP, cross-hole, OBS, OBC.
Package FD:
Model: Using package MODEL.
Type of waves: Complete elastic wave field.
Computations: Presently 2-D, without fluids.
Aquisition schemes: Surface seismics (land).
Package FORMS:
Computations:
Subroutines used by other program packages including
data input and output subroutines,
management and plotting of synthetic seismograms,
2-D and 3-D graphics including 3-D virtual reality
with VRML and GOCAD visualization,
manipulation and calculation with gridded data (data cubes),
programs for matrix and vector operations necessary for inversion,
other general-purpose seismic software.
2. Sample data for the program packages
Examples of the input data describing or approximating models
delivered by the consortium members or other typical models
will be prepared.
Examples of the input data to perform calculations in such
models will also be prepared.
3. Two-point ray tracing in complex isotropic 3-D structures
Properties of the projection of ray coordinates onto Cartesian
coordinates within individual ray histories will be studied. New
applications of ray histories to numerical algorithms will be proposed.
The two--point ray tracing code will further be tested and applied
to various models.
4. Paraxial Fresnel edge waves
Contributions of single diffractions to the Green function may be
approximated by paraxial expansions of Fresnel edge waves in
vicinities of boundary rays between ray histories.
5. Synthetic seismograms in 3-D isotropic complex structures
Methods to calculate synthetic seismograms in complex structures will
be studied, mutually compared and combined.
Emphasis will be put on numerical implementation of Gaussian beams.
6. Seismic wave propagation in weakly anisotropic inhomogeneous media
Derivation of anisotropic common--ray dynamic ray tracing equations.
Derivation of coupling ray theory from the elastodynamic equation
concentrated on the study of errors due to neglected terms.
Study of the inaccuracy of the coupling ray theory due to the
differences between the reference and exact rays.
Estimation of frequencies at which the quasi-isotropic common-ray
or anisotropic common-ray approximation of the coupling ray theory
should be replaced by the anisotropic ray theory, or by the coupling
ray theory along the anisotropic ray theory rays.
7. Applicability of the high-frequency asymptotics in the
vicinity of S-wave singularities
Anisotropic S-wave ray tracing in a vicinity of singularities
in inhomogeneous anisotropic media. Investigation of accuracy
of various high-frequency approximations (zeroth-order anisotropic
ray approximation, quasi-isotropic approximation, higher-order
anisotropic ray approximations, Gaussian beams, etc.) for S waves,
propagating in strongly or weakly anisotropic homogeneous or
inhomogeneous media in the vicinity of singularities, by comparison
with other more precise methods (exact solutions, finite differences,
etc.).
8. R/T coefficients for inhomogeneous waves in dissipative elastic media
Computation of R/T coefficients at structural interfaces and
thin transition layers between two isotropic dissipative media will be
tested further. The main attention will be devoted to
frequency-independent reference R/T coefficients.
9. Computation of ray-theory travel times, amplitudes
and other quantities at the nodes of 3-D grids
Algorithms of fast calculation of ray-theory travel times in dense
rectangular grids will be investigated further.
The accuracy and efficiency of the interpolation of ray-theory
travel times within ray cells in 3-D models will be studied further,
and the relevant numerical algorithms will be improved,
or new ones will be proposed.
If possible, attention will also be devoted to the interpolation
between different shot and receiver positions.
10. Accuracy of seismic modelling
The research will be concentrated mainly on the accuracy of travel-time
calculations,
on the accuracy of finite-difference modelling of seismic wave fields,
and on the accuracy of other modelling methods designed or studied
within the framework of the project.
The main attention will be devoted to the estimation of the
feasibility and costs of ray tracing,
and to the high-frequency validity of velocity models.
11. Lyapunov exponents and model smoothing
Construction and smoothing of velocity macro models will further
be studied, with emphasis on the application of Sobolev scalar
products and Lyapunov exponents.
Attention will be paid to a possible extension of the estimation
of Lyapunov exponents to smooth 3-D models and to 2-D models with
structural interfaces.
12. Seismic tomography
Development of theory, algorithms and programs applicable
in seismic travel-time tomography and inversion of the coherency
panels, with emphasis on the estimation of the accuracy of the
resulting model compared to the geological structure.
Further development of algorithms for seismic travel-time tomography in
anisotropic media, synthetic tests of possibility to discriminate
anisotropy from inhomogeneity. Extension to other types of data
acquisition, extension to shear wave data, study of effects of noise.
The theoretical investigation will also be aimed at the determination
of the resolution of the elastic parameters with respect to the
measurement geometry.
13. Local anisotropy parameter estimation from VSP measurements
Further development of algorithms for the local determination of
elastic parameters from the VSP measurements.
Incorporation of qS-wave data.
Testing various observation configurations to discriminate
between various anisotropy symmetries of studied media.
Inversion for parameters of TI media and for angles of axis
of symmetry in case of the inclined TI symmetry.
14. Decomposition of a wave field into Gabor wavelets
Discretized integral decomposition of a spatial wave field
at a fixed time or a time-dependent wavefield along a smooth
surface into Gabor wavelets in 2-D or 3-D.
The Gabor wavelets may be frequency-dependent and their
shape may be smoothly varying in space and time, which
requires much more general decomposition than the integral
Gabor transform.
The decomposition into Gabor wavelets may be useful for the
decomposition of a general wave field into Gaussian beams or packets.
15. Migrations
Resolution and accuracy of migrations will be studied. Attention will
be paid to the physical meaning of the migrated sections and to their
sensitivity to the velocity model, including its anisotropy.
Algorithm of the Gaussian-packet prestack depth migration
will be developed.
Gaussian packets should be very efficient and offer explicit
correspondence between the time and depth sections.
Attention will be paid to the optimization of the shape of
Gaussian packets.
16. Gaussian-packet "finite elements"
In the beginning of each time interval of a given duration,
the wave field may be decomposed into Gaussian packets
propagating throughout the corresponding time interval and
yielding the wave field in the beginning of the next time interval.
17. Finite-difference solutions of elastodynamic equations
The main attention will be paid to the hybrid DW--FD and Ray--FD methods.
18. Concluding remarks
In addition to this programme, we will certainly be responsive to
specific technical suggestions and recommendations of the consortium
members within the general framework of the project.
The research in most directions listed above will continue into the
future years of the project.
You may download PostScript file
prog02.ps (65 kB) with the Research Programme.
SW3D
- main page of consortium Seismic Waves in Complex 3-D Structures .