This is an old revision of the document!

# GP211 - Labs

Homework Resources
Program Documentation: Many of the assignments use Sergey Fomel's Fortran90 inversion library. Short explanations are provided for all main programs.
SEP I/O module: Sergey Fomel has written a Fortran90 Input/Output (I/O) module which allows you to conveniently read, write, and manipulate SEP-format datafiles.
Fortran90 FAQ: Although you are not expected to be an expert programmer, you will modify Fortran90 programs in nearly all assignments. If you have programming questions, contact the TA, or better yet, persue this list of Frequently Asked F90 Questions.
 Assignment List Note: all the code from the book, and subroutines you will link to in some of your programs are available in ~prof/gee/Src and ~prof/gee/Ops. I suggest setting up some links when you're working on the labs.
 Missing Data Estimation and Multiple Realizations Lab4.pdf Geostatistics has gained much recent attention for its ability to solve missing data problems not with a single answer, but with a suite of so-called equiprobable models, or models that fit the data equally well. The geostatistical approach gives better results in fluid flow simulation (for instance), and also provides a useful tool for multiple scenario analysis. Clapp (SEP-105) showed that multiple realizations can easily be incorporated into the existing SEP framework to solving missing data problems. This lab has two parts. In the first, you will modify a conjugate gradient algorithm to allow multiple realizations, and then test it on the 1-D inverse linear interpolation problem given in Chapter 3 of the book. In the second part, you will use your new solver to tackle the Wells not matching the seismic map problem of Chapter 3. You will improve the ``extension'' of the seismic data into empty regions of the model space, first by writing a program to de-trend the seismic data with a least squares best-fitting plane, and second, by using your new solver to come up with multiple realizations of the extended seismic data. 10/25 Lab4.tar
 Interval Velocity Estimation Lab5.pdf This lab will guide you through most of the steps - and pitfalls - of interval velocity estimation, perhaps the most important problem in exploration geophysics. You will apply the least squares velocity scan of Lab 2, and then pick RMS velocities using an autopicker and by hand. You will apply the multiple realizations concept to bound uncertainty on your picked RMS velocity. You will solve the Dix equation as a least squares problem to estimate interval velocity in 1-D and make some preliminary steps in the direction of 2-D. 10/29 Lab5.tar
 Debluring by inversion with regularization Lab6.pdf In this Lab, we test different regularization schemes for deblurring a text image with random noise. Regularization is very important in inversion when there are not enough data. Prior information can be included in the regularization to make the inversion converge fast. In this Lab, you will code the adjoint of the given operator. You will then apply two regularization operators to the linear inversion. In the final part, you will apply an edge-preserving regularization scheme to the non-linear inversion. 11/08 Lab6.tar
 Madagascar Lab7.pdf The Madagascar satellite data set provides images of a spreading ridge off the coast of Madagascar. This data set has two regions: the southern half is densely sampled and the northern half is sparsely sampled. The sparsely sampled region presents a missing data problem. In this exercise, you will apply basic fitting goals to image the dense area. Then, you will be asked to fill the missing data in the sparsely sampled region. Additionally, you will use prediction-error filters estimated on the dense southern half to fill data on the sparse half. Preconditioning will be used to speed convergence. The prediction-error filters will attempt to spread the texture of the spreading ridge to the sparse tracks. 11/15 Lab7.tar
 Model space residual Lab8.pdf Steepest descent minimizes the residual in the data space. Alternatively, we can try to minimize the residual in the model space. This gives opportunities of spatially variable step length. 11/21 Lab8.tar

 SEP Stanford Exploration Project GP Department of Geophysics SES School of Earth Sciences Stanford Stanford University

Page Maintainer: webmaster@sep.stanford.edu