WO1997026559A1 - Method for removing borehole rugosity noise from well log data - Google Patents

Method for removing borehole rugosity noise from well log data Download PDF

Info

Publication number
WO1997026559A1
WO1997026559A1 PCT/US1996/000449 US9600449W WO9726559A1 WO 1997026559 A1 WO1997026559 A1 WO 1997026559A1 US 9600449 W US9600449 W US 9600449W WO 9726559 A1 WO9726559 A1 WO 9726559A1
Authority
WO
WIPO (PCT)
Prior art keywords
log data
noise
caliper
wavenumber
amplitude spectrum
Prior art date
Application number
PCT/US1996/000449
Other languages
French (fr)
Inventor
Robert G. Keys
John A. Nieto
Keh Pann
Denis P. Schmitt
Original Assignee
Mobil Oil Corporation
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 Mobil Oil Corporation filed Critical Mobil Oil Corporation
Priority to PCT/US1996/000449 priority Critical patent/WO1997026559A1/en
Priority to EP96902668A priority patent/EP0991958A1/en
Publication of WO1997026559A1 publication Critical patent/WO1997026559A1/en
Priority to NO983262A priority patent/NO983262D0/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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48—Processing data
    • 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
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/30—Noise handling
    • G01V2210/32—Noise reduction

Definitions

  • the present invention pertains to noise removal from well log data and more particularly the removal of induced rugosity noise from well log data.
  • the method of the present invention eliminates sinusoidal noise without affecting the response of the formation. This means that the actual formation responses of the logging tools are recovered, and the logs can be used quantitatively.
  • the present invention discloses a robust, fast and easy to use processing module which removes "sinusoidal" noise while preserving quantitative data integrity.
  • the wavenumber is the reciprocal of wavelength.
  • Preserving the quantitative data integrity is accomplished by approximating the log with a least-squares cubic spline which retains the overall log character, ignoring the sinusoidal noise.
  • a zero mean stationary series is formed by subtracting the least-squares cubic spline from the data.
  • the remaining steps namely, spectral or Fourier analysis and filtering, are performed on the difference series. Recombining the filtered series with the spline restores the log data without the sinusoidal noise.
  • Figure IA is a graphical representation of a caliper log vs. depth and a gamma ray log vs. depth.
  • Figure IB is a graphical representation of a density log vs. depth and a neutron log vs. depth.
  • Figure 2A is a graphical representation of caliper spectra before and after filtering.
  • Figure 2B is a graphical representation of neutron spectra before and after filtering.
  • Figure 3 is a flow chart illustrating a method for removing borehole rugosity noise from well log data.
  • Figure 4A is a graphical representation of the log of Figure IA after filtering.
  • Figure 4B is a graphical representation of the log of Figure IB after filtering.
  • Figures 1A-2B and 4A-4B show the raw log as recorded with associated borehole "noise".
  • a distinctive "sinusoidal" borehole rugosity can be observed on the wireline caliper log identified as C in Figure IA.
  • This borehole rugosity can have an adverse effect on the response of logging devices which rely upon tool to borehole contact.
  • primary porosity and lithology identification tools such as density and neutron logging tools.
  • Figure IB shows clearly the degradation of the response of these devices.
  • FIGS. 2A and 2B illustrate the result of the Fourier analysis (amplitude spectrum) of the difference series on the wireline caliper and neutron logs. Note that the major peak for both devices occurs at the normalized wavenumber 0.26.
  • This data takes the form of a caliper log and another concomitant well log, such as a gamma ray log as described in the preferred embodiment.
  • the wavenumbers corresponding to the sinusoidal noise must be identified in the caliper log. This is accomplished by reducing the caliper log to a zero mean stationary series which contains the borehole rugosity noise.
  • the background trend in the caliper log is approximated with a piece wise cubic spline, in the least squares sense.
  • the spline nodes are chosen so that the cubic spline retains the trend information from the caliper log, but not the sinusoidal noise which has a large wavenumber. If the interval between spline nodes exceeds the reciprocal of twice the sinusoidal wavenumber, then the cubic spline will not contain the sinusoidal noise. The distance between nodes cannot be so large, however, that the cubic spline fails to accurately represent the low wavenumber components of the caliper log.
  • Standard methods for reducing the data to a zero mean series, such as removing the mean, will not preserve the low wavenumber information.
  • a zero mean stationary series is formed by subtracting the least-squares cubic spline from the caliper log.
  • This residual difference series contains only the high frequency information from the caliper log, including the sinusoidal noise.
  • the Fourier amplitude spectrum of the residual series from block 14 is computed. Sinusoidal borehole rugosity noise appears as peaks in this spectrum. Trend information resides in the low wavenumber components of the amplitude spectrum. If the background trend has been properly removed from the residual series, the Fourier amplitude spectrum will be zero in an interval about the zero wavenumber.
  • the amplitude spectrum of the residual series can be used to check the adequacy of the choice for spline nodes.
  • the amplitude spectrum of the residuals should contain no low wavenumber components, but it should have a large peak, and possibly harmonic peaks, indicating the presence of sinusoidal noise. If the spline node interval is too fine, the sinusoidal noise will be included in the cubic spline and the peak that corresponds to the sinusoidal noise will be absent from the residual amplitude spectrum. In the limit, when the spline node interval equals the sample interval for the well log data, the cubic spline will exactly coincide with the well log data and the residual amplitude spectrum will be identically zero with no peaks.
  • the residual amplitude spectrum will contain DC or low wavenumber components. These wavenumber components can be affected by the notch filter that is used to remove the sinusoidal noise. Consequently, combining the residual series with the cubic spline will corrupt low wavenumber trend information in the well log.
  • peaks that correspond to sinusoidal borehole rugosity noise are identified in the amplitude spectrum of the residuals, and their wavenumbers are saved for later use.
  • a notch filter is applied to the residual series to remove the sinusoidal noise.
  • a recursive (IIR) filter was used as a notch filter in tests of this method.
  • IIR recursive
  • y n x n - 2cos( ⁇ k 0 ) x n _ ! + x n _ 2 + 2 ⁇ cos(jrk 0 ) y n . ⁇ - ⁇ 2 y n _ 2
  • k 0 is the wavenumber to be eliminated
  • a is a design parameter that determines the location of the poles of the filter transfer function. Satisfactory results were obtained with a equal to 0.8 - 0.9. The time reverse of this filter was applied to achieve a zero phase filter response. It is understood that other types of notch filters can be applied.
  • the notch filter is designed to remove wavenumber components of the amplitude spectrum that were identified as sinusoidal borehole rugosity noise. Note that since the residual series contains no low wavenumber components, the low wavenumber trend information in the caliper log cannot be affected by this procedure. Thus all trend information is preserved in the caliper log data.
  • harmonics or other peaks are determined. Occasionally, after application of the filter, harmonic noise peaks will be present at multiples of the sinusoidal noise wavenumber. These peaks can be eliminated in a similar manner.
  • the amplitude spectrum of the filtered residual series is computed in block 28 and the procedure outlined in blocks 22 and 24 are repeated. Otherwise, the process proceeds to block 30.
  • the residual series is added to the least- squares spline determined from blocks 14-20. This restores the caliper log data, without the sinusoidal borehole rugosity noise.
  • the process steps at blocks 32-38 are to remove the background trends from other well logs that are affected by borehole rugosity. These steps duplicate the procedure of blocks 14-20 that was used to remove the background trend from the caliper log. However, these steps produce zero mean stationary residual series for each well log. Any borehole rugosity noises that are present in the well logs are also contained in the residual series.
  • the wavenumbers corresponding to the sinusoidal noise must be identified in the well log. This is accomplished by reducing the well log to a zero mean stationary series which contains the borehole rugosity noise.
  • the background trend in the well log is approximated with a piece wise cubic spline, in the least squares sense.
  • the spline nodes are chosen so that the cubic spline retains the trend information from the well log, but not the sinusoidal noise which has a large wavenumber. If the interval between spline nodes exceeds the reciprocal of twice the sinusoidal wavenumber, then the cubic spline will not contain the sinusoidal noise. The distance between nodes cannot be so large, however, that the cubic spline fails to accurately represent the low wavenumber components of the well log.
  • Standard methods for reducing the data to a zero mean series, such as removing the mean, will not preserve the low wavenumber information.
  • a zero mean stationary series is formed by subtracting the least-squares cubic spline from the well log. This residual difference series contains only the high frequency information from the well log, including the sinusoidal noise.
  • the Fourier amplitude spectrum of the residual series from block 32 is computed. Sinusoidal borehole rugosity noise appears as peaks in this spectrum. Trend information resides in the low wavenumber components of the amplitude spectrum. If the background trend has been properly removed from the residual series, the Fourier amplitude spectrum will be zero in an interval about the zero wavenumber.
  • the amplitude spectrum of the residual series can be used to check the adequacy of the choice for spline nodes.
  • the amplitude spectrum of the residuals from the well log data should contain no low wavenumber components, but it should have a large peak, and possibly harmonic peaks, indicating the presence of sinusoidal noise. If the spline node interval is too fine, the sinusoidal noise will be included in the cubic spline and the peak that corresponds to the sinusoidal noise will be absent from the residual amplitude spectrum. In the limit, when the spline node interval equals the sample interval for the well log data, the cubic spline will exactly coincide with the well log data and the residual amplitude spectrum will be identically zero with no peaks.
  • the residual amplitude spectrum will contain DC or low wavenumber components. These wavenumber components can be affected by the notch filter that is used to remove the sinusoidal noise. Consequently, combining the residual series with the cubic spline will corrupt low wavenumber trend information in the well log.
  • a potential borehole rugosity peak is identified in the amplitude spectrum of the well log residual series.
  • the potential borehole rugosity peak is compared to the borehole rugosity peaks identified with the caliper log.
  • Block 44 the result of the comparison of block 42 is evaluated. If there is a match, indicating that the peak is caused by a borehole rugosity, the process proceeds to block 46 where a notch filter is used to remove the noise. Otherwise, the process proceeds to Block 48.
  • the restored well log is obtained by combining the least squares cubic spline with the filtered residual series. Elimination of the sinusoidal noise means that the actual formation responses of the logging tools are recovered, and the logs can now be used quantitatively. Elimination of the sinusoidal noise exposes the underlying signal in the well log data without corrupting the low wavenumber trend.
  • results from blocks 30 and 52 may be displayed. This may be by any means currently in use in the art, such as a cathode ray tube, graphical printout etc..
  • Figures 4A and 4B show the restored well log obtained by combining the least squares cubic spline with the filtered difference series.
  • Caliper log C is designated by C after the process of the present invention has been performed.
  • Neutron log N is designated by N' and density log D by D' after the process of the present invention has been performed.
  • Elimination of the sinusoidal noise means that the actual formation response of the logging tools is recovered, and the logs can now be used quantitatively. As expected, the period of the rugosity dictates which logging tool responses are affected.
  • This procedure has an important log quality control application. This application is that the spectrum of the caliper log can be used to verify that peaks in the spectra of other logs are caused by borehole rugosity.

Abstract

A method for removing sinusoidal noise without affecting the response of the formation means that the actual formation responses of the logging tools are recovered, and the logs can be used quantitatively. Removal of sinusoidal noise from a log is accomplished in three steps. First, the log is reduced to a zero-mean, stationary series. Second, the wavenumber of the sinusoidal noise is identified by its peak in the Fourier amplitude spectrum (18). And third, the noise is removed by applying a zero-phase notch filter (46). In order to preserve the quantitative data integrity, the low wavenumber trend is kept. Preserving the quantitative data integrity is accomplished by approximating the log with a least-squares cubic spline which retains the overall log character, ignoring the sinusoidal noise. A zero mean stationary series is formed by subtracting the least-squares cubic spline from the data. The remaining steps, Fourier analysis and filtering are performed on the difference series. Recombining the filtered series with the spline restores the log data without the sinusoidal noise.

Description

JMETHOD FOR REMOVING BOREHOLE UGOSITY NOISE FROM WELL LOG DATA
The present invention pertains to noise removal from well log data and more particularly the removal of induced rugosity noise from well log data.
Under certain drilling conditions, often in an oil based mud environment, a distinctive "sinusoidal" borehole rugosity is observed on the wireline caliper log. This borehole rugosity can have an adverse effect on the response of logging devices which require borehole contact with the tool. Among such devices are the primary porosity and lithology identification tools, density and neutron logging tools.
The method of the present invention eliminates sinusoidal noise without affecting the response of the formation. This means that the actual formation responses of the logging tools are recovered, and the logs can be used quantitatively.
Under certain conditions, borehole rugosity can adversely affect the response of wireline logging tools. The present invention discloses a robust, fast and easy to use processing module which removes "sinusoidal" noise while preserving quantitative data integrity.
Removal of sinusoidal noise from a log is accomplished in three steps. First, the log is reduced to a zero-mean, stationary series. Second, the wavenumber of the sinusoidal noise is identified by its peak in the Fourier amplitude spectrum. And third, the noise is removed by applying a zero-phase notch filter.
In order to preserve the quantitative data integrity, the low wavenumber trend must be kept. For purposes of this discussion, the wavenumber is the reciprocal of wavelength. Preserving the quantitative data integrity is accomplished by approximating the log with a least-squares cubic spline which retains the overall log character, ignoring the sinusoidal noise. A zero mean stationary series is formed by subtracting the least-squares cubic spline from the data. The remaining steps, namely, spectral or Fourier analysis and filtering, are performed on the difference series. Recombining the filtered series with the spline restores the log data without the sinusoidal noise. An advantage of this procedure is that the results are independent of the depth interval used to perform the analysis.
Figure IA is a graphical representation of a caliper log vs. depth and a gamma ray log vs. depth.
Figure IB is a graphical representation of a density log vs. depth and a neutron log vs. depth.
Figure 2A is a graphical representation of caliper spectra before and after filtering.
Figure 2B is a graphical representation of neutron spectra before and after filtering.
Figure 3 is a flow chart illustrating a method for removing borehole rugosity noise from well log data.
Figure 4A is a graphical representation of the log of Figure IA after filtering.
Figure 4B is a graphical representation of the log of Figure IB after filtering.
The procedure of the present invention is illustrated in Figures 1A-2B and 4A-4B. Figures IA and IB show the raw log as recorded with associated borehole "noise". Under drilling conditions, as stated previously, often in an oil based mud environment, a distinctive "sinusoidal" borehole rugosity can be observed on the wireline caliper log identified as C in Figure IA. This borehole rugosity can have an adverse effect on the response of logging devices which rely upon tool to borehole contact. Among such devices are primary porosity and lithology identification tools, such as density and neutron logging tools. Figure IB shows clearly the degradation of the response of these devices. This degradation renders the log of the density tool, identified as D and the log of the neutron logging tool, identified as N, both shown in Figure IB, of qualitative use only. Figures 2A and 2B illustrate the result of the Fourier analysis (amplitude spectrum) of the difference series on the wireline caliper and neutron logs. Note that the major peak for both devices occurs at the normalized wavenumber 0.26.
Subsequent application of a notch filter to this series does not affect the low wavenumber information in the well log data. The curves in Figures 2A and 2B show the spectra of the notch filtered logs. Curve Cl illustrates the caliper log before filtering and C2 illustrates the caliper log after filtering. Similarly, curve Nl illustrates the neutron log before filtering and N2 illustrates the neutron log after filtering. The peak caused by the borehole rugosity has been eliminated. Occasionally, harmonic noise peaks will be present at multiples of the sinusoidal noise wavenumber (see caliper) . These peaks can be eliminated in a similar manner.
Referring now to Figure 3, in this embodiment, it is assumed that high frequency sinusoidal noise has been detected on the caliper log. At block 12 data is received. This data takes the form of a caliper log and another concomitant well log, such as a gamma ray log as described in the preferred embodiment.
At block 14, to remove sinusoidal borehole rugosity noise from a well log without corrupting low wavenumber trend information, the wavenumbers corresponding to the sinusoidal noise must be identified in the caliper log. This is accomplished by reducing the caliper log to a zero mean stationary series which contains the borehole rugosity noise.
To reduce the caliper log to a zero mean stationary series, the background trend in the caliper log is approximated with a piece wise cubic spline, in the least squares sense. The spline nodes are chosen so that the cubic spline retains the trend information from the caliper log, but not the sinusoidal noise which has a large wavenumber. If the interval between spline nodes exceeds the reciprocal of twice the sinusoidal wavenumber, then the cubic spline will not contain the sinusoidal noise. The distance between nodes cannot be so large, however, that the cubic spline fails to accurately represent the low wavenumber components of the caliper log.
Standard methods for reducing the data to a zero mean series, such as removing the mean, will not preserve the low wavenumber information.
At block 16, a zero mean stationary series is formed by subtracting the least-squares cubic spline from the caliper log. This residual difference series contains only the high frequency information from the caliper log, including the sinusoidal noise.
At block 18, the Fourier amplitude spectrum of the residual series from block 14 is computed. Sinusoidal borehole rugosity noise appears as peaks in this spectrum. Trend information resides in the low wavenumber components of the amplitude spectrum. If the background trend has been properly removed from the residual series, the Fourier amplitude spectrum will be zero in an interval about the zero wavenumber.
At block 20 the amplitude spectrum of the residual series can be used to check the adequacy of the choice for spline nodes. The amplitude spectrum of the residuals should contain no low wavenumber components, but it should have a large peak, and possibly harmonic peaks, indicating the presence of sinusoidal noise. If the spline node interval is too fine, the sinusoidal noise will be included in the cubic spline and the peak that corresponds to the sinusoidal noise will be absent from the residual amplitude spectrum. In the limit, when the spline node interval equals the sample interval for the well log data, the cubic spline will exactly coincide with the well log data and the residual amplitude spectrum will be identically zero with no peaks.
If the spline interval is too coarse, the residual amplitude spectrum will contain DC or low wavenumber components. These wavenumber components can be affected by the notch filter that is used to remove the sinusoidal noise. Consequently, combining the residual series with the cubic spline will corrupt low wavenumber trend information in the well log.
If the amplitude spectrum indicates that trend information remains in the residual series, then the interval between spline nodes must be reduced and the procedure in blocks 14-18 must be repeated. If the amplitude spectrum of the residual series shows that the trend information is removed, then the algorithm proceeds to block 22.
At block 22 peaks that correspond to sinusoidal borehole rugosity noise are identified in the amplitude spectrum of the residuals, and their wavenumbers are saved for later use.
At block 24 a notch filter is applied to the residual series to remove the sinusoidal noise. A recursive (IIR) filter was used as a notch filter in tests of this method. In the preferred embodiment the filter is given by:
yn = xn - 2cos(πk0) xn_! + xn_2 + 2αcos(jrk0) yn.α - α2yn_2
where k0 is the wavenumber to be eliminated, and a is a design parameter that determines the location of the poles of the filter transfer function. Satisfactory results were obtained with a equal to 0.8 - 0.9. The time reverse of this filter was applied to achieve a zero phase filter response. It is understood that other types of notch filters can be applied.
The notch filter is designed to remove wavenumber components of the amplitude spectrum that were identified as sinusoidal borehole rugosity noise. Note that since the residual series contains no low wavenumber components, the low wavenumber trend information in the caliper log cannot be affected by this procedure. Thus all trend information is preserved in the caliper log data.
At block 26 the presence of harmonics or other peaks is determined. Occasionally, after application of the filter, harmonic noise peaks will be present at multiples of the sinusoidal noise wavenumber. These peaks can be eliminated in a similar manner.
If harmonics or other peaks remain in the amplitude spectrum of the residual series, the amplitude spectrum of the filtered residual series is computed in block 28 and the procedure outlined in blocks 22 and 24 are repeated. Otherwise, the process proceeds to block 30.
At block 30 the residual series is added to the least- squares spline determined from blocks 14-20. This restores the caliper log data, without the sinusoidal borehole rugosity noise.
Having determined the wavenumber and its harmonics associated with the sinusoidal noise from the caliper log, the steps above are repeated on other logs which are known to be potentially affected by borehole rugosity. The period of the rugosity dictates which logging tool responses are affected. This procedure has an important log quality control application. The spectrum of the caliper log can be used to verify that peaks in the spectra of other logs are caused by borehole rugosity.
If the wavenumber or wavenumbers obtained from the caliper log are in agreement with the ones obtained from the other log, the sinusoidal noise present in the log is indeed due to borehole rugosity. Removal of the sinusoidal noise is then fully justified.
The process steps at blocks 32-38 are to remove the background trends from other well logs that are affected by borehole rugosity. These steps duplicate the procedure of blocks 14-20 that was used to remove the background trend from the caliper log. However, these steps produce zero mean stationary residual series for each well log. Any borehole rugosity noises that are present in the well logs are also contained in the residual series.
At block 32, to remove sinusoidal borehole rugosity noise from a well log without corrupting low wavenumber trend information, the wavenumbers corresponding to the sinusoidal noise must be identified in the well log. This is accomplished by reducing the well log to a zero mean stationary series which contains the borehole rugosity noise.
To reduce the well log to a zero mean stationary series, the background trend in the well log is approximated with a piece wise cubic spline, in the least squares sense. The spline nodes are chosen so that the cubic spline retains the trend information from the well log, but not the sinusoidal noise which has a large wavenumber. If the interval between spline nodes exceeds the reciprocal of twice the sinusoidal wavenumber, then the cubic spline will not contain the sinusoidal noise. The distance between nodes cannot be so large, however, that the cubic spline fails to accurately represent the low wavenumber components of the well log.
Standard methods for reducing the data to a zero mean series, such as removing the mean, will not preserve the low wavenumber information.
At block 34, a zero mean stationary series is formed by subtracting the least-squares cubic spline from the well log. This residual difference series contains only the high frequency information from the well log, including the sinusoidal noise.
At block 36, the Fourier amplitude spectrum of the residual series from block 32 is computed. Sinusoidal borehole rugosity noise appears as peaks in this spectrum. Trend information resides in the low wavenumber components of the amplitude spectrum. If the background trend has been properly removed from the residual series, the Fourier amplitude spectrum will be zero in an interval about the zero wavenumber.
At block 38 the amplitude spectrum of the residual series can be used to check the adequacy of the choice for spline nodes. As in the case with the caliper log data, the amplitude spectrum of the residuals from the well log data should contain no low wavenumber components, but it should have a large peak, and possibly harmonic peaks, indicating the presence of sinusoidal noise. If the spline node interval is too fine, the sinusoidal noise will be included in the cubic spline and the peak that corresponds to the sinusoidal noise will be absent from the residual amplitude spectrum. In the limit, when the spline node interval equals the sample interval for the well log data, the cubic spline will exactly coincide with the well log data and the residual amplitude spectrum will be identically zero with no peaks.
If the spline interval is too coarse, the residual amplitude spectrum will contain DC or low wavenumber components. These wavenumber components can be affected by the notch filter that is used to remove the sinusoidal noise. Consequently, combining the residual series with the cubic spline will corrupt low wavenumber trend information in the well log.
If the amplitude spectrum indicates that trend information remains in the residual series, then the interval between spline nodes must be reduced and the procedure in blocks 32-36 must be repeated. If the amplitude spectrum of the residual series shows that the trend information is removed, then the algorithm proceeds to block 40.
At block 40 a potential borehole rugosity peak is identified in the amplitude spectrum of the well log residual series.
At block 42 the potential borehole rugosity peak is compared to the borehole rugosity peaks identified with the caliper log.
At block 44, the result of the comparison of block 42 is evaluated. If there is a match, indicating that the peak is caused by a borehole rugosity, the process proceeds to block 46 where a notch filter is used to remove the noise. Otherwise, the process proceeds to Block 48.
At block 48 a determination is made whether there are other peaks or harmonics in the amplitude spectrum of the well log residuals. If there are, the amplitude spectrum of the filtered residuals is computed at block 50, and the process of blocks 40-46 are repeated. Otherwise, the process continues to block 52.
At block 52, the restored well log is obtained by combining the least squares cubic spline with the filtered residual series. Elimination of the sinusoidal noise means that the actual formation responses of the logging tools are recovered, and the logs can now be used quantitatively. Elimination of the sinusoidal noise exposes the underlying signal in the well log data without corrupting the low wavenumber trend.
Finally, at block 54, the results from blocks 30 and 52 may be displayed. This may be by any means currently in use in the art, such as a cathode ray tube, graphical printout etc..
Figures 4A and 4B show the restored well log obtained by combining the least squares cubic spline with the filtered difference series. Caliper log C is designated by C after the process of the present invention has been performed. Neutron log N is designated by N' and density log D by D' after the process of the present invention has been performed. Elimination of the sinusoidal noise means that the actual formation response of the logging tools is recovered, and the logs can now be used quantitatively. As expected, the period of the rugosity dictates which logging tool responses are affected. This procedure has an important log quality control application. This application is that the spectrum of the caliper log can be used to verify that peaks in the spectra of other logs are caused by borehole rugosity.

Claims

1. A method for removal of sinusoidal noise from a log comprising the steps of: receiving caliper log data and well log data having sinusoidal noise; reducing said caliper log data to a zero-mean, stationary series; identifying the wavenumber of said sinusoidal noise in said caliper log data by its peak in a Fourier amplitude spectrum of residuals after said reducing step; and removing said sinusoidal noise from said caliper log data by applying a filter.
2. The method according to Claim 1 wherein said step of reducing said caliper log data also includes the steps of: approximating the background trend in said caliper log data with a piece wise cubic spline; and subtracting said piece wise cubic spline from said caliper log data.
3. The method according to Claim 1 wherein said step of identifying the wavenumber of said sinusoidal noise in said caliper log data includes the steps of: computing the Fourier amplitude spectrum of a residual series in said caliper log data from said reducing step; determining wavenumber components of said amplitude spectrum; and correlating peaks in said amplitude spectrum to borehole rugosity noise.
. The method according to Claim 1 wherein said step of removing sinusoidal noise from said caliper log data includes the step of: applying a notch filter to said caliper log data.
5. The method according to Claim 4 wherein said step of applying a notch filter to said caliper log data includes the step of: using a zero-phase notch filter given by:
Yn - xn - 2cos(τrk0) X-,.! + xn_2 + 2acos(ffk0) y,^ - α2yn_2
where k0 is the wavenumber to be eliminated, and α is a design parameter that determines the location of the poles of the filter transfer function.
6. The method according to Claim 1 also including the steps of: reducing said well log data to a zero-mean, stationary series; identifying the wavenumber of said sinusoidal noise in said well log data by its peak in the Fourier amplitude spectrum; comparing said peak of said sinusoidal noise in said well log data with said peak of said sinusoidal noise in said caliper log data; and removing said sinusoidal noise from said well log data by applying a filter.
7. The method according to Claim 6 wherein said step of reducing said well log data step also includes the steps of: approximating the background trend in said well log data with a piece wise cubic spline; and subtracting said piece wise cubic spline from said well log data.
8. The method according to Claim 6 wherein said step of identifying the wavenumber of said sinusoidal noise in said well log data includes the steps of: computing the Fourier amplitude spectrum of a residual series in said well log data from said reducing step; determining wavenumber components of said amplitude spectrum; and correlating peaks in said amplitude spectrum to borehole rugosity noise.
9. The method according to claim 6 wherein said step of comparing includes the step of; determining whether said peaks in said well log data correspond to said peaks in said caliper log data.
10. The method according to Claim 6 wherein said step of removing sinusoidal noise from said well log data includes the step of: applying a notch filter to said well log data.
11. The method according to Claim 10 wherein said step of applying a notch filter to said caliper log data includes the step of: using a zero-phase notch filter given by:
yn = xn - 2cos(πk0) x^ + xn_2 + 2αcos(τrk0) y^ - α2yn_2
where k0 is the wavenumber to be eliminated, and α is a design parameter that determines the location of the poles of the filter transfer function.
12. A method for removal of sinusoidal noise from a log comprising the steps of: receiving caliper log data and well log data having sinusoidal noise; reducing said caliper log data to a zero-mean, stationary series; identifying the wavenumber of said sinusoidal noise by its peak in a Fourier amplitude spectrum of residuals after said reducing step; removing said sinusoidal noise from said caliper log data by applying a notch filter; reducing said well log data to a zero-mean, stationary series; identifying the wavenumber of said sinusoidal noise in said well log data by its peak in the Fourier amplitude spectrum; comparing said peak of said sinusoidal noise in said well log data with said peak of said sinusoidal noise in said caliper log data; and removing said sinusoidal noise from said well log data by applying a notch filter.
13. The method according to Claim 12 wherein said step of reducing said caliper log data also includes the steps of: approximating the background trend in said caliper log data with a piece wise cubic spline; and subtracting said piece wise cubic spline from said caliper log data.
14. The method according to Claim 12 wherein said step of identifying the wavenumber of said sinusoidal noise in said caliper log data includes the steps of: computing the Fourier amplitude spectrum of a residual series in said caliper log data from said reducing step; determining wavenumber components of said amplitude spectrum; and correlating peaks in said amplitude spectrum to borehole rugosity noise.
15. The method according to Claim 12 wherein said step of removing sinusoidal noise from said caliper log data by applying a notch filter includes the step of: applying a zero-phase notch filter given by:
yn = Xn ~ 2cos(τrk0) X-..J + xn_2 + 2ocos(π-k0) yn.: - αyn_2
where k0 is the wavenumber to be eliminated, and α is a design parameter that determines the location of the poles of the filter transfer function.
16. The method according to Claim 12 wherein said step of reducing said well log data step also includes the steps of: approximating the background trend in said well log data with a piece wise cubic spline; and subtracting said piece wise cubic spline from said well log data.
17. The method according to Claim 12 wherein said step of identifying the wavenumber of said sinusoidal noise in said well log data includes the steps of: computing the Fourier amplitude spectrum of a residual series in said well log data from said reducing step; determining wavenumber components of said amplitude spectrum; and correlating peaks in said amplitude spectrum to borehole rugosity noise.
18. The method according to claim 12 wherein said step of comparing includes the step of: determining whether said peaks in said well log data correspond to said peaks in said caliper log data.
19. The method according to Claim 12 wherein said step of removing sinusoidal noise from said well log data by applying a notch filter includes the step of: applying a zero-phase notch filter given by: yn = xn - 2cos (7tk0) xn_. + xn.2 + 2αcos(τrk0) yn.α - α2yn_2
where k0 is the wavenumber to be eliminated, and a is a design parameter that determines the location of the poles of the filter transfer function.
20. A method for removal of sinusoidal noise from a omprising the steps of: receiving caliper log data and well log data having sinusoidal noise; reducing said caliper log data to a zero-mean, stationary series by approximating the background trend in said caliper log data with a piece wise cubic spline and subtracting said piece wise cubic spline from said caliper log data; identifying the wavenumber of said sinusoidal noise by its peak in a Fourier amplitude spectrum of residuals after said reducing step by computing the Fourier amplitude spectrum of a residual series in said caliper log data from said reducing step, determining wavenumber components of said amplitude spectrum and correlating peaks in said amplitude spectrum to borehole rugosity noise; removing said sinusoidal noise from said caliper log data by applying a notch filter given by:
Yn = xn ~ 2cos(rk0) x^ + xn_2 + 2αcos(τrk0) y^ - α2yn_2
where k0 is the wavenumber to be eliminated, and a is a design parameter that determines the location of the poles of the filter transfer function; reducing said well log data to a zero-mean, stationary series by approximating the background trend in said well log data with a piece wise cubic spline and subtracting said piece wise cubic spline from said well log data; identifying the wavenumber of said sinusoidal noise in said well log data by its peak in the Fourier amplitude spectrum- by computing the Fourier amplitude spectrum of a residual series in said well log data from said reducing step, determining wavenumber components of said amplitude spectrum and correlating peaks in said amplitude spectrum to borehole rugosity noise; comparing said peak of said sinusoidal noise in said well log data with said peak of said sinusoidal noise in said caliper log data by determining whether said peaks in said well log data correspond to said peaks in said caliper log data; and removing said sinusoidal noise from said well log data by applying a notch filter given by:
yn = xn - 2cos(τrk0) x^ + xn_2 + 2αcos(τrk0) yn.j - α2yn_2
where k0 is the wavenumber to be eliminated, and α is a design parameter that determines the location of the poles of the filter transfer function.
PCT/US1996/000449 1996-01-16 1996-01-16 Method for removing borehole rugosity noise from well log data WO1997026559A1 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
PCT/US1996/000449 WO1997026559A1 (en) 1996-01-16 1996-01-16 Method for removing borehole rugosity noise from well log data
EP96902668A EP0991958A1 (en) 1996-01-16 1996-01-16 Method for removing borehole rugosity noise from well log data
NO983262A NO983262D0 (en) 1996-01-16 1998-07-15 Removal procedure for stay p.g.a. borehole rugosity from burn data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/US1996/000449 WO1997026559A1 (en) 1996-01-16 1996-01-16 Method for removing borehole rugosity noise from well log data

Publications (1)

Publication Number Publication Date
WO1997026559A1 true WO1997026559A1 (en) 1997-07-24

Family

ID=22254610

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US1996/000449 WO1997026559A1 (en) 1996-01-16 1996-01-16 Method for removing borehole rugosity noise from well log data

Country Status (2)

Country Link
EP (1) EP0991958A1 (en)
WO (1) WO1997026559A1 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2343254A (en) * 1998-08-25 2000-05-03 Schlumberger Holdings Removing periodic errors from well log data
GB2347753A (en) * 1999-03-08 2000-09-13 Baker Hughes Inc A method of focussing measurements assuming an inhomogeneous background signal
GB2392247A (en) * 2002-07-01 2004-02-25 Halliburton Energy Serv Inc Crosstalk reduction in downhole tools
WO2009014995A1 (en) * 2007-07-25 2009-01-29 Schlumberger Canada Limited Methods and apparatuses for formation tester data processing
WO2009118511A2 (en) * 2008-03-25 2009-10-01 Geco Technology B.V. Noise attenuation of seismic data

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4350887A (en) * 1980-11-05 1982-09-21 The United States Of America As Represented By The United States Department Of Energy Method of assaying uranium with prompt fission and thermal neutron borehole logging adjusted by borehole physical characteristics
US4809236A (en) * 1986-10-15 1989-02-28 Schlumberger Technology Corporation Method and apparatus for determining the magnitude of components of measurements made from inside a borehole
US4853903A (en) * 1988-10-19 1989-08-01 Mobil Oil Corporation Method and apparatus for removing sinusoidal noise from seismic data
US5010526A (en) * 1990-07-30 1991-04-23 Mobil Oil Corporation Method of digital filtering

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU555850B2 (en) * 1983-09-26 1986-10-09 Exploration Logging Inc. Noise subtraction filter
US4885722A (en) * 1988-01-19 1989-12-05 Mobil Oil Corporation Method for removing low-frequency noise from records with data impulse
US5479251A (en) * 1993-12-30 1995-12-26 Corning Incorporated Methods for determining optical properties of optical waveguide fibers using an optical time domain reflectometer

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4350887A (en) * 1980-11-05 1982-09-21 The United States Of America As Represented By The United States Department Of Energy Method of assaying uranium with prompt fission and thermal neutron borehole logging adjusted by borehole physical characteristics
US4809236A (en) * 1986-10-15 1989-02-28 Schlumberger Technology Corporation Method and apparatus for determining the magnitude of components of measurements made from inside a borehole
US4853903A (en) * 1988-10-19 1989-08-01 Mobil Oil Corporation Method and apparatus for removing sinusoidal noise from seismic data
US5010526A (en) * 1990-07-30 1991-04-23 Mobil Oil Corporation Method of digital filtering

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP0991958A4 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2343254A (en) * 1998-08-25 2000-05-03 Schlumberger Holdings Removing periodic errors from well log data
GB2343254B (en) * 1998-08-25 2001-05-30 Schlumberger Holdings Removing periodic errors from well log data
GB2347753A (en) * 1999-03-08 2000-09-13 Baker Hughes Inc A method of focussing measurements assuming an inhomogeneous background signal
GB2347753B (en) * 1999-03-08 2003-09-24 Baker Hughes Inc Electromagnetic induction well logging
AU2003204779A8 (en) * 2002-07-01 2010-04-29 Halliburton Energy Services, Inc. Method and apparatus employing phase cycling for reducing crosstalk in downhole tools
US6765384B2 (en) 2002-07-01 2004-07-20 Halliburton Energy Services, Inc. Method and apparatus employing phase cycling for reducing crosstalk in downhole tools
GB2392247B (en) * 2002-07-01 2006-02-15 Halliburton Energy Serv Inc Crosstalk reduction in downhole tools
AU2003204779B2 (en) * 2002-07-01 2007-09-20 Halliburton Energy Services, Inc. Method and apparatus employing phase cycling for reducing crosstalk in downhole tools
GB2392247A (en) * 2002-07-01 2004-02-25 Halliburton Energy Serv Inc Crosstalk reduction in downhole tools
AU2003204779B8 (en) * 2002-07-01 2010-04-29 Halliburton Energy Services, Inc. Method and apparatus employing phase cycling for reducing crosstalk in downhole tools
WO2009014995A1 (en) * 2007-07-25 2009-01-29 Schlumberger Canada Limited Methods and apparatuses for formation tester data processing
GB2464027A (en) * 2007-07-25 2010-04-07 Schlumberger Holdings Methods and apparatuses for formation tester data processing
GB2464027B (en) * 2007-07-25 2011-11-30 Schlumberger Holdings Methods and apparatuses for formation tester data processing
WO2009118511A2 (en) * 2008-03-25 2009-10-01 Geco Technology B.V. Noise attenuation of seismic data
WO2009118511A3 (en) * 2008-03-25 2009-11-26 Geco Technology B.V. Noise attenuation of seismic data
US8855934B2 (en) 2008-03-25 2014-10-07 Westerngeco L.L.C. Noise attenuation of seismic data

Also Published As

Publication number Publication date
EP0991958A4 (en) 2000-04-12
EP0991958A1 (en) 2000-04-12

Similar Documents

Publication Publication Date Title
US5579248A (en) Method for removing borehole rugosity noise from well log data
US7457422B2 (en) Method and implementation for detecting and characterizing audible transients in noise
CA2923888C (en) Dynamic threshold methods, systems, computer readable media, and program code for filtering noise and restoring attenuated high-frequency components of acoustic signals
GB2196167A (en) Encoded marking of a recording signal
NO992355D0 (en) Hearing-adapted quality assessment of audio signals that have undergone coding and decoding
WO1997026559A1 (en) Method for removing borehole rugosity noise from well log data
US4885722A (en) Method for removing low-frequency noise from records with data impulse
Evans Running median filters and a general despiker
US5672867A (en) Method for filtering gamma ray well logging tool response to enhance vertical detail while suppressing statistical noise
CA2172470A1 (en) Method for removing borehole rugosity noise from well log data
CN110007346A (en) Adaptive intermediate value very noisy based on wavelet decomposition interferes seismic data noise attenuation method
DE69737875T2 (en) Demodulation device for data-modulated carriers
CN107610055A (en) The noise measuring of Fourier transform spectrometer, interference pattern and suppressing method
US6049757A (en) Parametric modeling of well log data to remove periodic errors
US6853933B2 (en) Method of identifying spectral impulses for Rj Dj separation
Hargrove et al. Improvement of algorithm for peak detection in automatic gas chromatography-mass spectrometry data processing
Liu et al. Analyzing engineering surface texture using wavelet filter
US6807488B2 (en) Method of identification of non-primary events in seismic data
CN112764109A (en) Method and device for separating and extracting dipole transverse wave reflected wave
CN112925024B (en) Square wave pressing method for seismic recording
Nieto et al. Removal of borehole induced noise from well logs
CN115236753A (en) Compact sandstone reservoir NMR echo data variation modal domain adaptive filtering method
CN110779908A (en) Method and system for determining signal-to-noise ratio of Raman spectrum
JP3293713B2 (en) Discrimination method of water leak sound by high precision frequency analysis method
CN117395110A (en) Open loop demodulation method and device for short frame burst communication with low signal-to-noise ratio

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): NO

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE CH DE DK ES FR GB GR IE IT LU MC NL PT SE

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
WWE Wipo information: entry into national phase

Ref document number: 1996902668

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 1996902668

Country of ref document: EP