![next](icons/next.png) |
![up](icons/up.png) |
![previous](icons/previous.png) |
![](icons/left.png) | Wave-equation traveltime tomography by global optimization | ![](icons/right.png) |
![[pdf]](icons/pdf.png) |
Next: Synthetic example
Up: Almomin: Tomography
Previous: Introduction
The objective function of FWI can be written as follows:
![$\displaystyle \mathbf J_{FWI}(\mathbf v) = \sum_{\mathbf x_s} \sum_{\mathbf x_g...
...athbf x_s;\mathbf v)-\mathbf d_{obs}(t,\mathbf x_g,\mathbf x_s)\vert\vert _2^2,$](img2.png) |
(1) |
where
and
are the source and receiver locations,
is the modeled data with velocity
, and
is the observed data. By setting the first derivative of equation (1) around the velocity
to zero, the velocity update can be expressed as follows:
![$\displaystyle \Delta \mathbf v = \frac{-s}{\mathbf v_0^3}\sum_{\mathbf x_s} \su...
...\mathbf x_g,\mathbf x_s;\mathbf v)-\mathbf d_{obs}(t,\mathbf x_g,\mathbf x_s)),$](img9.png) |
(2) |
where
is the step size,
is the source signature, and
and
are the forward wave propagation operator and its adjoint, respectively.
The objective function wave-equation traveltime tomography can be written as follows:
![$\displaystyle \mathbf J_{\Delta \tau}(\mathbf v) = \sum_{\mathbf x_s} \sum_{\mathbf x_g} \vert\vert\Delta \tau(\mathbf x_g,\mathbf x_s;\mathbf v)\vert\vert _2^2,$](img14.png) |
(3) |
where
is the lag of the maximum cross-correlation between the observed data and the data modeled by a velocity model
. Again, the first derivative of equation (3) around the lags
is set to zero to get the velocity update, which can be expressed as follows:
![$\displaystyle \Delta \mathbf v = \frac{s}{\mathbf v_0^3}\sum_{\mathbf x_s} \sum...
...c{\partial}{\partial t} \mathbf d_{obs}(t+\Delta \tau,\mathbf x_g,\mathbf x_s),$](img16.png) |
(4) |
where
is defined as follows:
![$\displaystyle \xi = \sum_{t} \frac{\partial}{\partial t} \mathbf d_{cal}(t,\mat...
...c{\partial}{\partial t} \mathbf d_{obs}(t+\Delta \tau,\mathbf x_g,\mathbf x_s),$](img18.png) |
(5) |
By examining equations (2) and (4), it can be shown that (WT) can handle much larger velocity errors than (FWI).
Now, I cast the picking procedure of the lags
as a global optimization problem with an objective function as follows:
![$\displaystyle \mathbf C(\Delta \tau) = \sum_{\mathbf x_s} \sum_{\mathbf x_g} f (\mathbf A \Delta \tau(\mathbf y_g,\mathbf y_s)) ,$](img19.png) |
(6) |
where
and
are a sparse representation of the source and receiver locations,
is a bicubic spline interpolation operator that maps the sparse coordinates
and
to the original coordinates
and
, and
evaluate the correlation value at
.
The goal of the global optimization is to maximize the function described by equation (6), which is to maximize the stacking power along the interpolated spline surface. The searching procedure is a simulated annealing algorithm, which varies the spline points along the time axis in a stochastic sense until a satisfying solution is reached. In the following section, I show the results of using such global scheme to pick the correlation lags.
|
---|
velali
Figure 1. The true velocity model used to create the data.
|
---|
|
---|
![next](icons/next.png) |
![up](icons/up.png) |
![previous](icons/previous.png) |
![](icons/left.png) | Wave-equation traveltime tomography by global optimization | ![](icons/right.png) |
![[pdf]](icons/pdf.png) |
Next: Synthetic example
Up: Almomin: Tomography
Previous: Introduction
2010-05-19