The exercises in light gray color are extra credit and are not required.
Reading: Claerbout, 1992, Chapter 2 and Chapter 3.
cd [folder where you put jrg.jar] java -cp jrg.jar ZPlane &Read about the Zplane program in Appendix E of the text. In this new version, point at the Z-plane field to see the complex coordinates, omega, and rho of the place you are pointing at. To place a Zero, point where you want the zero to be and click the mouse. To place a Pole, hold down the Option key when you click. To remove a Pole or Zero, point at the root you want to remove, hold down the Shift key, and click the mouse on it. You can also point and type the Z or P key to place Zeroes and Poles, and point and type the X key to remove them. You can also pick up and drag roots. To put a root on an axis, you can point just outside the axis.
Type L to load a data set, see its summed spectrum, and filter it. Download this PVIZplaneDataExample.zip ZIP archive, which contains the reflection record shown by the original Zplane program, as a JRG Pack and as SEG-Y. Unzip it. When you press L, you see a blank, zero data record. Click on it and press Command-L (Control-L on Windows) to load a JRG Pack. You can navigate to where you unzipped the archive, open the PVIZplaneDataExample folder that contains the JRG Pack, and double-click on any file inside. The data will open much as it appeared in the original Zplane interface. Select the example data set with the choice tool near the bottom of the ZPlane window.
Turn in the answers to the following questions:
Deconvolve a single-channel seismic-reflection data set recorded in Lake Washington, that crosses the Seattle fault several times. Use both recorded source wavelets and spectral whitening in your trials.
Invert single-channel reflection data from Lake Washington for the reflectivity sequence, given data sets recorded at different ranges of frequency over the same location. Develop a transfer function between the reflectivity and other geophysical measurements of the profile.
Construct a minimum-phase multiple-notch filter for a data set contaminated by 60 Hz noise and harmonics. Filter the data set, and describe artifacts, confidence in the content of the filtered spectrum, and any phase changes in the data.
Implement a flexible time-varying, minimum-phase Butterworth filter that will accomplish band pass or reject operations on aritrarily-sized seismic data. Test it on a real data set.
Investigate the use of the slant-stack operator for estimating missing data. Implement a moving-window routine for computing local slantstacks, and test it by filling gaps in stacks, migrations, shot records, and images.
Investigate the properties of a tree-ring data set. Put quantitative bounds on the confidence of cross-correlations between overlapping data, and on the periodicities of climatic cycles.
Develop a method for computing a spectrum from unequally-spaced data. Explore the advantages and disadvantages of linear and spline interpolation, integrations, and model fitting. Speculate on 2-d spectra from randomly-spaced 2-d data.