Method and apparatus for removing tube wave interference from optical fiber acoustic wave sensing seismic data

Patent metadata

Field	Value
Publication number	US-12326533-B2
Application number	US-202117797708-A
Country	US
Kind code	B2
Filing date	Jan 28, 2021
Priority date	Feb 5, 2020
Publication date	Jun 10, 2025
Grant date	Jun 10, 2025

Abstract

A method for removing tube wave interference from optical fiber acoustic wave sensing seismic data, including: acquiring seismic wavefield data which contains a tube wave and is collected by an optical fiber acoustic wave sensing instrument; calculating a root-mean-square amplitude of the waveform data cut on the seismic trace as an amplitude normalization factor; performing normalization processing on the amplitude value; performing de-tail mean filtering processing on the normalized amplitude value along the travel time of the tube wave, to obtain a predicted amplitude value; performing tube wave interference removal processing on each seismic trace, and performing inverse normalization processing to obtain the seismic wavefield data without tube wave interference. The method effectively suppresses the tube wave interference in the optical fiber acoustic wave sensing seismic data. An apparatus for removing tube wave interference from optical fiber acoustic wave sensing seismic data, and a computer device are further provided.

First claim

The invention claimed is:

1. A method for removing tube wave interference from optical fiber acoustic wave sensing seismic data, the method comprising:
suspending optical fibers in a well, the optical fibers being connected to an optical fiber acoustic wave sensing instrument;
exciting a seismic wavefield by an explosion source or artificial vibroseis on the ground or in the well;
collecting, by the optical fiber acoustic wave sensing instrument, seismic wavefield data which contains a tube wave, wherein the seismic wavefield data comprises amplitude values at each sampling time point of seismic traces excited by a single shot, and the tube wave has an amplitude value that is inversely proportional to a degree of coupling between the optical fibers and a borehole wall of the well;
acquiring, by a processor, the seismic wavefield data;
downwardly cutting, by the processor, waveform data with a preset time window length on each seismic trace along travel time of the tube wave, and calculating, by the processor, a root-mean-square amplitude of the waveform data cut on the seismic trace as an amplitude normalization factor of the seismic trace;
performing, by the processor, normalization processing on the amplitude value of each seismic trace at each sampling time point based on the amplitude normalization factor of the seismic trace, to obtain a normalized amplitude value of the seismic trace at the sampling time point;
performing, by the processor, de-tail mean filtering processing on the normalized amplitude value of each seismic trace at each sampling time point along the travel time of the tube wave, to obtain a predicted amplitude value of the tube wave at the sampling time point;
performing, by the processor, tube wave interference removal processing on each seismic trace at each sampling time point based on the predicted amplitude value of the tube wave at the sampling time point, and performing, by the processor, inverse normalization processing on the seismic trace at the sampling time point using the amplitude normalization factor of the seismic trace, so as to obtain the seismic wavefield data without tube wave interference; and
performing seismic exploration using the seismic wavefield data without tube wave interference.

2. The method according to claim 1 , wherein performing the de-tail mean filtering processing on the normalized amplitude value of each seismic trace at each sampling time point along the travel time of the tube wave, to obtain a predicted amplitude value of the tube wave at the sampling time point, comprises:
sorting the normalized amplitude values of the seismic traces at each sampling time point along a direction of the travel time of the tube wave, removing the first n normalized amplitude values and the last n normalized amplitude values, and performing a mean value calculation on the remaining normalized amplitude values to obtain the predicted amplitude value of the tube wave at the sampling time point, where n is a positive integer.

3. The method according to claim 2 , wherein a value of n is less than half of a total number of the seismic traces involved in filtering.

4. The method according to claim 1 , wherein:
performing the normalization processing on the amplitude value of each seismic trace at each sampling time point based on the amplitude normalization factor of the seismic trace, to obtain the normalized amplitude value of the seismic trace at the sampling time point, comprises: dividing the amplitude value of each seismic trace at each sampling time point by the amplitude normalization factor of the seismic trace, to obtain the normalized amplitude value at the sampling time point of the seismic trace; and
performing the tube wave interference removal processing on each seismic trace at each sampling time point based on the predicted amplitude value of the tube wave at the sampling time point, and performing the inverse normalization processing on the seismic trace at the sampling time point using the amplitude normalization factor of the seismic trace to obtain the seismic wavefield data without the tube wave interference, comprises: subtracting the predicted amplitude value of the tube wave at each sampling time point from the normalized amplitude value of each seismic trace at the sampling time point to obtain a result, and multiplying the result by the amplitude normalization factor of the seismic trace to obtain the seismic wavefield data without the tube wave interference.

5. The method according to claim 1 , wherein the preset time window length is 2 to 10 times of a dominant period of the tube wave.

6. An apparatus comprising:
optical fibers for being suspended in a well;
an optical fiber acoustic wave sensing instrument connected to the optical fibers and configured to collect seismic wavefield data which contains a tube wave, wherein the seismic wavefield data comprises amplitude values at each sampling time point of seismic traces excited by a single shot, the tube wave has an amplitude value that is inversely proportional to a degree of coupling between the optical fibers and a borehole wall of the well, and the seismic wavefield is excited by an explosion source or artificial vibroseis on the ground or in the well;
a seismic wavefield data acquisition unit configured to acquire the seismic wavefield data;
a normalization factor determination unit configured to downwardly cut waveform data with a preset time window length on each seismic trace along travel time of the tube wave, and calculate a root-mean-square amplitude of the waveform data cut on the seismic trace as an amplitude normalization factor of the seismic trace;
a normalization processing unit configured to perform normalization processing on the amplitude value of each seismic trace at each sampling time point based on the amplitude normalization factor of the seismic trace, to obtain a normalized amplitude value of the seismic trace at the sampling time point;
a tube wave amplitude prediction unit configured to perform de-tail mean filtering processing on the normalized amplitude value of each seismic trace at each sampling time point along the travel time of the tube wave, to obtain a predicted amplitude value of the tube wave at the sampling time point; and
a tube wave interference removal unit configured to perform tube wave interference removal processing on each seismic trace at each sampling time point based on the predicted amplitude value of the tube wave at the sampling time point, and perform inverse normalization processing on the seismic trace at the sampling time point using the amplitude normalization factor of the seismic trace to obtain the seismic wavefield data without tube wave interference for seismic exploration.

7. The apparatus according to claim 6 , wherein the tube wave amplitude prediction unit is further configured to sort the normalized amplitude values of the seismic traces at each sampling time point along a direction of the travel time of the tube wave, remove the first n normalized amplitude values and the last n normalized amplitude values, and perform a mean value calculation on the remaining normalized amplitude values to obtain a predicted amplitude value of the tube wave at the sampling time point, where n is a positive integer.

8. The apparatus according to claim 6 , wherein the normalization factor determination unit is further configured to divide the amplitude value of each seismic trace at each sampling time point by the amplitude normalization factor of the seismic trace, to obtain the normalized amplitude value of the seismic trace at the sampling time point; and the tube wave interference removal unit is further configured to subtract the predicted amplitude value of the tube wave at e

Assignees

• China Nat Petroleum Corp
• Bgp Inc China Nat Pet Corp

Inventors

• Li Yanpeng
• Li Fei

Classifications

• G01V1/18

  Receiving elements, e.g. seismometer, geophone {or torque detectors, for localised single point measurements} · CPC title

• G01V2210/324

  Filtering · CPC title

• G01V1/50

  Analysing data · CPC title

• G01V2210/63

  Seismic attributes, e.g. amplitude, polarity, instant phase · CPC title

• G01V1/307

  for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity · CPC title