US4467425A
Expired / Lapsed
US
⏱️ 0.0 years remaining
Deconvolution filter for induction log processing
Assignee: Schlumberger Technology Corp
•
Filing: 1982-01-12
•
Granted: 1984-08-21
•
Expires: 2002-01-12
Abstract
Abstract A method and system for processing induction measurements of sub-surface formations taken by an induction logging system at various depths in a borehole is disclosed. The method reduces the unwanted contributions in each measurement from induction currents flowing in the formations spaced apart from each measurement depth. The spatial domain sonde response function obtained at zero conductivity is Fourier transformed and truncated at a spatial frequency less than the spatial frequency at which the transformed function first goes to zero. A Kaiser window transfer function is then selected and the inverse Fourier transform of the ratio of the Kaiser window function to the transformed truncated sonde response function is formed to obtain a deconvolution filter function which when convolved with the sonde response function reduces the non-zero sidelobes in the system response function and thus reduces the shoulder effect.
Inventors
- 👤 Richard T. Schaefer
- 👤 Clinton H. Dutcher, Jr.
Bibliographic Details
Application Number
—
Jurisdiction
US
Data Source
google_patents
Last Synchronized
2026-10-28 07:58
CPC Classifications
G; G01; G01V; G01V3/00; G01V3/18; G01V3/26; G01V3/28
Claims Analysis
Total Claims: 27
Claim 1
Independent Claim
208 words
A method of processing induction measurements of sub-surface formations taken by an induction logging system at various depths in a borehole in the earth, the logging system having a sonde response function which varies with the conductivity of the sub-surface formations being investigated, the spatial domain sonde response function characterized by a main lobe and side-lobes, said method reducing unwanted contributions in each measurement from formation currents flowing in formations spaced apart from each measurement depth, said method comprising the steps of: (a) windowing the sonde response function obtained at zero conductivity by truncating its Fourier transform at a frequency less than the frequency at which the transformed function first goes to zero; (b) selecting a target transfer function containing frequencies within the window; (c) generating a spatial domain filter response function from the truncated transformed zero conductivity sonde response function and the target transfer function which when convolved with the zero conductivity spatial domain sonde response function results in a system response function having minimum sidelobes; and (d) convolving the spatial domain filter response function with the log measurements to obtain a processed log in which the unwanted contributions in each measurement from formation currents flowing in formations spaced apart to each measurement depth point are reduced.
Claim 2
Independent Claim
208 words
A method of processing induction measurements of sub-surface formations taken by an induction logging system at various depths in a borehole in the earth, the logging system having a sonde response function which varies with the conductivity of the sub-surface formations being investigated, the spatial domain sonde response function characterized by a main lobe and side-lobes, said method reducing unwanted contributions in each measurement from formation currents flowing in formations spaced apart from each measurement depth, said method comprising the steps of: (a) windowing the sonde response function obtained at zero conductivity by truncating its Fourier transform at a frequency less than the frequency at which the transformed function first goes to zero; (b) selecting a target transfer function containing frequencies within the window; (c) generating a spatial domain filter response function from the truncated transformed zero conductivity sonde response function and the target transfer function which when convolved with the zero conductivity spatial domain sonde response function results in a system response function having minimum sidelobes; and (d) convolving the spatial domain filter response function with the log measurements to obtain a processed log in which the unwanted contributions in each measurement from formation currents flowing in formations spaced apart to each measurement depth point are reduced.
Claim 3
Dependent Claim
41 words
The method of claim 1 wherein the step of selecting a target transfer function involves selecting a target transfer function so that the spatial domain filter response function produces minimum ripple in the processed log for step changes in formation conductivity.
Claim 4
Independent Claim
30 words
The method of claims 1 or 2 wherein the target transfer function is a Kaiser window function containing the frequencies remaining in the truncated transformed zero conductivity sonde response function.
Claim 5
Dependent Claim
31 words
The method of claim 1 wherein the spatial domain filter response function is determined from the ratio of the target transfer function to the truncated transformed zero conductivity sonde response function.
Claim 6
Independent Claim
197 words
A method of processing induction measurements of sub-surface formation materials taken by an induction logging system at various depths in a borehole in the earth, the system having a spatial domain sonde response function which varies with the conductivity of the sub-surface formations being investigated and is characterized by a main lobe and non-zero sidelobes, the method reducing unwanted contributions in each measurement from formation current flowing in spaced apart formations as a result of the non-zero sidelobes in the sonde response function, the method comprising the steps of: (a) Fourier transforming the spatial domain sonde response function obtained at zero conductivity; (b) truncating the transformed sonde response function at a frequency less than the frequency at which the transformed function first goes to zero; (c) determining a target transfer function for the frequencies remaining in the truncated transformed function; (d) determining a spatial domain filter response function from the ratio of the target transfer function to the transformed truncated sonde response function; and (e) convolving the spatial domain filter response function obtained in step (d) with the log measurements to remove from each log measurement the unwanted contributions from formation currents flowing in spaced apart formations.
Claim 7
Independent Claim
197 words
A method of processing induction measurements of sub-surface formation materials taken by an induction logging system at various depths in a borehole in the earth, the system having a spatial domain sonde response function which varies with the conductivity of the sub-surface formations being investigated and is characterized by a main lobe and non-zero sidelobes, the method reducing unwanted contributions in each measurement from formation current flowing in spaced apart formations as a result of the non-zero sidelobes in the sonde response function, the method comprising the steps of: (a) Fourier transforming the spatial domain sonde response function obtained at zero conductivity; (b) truncating the transformed sonde response function at a frequency less than the frequency at which the transformed function first goes to zero; (c) determining a target transfer function for the frequencies remaining in the truncated transformed function; (d) determining a spatial domain filter response function from the ratio of the target transfer function to the transformed truncated sonde response function; and (e) convolving the spatial domain filter response function obtained in step (d) with the log measurements to remove from each log measurement the unwanted contributions from formation currents flowing in spaced apart formations.
Claim 8
Dependent Claim
43 words
The method of claim 5 wherein the target transfer function is selected as that the spatial domain filter response function produces a minimum ripple in the final processed a log measurements in response to step changes in conductivity of the formations being investigated.
Claim 9
Independent Claim
30 words
The method of claims 5 or 6 wherein the target transfer function is a Kaiser window function containing the frequencies remaining in the truncated transformed zero conductivity sonde response function.
Claim 10
Independent Claim
212 words
A method of processing induction measurements of sub-surface formations taken by an induction logging system at various depths in a borehole in the earth, the logging system having a spatial domain sonde response function which varies with the conductivity of the sub-surface formations being investigated and is characterized by a main lobe and non-zero sidelobes, said method reducing both the unwanted contributions in each measurement from formation currents flowing in formations spaced apart from each measurement depth and the effects of variations in the sonde response function with formation conductivity, the unwanted contributions in each measurement from spaced apart formations resulting from non-zero sidelobes in the sonde response function, said method comprising the steps of: (a) determining for each of a plurality of formation conductivity ranges a deconvolution filter response function based on the sonde response function where each filter function, when convolved with the measurements of formations having conductivities in the range of that filter, reduces the unwanted contributions in the measurements from formation currents flowing in formations spaced apart from the depth at which a measurement was made; and (b) processing each measurement by convolving the measurements with a deconvolution filter response function selected from among the plurality of filter functions according to the magnitude of the measurements being convolved.
Claim 11
Independent Claim
212 words
A method of processing induction measurements of sub-surface formations taken by an induction logging system at various depths in a borehole in the earth, the logging system having a spatial domain sonde response function which varies with the conductivity of the sub-surface formations being investigated and is characterized by a main lobe and non-zero sidelobes, said method reducing both the unwanted contributions in each measurement from formation currents flowing in formations spaced apart from each measurement depth and the effects of variations in the sonde response function with formation conductivity, the unwanted contributions in each measurement from spaced apart formations resulting from non-zero sidelobes in the sonde response function, said method comprising the steps of: (a) determining for each of a plurality of formation conductivity ranges a deconvolution filter response function based on the sonde response function where each filter function, when convolved with the measurements of formations having conductivities in the range of that filter, reduces the unwanted contributions in the measurements from formation currents flowing in formations spaced apart from the depth at which a measurement was made; and (b) processing each measurement by convolving the measurements with a deconvolution filter response function selected from among the plurality of filter functions according to the magnitude of the measurements being convolved.
Claim 12
Dependent Claim
108 words
The method of claim 8 wherein the step of determining a deconvolution filter response function for a given formation conductivity range comprises the steps of: (a) windowing the sonde response function for the given conductivity range by truncating its Fourier transform at a frequency less than the frequency at which the transformed function first goes to zero; (b) selecting a target transfer function containing the window frequencies; and (c) generating a spatial domain filter response function from the truncated transformed sonde response function and the target transfer functions which when convolved with the given conductivity range sonde response function results in a system response function having reduced sidelobes.
Claim 13
Dependent Claim
46 words
The method of claim 9 wherein the step of selecting a target transfer function includes selecting a transfer function so that the spatial domain filter transfer function produces minimum ripple in the processed log in response to step changes in conductivity of the formations being investigated.
Claim 14
Independent Claim
28 words
The method of claims 9 or 10 wherein the target transfer function is a Kaiser window function containing the frequencies remaining in the truncated transformed sonde response function.
Claim 15
Dependent Claim
29 words
The method of claim 9 wherein the spatial domain filter response function is determined from the ratio of the target transfer function to the truncated transformed sonde response function.
Claim 16
Independent Claim
245 words
A method of processing induction measurements of sub-surface formations taken by an induction logging system at various depths in a borehole drilled into the earth, the logging system having a spatial domain sonde response function which varies with the conductivity of the sub-surface formations being investigated, said method reducing the effects of variations in the sonde response function with formation conductivity, said method comprising the steps of: (a) determining at each of a plurality of conductivity values a deconvolution filter response function based on the sonde response function at each value of conductivity which when convolved with the measurements reduce the unwanted contributions in each measurement to formation currents flowing in formations spaced apart from the measurement depth; (b) determining a digital filter for implementing the deconvolution filter function of step (a), said digital filter having a plurality of coefficients; (c) implementing each deconvolution filter by determining the coefficients of said digital filter; (d) curve fitting a function for each coefficient to the plurality of values for that coefficient obtained in step (c) to obtain coefficient functions which vary as a function of a conductivity variable; and (e) processing the measurements to reduce the effects of variations in the logging system transfer function with formation conductivity by convolving the measurements with a deconvolution filter implemented by said digital filter in which the coefficients of the convolved digital filter are determined by solving the coefficient functions for each coefficient using the measurements as the conductivity variable.
Claim 17
Independent Claim
245 words
A method of processing induction measurements of sub-surface formations taken by an induction logging system at various depths in a borehole drilled into the earth, the logging system having a spatial domain sonde response function which varies with the conductivity of the sub-surface formations being investigated, said method reducing the effects of variations in the sonde response function with formation conductivity, said method comprising the steps of: (a) determining at each of a plurality of conductivity values a deconvolution filter response function based on the sonde response function at each value of conductivity which when convolved with the measurements reduce the unwanted contributions in each measurement to formation currents flowing in formations spaced apart from the measurement depth; (b) determining a digital filter for implementing the deconvolution filter function of step (a), said digital filter having a plurality of coefficients; (c) implementing each deconvolution filter by determining the coefficients of said digital filter; (d) curve fitting a function for each coefficient to the plurality of values for that coefficient obtained in step (c) to obtain coefficient functions which vary as a function of a conductivity variable; and (e) processing the measurements to reduce the effects of variations in the logging system transfer function with formation conductivity by convolving the measurements with a deconvolution filter implemented by said digital filter in which the coefficients of the convolved digital filter are determined by solving the coefficient functions for each coefficient using the measurements as the conductivity variable.
Claim 18
Dependent Claim
129 words
The method of claim 13 wherein the step of determining a deconvolution filter response function at a given conductivity value comprises the steps of: (a) windowing the spatial domain sonde response function for the given conductivity value by truncating its Fourier transform at a frequency less than the frequency at which the transformed function first goes to zero; (b) selecting a target transfer function containing the window frequencies; and (c) generating a spatial domain filter response function from the truncated transformed sonde response function and the target transfer function which when convolved with the measurements reduces the unwanted contributions in each measurement from formation currents flowing in formations spaced apart from the depth at which the measurement was made, the filter response function comprising the deconvolution filter response function.
Claim 19
Dependent Claim
52 words
The method of claim 14 wherein the step of selecting a target transfer function involves selecting a transfer function so the spatial domain filter response function when convolved with the measurements results in minimum ripple in the processed log in response to a step change in conductivity in the formations being investigated.
Claim 20
Independent Claim
28 words
The method of claims 14 or 15 wherein the target transfer function is a Kaiser window function containing the frequencies remaining in the truncated transformed sonde response function.
Claim 21
Dependent Claim
29 words
The method of claim 14 wherein the spatial domain filter response function is determined from the ratio of the target transfer function to the truncated transformed sonde response function.
Claim 22
Independent Claim
211 words
A method of processing induction measurements of sub-surface formations taken by an induction logging system at various depths in a borehole in the earth, the logging system having a spatial domain sonde response function g(z, σ F ) at each value of formation conductivity σ F (z) where g(z, σ F ) varies with the σ F of the sub-surface formation being investigated, and is further characterized by a main lobe and non-zero sidelobes, said method reducing the unwanted contributions in the measurements from formation currents flowing in formations spaced apart from each measurement depth, said method comprising the steps of: (a) obtaining a deconvolution filter h(z) from the inverse Fourier transform of H(ω) where, ##EQU38## T(ω) comprising a target transfer function for frequency ω≦ω c so that when h(z) is convolved with g(z, σ F ) results in a system response function having minimum sidelobes, G(ω) comprising the Fourier transform of g(z, σ F ) at σ F (z)=0, and ω c comprising an upper frequency less than the frequency at which G(ω) first goes to zero; and (b) convolving h(z) with the measurements to obtain a processed log in which the unwanted contributions in each measurement from formation currents flowing in spaced apart formations to each measurement are reduced.
Claim 23
Independent Claim
211 words
A method of processing induction measurements of sub-surface formations taken by an induction logging system at various depths in a borehole in the earth, the logging system having a spatial domain sonde response function g(z, σ F ) at each value of formation conductivity σ F (z) where g(z, σ F ) varies with the σ F of the sub-surface formation being investigated, and is further characterized by a main lobe and non-zero sidelobes, said method reducing the unwanted contributions in the measurements from formation currents flowing in formations spaced apart from each measurement depth, said method comprising the steps of: (a) obtaining a deconvolution filter h(z) from the inverse Fourier transform of H(ω) where, ##EQU38## T(ω) comprising a target transfer function for frequency ω≦ω c so that when h(z) is convolved with g(z, σ F ) results in a system response function having minimum sidelobes, G(ω) comprising the Fourier transform of g(z, σ F ) at σ F (z)=0, and ω c comprising an upper frequency less than the frequency at which G(ω) first goes to zero; and (b) convolving h(z) with the measurements to obtain a processed log in which the unwanted contributions in each measurement from formation currents flowing in spaced apart formations to each measurement are reduced.
Claim 24
Dependent Claim
43 words
The method of claim 18 wherein T(ω) is a Kaiser window function which when its spatial domain inverse transform h(z) is convolved with a log having a step change in conductivity, a processed log having minimum ripple in response to the step results.
Claim 25
Dependent Claim
308 words
An induction logging system for generating an induction log from conductivity measurements of sub-surface formations taken by a an induction logging tool at various depths in a borehole in the earth, the logging tool having a sonde which has a spatial domain sonde response function which varies with the conductivity of the sub-surface formations being investigated and is characterized by a main lobe and non-zero sidelobes, said system reducing the unwanted contributions in each measurement from formation currents flowing in formations spaced apart from each measurement depth (shoulder effect), the system comprising: (a) a means for positioning the logging tool at various depths in the borehole, said positioning means including a wireline cable for transmitting signals to and from the tool, the tool including a means for generating an in-phase component for each conductivity measurement where the sonde response function transforms the formation conductivity distribution into the in-phase components; (b) a processing unit responsive to the conductivity measurements for processing the measurements to reduce shoulder effect, said processing unit including, a deconvolution filtering means responsive to the in-phase component measurements for generating deconvolved conductivity measurements in which shoulder effect has been reduced, the filter response function of said deconvolution filtering means derived from the sonde response function obtained at zero formation conductivity by truncating the Fourier transform of the zero conductivity sonde response function at a frequency less than the frequency at which the transformed function first goes to zero and selecting a target transfer function so that the inverse Fourier transform of the ratio of the target transfer function to the truncated transformed zero conductivity sonde response function, when convolved with the sonde response function, produces a system response function having minimum sidelobes; and (c) a recorder responsive to said processing unit and said positioning means for recording the deconvolved measurement as the log of the system.
Claim 26
Dependent Claim
308 words
An induction logging system for generating an induction log from conductivity measurements of sub-surface formations taken by a an induction logging tool at various depths in a borehole in the earth, the logging tool having a sonde which has a spatial domain sonde response function which varies with the conductivity of the sub-surface formations being investigated and is characterized by a main lobe and non-zero sidelobes, said system reducing the unwanted contributions in each measurement from formation currents flowing in formations spaced apart from each measurement depth (shoulder effect), the system comprising: (a) a means for positioning the logging tool at various depths in the borehole, said positioning means including a wireline cable for transmitting signals to and from the tool, the tool including a means for generating an in-phase component for each conductivity measurement where the sonde response function transforms the formation conductivity distribution into the in-phase components; (b) a processing unit responsive to the conductivity measurements for processing the measurements to reduce shoulder effect, said processing unit including, a deconvolution filtering means responsive to the in-phase component measurements for generating deconvolved conductivity measurements in which shoulder effect has been reduced, the filter response function of said deconvolution filtering means derived from the sonde response function obtained at zero formation conductivity by truncating the Fourier transform of the zero conductivity sonde response function at a frequency less than the frequency at which the transformed function first goes to zero and selecting a target transfer function so that the inverse Fourier transform of the ratio of the target transfer function to the truncated transformed zero conductivity sonde response function, when convolved with the sonde response function, produces a system response function having minimum sidelobes; and (c) a recorder responsive to said processing unit and said positioning means for recording the deconvolved measurement as the log of the system.
Claim 27
Dependent Claim
16 words
The logging system of claim 20 wherein said processing unit is a general purpose digital computer.
References Cited (Backward Citations)
US3340464A
Cited
US2582314A
Cited
US2790138A
Cited
US3067383A
Cited
US3166709A
Cited
US3056917A
Cited
US3147429A
Cited
US3179879A
Cited
US3226633A
Cited
US3750009A
Cited
US4314339A
Cited
USPTO Annuity & Maintenance Windows
| Annuity Stage | Window Opens | Payment Deadline | Surcharge Period Until | Estimated Fee | Status |
|---|---|---|---|---|---|
| 3.5 Years (Year 4) | 1987-08-21 | 1988-02-20 | 1988-08-21 | $2000.00 | Window Closed / Lapsed |
| 7.5 Years (Year 8) | 1991-08-21 | 1992-02-20 | 1992-08-21 | $3760.00 | Window Closed / Lapsed |
| 11.5 Years (Year 12) | 1995-08-21 | 1996-02-20 | 1996-08-21 | $7700.00 | Window Closed / Lapsed |