I have set up a least-squares optimization problem in which I minimize the difference between data and its reciprocal counterpart. Analyzing the solution of that problem shows that none of the eigenvectors has a pure DC component. Thus the absolute radiation pattern cannot be determined. However, relative changes are perfectly resolvable within a cable length. This relativity requires us to specify a reference with respect to which the prediction error filters are computed. There are a multitude of choices. One possibility would be to average the sources over some range and establish the average as a reference source. The disadvantage of this method is that a relatively noisy (bad) common source gather biases the estimation procedure unfavourably. Another choice would be to design or give some training data as reference, but that would have the effect of wavelet shaping. Instead, I took the first source as a reference source, all other filters show differences with respect to that source. My initial guess is a single spike at the zero lag filter location. In the minimization process the filters are not constrained at all (except for being specified as to the number of filter coefficients). Once I decided on the filter shape (1-D, short), I tested the following objective functions governing the minimization process.

First, we find the global minimum of the squared difference for all the components at the same time. This requires simultaneous accessibility of all 9 component data. For each iteration that data must be read once.

(5) |

*E _{1}* represents the complete
misfit of filtered data samples;

For pressure data, the index range for *l* and *k*
degenerates to just being unity.
For vector wave field data these indices
indicate the type of component (vertical, horizontal, inline, crossline).
To each ensemble a filter pair *s*_{i}^{kl} and *s*_{j}^{lk} is applied.
In equation (5)
*s* is estimated, for each source component pair, as a short function in
time.

Instead of equalizing all components simultaneously, we can treat
each reciprocal component pair separately. The advantage is smaller
data sets to be worked on. Equation (6) gives
the appropriate objective function. The filter *s*^{kl} is then a
single scalar value for each separate minimization process.
*E _{2}*

(6) |

A third minimization alternative, however, is to find a global minimum for a filter that consists of a scalar array multiplied by a source location filter, as follows

(7) |

11/18/1997