Page 224 - Acquisition and Processing of Marine Seismic Data
P. 224

4.3 CONVOLUTION                                215

           FIG. 4.2  Schematical illustration of autocorrela-
           tion calculation of a time series x(t)consisting offour
           discrete elements. Mutual elements of the series in
           yellow boxes are multiplied and the results are
           summed. τ is the autocorrelation lag.






           Wiener-Khintchine  Theorem    and  can   be  integration. If x(t) and y(t) functions are of differ-
           expressed as                                 ent time lengths, then the length of R xy (τ) will
                                                        equal to the length of x(t). Calculation of
                                                  2
              ℑ R xx τðÞg ¼ R xx ωðÞ ¼ X ωðÞ   X ωðÞ ¼ X ωðÞj
               f
                                             j
                                                        crosscorrelation of x(t) and y(t) time series in
                                                  (4.3)  time domain corresponds to the multiplying of
                                                        amplitude spectra of both functions in fre-
           Eq. (4.3) indicates that taking the autocorrelation
                                                        quency domain, and the phase spectra of both
           of an x(t) function in time domain corresponds
                                                        functions are subtracted. Order of the functions
           to calculating the power spectrum of x(t) in fre-
                                                        in crosscorrelation calculation is important, and
           quency domain.
                                                        R xy (τ) ¼ R yx ( τ).
              Fig.4.2schematicallyshowstheautocorrelation
           calculation of discrete time series. Corresponding  If both functions are similar, then the
           elements of the input series are multiplied and the  crosscorrelation calculation produces a large
           results are summed (τ ¼ 0). Then one of the series  positive number, while smaller correlation
           islaggedonetimesample(τ ¼ 1)andtheprocessis  values are obtained if the functions are not alike.
           repeated, and calculations continue until the last  Inpractice,twotimeserieshavesimilaritiesifone
                                                        is the time-lagged version of another. Hence,
           sample of the series is reached.
                                                        crosscorrelation actually indicates the amount
                                                        of time shift of a function relative to other, and
                 4.2 CROSSCORRELATION                   is used to obtain deconvolution operator in pre-
                                                        dictive deconvolution (Section 6.4). Fig. 4.3 sche-
                                                        matically shows the crosscorrelation calculation
              Crosscorrelation is used to obtain the degree
           of the similarity between two different time  of discrete time series. Calculations are per-
           series. The crosscorrelation of x(t) and y(t) series  formed as for the autocorrelation calculations.
           is expressed analytically as
                      Z ∞             Z ∞                         4.3 CONVOLUTION

                                         ð
                            ð
              R xy τðÞ ¼  xtðÞyt + τÞdt ¼  xt τÞytðÞdt
                      ∞               ∞                    Convolution is the mathematical process
                      Z ∞             Z ∞               betweeninputandoutputofalinearsystem.That
                                         ð
              R yx τðÞ ¼  ytðÞxt + τÞdt ¼  yt τÞxtðÞdt  is, there is a convolution process between input
                            ð
                      ∞               ∞                 and output of a stationary system. For x(t) and
                                                  (4.4)  y(t) time series, it is analytically expressed as
                                                                          Z ∞
           where τ is the time lag. According to Eq. (4.4),
           y(t) analytical function is lagged by an amount          ftðÞ ¼   x τðÞyt τÞdτ      (4.5)
                                                                                 ð
           of τ, multiplied by x(t), and summed up via                    ∞
   219   220   221   222   223   224   225   226   227   228   229