WO2012125161A1 - Acoustic signal processing using model-based adaptive filtering - Google Patents
Acoustic signal processing using model-based adaptive filtering Download PDFInfo
- Publication number
- WO2012125161A1 WO2012125161A1 PCT/US2011/028516 US2011028516W WO2012125161A1 WO 2012125161 A1 WO2012125161 A1 WO 2012125161A1 US 2011028516 W US2011028516 W US 2011028516W WO 2012125161 A1 WO2012125161 A1 WO 2012125161A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- slowness
- parameter
- formation
- frequency filter
- determining
- Prior art date
Links
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 60
- 238000012545 processing Methods 0.000 title claims abstract description 35
- 238000001914 filtration Methods 0.000 title claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 87
- 238000000034 method Methods 0.000 claims abstract description 42
- 239000012530 fluid Substances 0.000 claims description 12
- 230000006835 compression Effects 0.000 claims description 7
- 238000007906 compression Methods 0.000 claims description 7
- 230000003247 decreasing effect Effects 0.000 claims 4
- 238000005755 formation reaction Methods 0.000 description 68
- 230000005284 excitation Effects 0.000 description 62
- 230000008569 process Effects 0.000 description 10
- 238000005553 drilling Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 238000003672 processing method Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 239000006185 dispersion Substances 0.000 description 3
- 238000012805 post-processing Methods 0.000 description 3
- 238000003908 quality control method Methods 0.000 description 3
- 230000004075 alteration Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000010304 firing Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000005404 monopole Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
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. for interpretation or for event detection
- 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/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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
- G01V2210/21—Frequency-domain filtering, e.g. band pass
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
Definitions
- Acoustic waveform processing is frequently used to estimate properties of a geological formation.
- Acoustic waveform data may be collected using an acoustic logging tool.
- the acoustic logging tool typically includes an acoustic source (transmitter), and a set of receivers that are spaced several inches or feet apart.
- An acoustic signal is transmitted by the acoustic source, and received at the receivers of the tool, which are spaced apart from the acoustic source. Measurements are repeated every few inches as the tool is drawn up (or down) the borehole.
- the acoustic signal from the acoustic source travels through the formation adjacent to the borehole to the receiver array. The arrival times and perhaps other characteristics of the received signals are recorded.
- compressional wave P-wave
- shear wave S-wave
- Stoneley wave arrivals are observed by the receivers and are processed into an acoustic waveform data format.
- the acoustic waveform data can be processed to determine other formation characteristics such as formation slowness (which is the inverse of acoustic speed), from which pore pressure, porosity, and other formation properties can be determined.
- Figure 1 A illustrates acquiring acoustic waveform data in a wireline environment.
- Figure IB illustrates acquiring acoustic waveform data in a logging while drilling environment.
- Figures 2A-2C illustrate dispersion slowness behavior and the corresponding formation excitation from a dipole source.
- Figure 3 is a chart illustrating the steps for determining the polynomial coefficients used for an adaptive filter.
- Figure 4 is a chart illustrating the steps of a real-time processing embodiment of the acoustic waveform processing method disclosed herein.
- Figure 5 is a chart illustrating the steps of a real-time processing embodiment with adaptive waveform firing of the acoustic waveform signal processing method disclosed herein.
- Figure 6 is a chart illustrating the processing steps of a post-processing embodiment of the acoustic waveform processing method disclosed herein.
- Figures 7A-7H illustrate results from analyzing nearly 2,376 different models with variation of borehole size, borehole fluid properties, and all possible earth formations.
- Figures 8A illustrates the polynomials fitting the minimum excitation and peak excitation frequencies of various models of formation compressional slowness.
- Figures 8B illustrates the polynomials fitting the minimum excitation and peak excitation frequencies of various models of formation shear slowness.
- Figure 9A is a graph of the lower bound and the upper bound of an adaptive frequency filter for the compressional slowness.
- Figure 9B is a graph of the lower bound and the upper bound of an adaptive frequency filter for the shear slowness.
- Figure 1 OA illustrates a fitting surface for the lower bound of an adaptive frequency filter.
- Figure 1 OB illustrates a fitting surface for the upper bound of an adaptive frequency filter.
- a drilling platform supports a derrick 102 having a traveling block for raising and lowering an acoustic logging tool 104.
- the acoustic logging tool 104 has one or more acoustic sources (“transmitters”) that are spaced apart from one or more receivers.
- the acoustic logging tool 104 is connected by a wire 106 to a computer processing system positioned above ground in a mobile unit 108.
- the computer processing system is capable of sending commands to the acoustic logging tool 104, and receiving data from the acoustic logging tool 104.
- acoustic waveform data received from the acoustic logging tool 104 is processed in real-time by the computer processing system. In another embodiment, the acoustic waveform data received from the acoustic logging tool 104 is received by the computer processing system and stored for further processing at a remote computing center 1 10.
- FIG. IB an illustrative logging while drilling (“LWD”) environment is shown.
- a drilling platform supports a derrick 112 having a traveling block for raising and lowering a an acoustic logging tool 114 that is on a drill string 1 16.
- the acoustic logging tool 114 has one or more acoustic sources (“transmitters”) that are spaced apart from one or more receivers.
- the acoustic logging tool 114 has a computer processing system connected to the transmitters and receivers. In one embodiment, the computer processing system is capable of performing real- time processing of acoustic waveform data acquired by the receivers.
- the computer processing system stores the acoustic waveform data in a memory for later retrieval and processing at a remote computing center 1 10.
- the computer processing system may be capable of sending and receiving data to a computer processing system at the surface or a remote computing center 110 using mud pulses, or other data transmission techniques. Data may also be retrieved from the acoustic logging tool 1 14 when it is brought to the surface. The data may be transmitted to a remote computing center 1 10 for further processing.
- an acoustic wave propagates along a borehole depends upon the type of source excitation generated by the acoustic source.
- an explosive acoustic source at the center of the borehole produces a wave that propagates in a monopole or hoop mode.
- a dipole acoustical source produces a wave that propagates in a flexural mode along the borehole.
- a quadrupole acoustic source produces a wave that propagates in a screw mode along the borehole.
- Each mode exhibits different characteristics that must be considered when using acoustic waveform data to calculate formation properties, particularly when the formation velocities are slower than the borehole fluid.
- all of the excitation modes are dispersive, which must be considered in order to calculate formation properties accurately.
- STC processing can be used to separate compressional, shear, flexural, and Stoneley wave arrivals.
- the measured signal is time window "filtered” and stacked, and a semblance function is computed.
- the semblance function relates the presence or absence of an arrival with a particular assumed slowness and particular assumed arrival time. If the assumed slowness and arrival time do not coincide with that of the measured arrival, the semblance takes on a smaller value. Consequently, arrivals in the received waveforms manifest themselves as local peaks in a plot of semblance versus slowness and arrival time.
- the dispersive nature of the acoustic waves result in estimates of formation slowness that are far from the actual formation slowness.
- Figures 2A-2C illustrate dispersion slowness behavior and the corresponding formation excitation from a dipole source.
- Figure 2A illustrates a model of a fast formation where the compression slowness (DTC) is 50 microseconds per foot, and the shear slowness (DTS) is 85 microseconds per foot.
- DTC compression slowness
- DTS shear slowness
- a formation is "fast" when the formation's shear velocity is greater than the acoustic velocity (which is the inverse of slowness) of the fluid in the borehole.
- Three boreholes were excited using an acoustic logging tool with a dipole acoustic source at various frequencies at constant power. The slowness at each frequency was then estimated using a theoretical modal analysis of these three boreholes with appropriate parameters.
- Lines 202, 206, and 210 show the estimated slowness for the 16-inch, 9-inch, and 6-inch boreholes respectively.
- the excitation energy of the formation at each frequency was also calculated from the same theoretical modal analysis.
- Lines 204, 208, and 212 display the relative energy levels detected in the formation from the 16-inch, 9-inch, and 6-inch boreholes respectively (energy scale not shown). From this data, the minimum excitation frequency and the peak excitation frequency are determined for each borehole.
- the minimum excitation frequency is the lowest frequency at which an excitation can be detected in the formation using an acoustic source at a given power. Typically, excitation in the formation can be detected at a lower frequency only if the power of the acoustic source is significantly increased.
- the peak excitation frequency is the frequency at which the most energy can be anticipated from the formation to the acoustic logging tool.
- the minimum excitation frequency was observed at approximately 1.8 kHz, and the maximum excitation frequency was observed at approximately 2.2 kHz.
- the minimum excitation frequency was observed at approximately 3 kHz, and the maximum excitation frequency was observed at approximately 4 kHz.
- the minimum excitation frequency was observed at approximately 4.7 kHz, and the maximum excitation frequency was observed at approximately 6 kHz.
- Figure 2B illustrates a model of a medium formation where the compression slowness (DTC) is 100 microseconds per foot, and the shear slowness (DTS) is 215 microseconds per foot.
- Lines 220, 224, and 228 show the calculated formation slowness for the 16-inch, 9-inch, and 6- inch boreholes respectively.
- Lines 222, 226, and 230 display the relative excitation levels detected in the formation from the 16-inch, 9-inch, and 6-inch boreholes respectively (energy scale not shown). For the 16-inch borehole, the minimum excitation frequency was observed at approximately 1 kHz, and the maximum excitation frequency was observed at approximately 2 kHz.
- the minimum excitation frequency was observed at approximately 1.5 kHz, and the maximum excitation frequency was observed at approximately 4 kHz.
- the minimum excitation frequency was observed at approximately 2 kHz, and the maximum excitation frequency was observed at approximately 6.5 kHz.
- Figure 2C illustrates a model of a slow formation where the compression slowness (DTC) is 170 microseconds per foot, and the shear slowness (DTS) is 700 microseconds per foot.
- DTC compression slowness
- DTS shear slowness
- a formation is "slow" when the formation's shear velocity is lower than the acoustic velocity of the fluid in the borehole.
- Lines 240, 244, and 248 show the calculated slowness for the 16-inch, 9- inch, and 6-inch boreholes respectively.
- Lines 242, 246, and 250 display the relative excitation levels detected in the formation from the 16-inch, 9-inch, and 6-inch boreholes respectively (energy scale not shown).
- the minimum excitation frequency was observed at approximately 0.4 kHz, and the maximum excitation frequency was observed at approximately 0.9 kHz.
- the minimum excitation frequency was observed at approximately 0.5 kHz, and the maximum excitation frequency was observed at approximately 1.2 kHz.
- the minimum excitation frequency was observed at approximately 0.7 kHz, and the maximum excitation frequency was observed at approximately 1.8 kHz.
- Formation slowness calculations using slowness-time coherence methods are frequently performed at higher acoustic frequencies because the signal strength is greater.
- these estimates can overstate the slowness of a formation.
- theoretical models have been used to correct for the dispersive slowness behavior exhibited at higher acoustic frequencies.
- Theoretical models require accurate information about the diameter of the borehole and properties of the borehole fluid, which may vary at different depths, as well as downhole tool properties and the calculated formation slowness.
- the acoustic waveform processing method disclosed herein uses adaptive filtering based upon various models to more accurately estimate formation slowness.
- the foregoing figures illustrate two principles that are exploited by the method disclosed herein. First, when formation slowness is calculated at the minimum excitation frequency, it closely estimates the actual formation slowness. Second, when formation slowness is calculated at the peak excitation frequency or higher, the calculated formation slowness is significantly greater than the actual formation slowness.
- FIG. 3 a chart illustrating the steps for determining the polynomial coefficients used for an adaptive filter is shown.
- Data describing the types of downhole tools 302, borehole diameters 304, fluid properties (velocities, densities, or both) 306, and formation properties (velocities, densities, or both) 308 may be used as parameters for modal analysis 310.
- step 315 two sets of data plots are created. In the first plot, the minimum excitation frequencies are plotted with the corresponding formation slowness (the inverse of the formation's acoustic velocity) and a third parameter. In the second plot, the peak excitation frequencies are plotted with the corresponding formation slowness and the third parameter.
- fit appropriate polynomials are determined for each data plot.
- the polynomials may be fitted to the respective curves or surfaces using the method of least squares.
- the fit appropriate polynomial determined for the data plot using the minimum excitation frequencies is used to establish a lower bound of the adaptive filter.
- the fit appropriate polynomial determined for the data plot using the peak excitation frequencies is used to establish an upper bound of the adaptive filter.
- the coefficients of the polynomials are used to generate the adaptive filter band at each depth from the observed slowness.
- a default filter band 401 and basic input parameters 405 are supplied.
- a default filter band of 1 kHz to 6 kHz may be used.
- a default filter band of 750 Hz to 4 kHz may be used.
- a default filter band of 500 Hz to 2 kHz may be used.
- the basic input parameters 405 may include the slowness window of processing, tool types, borehole diameters, and fluid properties.
- the tool type parameter may indicate whether a wireline or logging-while-drilling tool is being used.
- Wireline tools have less of an effect on formation property measurements than logging-while-drilling tools because logging-while-drilling tools have a larger metal body presence in the borehole.
- the borehole diameters may be measured by a caliper tool.
- the fluid properties may be calculated (for example, the mud weight may be used to estimate the mud velocity), or may be measured by a sensor.
- the default filter band 401 is used to initialize the frequency filter 410.
- the basic input parameters 405 are used to initialize the point-to-point slowness-time coherence process 415 and the adaptive filter band estimator 420.
- streaming acoustic waveform data 425 acquired from an acoustic logging tool such as the acoustic logging tool 104 shown in figure 1A or the acoustic logging tool 114 shown in figure IB, is fed into the frequency filter 410. If a data point from the streaming acoustic waveform data 425 is above the lower bound frequency of the current adaptive frequency filter in use by the frequency filter 410 and below the upper bound of the current adaptive frequency filter in use by the frequency filter 410, then the data point is fed into the point-to-point slowness-time coherence process 415. The slowness-time coherence process 415 calculates the slowness from the data point.
- slowness peaks may be noted, and other measures to ensure quality control may be calculated.
- a semblance function is calculated in the slowness-time coherence process 415 .
- the semblance value ranges from 0, indicating no correlation, to 1, indicating perfect correlation.
- the semblance value may be used to determine whether to ignore a particular data point or not. For example, if the semblance value for a particular data point is below a certain threshold, such as 0.25, then that particular data point may not be used by the up-log/down-log slowness trend evaluator 440 because the data point is considered an outlier from the slowness trend.
- the resulting calculations may then be logged in the slowness logs 435.
- the estimated slowness is fed into the up-log/down-log slowness trend evaluator 440, which tracks whether the estimated slowness of the formation is trending up or trending down.
- the detected trend is then provided to the adaptive filter band estimator 420.
- the adaptive filter band estimator uses the supplied basic input parameters 405, polynomial coefficients 445 calculated from the modal analysis, and the estimated formation slowness trends from the up-log/down-log slowness trend evaluator 440 to adjust the frequency filter that is to be applied to the next set of data points received from the streaming acoustic waveform data 425.
- FIG. 5 a chart illustrating the steps of a real-time processing embodiment with adaptive waveform firing of the acoustic waveform signal processing method disclosed herein is shown.
- a default filter band 401 and basic input parameters 405 are supplied.
- the default filter band 401 is used to initialize the frequency filter 410.
- the basic input parameters 405 are used to initialize the point-to-point slowness-time coherence process 415 and the adaptive filter band estimator 420.
- streaming acoustic waveform data 425 acquired from an acoustic logging tool such as the acoustic logging tool 104 shown in figure 1A or the acoustic logging tool 114 shown in figure IB, is fed into the frequency filter 410. If a data point from the streaming acoustic waveform data 425 is above the lower bound of the current adaptive frequency filter in use by the frequency filter 410 and below the upper bound of the current adaptive frequency filter in use by the frequency filter 410, then the data point is fed into the point-to-point slowness- time coherence process 415. The slowness-time coherence process 415 calculates the slowness from the data point.
- slowness peaks may be noted, and other calculations, such as the semblance value, to ensure quality control may be calculated.
- the resulting calculations may then be logged in the slowness logs 435.
- the estimated slowness is fed into the up-log/down-log slowness trend evaluator 440, which tracks whether the estimated slowness of the formation is trending up or trending down.
- the detected trend is then provided to the adaptive filter band estimator 420.
- the adaptive filter band estimator uses the supplied basic input parameters 405, polynomial coefficients 445 calculated from the modal analysis, and the estimated formation slowness trends from the up-log/down-log slowness trend evaluator 440 to adjust the frequency filter that is to be applied to the next set of data points received from the streaming acoustic waveform data 425.
- the updated filter band from the adaptive filter band estimator 420 is provided to a selection frequency band module 505.
- the selection frequency band module 505 uses the updated filter band to select the frequencies at which the acoustic source will excite the formation.
- the selected frequency band is sent to the dynamic wavelet builder 510 that creates a digital wavelet containing the selected frequency band 505 and appropriate amplitude to excite the formation.
- the digital-to-analog converter (D/A) 515 transforms the digital wavelet received from the dynamic wavelet builder 510 into an analog signal for input into the linear driver 520.
- the linear driver 520 is a preamplifier used to amplify the wavelet so that there will be no distortion of the signal when it is used as the input to the power amplifier 525.
- the power amplifier 525 amplifies the received wavelet so that it can drive the acoustic transmitter 530.
- the acoustic transmitter 530 is used to excite the formation. For dipole mode generation, the acoustic transmitter 530 may be a PZT bender bar.
- FIG. 6 a chart illustrating the processing steps of a post-processing embodiment of the acoustic waveform processing method disclosed herein is shown.
- all of the acoustic waveform data and logs are retrieved from an acoustic logging tool, such as the acoustic logging tool 114 shown in figure IB, or may be retrieved from a computer processing system 108 which may have acquired the data live from the acoustic logging tool 104 shown in figure 1A.
- a default filter band 401 and basic input parameters 405 are supplied.
- the basic input parameters 405 may be retrieved from logs.
- the default filter band 401 is used to initialize the frequency filter 410.
- the basic input parameters 405 are used to initialize the point-to-point slowness-time coherence process 415 and the adaptive filter band estimator 420.
- stored acoustic waveform data 625 is fed into the frequency filter 410. If a data point from the stored acoustic waveform data 625 is above the lower bound of the current adaptive frequency filter in use by the frequency filter 410 and below the upper bound of the current adaptive frequency filter in use by the frequency filter 410, then the data point is fed into the point-to-point slowness-time coherence process 415.
- the slowness-time coherence process 415 calculates the slowness from the data point.
- slowness peaks may be noted, and other measures to ensure quality control, such as the semblance value, may be calculated.
- the resulting calculations may then be logged in slowness logs 435.
- the estimated slowness is fed into a module that tracks the average slowness within a depth range 640.
- the updated average slowness for the current depth range is then provided to the adaptive filter band estimator 420.
- the adaptive filter band estimator uses the supplied basic input parameters 405, polynomial coefficients 445 calculated from the modal analysis, and the average slowness within the current depth range, supplied from module 640, to adjust the frequency filter that is to be applied to the next set of data points received from the stored acoustic waveform data 625.
- Figures 7A-7H illustrate results from analyzing nearly 2,376 different models.
- the sample data included data acquired from boreholes with diameters ranging from five to sixteen inches, borehole fluid slowness ranging from 165 to 250 microseconds per foot, compressional slowness ranging from 40 to 210 microseconds per feet, and shear slowness ranging from 60 to 800 microseconds per foot.
- Figure 7A is a data plot of the minimum excitation frequencies and peak excitation frequencies observed in formations of varying slowness.
- Figure 7B is a data plot of the minimum excitation frequencies and peak excitation frequencies with the observed relative excitation energies.
- Figure 7C is a data plot showing the lower and upper limits of the observed minimum excitation frequencies, and the lower and upper limits of the observed peak excitation frequencies.
- Figure 7D is a data plot showing the observed shear slowness at the peak excitation energy levels.
- Figure 7E is a data plot showing the calculated formation compression slowness at the peak excitation frequencies.
- Figure 7F is a data plot showing the calculated formation shear slowness at peak excitation and the formation shear slowness.
- Figure 7G is a data plot showing the minimum and peak excitation frequencies and the observed fluid slowness.
- Figure 7H is a data plot showing the minimum and peak excitation frequencies and the observed borehole diameter
- Figures 8A illustrates the polynomials fitting the minimum excitation and peak excitation frequencies of various models of formation compressional slowness.
- the fitting polynomial for the lower bound of the adaptive frequency filter is based on the lowest minimum excitation frequencies from the model.
- the fitted polynomials were found using the least squares method.
- Figure 9A is a graph of the lower bound and the upper bound of an adaptive frequency filter for the compressional slowness.
- Figures 8B illustrates the polynomials fitting the minimum excitation and peak excitation frequencies of various models of formation shear slowness.
- the fitting polynomial for the lower bound is based on the lowest minimum excitation frequencies from the model.
- the fitted polynomials were found using the least squares method.
- Figure 9B is a graph of the lower bound and the upper bound of an adaptive frequency filter for the shear slowness.
- Figure 1 OA illustrates a fitting surface for the lower bound of an adaptive frequency filter.
- the minimum excitation frequencies are plotted with the corresponding formation shear slowness and the borehole diameter.
- Figure 10B illustrates a fitting surface for the upper bound of an adaptive frequency filter.
- the peak excitation frequencies are plotted with the corresponding formation shear slowness and the borehole diameter.
- a polynomial function can be found that fits each surface, and the resulting functions can be used for an adaptive frequency filter.
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Description
Claims
Priority Applications (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
MX2013010518A MX2013010518A (en) | 2011-03-15 | 2011-03-15 | Acoustic signal processing using model-based adaptive filtering. |
US14/004,242 US9354342B2 (en) | 2011-03-15 | 2011-03-15 | Acoustic signal processing using model-based adaptive filtering |
PCT/US2011/028516 WO2012125161A1 (en) | 2011-03-15 | 2011-03-15 | Acoustic signal processing using model-based adaptive filtering |
CA2829764A CA2829764A1 (en) | 2011-03-15 | 2011-03-15 | Acoustic signal processing using model-based adaptive filtering |
EP11712092.3A EP2686708A1 (en) | 2011-03-15 | 2011-03-15 | Acoustic signal processing using model-based adaptive filtering |
BR112013023475A BR112013023475A2 (en) | 2011-03-15 | 2011-03-15 | method for processing acoustic waveform data. |
AU2011362626A AU2011362626B2 (en) | 2011-03-15 | 2011-03-15 | Acoustic signal processing using model-based adaptive filtering |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/US2011/028516 WO2012125161A1 (en) | 2011-03-15 | 2011-03-15 | Acoustic signal processing using model-based adaptive filtering |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2012125161A1 true WO2012125161A1 (en) | 2012-09-20 |
Family
ID=44561133
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2011/028516 WO2012125161A1 (en) | 2011-03-15 | 2011-03-15 | Acoustic signal processing using model-based adaptive filtering |
Country Status (7)
Country | Link |
---|---|
US (1) | US9354342B2 (en) |
EP (1) | EP2686708A1 (en) |
AU (1) | AU2011362626B2 (en) |
BR (1) | BR112013023475A2 (en) |
CA (1) | CA2829764A1 (en) |
MX (1) | MX2013010518A (en) |
WO (1) | WO2012125161A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11550073B2 (en) | 2016-10-25 | 2023-01-10 | Halliburton Energy Services, Inc. | Enhanced-resolution rock formation body wave slowness determination from borehole guided waves |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2012396820A1 (en) * | 2012-12-11 | 2015-05-14 | Halliburton Energy Services, Inc. | Method and system for direct slowness determination of dispersive waves in a wellbore environment |
WO2017172805A1 (en) * | 2016-04-01 | 2017-10-05 | Halliburton Energy Services, Inc. | Automatic slowness-frequency range determination for advanced borehole sonic data processing |
DK3465286T3 (en) | 2016-05-25 | 2022-10-31 | Geoquest Systems Bv | Elastic parameter estimation |
US11435493B2 (en) * | 2017-06-15 | 2022-09-06 | Halliburton Energy Services, Inc. | Enhanced waveform analysis for target modes of borehole waves |
NO20210652A1 (en) * | 2018-12-28 | 2021-05-20 | Halliburton Energy Services Inc | Subsurface wave slowness prediction system |
CN112346038B (en) * | 2020-10-21 | 2022-07-19 | 自然资源部第三海洋研究所 | Sound channel response estimation method based on shallow sea explosion sound propagation signal |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040122595A1 (en) * | 2002-12-23 | 2004-06-24 | Henri-Pierre Valero | Methods and apparatus for processing acoustic waveforms received in a borehole |
WO2005116693A1 (en) * | 2004-05-17 | 2005-12-08 | Schlumberger Technology B.V. | Methods for processing dispersive acoustic waveforms |
US20090168595A1 (en) * | 2007-12-27 | 2009-07-02 | Schlumberger Technology Corporation | Method for gas zone detection using sonic wave attributes |
WO2010040087A2 (en) * | 2008-10-03 | 2010-04-08 | Baker Hughes Incorporated | A novel curve-fitting technique for determining dispersion characteristics of guided elastic waves |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5587966A (en) * | 1994-10-13 | 1996-12-24 | Schlumberger Technology Corporation | Sonic well logging methods and apparatus for processing flexural wave in optimal frequency band |
US6192316B1 (en) * | 1999-05-26 | 2001-02-20 | Atlantic Richfield Company | Fracture discrimination using borehole frequency response of stoneley waves |
US6510104B1 (en) * | 2000-06-07 | 2003-01-21 | Schlumberger Technology Corporation | Acoustic frequency selection in acoustic logging tools |
US6748329B2 (en) * | 2000-12-08 | 2004-06-08 | Halliburton Energy Services, Inc. | Acoustic signal processing method using array coherency |
US8964503B2 (en) * | 2009-04-28 | 2015-02-24 | Baker Hughes Incorporated | Petrophysics-guided processing of LWD acoustic data |
-
2011
- 2011-03-15 CA CA2829764A patent/CA2829764A1/en not_active Abandoned
- 2011-03-15 AU AU2011362626A patent/AU2011362626B2/en not_active Ceased
- 2011-03-15 MX MX2013010518A patent/MX2013010518A/en active IP Right Grant
- 2011-03-15 EP EP11712092.3A patent/EP2686708A1/en not_active Withdrawn
- 2011-03-15 US US14/004,242 patent/US9354342B2/en active Active
- 2011-03-15 BR BR112013023475A patent/BR112013023475A2/en not_active Application Discontinuation
- 2011-03-15 WO PCT/US2011/028516 patent/WO2012125161A1/en active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040122595A1 (en) * | 2002-12-23 | 2004-06-24 | Henri-Pierre Valero | Methods and apparatus for processing acoustic waveforms received in a borehole |
WO2005116693A1 (en) * | 2004-05-17 | 2005-12-08 | Schlumberger Technology B.V. | Methods for processing dispersive acoustic waveforms |
US20090168595A1 (en) * | 2007-12-27 | 2009-07-02 | Schlumberger Technology Corporation | Method for gas zone detection using sonic wave attributes |
WO2010040087A2 (en) * | 2008-10-03 | 2010-04-08 | Baker Hughes Incorporated | A novel curve-fitting technique for determining dispersion characteristics of guided elastic waves |
Non-Patent Citations (1)
Title |
---|
KIMBALL C V: "Shear slowness measurement by dispersive processing of the borehole flexural mode", GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS, US, vol. 63, no. 2, 1 March 1998 (1998-03-01), pages 337 - 344, XP002284422, ISSN: 0016-8033, DOI: 10.1190/1.1444333 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11550073B2 (en) | 2016-10-25 | 2023-01-10 | Halliburton Energy Services, Inc. | Enhanced-resolution rock formation body wave slowness determination from borehole guided waves |
Also Published As
Publication number | Publication date |
---|---|
US9354342B2 (en) | 2016-05-31 |
AU2011362626B2 (en) | 2014-08-21 |
AU2011362626A1 (en) | 2013-10-17 |
CA2829764A1 (en) | 2012-09-20 |
US20140005946A1 (en) | 2014-01-02 |
BR112013023475A2 (en) | 2016-12-06 |
MX2013010518A (en) | 2013-12-06 |
EP2686708A1 (en) | 2014-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9354342B2 (en) | Acoustic signal processing using model-based adaptive filtering | |
US6552962B1 (en) | Method and apparatus for acoustic logging | |
RU2369884C2 (en) | Facility and methods of meausurement of interval time for drilling agent in borehole | |
AU2005325123B2 (en) | Method and apparatus for estimating formation slowness | |
US9529105B2 (en) | Acoustic logging systems and methods employing multi-mode inversion for anisotropy and shear slowness | |
US10353090B2 (en) | Time-based processing of broadband borehole acoustic data | |
US20050254343A1 (en) | Methods for processing dispersive acoustic waveforms | |
US6510104B1 (en) | Acoustic frequency selection in acoustic logging tools | |
US7660195B2 (en) | Attenuation mapping apparatus, systems, and methods | |
US11360232B2 (en) | Mapping wave slowness using multi-mode semblance processing techniques | |
WO2009090465A2 (en) | Method for permeable zone detection | |
US7739049B2 (en) | Method and apparatus for multi-mode signal processing | |
US20050226098A1 (en) | Dynamic acoustic logging using a feedback loop | |
WO2020106287A1 (en) | Enhanced anisotropy analysis with multi-component dipole sonic data | |
WO2008142500A2 (en) | Methods and systems for processing acoustic waveform data | |
NO20210922A1 (en) | Waveform processing utilizing an amplitude adaptive data mask | |
WO2015153219A1 (en) | Formation measurements using flexural modes of guided waves | |
Market et al. | Processing and quality control of LWD dipole sonic measurements | |
US20180031722A1 (en) | Systems and methods employing a menu-based graphical user interface (gui) to derive a shear slowness log | |
EP3850188B1 (en) | Method for estimating a pore pressure value in geological formations to be drilled by a drilling apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 11712092 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2829764 Country of ref document: CA |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14004242 Country of ref document: US |
|
WWE | Wipo information: entry into national phase |
Ref document number: MX/A/2013/010518 Country of ref document: MX |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2011712092 Country of ref document: EP |
|
ENP | Entry into the national phase |
Ref document number: 2011362626 Country of ref document: AU Date of ref document: 20110315 Kind code of ref document: A |
|
REG | Reference to national code |
Ref country code: BR Ref legal event code: B01A Ref document number: 112013023475 Country of ref document: BR |
|
ENP | Entry into the national phase |
Ref document number: 112013023475 Country of ref document: BR Kind code of ref document: A2 Effective date: 20130912 |