WO1997044685A1 - Method for attenuation of reverberations using a pressure-velocity bottom cable - Google Patents

Method for attenuation of reverberations using a pressure-velocity bottom cable Download PDF

Info

Publication number
WO1997044685A1
WO1997044685A1 PCT/IB1997/000435 IB9700435W WO9744685A1 WO 1997044685 A1 WO1997044685 A1 WO 1997044685A1 IB 9700435 W IB9700435 W IB 9700435W WO 9744685 A1 WO9744685 A1 WO 9744685A1
Authority
WO
WIPO (PCT)
Prior art keywords
result
data
velocity
sensor
pressure
Prior art date
Application number
PCT/IB1997/000435
Other languages
French (fr)
Inventor
Nicolae Moldoveanu
Original Assignee
Schlumberger Holdings Limited
Schlumberger Canada Limited
Services Petroliers Schlumberger
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Schlumberger Holdings Limited, Schlumberger Canada Limited, Services Petroliers Schlumberger filed Critical Schlumberger Holdings Limited
Priority to GB9822606A priority Critical patent/GB2326945B/en
Priority to CA002255631A priority patent/CA2255631C/en
Priority to AU23047/97A priority patent/AU710826B2/en
Priority to EP97915653A priority patent/EP0900393B1/en
Publication of WO1997044685A1 publication Critical patent/WO1997044685A1/en
Priority to NO19985370A priority patent/NO327298B1/en

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/364—Seismic filtering
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • G01V1/3808—Seismic data acquisition, e.g. survey design
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/56—De-ghosting; Reverberation compensation

Definitions

  • This invention relates to seismic profiling of formations located under bodies of water using pressure-velocity bottom cables.
  • Marine seismic exploration is an extremely important tool for location of off ⁇ shore reserves.
  • One known procedure for marine seismic exploration involves use of an ocean bottom cable as illustrated in prior art Fig. 1.
  • Surveys using ocean bottom cables are typically employed in areas populated with numerous obstacles, such as drilling and production platforms.
  • several miles of bottom cables 11 (only one shown in the Fig. 1) are deployed along the sea floor 13 by vessel 15.
  • multiple cables 11 are deployed in parallel, as shown in Fig. 2.
  • the bottom cable 11 is provided with a plurality of sensor pairs 17 placed at regular intervals along the cable, each sensor pair 17 containing a pressure sensor (e.g. , hydrophone) and particle velocity sensor (e.g., a geophone).
  • a pressure sensor e.g. , hydrophone
  • particle velocity sensor e.g., a geophone
  • Acoustic energy is generated in the vicinity of the cable using an air gun array or a marine vibrator array 19.
  • the source wavelet travels downward through the earth and is partially reflected by subsurface layers (formation 21 in Fig. 1) that present an acoustic impedance contrast.
  • the primary reflected wavelet 23 travels upwardly from the subsurface layer, and the pressure waves generated by the upward-traveling reflection are detected by the sensor pairs 17.
  • Seismic exploration using ocean bottom cables is complicated by secondary waves such as wave 25, known as “ghosts,” that are received by the sensors pairs 17 as downward-traveling reflections after reflecting off the air/ water boundary at the surface 29.
  • the air/water boundary is an efficient reflector, and thus the ghosts are significant in amplitude and are difficult to differentiate from the primary waves.
  • These ghosts adversely affect the data obtained during the exploration by attenuating certain frequencies.
  • the spectral ghost notches fall in the seismic frequency band and drastically affect the seismic resolution. Resolution is further complicated by multiple reflection waves and water layer reverberations such as wave 27.
  • Fig. 3a illustrates amplitude versus time hydrophone response at the water bottom.
  • the hydrophone response is defined as T.
  • the hydrophone response is defined as T.
  • the geophone amplitude versus time response is shown in Fig. 3b.
  • the amplitude at t ⁇ is MT, where M is a sensitivity scaling factor that depends upon the particular type of sensors used.
  • M is a sensitivity scaling factor that depends upon the particular type of sensors used.
  • the geophone response is (l-r)MT.
  • the response is -r(l-r)MT, and at time 14 is r2(l-r)MT.
  • the attenuation of the reverberations can be achieved by adding the hydrophone and the geophone signals together after the signals have been suitably scaled.
  • the scale factor S can vary among different sensor pair locations on the same cable.
  • U.S. Patent No.5, 235,554 describes a method where a calibration survey is used to estimate the water bottom reflection coefficient. In such a calibration survey, a low energy source is fired over each sensor pair, and the scale is determined from the ratio of the peaks of the first arrivals of the hydrophone and geophone signals. Collection of this survey data requires additional time and cost over and above the data acquisition phase of the survey.
  • U.S. Patent No.5, 235,554 describes a method where a calibration survey is used to estimate the water bottom reflection coefficient. In such a calibration survey, a low energy source is fired over each sensor pair, and the scale is determined from the ratio of the peaks of the first arrivals of the hydrophone and geophone signals. Collection of this survey data requires additional time and cost over and above the data acquisition phase of the survey.
  • Patent No.5, 396,472 describes a method to derive the water bottom reflection coefficient that eliminates the need for a separate calibration survey, but involves complex mathematics including summing the pressure and velocity signals, multiplying the results by the inverse Backus operator, and then solving for the water bottom reflectivity r using an optimization algorithm.
  • U.S. Patent No. 5,365,492 describes a method wherein the hydrophone signal is used first to adaptively remove noise from the geophone signal, and then the cleaned geophone signals are scaled by a scaling factor and added to the hydrophone signals. The resulting signal is then auto-correlated and the relative amplitude of the first auto-correlation's function-side lobe is measured. The optimum scale factor for the geophone data is then found by optimizing the value of the scale factor with respect to the first side-lobe amplitude of the auto-correlation.
  • the invention relates to a method of generating a seismic profile using an ocean bottom cable having at least one sensor pair comprising a pressure sensor and a velocity sensor, comprising generating a seismic wave for reflection from earth strata; collecting pressure data with the pressure sensor and velocity data with the velocity sensor; multiplying the pressure data with the absolute value of the velocity data to produce a first result; multiplying the velocity data with the absolute value of the pressure data to produce a second result; summing the first result and the second result to create a third result; dividing the third result by a factor of 2 to obtain a fourth result; dividing the fourth result by a sensitivity scaling factor; recording a positive/negative sign for the fourth result; taking the square root of the absolute value of the fourth result to produce a fifth result; replacing the positive/negative sign of the fourth result into the fifth result to produce a sixth result; and incorporating the sixth result into the seismic profile.
  • the invention in another aspect, relates to a method of generating a seismic profile using an ocean bottom cable having a first sensor pair and a second sensor pair, each of the first sensor pair and the second sensor pair comprising a pressure sensor and a velocity sensor, comprising generating a seismic wave for reflection from earth strata; collecting pressure data with the pressure sensors and velocity data with the velocity sensors; multiplying the pressure data from the first sensor pair with the absolute value of the velocity data from the first sensor pair to produce a first result; multiplying the velocity data from the first sensor pair with the absolute value of the pressure data from the first sensor pair to produce a second result; summing the first result and the second result to create a third result; dividing the third result by a factor of 2 to obtain a fourth result; dividing the fourth result by a sensitivity scaling factor; recording a positive/negative sign for the fourth result; taking the square root of the absolute value of the fourth result to produce a fifth result; replacing the positive/negative sign of the fourth result into the fifth result to produce a sixth result; multiplying the pressure
  • the invention in another aspect, relates to a method of attenuating reverberations in seismic data, the seismic data including pressure data and velocity data containing reverberations, comprising reading the pressure data from at least one pressure sensor; reading the velocity data from at least one velocity sensor; and processing the pressure data and the velocity data to produce composite data wherein the reverberations have been attenuated, comprising the steps of: multiplying the pressure data with the absolute value of the velocity data to produce a first result; multiplying the velocity data with the absolute value of the pressure data to produce a second result; summing the first result and the second result to create a third result; dividing the third result by a factor of 2 to obtain a fourth result; dividing the fourth result by a sensitivity scaling factor; recording a positive/negative sign for the fourth result; taking the square root of the absolute value of the fourth result to produce a fifth result; and replacing the positive/negative sign of the fourth result into the fifth result to produce said composite data.
  • Fig. 1 is an illustration of a conventional ocean bottom cable technique
  • Fig. 2 is an illustration of an array of ocean bottom cables parallel to each other;
  • Figs. 3a and 3b are graphs illustrating hydrophone and geophone response at the water bottom
  • Fig. 4 is a flow chart illustrating a method in accordance with an embodiment of the invention.
  • Figs. 5a-5c are graphs illustrating results obtained from the responses shown in Figs. 3a and 3b using a method in accordance with the invention
  • Fig. 6 is a graph of hydrophone-geophone primary reflection data with reverberations
  • Fig. 7 is a graph of the data of Fig. 6 after application of a method in accordance with the invention.
  • the invention allows unwanted reverberations in ocean bottom cable (OBC) data to be attenuated without the need for determining water bottom reflectivity and without requiring complex mathematical operations.
  • Data may be collected by any known method including, e.g., the method shown in prior art Figs. 1 and 2.
  • Fig. 4 is a flow chart of the steps employed to attenuate these reverberations in accordance with an embodiment of the invention.
  • a seismic wave (51) is created and data are recorded with pressure and particle velocity sensors (53). Once the data are collected, the pressure response signal is multiplied by the absolute value of the velocity signal (55) at each time t, producing a result as shown in Fig. 5a.
  • the combined primary wave response signal at time tj now has an amplitude of MT 2 ; at t2, the amplitude of the signal is -M(l-r 2 )T 2 ; at t ⁇ , the amplitude of the signal is Mr 2 (l-r 2 )T 2 ; and finally, at .4, the amplitude is equal to -Mr4(l-r 2 )T 2 (additional reverberations are not shown).
  • Combined primary wave response signal at time t ⁇ has an amplitude of
  • the amplitude of the combined response signal for the first reverberation is equal to
  • Figs. 5a and 5b the multiplication process has made all of the reverberations have identical amplitudes, but opposite polarities as shown by their positive/negative signs. Adding these signals at each time interval (59) produces a graph with only one response signal located at t ⁇ with a magnitude of 2MT 2 . This result is shown in Fig. 5c. This resultant signal is a function of the primary signal alone; all reverberations have been cancelled out by their respective counterparts. This signal is then divided by two (61) and by the sensitivity scaling factor M (63). The sensitivity scaling factor M is dependent only on the type of sensors used in the bottom cable and this is specified by the manufacturer of the ocean bottom cable.
  • Fig. 6 illustrates synthetic OBC data containing six clusters of data (73, 75, 77, 79, and 81). Each cluster includes one or more primary reflections 83 and a plurality of reverberations 85. It is clear from the data shown in Fig. 6 that the reverberations 85 are significant in amplitude and will seriously distort any seismic profile produced using this data if they are not successfully attenuated.
  • Fig. 7 illustrates the data of Fig. 6 after application of the method described above and illustrated in Fig. 4. As is clearly shown, only the primary waves 83 are remaining in the data of Fig. 7. Thus, much higher resolution can be achieved in the seismic profile, without necessitating expensive operations to determine a scaling factor for each receiver position in the survey.
  • the processing of the seismic data is straightforward.
  • the processing sequence by which such data as the data shown in Fig. 7 are converted in a final seismic profile that can be geologically interpretated is well known in the art and thus will not be described in detail herein. For instance, this can include deconvolution, static corrections, velocity determination, normal moveout corrections, dip moveout corrections, stack and migration.
  • Detailed description of seismic data processing can be found in any of numerous texts on the subject, including for instance, "Seismic data processing" by Oz Yilmaz (Society of Exploration Geophysicists, Tulsa, Oklahoma, 1987).

Abstract

A method is described that attenuates reverberations in ocean bottom cable data used to generate a seismic profile. First, a seismic wave is generated and recorded by pressure and particle velocity sensor pairs. Adding the product of the pressure data times the absolute value of the velocity data and the product of the velocity-data times the absolute value of the pressure data eliminates-reverberations, leaving a combined signal that is a function of only the primary wave response signal. Dividing this combined signal by 2 and a sensitivity scaling factor, and then taking the square root of the result produces the magnitude of the primary signal. Replacing the positive/negative sign of the combined signal that existed before taking the absolute value provides the proper polarity. The data are then incorporated into a seismic profile.

Description

METHOD FOR ATTENUATION OF REVERBERATIONS USING A PRESSURE- VELOCITY BOTTOM CABLE
Field of the Invention
This invention relates to seismic profiling of formations located under bodies of water using pressure-velocity bottom cables.
Background of the Invention
Marine seismic exploration is an extremely important tool for location of off¬ shore reserves. One known procedure for marine seismic exploration involves use of an ocean bottom cable as illustrated in prior art Fig. 1. Surveys using ocean bottom cables are typically employed in areas populated with numerous obstacles, such as drilling and production platforms. In this technique, several miles of bottom cables 11 (only one shown in the Fig. 1) are deployed along the sea floor 13 by vessel 15. Usually, multiple cables 11 are deployed in parallel, as shown in Fig. 2. The bottom cable 11 is provided with a plurality of sensor pairs 17 placed at regular intervals along the cable, each sensor pair 17 containing a pressure sensor (e.g. , hydrophone) and particle velocity sensor (e.g., a geophone). Acoustic energy is generated in the vicinity of the cable using an air gun array or a marine vibrator array 19. The source wavelet travels downward through the earth and is partially reflected by subsurface layers (formation 21 in Fig. 1) that present an acoustic impedance contrast. The primary reflected wavelet 23 travels upwardly from the subsurface layer, and the pressure waves generated by the upward-traveling reflection are detected by the sensor pairs 17.
Seismic exploration using ocean bottom cables is complicated by secondary waves such as wave 25, known as "ghosts," that are received by the sensors pairs 17 as downward-traveling reflections after reflecting off the air/ water boundary at the surface 29. The air/water boundary is an efficient reflector, and thus the ghosts are significant in amplitude and are difficult to differentiate from the primary waves. These ghosts adversely affect the data obtained during the exploration by attenuating certain frequencies. In addition, when the water depth is large, the spectral ghost notches fall in the seismic frequency band and drastically affect the seismic resolution. Resolution is further complicated by multiple reflection waves and water layer reverberations such as wave 27.
The purpose of using both hydrophones and geophones in the ocean bottom cable is to capitalize on the differences between these two types of sensors to attenuate the downgoing waves which include the ghosts and the water layer reverberations. Their responses to the primary reflections are in phase, but are 180" out of phase to the ghosts and to the reverberations.
Fig. 3a illustrates amplitude versus time hydrophone response at the water bottom. For the hydrophone, at time l\ (the primary wave), the hydrophone response is defined as T. At time t2 (the first water layer reverberation), the hydrophone response is
-(l+r)T; at time 13 (the second reverberation) it is r(l-fr)T; and at time t4(the third reverberation) , it is -r2(l+r)T, where r is between 0 and 1 and is the water bottom reflectivity at each receiver position. Additional reverberations continue to decrease in amplitude.
The geophone amplitude versus time response is shown in Fig. 3b. The amplitude at t\ is MT, where M is a sensitivity scaling factor that depends upon the particular type of sensors used. At time t2, the geophone response is (l-r)MT. At time 131 the response is -r(l-r)MT, and at time 14 is r2(l-r)MT.
It is apparent from Figs. 3a and 3b that while the primary geophone and hydrophone responses are in phase, the responses to the water layer reverberations are 180° out of phase. Thus, the attenuation of the reverberations can be achieved by adding the hydrophone and the geophone signals together after the signals have been suitably scaled. Theoretically, the scale factor S=(l+r)/(l-r), where r is the water bottom reflectivity, should be applied to the geophone data. Determination of the water bottom reflectivity coefficient r depends upon the acoustic impedance of the bottom material. Thus, the scale factor S can vary among different sensor pair locations on the same cable.
There are several known methods for deriving the scaling factors for geophone signals. U.S. Patent No.5, 235,554 describes a method where a calibration survey is used to estimate the water bottom reflection coefficient. In such a calibration survey, a low energy source is fired over each sensor pair, and the scale is determined from the ratio of the peaks of the first arrivals of the hydrophone and geophone signals. Collection of this survey data requires additional time and cost over and above the data acquisition phase of the survey. U.S. Patent No.5, 396,472 describes a method to derive the water bottom reflection coefficient that eliminates the need for a separate calibration survey, but involves complex mathematics including summing the pressure and velocity signals, multiplying the results by the inverse Backus operator, and then solving for the water bottom reflectivity r using an optimization algorithm. U.S. Patent No. 5,365,492 describes a method wherein the hydrophone signal is used first to adaptively remove noise from the geophone signal, and then the cleaned geophone signals are scaled by a scaling factor and added to the hydrophone signals. The resulting signal is then auto-correlated and the relative amplitude of the first auto-correlation's function-side lobe is measured. The optimum scale factor for the geophone data is then found by optimizing the value of the scale factor with respect to the first side-lobe amplitude of the auto-correlation.
Summary of the Invention
In one aspect, the invention relates to a method of generating a seismic profile using an ocean bottom cable having at least one sensor pair comprising a pressure sensor and a velocity sensor, comprising generating a seismic wave for reflection from earth strata; collecting pressure data with the pressure sensor and velocity data with the velocity sensor; multiplying the pressure data with the absolute value of the velocity data to produce a first result; multiplying the velocity data with the absolute value of the pressure data to produce a second result; summing the first result and the second result to create a third result; dividing the third result by a factor of 2 to obtain a fourth result; dividing the fourth result by a sensitivity scaling factor; recording a positive/negative sign for the fourth result; taking the square root of the absolute value of the fourth result to produce a fifth result; replacing the positive/negative sign of the fourth result into the fifth result to produce a sixth result; and incorporating the sixth result into the seismic profile.
In another aspect, the invention relates to a method of generating a seismic profile using an ocean bottom cable having a first sensor pair and a second sensor pair, each of the first sensor pair and the second sensor pair comprising a pressure sensor and a velocity sensor, comprising generating a seismic wave for reflection from earth strata; collecting pressure data with the pressure sensors and velocity data with the velocity sensors; multiplying the pressure data from the first sensor pair with the absolute value of the velocity data from the first sensor pair to produce a first result; multiplying the velocity data from the first sensor pair with the absolute value of the pressure data from the first sensor pair to produce a second result; summing the first result and the second result to create a third result; dividing the third result by a factor of 2 to obtain a fourth result; dividing the fourth result by a sensitivity scaling factor; recording a positive/negative sign for the fourth result; taking the square root of the absolute value of the fourth result to produce a fifth result; replacing the positive/negative sign of the fourth result into the fifth result to produce a sixth result; multiplying the pressure data from the second sensor pair with the absolute value of the velocity data from the second sensor pair to produce a seventh result; multiplying the velocity data from the second sensor pair with the absolute value of the pressure data from the second sensor pair to produce an eighth result; summing the seventh result and the eighth result to create a ninth result; dividing the ninth result by a factor of 2 to obtain a tenth result; dividing the tenth result by a sensitivity scaling factor; recording a positive/negative sign for the tenth result; taking the square root of the absolute value of the tenth result to produce a eleventh result; replacing the positive/negative sign of the tenth result into the eleventh result to produce a twelfth result; and incorporating the sixth result and the twelfth result into the seismic profile.
In another aspect, the invention relates to a method of attenuating reverberations in seismic data, the seismic data including pressure data and velocity data containing reverberations, comprising reading the pressure data from at least one pressure sensor; reading the velocity data from at least one velocity sensor; and processing the pressure data and the velocity data to produce composite data wherein the reverberations have been attenuated, comprising the steps of: multiplying the pressure data with the absolute value of the velocity data to produce a first result; multiplying the velocity data with the absolute value of the pressure data to produce a second result; summing the first result and the second result to create a third result; dividing the third result by a factor of 2 to obtain a fourth result; dividing the fourth result by a sensitivity scaling factor; recording a positive/negative sign for the fourth result; taking the square root of the absolute value of the fourth result to produce a fifth result; and replacing the positive/negative sign of the fourth result into the fifth result to produce said composite data.
Brief Description of the Drawings
Fig. 1 is an illustration of a conventional ocean bottom cable technique;
Fig. 2 is an illustration of an array of ocean bottom cables parallel to each other;
Figs. 3a and 3b are graphs illustrating hydrophone and geophone response at the water bottom;
Fig. 4 is a flow chart illustrating a method in accordance with an embodiment of the invention;
Figs. 5a-5c are graphs illustrating results obtained from the responses shown in Figs. 3a and 3b using a method in accordance with the invention;
Fig. 6 is a graph of hydrophone-geophone primary reflection data with reverberations; and Fig. 7 is a graph of the data of Fig. 6 after application of a method in accordance with the invention.
Description of the Preferred Embodiments
Preferred embodiments of the invention will now be described with reference to the accompanying drawings.
The invention allows unwanted reverberations in ocean bottom cable (OBC) data to be attenuated without the need for determining water bottom reflectivity and without requiring complex mathematical operations. Data may be collected by any known method including, e.g., the method shown in prior art Figs. 1 and 2.
Fig. 4 is a flow chart of the steps employed to attenuate these reverberations in accordance with an embodiment of the invention. First, a seismic wave (51) is created and data are recorded with pressure and particle velocity sensors (53). Once the data are collected, the pressure response signal is multiplied by the absolute value of the velocity signal (55) at each time t, producing a result as shown in Fig. 5a. The combined primary wave response signal at time tj now has an amplitude of MT2; at t2, the amplitude of the signal is -M(l-r2)T2; at t^, the amplitude of the signal is Mr2(l-r2)T2; and finally, at .4, the amplitude is equal to -Mr4(l-r2)T2 (additional reverberations are not shown).
Similarly, the response signal from the velocity sensor is multiplied by the absolute value of the pressure signal (57) at each value of time t. Fig. 5b shows the result of this operation. Combined primary wave response signal at time t\ has an amplitude of
MT2, identical to the amplitude of the primary wave response signal in Fig. 5a. At time t2, the amplitude of the combined response signal for the first reverberation is equal to
M(l-r2)T2; at 13, the amplitude of the combined response signal is now -Mr2(l-r2)T2; at t4, the amplitude is Mr4(l-r2)T2.
As can be seen from Figs. 5a and 5b, the multiplication process has made all of the reverberations have identical amplitudes, but opposite polarities as shown by their positive/negative signs. Adding these signals at each time interval (59) produces a graph with only one response signal located at t\ with a magnitude of 2MT2. This result is shown in Fig. 5c. This resultant signal is a function of the primary signal alone; all reverberations have been cancelled out by their respective counterparts. This signal is then divided by two (61) and by the sensitivity scaling factor M (63). The sensitivity scaling factor M is dependent only on the type of sensors used in the bottom cable and this is specified by the manufacturer of the ocean bottom cable. Taking the square root of the absolute value of this last result (65) produces the true amplitude of the primary signal. Replacing the sign that existed before taking the absolute value in the previous step (67) provides the proper polarity for the signal. This result is then incorporated into a seismic profile (69).
Fig. 6 illustrates synthetic OBC data containing six clusters of data (73, 75, 77, 79, and 81). Each cluster includes one or more primary reflections 83 and a plurality of reverberations 85. It is clear from the data shown in Fig. 6 that the reverberations 85 are significant in amplitude and will seriously distort any seismic profile produced using this data if they are not successfully attenuated.
Fig. 7 illustrates the data of Fig. 6 after application of the method described above and illustrated in Fig. 4. As is clearly shown, only the primary waves 83 are remaining in the data of Fig. 7. Thus, much higher resolution can be achieved in the seismic profile, without necessitating expensive operations to determine a scaling factor for each receiver position in the survey. Once this data is produced with the reverberations attenuated, the processing of the seismic data is straightforward. The processing sequence by which such data as the data shown in Fig. 7 are converted in a final seismic profile that can be geologically interpretated is well known in the art and thus will not be described in detail herein. For instance, this can include deconvolution, static corrections, velocity determination, normal moveout corrections, dip moveout corrections, stack and migration. Detailed description of seismic data processing can be found in any of numerous texts on the subject, including for instance, "Seismic data processing" by Oz Yilmaz (Society of Exploration Geophysicists, Tulsa, Oklahoma, 1987).
Various embodiments of the invention have been shown and described. However, the invention is not so limited, but rather is limited only the scope of the amended claims.

Claims

1. A method of generating a seismic profile using an ocean bottom cable having at least one sensor pair comprising a pressure sensor and a velocity sensor, comprising:
generating a seismic wave for reflection from earth strata;
collecting pressure data with the pressure sensor and velocity data with the velocity sensor;
multiplying the pressure data with the absolute value of the velocity data to produce a first result;
multiplying the velocity data with the absolute value of the pressure data to produce a second result;
summing the first result and the second result to create a third result;
dividing the third result by a factor of 2 to obtain a fourth result;
dividing the fourth result by a sensitivity scaling factor;
recording a positive/negative sign for the fourth result;
taking the square root of the absolute value of the fourth result to produce a fifth result;
replacing the positive/negative sign of the fourth result into the fifth result to produce a sixth result; and
incorporating the sixth result into the seismic profile.
2. The method of claim 1 , wherein the velocity sensor is a geophone.
3. The method of claim 1 or claim 2, wherein the pressure sensor is a hydrophone.
4. A method of generating a seismic profile using an ocean bottom cable having a first sensor pair and a second sensor pair, each of the first sensor pair and the second sensor pair comprising a pressure sensor and a velocity sensor, comprising:
generating a seismic wave for reflection from earth strata; collecting pressure data with the pressure sensors and velocity data with the velocity sensors;
multiplying the pressure data from the first sensor pair with the absolute value of the velocity data from the first sensor pair to produce a first result;
multiplying the velocity data from the first sensor pair with the absolute value of the pressure data from the first sensor pair to produce a second result;
summing the first result and the second result to create a third result;
dividing the third result by a factor of 2 to obtain a fourth result;
dividing the fourth result by a sensitivity scaling factor;
recording a positive/ negative sign for the fourth result;
taking the square root of the absolute value of the fourth result to produce a fifth result;
replacing the positive/negative sign of the fourth result into the fifth result to produce a sixth result;
multiplying the pressure data from the second sensor pair with the absolute value of the velocity data from the second sensor pair to produce a seventh result;
multiplying the velocity data from the second sensor pair with the absolute value of the pressure data from the second sensor pair to produce an eighth result;
summing the seventh result and the eighth result to create a ninth result;
dividing the ninth result by a factor of 2 to obtain a tenth result;
dividing the tenth result by a sensitivity scaling factor;
recording a positive/negative sign for the tenth result;
taking the square root of the absolute value of the tenth result to produce a eleventh result;
replacing the positive/negative sign of the tenth result into the eleventh result to produce a twelfth result; and
incorporating the sixth result and the twelfth result into the seismic profile.
5. The method of claim 4, further comprising using a second ocean bottom cable having a first sensor pair and a second sensor pair, each of the first sensor pair and the second sensor pair comprising a pressure sensor and a velocity sensor.
6. The method of claim 4 or claim 5, wherein each velocity sensor is a geophone.
7. The method of any one of claims 5 to 6, wherein the pressure sensor is a hydrophone.
8. A method of attenuating reverberations in seismic data, the seismic data including pressure data and velocity data containing reverberations, comprising:
reading the pressure data from at least one pressure sensor;
reading the velocity data from at least one velocity sensor;
processing the pressure data and the velocity data to produce composite data wherein the reverberations have been reduced, comprising the steps of:
multiplying the pressure data with the absolute value of the velocity data to produce a first result;
multiplying the velocity data with the absolute value of the pressure data to produce a second result;
summing the first result and the second result to create a third result;
dividing the third result by a factor of 2 to obtain a fourth result;
dividing the fourth result by a sensitivity scaling factor;
recording a positive/negative sign for the fourth result;
taking the square root of the absolute value of the fourth result to produce a fifth result; and
replacing the positive/negative sign of the fourth result into the fifth result to produce a sixth result.
9. The method of claim 8, wherein the or each velocity sensor is a geophone
10. The method of claim 8 or claim 9, wherein the or each pressure sensor is a hydrophone.
PCT/IB1997/000435 1996-05-20 1997-04-24 Method for attenuation of reverberations using a pressure-velocity bottom cable WO1997044685A1 (en)

Priority Applications (5)

Application Number Priority Date Filing Date Title
GB9822606A GB2326945B (en) 1996-05-20 1997-04-24 Method for attenuation of reverberations using a pressure-velocity bottom cable
CA002255631A CA2255631C (en) 1996-05-20 1997-04-24 Method for attenuation of reverberations using a pressure-velocity bottom cable
AU23047/97A AU710826B2 (en) 1996-05-20 1997-04-24 Method for attenuation of reverberations using a pressure-velocity bottom cable
EP97915653A EP0900393B1 (en) 1996-05-20 1997-04-24 Method for attenuation of reverberations using a pressure-velocity bottom cable
NO19985370A NO327298B1 (en) 1996-05-20 1998-11-18 Method of attenuating reverberation by means of a pressure-velocity bottom cable

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US08/650,325 US5621700A (en) 1996-05-20 1996-05-20 Method for attenuation of reverberations using a pressure-velocity bottom cable
US08/650,325 1996-05-20

Publications (1)

Publication Number Publication Date
WO1997044685A1 true WO1997044685A1 (en) 1997-11-27

Family

ID=24608423

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/IB1997/000435 WO1997044685A1 (en) 1996-05-20 1997-04-24 Method for attenuation of reverberations using a pressure-velocity bottom cable

Country Status (7)

Country Link
US (1) US5621700A (en)
EP (1) EP0900393B1 (en)
AU (1) AU710826B2 (en)
CA (1) CA2255631C (en)
GB (1) GB2326945B (en)
NO (1) NO327298B1 (en)
WO (1) WO1997044685A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6775618B1 (en) 1999-03-22 2004-08-10 Schlumberger Technology Corporation Method and system for reducing effects of sea surface ghost contamination in seismic data

Families Citing this family (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2743896B1 (en) * 1996-01-23 1998-04-10 Geophysique Cie Gle METHOD FOR PROCESSING FOR CALIBRATION OF A PAIR OF HYDROPHONE/GEOPHONE SENSORS AND METHOD FOR SEISMIC PROSPECTING IMPLEMENTING THIS PROCESSING
FR2743897B1 (en) * 1996-01-23 1998-04-10 Geophysique Cie Gle PROCESS FOR MARINE SEISMIC PROSPECTION USING A COUPLE OF HYDROPHONE AND GEOPHONE SENSORS
US5774417A (en) * 1996-10-25 1998-06-30 Atlantic Richfield Company Amplitude and phase compensation in dual-sensor ocean bottom cable seismic data processing
US6021092A (en) * 1997-06-09 2000-02-01 Western Atlas International, Inc. Method for deriving surface consistent reflectivity map from dual sensor seismic data
GB9800741D0 (en) * 1998-01-15 1998-03-11 Geco As Multiple attenuation of multi-component sea-bottom data
US5991238A (en) * 1998-06-09 1999-11-23 Western Atlas International, Inc. Weighted backus filter method of combining dual sensor traces
US6246637B1 (en) * 1998-09-28 2001-06-12 Pgs Tensor, Inc. Method and system for combining three component seismic data
CA2431251A1 (en) * 2000-12-07 2002-06-13 Schlumberger Canada Limited A method of and apparatus for processing seismic data
AU2002310037B2 (en) 2001-05-25 2006-10-05 Exxonmobil Upstream Research Company Multiple suppression for ocean bottom seismic data
GB2384053B (en) * 2002-01-11 2005-09-14 Westerngeco Ltd A method of and apparatus for processing seismic data
US7239577B2 (en) * 2002-08-30 2007-07-03 Pgs Americas, Inc. Apparatus and methods for multicomponent marine geophysical data gathering
US6704244B1 (en) * 2002-10-08 2004-03-09 Pgs Americas, Inc. Method for acquiring and processing seismic survey data using ocean bottom cables and streamers
US7123543B2 (en) * 2003-07-16 2006-10-17 Pgs Americas, Inc. Method for seismic exploration utilizing motion sensor and pressure sensor data
GB2405473B (en) * 2003-08-23 2005-10-05 Westerngeco Ltd Multiple attenuation method
US6943729B2 (en) * 2003-10-01 2005-09-13 S5 Wireless, Inc. Method and system for time difference of arrival (TDOA) location services
US7359283B2 (en) * 2004-03-03 2008-04-15 Pgs Americas, Inc. System for combining signals of pressure sensors and particle motion sensors in marine seismic streamers
US7139217B2 (en) * 2004-05-27 2006-11-21 Pgs Americas, Inc. Water bottom cable seismic survey cable and system
US20060133202A1 (en) * 2004-12-22 2006-06-22 Tenghamn Stig R L Motion sensors in a marine seismic streamer
US7477992B2 (en) * 2005-02-18 2009-01-13 Exxonmobil Upstream Research Company Method for combining seismic data sets
US7768869B2 (en) * 2005-05-05 2010-08-03 Pgs Americas, Inc. Method for deghosting and water layer multiple reflection attenuation in marine seismic data
US7167413B1 (en) 2006-05-01 2007-01-23 Input/Output Towed streamer deghosting
US7599798B2 (en) * 2006-09-11 2009-10-06 Westerngeco L.L.C. Migrating composite seismic response data to produce a representation of a seismic volume
US7676327B2 (en) 2007-04-26 2010-03-09 Westerngeco L.L.C. Method for optimal wave field separation
US8553490B2 (en) * 2007-11-09 2013-10-08 Pgs Geophysical As Array grouping of seismic sensors in a marine streamer for optimum noise attenuation
US7675812B2 (en) * 2008-06-30 2010-03-09 Pgs Geophysical As Method for attenuation of multiple reflections in seismic data
US9110191B2 (en) * 2009-03-30 2015-08-18 Westerngeco L.L.C. Multiple attenuation for ocean-bottom seismic data
EP2526444B1 (en) * 2010-01-22 2020-09-09 ION Geophysical Corporation Seismic system with ghost and motion rejection
US20110310698A1 (en) 2010-06-21 2011-12-22 Sercel, Inc. Dual Axis Geophones For Pressure/Velocity Sensing Streamers Forming a Triple Component Streamer
US8982662B2 (en) 2010-09-02 2015-03-17 Ion Geophysical Corporation Multi-component, acoustic-wave sensor and methods
US9335430B2 (en) 2012-07-18 2016-05-10 Pgs Geophysical As Wave field separation by mixed domain inversion
EP2962131A2 (en) 2013-03-14 2016-01-06 ION Geophysical Corporation Seismic sensor devices, systems, and methods including noise filtering

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5235554A (en) * 1991-03-11 1993-08-10 Halliburton Geophysical Services, Inc. Method for correcting impulse response differences of hydrophones and geophones as well as geophone coupling to the water-bottom in dual-sensor, bottom-cable seismic operations
US5365492A (en) * 1993-08-04 1994-11-15 Western Atlas International, Inc. Method for reverberation suppression
US5396472A (en) * 1993-09-24 1995-03-07 Western Atlas International Method for deriving water bottom reflectivity in dual sensor seismic surveys

Family Cites Families (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2757356A (en) * 1954-01-08 1956-07-31 Texas Instruments Inc Method and apparatus for canceling reverberations in water layers
US2992694A (en) * 1955-01-17 1961-07-18 Socony Mobil Oil Co Inc Time delayed directional seismograph
US3943484A (en) * 1973-11-26 1976-03-09 Avance Oil & Gas Company, Inc. Method of attenuating unwanted seismic reflections in underwater seismic exploration
US4486865A (en) * 1980-09-02 1984-12-04 Mobil Oil Corporation Pressure and velocity detectors for seismic exploration
US4405036A (en) * 1980-12-11 1983-09-20 Marathon Oil Company Seafloor velocity and amplitude measurement apparatus method
US4516227A (en) * 1981-12-04 1985-05-07 Marathon Oil Company Subocean bottom explosive seismic system
US4520467A (en) * 1982-03-18 1985-05-28 Shell Oil Company Marine seismic system
US4558437A (en) * 1982-04-05 1985-12-10 Marathon Oil Company Seafloor velocity and amplitude measurement apparatus and method therefor
US4648080A (en) * 1984-06-15 1987-03-03 Western Geophysical Company Method for determining the far field signature of a marine seismic source from near-field measurements
US4992992A (en) * 1988-10-21 1991-02-12 Western Atlas International, Inc. Processing for seismic data from slanted cable
US5253217A (en) * 1989-04-14 1993-10-12 Atlantic Richfield Company Method for seismic exploration including compensation for near surface effects
US4937793A (en) * 1989-05-30 1990-06-26 Halliburton Geophysical Services, Inc. Processing method for marine seismic surveying utilizing dual streamers
US4935903A (en) * 1989-05-30 1990-06-19 Halliburton Geophysical Services, Inc. Reinforcement of surface seismic wavefields
US4958328A (en) * 1989-07-24 1990-09-18 Texaco Inc. Marine walkaway vertical seismic profiling
EG19158A (en) * 1989-08-25 1996-02-29 Halliburton Geophys Service System for attenuation of water-column reverberation
US5051961A (en) * 1989-10-26 1991-09-24 Atlantic Richfield Company Method and apparatus for seismic survey including using vertical gradient estimation to separate downgoing seismic wavefield
NO168557C (en) * 1989-10-26 1992-03-04 Norske Stats Oljeselskap SEISMIC DEVICE.
US5163028A (en) * 1991-09-27 1992-11-10 Halliburton Geophysical Services, Inc. Method for correcting impulse response differences of hydrophones and geophones as well as geophone coupling to the water-bottom in dual-sensor, bottom-cable seismic operations

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5235554A (en) * 1991-03-11 1993-08-10 Halliburton Geophysical Services, Inc. Method for correcting impulse response differences of hydrophones and geophones as well as geophone coupling to the water-bottom in dual-sensor, bottom-cable seismic operations
US5365492A (en) * 1993-08-04 1994-11-15 Western Atlas International, Inc. Method for reverberation suppression
US5396472A (en) * 1993-09-24 1995-03-07 Western Atlas International Method for deriving water bottom reflectivity in dual sensor seismic surveys

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6775618B1 (en) 1999-03-22 2004-08-10 Schlumberger Technology Corporation Method and system for reducing effects of sea surface ghost contamination in seismic data
USRE41656E1 (en) 1999-03-22 2010-09-07 Westerngeco L. L. C. Method and system for reducing effects of sea surface ghost contamination in seismic data
USRE43188E1 (en) 1999-03-22 2012-02-14 Schlumberger Technology Corporation Method and system for reducing effects of sea surface ghost contamination in seismic data

Also Published As

Publication number Publication date
AU2304797A (en) 1997-12-09
EP0900393A1 (en) 1999-03-10
US5621700A (en) 1997-04-15
EP0900393B1 (en) 2002-01-30
NO985370L (en) 1999-01-20
AU710826B2 (en) 1999-09-30
CA2255631A1 (en) 1997-11-27
GB2326945B (en) 2000-11-22
GB2326945A (en) 1999-01-06
NO327298B1 (en) 2009-06-02
GB9822606D0 (en) 1998-12-09
CA2255631C (en) 2000-09-19
NO985370D0 (en) 1998-11-18

Similar Documents

Publication Publication Date Title
US5621700A (en) Method for attenuation of reverberations using a pressure-velocity bottom cable
US6894948B2 (en) Method for processing dual sensor seismic data to attenuate noise
US5173880A (en) Method of generating seismic wavelets using seismic range equation
US5524100A (en) Method for deriving water bottom reflectivity in dual sensor seismic surveys
US4752916A (en) Method and system for removing the effect of the source wavelet from seismic data
RU2116657C1 (en) Method for processing of seismic data in order to suppress multiple reflections
US5995905A (en) Source signature determination and multiple reflection reduction
US6654693B2 (en) Angle dependent surface multiple attenuation for two-component marine bottom sensor data
AU2010201504B2 (en) Method for calculation of seismic attributes from seismic signals
EP0414344A2 (en) Marine seismic reflection geophysical surveying
US5991238A (en) Weighted backus filter method of combining dual sensor traces
AU723276B2 (en) Method for measuring the water bottom reflectivity
US5963507A (en) Method for improving the accuracy of ocean bottom reflectivity estimations using the inverse backus filter
US5442591A (en) Method for adaptively suppressing noise transients in summed co-sensor seismic recordings
WO1997037246A1 (en) Bottom-referenced seismic cable with vertical hydrophone arrays
WO2001020363A2 (en) Amplitude spectra estimation
GB2349697A (en) Determining water bottom reflectivity in marine seismic survey
Robinson REMOTE DETECTION USING DUAL SENSORS
Robinson Applied Seismology
Fatti Geophysics in oil exploration

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AL AM AT AU AZ BB BG BR BY CA CH CN CZ DE DK EE ES FI GB GE HU IL IS JP KE KG KP KR KZ LK LR LS LT LU LV MD MG MK MN MW MX NO NZ PL PT RO RU SD SE SG SI SK TJ TM TR TT UA UG UZ VN AM AZ BY KG KZ MD RU TJ TM

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH KE LS MW SD SZ UG AT BE CH DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA

DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
121 Ep: the epo has been informed by wipo that ep was designated in this application
ENP Entry into the national phase

Ref country code: GB

Ref document number: 9822606

Kind code of ref document: A

Format of ref document f/p: F

WWE Wipo information: entry into national phase

Ref document number: 1997915653

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 2255631

Country of ref document: CA

Ref country code: CA

Ref document number: 2255631

Kind code of ref document: A

Format of ref document f/p: F

NENP Non-entry into the national phase

Ref country code: JP

Ref document number: 97541891

Format of ref document f/p: F

WWP Wipo information: published in national office

Ref document number: 1997915653

Country of ref document: EP

REG Reference to national code

Ref country code: DE

Ref legal event code: 8642

WWG Wipo information: grant in national office

Ref document number: 1997915653

Country of ref document: EP