Figure 7 Hyperbolic Radon transform model domain of shot 08 from the Yilmaz data collection. 20 iterations of least-squares inversion were performed. 14 coordinates selected by LA are overlain. Without inclusion of the inappropriate high slowness values on the right side of the plot, the remaining picks remain stable.
I introduce Lloyd's algorithm as a tool to optimally represent the statistics of the model space from incomplete inversions. The modified algorithm selects coordinates with high amplitude surrounded by substantial energy. Thus isolated, powerful outliers are neglected. By optimally parsing potentially large multi-dimensional model spaces, the algorithm can cut short costly inversion iterations and focus an interpreter's attention to important locations within potentially large model domains. The algorithm returns stable solutions even in the presence of substantial noise.
The algorithm is very simple, easy to modify, has few parameters, and very fast. Using the algorithm depends on parameters being correlable. For multidimensional cases, uncorrelated parameters can simply be concatenated to an existing axis. Thus hypercubes of correlable and uncorrelable parameters can be evaluated simultaneously.
The next step in evaluating the effectiveness of using the algorithm to select optimal parameters would be to migrate data with a velocity model derived from RMS velocities selected by LA. This could potentially dovetail with the velocity uncertainty analysis presented in ().
Planewave decomposition of passive data to characterize non-obvious sources does not work. Until traces have been correlated, analysis transforms will simply redistribute the random character of the raw data. Unfortunately, correlating the wavefield destroys all the unique character of individual sources including timing, waveform, location and much of the spectral content information.
Thanks to Bob Clapp for the introduction of LA to the group and good discussions with Jeff Shragge on potential uses.