Previous: Symmetrizing the wave field
Up: Multi-component Source Equalization
Next: Results
Previous Page: Symmetric data acquisition
Next Page: Results

The optimization problem

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 spatial DC component. That DC component corresponds to zero eigenvalues. 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.

represents the complete misfit of filtered data samples; The indices and range over all different source and reciprocal source locations. The indices and (source and reciprocal source components) both stay within the range . The values of and may describe source location in 2-D or 3-D. Convolution is denoted by . For each source location, the general shape of in is

For pressure data, the index range for and 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 and is applied. In equation 5 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 is then a single scalar value for each separate minimization process. represents the objective function to be minimized.