MULTIPLE SUPPRESSION FOR OCEAN BOTTOM SEISMIC DATA

The invention is a method for removing trapped water bottom multiples (figure 2, item 19), receiver side peg-leg multiples (13), and source side peg-leg multiples (15) from dual sensor (figure 1, item 6), OBC (7) data, where the data includes both pressure signals and velocity signals. The pressure and velocity signals are compared to determine any polarity reversals between them (figure 3, item 21). Polarity reversals are used to identify and separate up-going and down-going wavefields in the pressure and velocity signals (22). A matching filter is applied to a portion of the velocity signal where polarity reversals exist (24). The down-going wavefield is then estimated by calculating the difference between the portion of the velocity signal where polarity reversals exist and the portion of the pressure signal where polarity reversals exist and applying a scaling factor to the result (26). An attenuated up-going pressure wavefield is then determined by combining the estimated down-going wavefield and the pressure signal (28). L'invention concerne un procédé d'élimination des réflexions multiples piégées dans la colonne d'eau, des réflexions multiples subocéaniques côté récepteur et des réflexions multiples subocéaniques côté source, présentes dans des données obtenues avec des câbles de fonds à capteurs doubles, ces données comprenant à la fois des signaux de pression et des signaux de vitesse. Les signaux de vitesse et de pression sont comparés afin de déterminer toute inversion de polarité entre eux. Les inversions de polarité sont utilisées afin d'identifier et de séparer des champs d'ondes montants et descendants dans les signaux de pression et de vitesse. Un filtre d'adaptation est appliqué à la portion du signal de vitesse dans laquelle existent des inversions de polarité. Le champ d'onde descendant est alors estimé par calcul de la différence entre la portion du signal de vitesse où les inversions de polarité existent et la portion du signal de pression où les inversions de polarité existent, puis par application d'un facteur d'échelle au résultat. On détermine alors un champ d'onde de pression montant par combinaison entre le champ d'onde descendant estimé et le signal de pression.