\prefacesection{Abstract} 

Wave-equation migration using one-way wavefield extrapolation operators is commonly used in industry to generate images of complex geologic structure from 3D seismic data.   By design, most conventional wave-equation approaches restrict propagation to downward continuation, where wavefields are recursively extrapolated to depth on Cartesian meshes.  In practice, this approach is limited in high-angle accuracy and is restricted to down-going waves, which precludes the use of some steep dip and all turning wave components important for imaging targets in such areas as steep salt body flanks.  
%
\par
%
This thesis discusses a strategy for improving wavefield extrapolation based on extending wavefield propagation to generalized coordinate system geometries that are more conformal to the wavefield propagation direction and permit imaging with turning waves.  Wavefield propagation in non-Cartesian coordinates requires properly specifying the Laplacian operator in the governing Helmholtz equation.  By employing differential geometry theory, I demonstrate how generalized a Riemannian wavefield extrapolation (RWE) procedure can be developed for any 3D non-orthogonal coordinate system, including those constructed by smoothing ray-based coordinate meshes formed from a suite of traced rays.  I present 2D and 3D generalized RWE propagation examples illustrating the improved steep-dip propagation afforded by the coordinate transformation. 
%
\par
%
One consequence of using non-Cartesian coordinates, though, is that the corresponding 3D extrapolation operators have up to 10 non-stationary coefficients, which can lead to imposing (and limiting) computer memory constraints for realistic 3D applications.  To circumvent this difficulty, I apply the generalized RWE theory to analytic coordinate systems, rather than numerically generated meshes.  Analytic coordinates offer the advantage of having straightforward analytic dispersion relationships and easy-to-implement extrapolation operators that add little computational overhead.  In particular, I demonstrate that the dispersion relationship for 2D elliptical geometry introduces only an effective velocity model stretch, permitting the use of existing high-order Cartesian extrapolators.   The results of elliptical coordinate shot-profile migration tests demonstrate the improvements in steep dip reflector imaging facilitated by the coordinate system transformation approach.
%
\par
%
I extend the analytic coordinate system approach to 3D geometries using tilted elliptical-cylindrical (TEC) meshes.  I demonstrate that propagation in a TEC coordinate system is equivalent to wavefield extrapolation in elliptically anisotropic media, which is easily handled by existing industry practice.  TEC coordinates also allow steep dip propagation in both the inline and cross-line directions by virtue of the associated elliptical and tilting Cartesian geometries.  Observing that a TEC coordinate system conforms closely to the shape of a line-source impulse response, I develop an TEC-coordinate, inline-delay-source migration strategy that enables the efficient migration of individual sail-line data.  I argue that this strategy is more robust than 3D plane-wave migration because of the reduced migration aperture requirements and, commonly, a lower number of total migration runs.  Synthetic imaging tests on a 3D wide-azimuth data set demonstrate the imaging advantages offered by the TEC coordinate transformation, especially in the cross-line direction.  Field data tests on a Gulf of Mexico data set similarly indicate the advantage of TEC coordinates.
%; though, this application is complicated by anisotropic sediments unaccounted for in the generalized RWE theory.

% rev. 9/18/07
%\footer{SEP--131}
\prefacesection{Preface}

\noindent
The electronic version of this report\footnote{http://sepwww.stanford.edu/private/docs/sep137} makes the included programs and applications available to the reader.  The markings {\bf ER}, {\bf CR}, and {\bf NR} are promises by the author about the reproducibility of each figure result. Reproducibility is a way of organizing computational research that allows both the author and the reader of a publication to verify the reported results.  Reproducibility facilitates the transfer  of knowledge within SEP and between SEP and its sponsors.
\begin{description}
\item {\bf ER} 
        denotes Easily Reproducible
        and are the results of processing described in the paper.
        The author claims that 
        you can reproduce such a figure from
        the programs, parameters, and makefiles
        included in the electronic document. 
        The data must either be included in the electronic distribution,
        be easily available to all researchers (e.g., SEG-EAGE data sets), or  
        be available in the SEP data 
  library\footnote{http://sepwww.stanford.edu/public/docs/sepdatalib/toc\_html}.
        We assume you have a UNIX workstation
        with Fortran, Fortran90, C, X-Windows system and the software
        downloadable from our website 
        (SEP makerules, SEPScons, SEPlib, and the SEP latex package), or other 
        free software such as SU.
        Before the publication of the electronic document, 
        someone other than the author tests the author's claim by 
        destroying and rebuilding all ER figures. 
        Some ER figures may not be reproducible by outsiders because
        they depend on data sets that are too large to distribute,
        or data that we do not have permission to redistribute but
        are in the SEP data library, or that the rules depend on commercial packages such as Matlab or Mathematica.
\item {\bf CR }
        denotes Conditional Reproducibility.
        The author certifies that the commands are in place to reproduce
        the figure if certain resources are available.
	    The primary reasons for the CR
        designation is that the processing requires 20 minutes or more. 
\item {\bf NR }
        denotes Non-Reproducible figures.
        SEP discourages authors from flagging their figures
        as NR except for figures that are used solely for 
        motivation, comparison, or illustration of the theory,
        such as:  artist drawings, scannings, or figures taken from SEP reports 
        not by the authors or from non-SEP publications.  Some 3D synthetic data examples associated with this report are classified as NR because they were generated externally to SEP and only the results were released by agreement.
\end{description}
Our testing is currently limited to LINUX 2.6 (using the Intel Fortran90 compiler) and the SEPlib-6.4.6 distribution, but the code should be portable to other architectures.  Reader's suggestions are welcome.  For more information on reproducing SEP's  electronic documents, please 
visit {\it http://sepwww.stanford.edu/research/redoc/}. 

\cleardoublepage
