Abstract: An embodiment of the invention relates to a method of processing seismic data comprising a plurality of traces, each trace having a source position, a receiver position and a midpoint position, each having coordinates in an inline direction and a crossline direction. The invention uses a data set comprising traces having one of the source position, the receiver position or the midpoint position having coordinates which are invariant in one of the inline or the crossline directions. A non-iterative surface related multiple attenuation algorithm is applied to the data set which assumes substantial non-variation of depth of reflectors in the same one of the inline or the crossline directions.

Abstract: An embodiment of the invention relates to a method of processing seismic data comprising a plurality of traces, each trace having a source position, a receiver position and a midpoint position, each having coordinates in an inline direction and a crossline direction. The invention uses a data set comprising traces having one of the source position, the receiver position or the midpoint position having coordinates which are invariant in one of the inline or the crossline directions. A non-iterative surface related multiple attenuation algorithm is applied to the data set which assumes substantial non-variation of depth of reflectors in the same one of the inline or the crossline directions.

Abstract: A mehtod of determinig an integral of a function over an n-dimensional intergratio domain from the values of the function at a plurality of discrete data points within the integration domain, comprises partitioning the integration domain into a plurality of k-dimensional simplexes, where k<=n. The function is integrated over each simplex, and the results are summed. The method may be applied in the processing of geophysical data such as, for example, seismic data, where the function is a function of at least the geophysical data. The processed data may then be used to obtain information about the physical properties of the earth's interior.

Abstract: A number of seismic stacks are precomputed (20) for known velocity fields. The velocity fields are chosen to span the range of velocities of interest. The stacks are then arranged (21) in the 3D memory of a graphics computer (10–14) using time and position as first dimensions and the index of the velocity field as the last dimension. In such 3D space, any velocity field to be used for stacking appears as a surface (S) within a volume. Projecting the seismic stacks onto that surface provides the seismic line stacked for the velocity field of interest.

Abstract: A method of determining an integral of a function over an n-dimensional integration domain from the values of the function at a plurality of discrete data points within the integration domain, comprises partitioning the integration domain into a plurality of k-dimensional simplexes, where k?n. The function is integrated over each simplex, and the results are summed. The method may be applied in the processing of geophysical data such as, for example, seismic data, where the function is a function of at least the geophysical data. The processed data may then be used to obtain information about the physical properties of the earth's interior.

Abstract: A number of seismic stacks are precomputed (20) for known velocity fields. The velocity fields are chosen to span the range of velocities of interest. The stacks are then arranged (21) in the 3D memory of a graphics computer (10-14) using time and position as first dimensions and the index of the velocity field as the last dimension. In such 3D space, any velocity field to be used for stacking appears as a surface (S) within a volume. Projecting the seismic stacks onto that surface provides the seismic line stacked for the velocity field of interest.