Frequency domain seismic forward modelling: A tool for waveform ...
Frequency domain seismic forward modelling: A tool for waveform ...
Frequency domain seismic forward modelling: A tool for waveform ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
that may fail to adequately modelthewave<strong>for</strong>ms in complex media. Asymptotic ray<br />
theory assumes a high-frequency wave behaviour; this puts certain constraints on the<br />
model complexity asa function of the lowest wavelength. If velocity discontinuities<br />
are reached, explicit boundary conditions must be applied in order to divide the ray<br />
into reected and transmitted (i.e., refracted) rays, each of which are further traced<br />
through the model. The high frequency restriction limits the use of the technique<br />
to simple models with relatively few data phases, usually specied in advance. The<br />
second group of <strong>modelling</strong> methods comprise the numerical methods based on partial<br />
dierential or integro-dierential wave equations, without the use of a high frequency<br />
approximation. These methods are usually <strong>for</strong>mulated as nite dierence or nite<br />
element problems. Such wave equation methods equation guarantee the simulation<br />
of all possible phases (within the assumptions built into the initial wave equation).<br />
The generation of mode conversions, reections and refractions is not determined<br />
by the choice of input parameters (as in asymptotic ray theory), but is instead<br />
an integral feature of the <strong>modelling</strong> itself. [An exception to this are the numerical<br />
methods of (Madariaga, 1984), based on matrix propagator methods. However these<br />
methods are usually only available <strong>for</strong> 1-D models]. As a result, relating the phases<br />
in the <strong>seismic</strong> record to individual features in the model may not be straight<strong><strong>for</strong>ward</strong><br />
in complex models.<br />
Wave equation methods can be further sub-divided intoanumber of classes,<br />
depending on the <strong>domain</strong> in which the initial wave equation is solved.<br />
Possible<br />
choices of <strong>domain</strong> include any combination of time/frequency, space/wavenumber<br />
or other <strong>domain</strong>s, such as the , p trans<strong>for</strong>m <strong>domain</strong>. Each <strong>domain</strong> has its own<br />
advantages and disadvantages. For 2-D earth models, time <strong>domain</strong> methods have<br />
dominated the literature.<br />
In contrast, this thesis will be largely concerned with<br />
numerical <strong>modelling</strong> in the frequency-space <strong>domain</strong>. The primary reason <strong>for</strong> this is<br />
that the <strong>modelling</strong> algorithm is tightly coupled to a <strong>for</strong>mal method <strong>for</strong> the automatic,<br />
frequency-space <strong>domain</strong> inversion of <strong>seismic</strong> wave<strong>for</strong>m data. The results obtained by<br />
17