Abstract: Described is a finite-difference methodology for efficiently computing the response from a model subject to changes within sub-volumes with the sub-volume enclosed with an extrapolation surface within an injection boundary using Green's functions to update the injection boundary and transmitting an updated wavefield from the injection boundary back into the altered sub-volume to include higher order reflections from outside of the altered sub-space.

Abstract: A method of processing seismic data that includes including first and second modes of seismic energy where the second mode has been generated by partial mode conversion of the first mode at a boundary face of a layer of the seabed includes the step of cross-correlating a trace acquired at a receiver and including events corresponding to the first mode with a trace acquired at the same receiver and including events corresponding to the second mode. An event in the cross-correlated data that corresponds to partial mode conversion is identified, and the amplitude of this event is normalized, for example relative to the amplitude of the peak in the cross-correlogram at zero time delay. Information about the effects of the static shift produced by the layer and/or about vector infidelity can be arrived from normalized cross-correlograms for receivers in a receiver array.

Abstract: Method for computing one or more timing errors that occur in one or more multiples predicted by a multiple prediction algorithm. The method includes the steps of generating one or more actual three-dimensional primary travel times and one or more actual three-dimensional multiple travel times, applying the multiple prediction algorithm to the actual three-dimensional primary travel times to generate one or more travel times for the multiples predicted by the multiple prediction algorithm, and subtracting the actual three-dimensional multiple travel times from the travel times for the multiples predicted by the multiple prediction algorithm.

Abstract: Method and apparatus for predicting surface multiples, which includes (a) selecting a target subsurface line (SSL); (b) selecting an input SSL within an aperture of the target SSL; (c) selecting a point on a line twice the distance between the input SSL and the target SSL, the point corresponding to a potential downward reflection point of the surface multiples for a trace; (d) generating a potential surface multiple for the trace corresponding to the point; (e) repeating steps (c) through (d) for each point on the line to generate an inline of potential surface multiples corresponding to each point on the line; (f) repeating steps (b) through (e) for each input SSL within the aperture of the target SSL to generate potential surface multiples for the trace corresponding to each input SSL within the aperture; and (g) adding the potential surface multiples to generate a surface multiple for the trace.

Abstract: A method of determining seismic properties of a layer of the seabed, in particular a surface or near-surface layer (5), comprises directing seismic energy propagating in a first mode at a boundary face of the layer so as to cause partial mode conversion of the seismic energy at the boundary face. For example partial mode conversion may occur when seismic energy propagates upwards through the interface between a surface or near-surface layer (5) and the basement (6), owing to the difference in seismic properties between the surface or near-surface layer (5) and the basement (6). In the invention, the two modes of seismic energy—that is the initial mode and the mode generated by mode conversions at the interface—are received at a receiver. The difference in travel time of the two modes between the interface and the receiver is determined from the seismic data acquired by the receiver.

Abstract: A method of controlling a streamer positioning device (18) configured to be attached to a marine seismic streamer (12) and towed by seismic survey vessel (10) and having a wing and a wing motor for changing the orientation of the wing. The method includes the steps of: obtaining an estimated velocity of the streamer positioning device, calculating a desired change in the orientation of the wing using the estimated velocity of the streamer positioning device, and actuating the wing motor to produce the desired change in the orientation of the wing. The invention also involves an apparatus for controlling a streamer positioning device including means for obtaining an estimated velocity of the streamer positioning device, means for calculating a desired change in the orientation of the wing using the estimated velocity of the streamer positioning device, and means for actuating the wing motor to produce the desired change in the orientation of the wing.

Abstract: A deflector for offsetting a streamer being towed by its lead-in by a seismic survey vessel comprises a wing-shaped body which is coupled to the lead-in and shaped to produce a sideways force which urges the lead-in laterally with respect to the direction of movement of the vessel. A boom extends rearwardly from the wing-shaped body, and forms part of an arrangement by which the angle of the wing-shaped body in the water, and therefore the sideways force produced by it, can be varied. The upper end of the wing-shaped body is directly connected to the underside of an elongate float.

Abstract: A method of controlling a streamer positioning device configured to be attached to a marine seismic streamer and towed by a seismic survey vessel and having a wing and a wing motor for changing the orientation of the wing. The method includes the steps of: obtaining an estimated velocity of the streamer positioning device, calculating a desired change in the orientation of the wing using the estimated velocity of the streamer positioning device, and actuating the wing motor to produce the desired change in the orientation of the wing. The invention also involves an apparatus for controlling a streamer positioning device including means for obtaining an estimated velocity of the streamer positioning device, means for calculating a desired change in the orientation of the wing using the estimated velocity of the streamer positioning device, and means for actuating the wing motor to produce the desired change in the orientation of the wing.

Abstract: A method of performing a seismic survey of a hydrocarbon reservoir in the earth formations beneath a body of water includes deploying a seismic cable from a drum carried by a remotely operated vehicle on the seabed. The cable is deployed into a lined trench, which is formed either concurrently with cable deployment or during a previous survey, to ensure good repeatability of successive surveys of the reservoir, in order to enable changes in characteristics of the reservoir, eg due to depletion, to be monitored.

Abstract: A method of de-ghosting seismic data acquired in a marine seismic survey comprises estimating ?zp (where p is the fluid pressure) from a time-derivative of p and a derivative of p in a horizontal direction such as the x-direction. The time-derivative of p and the derivative of p in the horizontal direction can be readily obtained from the acquired seismic data, and are used to provide an improved estimate of ?zp which, in general, cannot be directly determined from the acquired seismic data. Once ?zp has been estimated it can be used to de-ghost the acquired seismic data.

Abstract: The invention provides a method of determining the orientation of a seismic receiver (4, 6) from seismic data acquired at the receiver (4, 6). The determined orientation can be taken into account in subsequent analysis of the seismic data. This avoids inaccuracies that can occur if the receiver orientation is estimated. In a preferred embodiment the horizontal spatial derivatives of the pressure measured at the receiver are used in the determination of the orientation of the receiver. These may be calculated on either the source side or the receiver side. These are combined with horizontal components of the particle displacement, velocity or acceleration (or a higher time-derivative of the particle displacement). Alternatively, horizontal spatial derivatives of the particle displacement measured at the receiver may be used in place of the horizontal spatial derivatives of the pressure.

Abstract: A method of performing a seismic survey of a hydrocarbon reservoir in the earth formations beneath a body of water includes deploying a seismic cable from a drum carried by a remotely operated vehicle on the seabed. The cable is deployed into a lined trench, which is formed either concurrently with cable deployment or during a previous survey, to ensure good repeatability of successive surveys of the reservoir, in order to enable changes in characteristics of the reservoir, eg due to depletion, to be monitored.

Abstract: A method of processing data sequences obtained at substantially the same time at locations spatially separated from one another is disclosed. A sequence indicative of a measure of spread of the data sequences or of selected data sequences is determined. A common signal data sequence may also be determined from the (selected) data sequences. The measure of spread of the data sequences is preferably normalized, for example relative to the absolute value of the common signal data sequence. The measure of spread is a measure of the noise to signal ratio in the initially-obtained data sequences. It may be used to control the parameters of other processing steps performed on the data sequences. Alternatively or additionally, the measure of spread may be output, for example for monitoring by an observer. The method can be applied to the processing of seismic data obtained using a single sensor seismic data acquisition system.

Abstract: The present invention provides a method and apparatus for analyzing seismic data. The method includes determining a portion of the seismic data recorded substantially before the arrival of a seismic signal and adding a test signal to the portion of the seismic data recorded substantially before the arrival of the seismic signal.

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 for performing amplitude variation with offset (AVO) analysis of a plurality of seismic data traces. The method includes fitting at least a two-term AVO equation to at least three seismic data traces having small angles of incidence using a curve fitting technique to generate an AVO intercept and an AVO gradient, computing a plurality of synthetic seismic data traces using the AVO intercept and the AVO gradient, subtracting the synthetic seismic traces from the plurality of seismic data traces to generate a plurality of higher-order seismic data traces characterized by a residual AVO equation, and fitting the residual AVO equation to the higher order seismic data traces having large angles of incidence using the curve fitting technique to generate a higher order AVO attribute.

Abstract: A staggered vertical marine seismic source contains upper and lower arrays (10, 11) of emitters of seismic energy (S11, S12; S21, S22). The upper array (10) is horizontally displaced relative to the lower array (11). The source is used in a marine seismic surveying arrangement that has means for moving the source and at least one seismic receiver. In use, the source is moved through the water in a direction parallel to the direction in which the two arrays are displaced. The arrays (10, 11) are fired sequentially, and the time delay between the firing of the first-fired array and firing of the second-fired array is chosen such that each seismic emitter in one array is fired at the same x- and y-co-ordinates as the corresponding emitter in the other array. The seismic wavefields generated by firing the two arrays are thus generated at the same x- and y-co-ordinates, but at different depths.

Abstract: A method of generating a semblance panel comprises summing two or more gathers of traces. The dip of the reflector is taken into account in the summation process, and this prevents semblance peaks from becoming smeared during the summation process and so allows a greater number of gathers to be used to generate the semblance panel. In an embodiment of the invention the dip is determined from the seismic data. Alternatively, the reflector dip used in the summation process may be obtained from pre-existing data acquired at the survey location, or the reflector dip may be already known. The invention can be applied to seismic data containing events from more than one reflector.

Abstract: A method of processing seismic data using a seismic energy propagation model of the subsurface is disclosed. The method assigns seismic source, seismic receiver, and reflection point locations to the propagation model; identifies a plurality of alternative raypaths consistent with the propagation model that originate at said seismic source location, reflect at said reflection point location, and terminate at said seismic receiver location; selects a raypath from the plurality of alternative raypaths having a shortest ray length, and utilizes the selected raypath in subsequent seismic processing.

Abstract: The present invention provides marine vibroseis methods and apparatus that operate in an acquisition time frame significantly shorter than current marine vibroseis acquisition methods and comparable to marine airgun seismic source acquisition, while providing data seismic energy comparable or superior to marine airgun seismic sources, by apportioning a desired vibroseis bandwidth signal over a plurality or multiplicity of vibroseis projectors (sources).