Next: 1-D Model Regularization: Discontinuity Up: Brown and Clapp: Integrated Previous: Introduction

# Least Squares Formulation

The fundamental data unit of this paper is the time/depth pair'', which is quite simply the traveltime of seismic waves to a specified depth, along an assumed vertical raypath. We denote time depth pairs by (), indexed by p, the pair'' index. The output velocity function is linear within layers. We denote layer boundaries, which are independent of the time/depth pairs, by tl, indexed by l, the layer'' index.

We begin by deriving the 1-D data residual - the depth error between the depth component of a time/depth pair () and its time component () after vertical stretch with the (unknown) interval velocity:
 (1)
For implementation purposes, we break the integral into the sum of integrals between neighboring time/depth pairs ():
 (2)
We assume that the interval velocity in layer l is linear,
 (3)
so the integral in equation (2) has a closed form. To obtain a correspondence between time/depth pairs and layer boundaries, note that, given a time/depth pair, we can always determine in which layer it resides. In other words, we can unambiguously write l as a function of p, l[p]. Now we can evaluate the integrals of equation (2):
 (4)
Equation (4) defines the misfit for a single time/depth pair, as a function of the model parameters. Now pack the individual misfits from equation (4) into a residual vector, :
 (5)
The elements of vector are the time/depth pair depth values, is the summation operator suggested by equation (4). and is the unknown vector of intercept and slope parameters. The primary goal of least squares optimization is to minimize the mean squared error of the data residual, hence the familiar fitting goal () notation.

Next: 1-D Model Regularization: Discontinuity Up: Brown and Clapp: Integrated Previous: Introduction
Stanford Exploration Project
4/29/2001