EP3427087B1 - Zeitumgekehrte nichtlineare akustik für den bohrlochdruckmessungen - Google Patents

Zeitumgekehrte nichtlineare akustik für den bohrlochdruckmessungen Download PDF

Info

Publication number
EP3427087B1
EP3427087B1 EP17764114.9A EP17764114A EP3427087B1 EP 3427087 B1 EP3427087 B1 EP 3427087B1 EP 17764114 A EP17764114 A EP 17764114A EP 3427087 B1 EP3427087 B1 EP 3427087B1
Authority
EP
European Patent Office
Prior art keywords
signal
transceiver
time
acoustic signal
reversed
Prior art date
Legal status (The legal status 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 status listed.)
Active
Application number
EP17764114.9A
Other languages
English (en)
French (fr)
Other versions
EP3427087A1 (de
EP3427087A4 (de
Inventor
Paul A. Johnson
Timothy J. ULRICH II
Pierre-Yves Le Bas
Robert A. Guyer
Harvey E. Goodman
Marcel C. REMILLIEUX
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chevron USA Inc
Triad National Security LLC
Original Assignee
Chevron USA Inc
Triad National Security LLC
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 Chevron USA Inc, Triad National Security LLC filed Critical Chevron USA Inc
Publication of EP3427087A1 publication Critical patent/EP3427087A1/de
Publication of EP3427087A4 publication Critical patent/EP3427087A4/de
Application granted granted Critical
Publication of EP3427087B1 publication Critical patent/EP3427087B1/de
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/003—Seismic data acquisition in general, e.g. survey design
    • G01V1/005—Seismic data acquisition in general, e.g. survey design with exploration systems emitting special signals, e.g. frequency swept signals, pulse sequences or slip sweep arrangements
    • 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/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • 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
    • G01V2200/00—Details of seismic or acoustic prospecting or detecting in general
    • G01V2200/10—Miscellaneous details
    • G01V2200/16—Measure-while-drilling or logging-while-drilling
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/10—Aspects of acoustic signal generation or detection
    • G01V2210/12—Signal generation
    • G01V2210/121—Active source
    • G01V2210/1216—Drilling-related
    • 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
    • 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/624—Reservoir parameters
    • G01V2210/6248—Pore pressure
    • 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase

Definitions

  • Sedimentary rocks and man-made materials like concrete may be described as a network of mesoscopic-sized "hard” elements (e.g., grains with characteristic lengths ranging from tens to hundreds of microns) embedded in a "soft" bond system (e.g., cement between grains, pore space, fluid).
  • a "soft” bond system e.g., cement between grains, pore space, fluid.
  • Such systems belong to a wider class of materials referred to as Nonlinear Mesoscopic Elastic Materials (NMEMs).
  • NMEMs Nonlinear Mesoscopic Elastic Materials
  • the microscopic-sized imperfections at the interfaces between "hard” and “soft” subsystems are believed to be responsible for a number of interesting properties related to nonlinear and nonequilibrium dynamics, including the dependence of elastic parameters and attenuation on strain amplitude, slow dynamics, and hysteresis with end-point memory. Understanding and predicting these properties is basic for numerous applications including oil exploration.
  • drilling and/or extraction processes can be hampered and/or disrupted.
  • Pore pressures are the fluid pressures in the pore spaces in the porous formations.
  • the pore pressure gradient is used in drilling for determining mud weight, which is selected based on pore pressure gradient, wellbore stability and fracture gradient prior to setting and cementing a casing.
  • the drilling fluid is applied in the form of mud pressure to support the wellbore walls for preventing influx and wellbore collapse during drilling.
  • WO 2013/074112 A1 which relates to analyzing a subterranean formation.
  • a first signal is transmitted from a transmitter to the formation and a second signal which is a reflection of the first signal is received.
  • a third signal which is the second signal reversed in time, is then transmitted to the formation.
  • a fourth signal which is a reflection of the third signal from the formation is then received and monitored.
  • US 2012/075951 A1 which relates to a device and method for imaging non-linear and linear properties of formations surrounding a borehole.
  • WO 02/31538 A1 discloses determining pore pressure using a logging tool transmitting and receiving an acoustic signal. Pore pressure is estimated based on a frequency analysis of the received signal.
  • Benefits and advantages of embodiments of the present invention include, but are not limited to, providing an apparatus and method for determining the existence of and the distance to a down hole over-pressured region in advance of a drilling bit, using a combination of time reversal and elastic nonlinearity.
  • Overpressure rock has a signature elastic response that can be detected by combining Time Reversal techniques with Elastic Nonlinearity in a technique which is known as Time Reversal Nonlinear Elastic Wave Spectroscopy (TR NEWS).
  • TR NEWS Time Reversal Nonlinear Elastic Wave Spectroscopy
  • the nonlinear elastic wave response is directly related to the effective pressure (hydrostatic load minus the pore pressure).
  • Time reversal is a method for focusing acoustic waves such that large wave amplitudes are obtained in a localized region of space.
  • harmonics may be generated (and sum and difference frequencies if two waves are present). These harmonic frequencies are detected at the focus and, as will be discussed in more detail below, changes in the amplitude of the detected harmonics indicate that high pressure may be present.
  • Nonlinear materials exhibit a nonlinear stress-strain relation which can be probed by acoustic waves, leading to pressure-specific acoustic signatures. Harmonics of the incident acoustic frequencies are created when the acoustic waves are focused.
  • K is the linear stiffness constant
  • ⁇ is the strain
  • ⁇ ⁇ is the strain amplitude
  • ⁇ denotes the partial derivative with respect to time
  • sign is a function returning the sign (positive or negative) of the argument
  • ⁇ and ⁇ are combinations of third- and fourth-order elastic constants representing the acoustoelasticity (quadratic and cubic classical nonlinearity)
  • the parameter ⁇ relates to the strength of the hysteresis, according to the Preisach- Mayergoyz model of elasticity.
  • the parameters ⁇ , ⁇ , and ⁇ may be obtained from the time reversal signal, with ⁇ being obtained from the velocity change of the focused signal as a function of strain amplitude.
  • the velocity change may be also measured using cross correlation or another standard technique on a low amplitude (linear) wave at the time reversal focus, and the progressive delays caused by using progressively larger amplitude excitation waves.
  • Cross correlation is a commonly applied method for measuring time delays between a reference signal and a signal that has experienced a velocity change.
  • ⁇ is obtained from the amplitude dependence of the second harmonic of a pulsed pure sinusoid or the amplitude dependence of sum ( ⁇ 1 + ⁇ 2 ) and difference ( ⁇ 1 - ⁇ 2 ) frequencies if two waves are employed. See, also, TenCate, J.A. et al. (1996) "Laboratory Study Of Linear And Nonlinear Elastic Pulse Propagation In Sandstone," J. Acoust. Soc. Am. 100(3), 1383-1391 .
  • ⁇ is obtained from the amplitude dependence of the third harmonic of the fundamental drive amplitude at small, but still nonlinear amplitudes and, in general, can be ignored. At larger amplitudes, however, ⁇ dominates and ⁇ becomes overwhelmed and can be ignored.
  • the second derivative of u with respect to t is the particle acceleration measured in the frequency domain, f is the wave fundamental frequency, and ⁇ is the strain measured at frequency f in the focal region as the signal source amplitude is increased.
  • alpha can be obtained from the third harmonic amplitude also when wave amplitudes are large.
  • beta and delta are shown.
  • L is the wavelength of the fundamental frequency divided by two, equivalent to the radius of the focal region
  • the second derivative of u with respect to time, 3f is the third harmonic acceleration amplitude
  • the second derivative of u with respect to time, 2f is the second harmonic acceleration amplitude
  • the second derivative of u with respect to time, 1f is the fundamental harmonic acceleration amplitude
  • ⁇ 2 ⁇ f , where f is the fundamental frequency.
  • FIGURE 1A is a graph of ⁇ as a function of pore pressure
  • FIG. 1B is a graph of ⁇ as a function of pore pressure.
  • Time reversal permits the generation of focused, intense (non-damaging) sound in a region to induce local nonlinearities if high pressure is present, by taking advantage of the above relation for u 2f , thereby permitting detection and imaging of overpressure regions.
  • waves may be introduced into a specimen using a piezoelectric transducer. The waves are recorded on another transducer located elsewhere on the sample surface. The recorded waves are then reversed in time, and emitted from the detecting transducers, where they follow their forward wave paths backwards-in-space, and coalesce, focusing at the original source transducer, since the elastic wave equation is symmetric with respect to time. That is, the wave equation may be evaluated either forward or backward in time, the physics being identical.
  • Amplitudes at the time-reversed focus are large due to conservation of energy, since all of the energy contained in the long-duration scattered-signal is collapsed onto the focal point in space and time. Since wave amplitudes are largest at the focus, the local response may be nonlinear, but only at the focus.
  • FIG. 2A a schematic representation of an embodiment of the present apparatus, 10, for practicing look-ahead drilling based on TR NEWS is illustrated.
  • Other signal types may be employed, including pulsed signals.
  • transceiver mount 24 may be placed between just behind drilling bit 14 and 10 - 20 m therefrom.
  • Apparatus, 25, is employed to change the direction of drilling bit 14 in response to pore pressure measurements in accordance with embodiments of the present invention.
  • the detected signals are digitally reversed-in-time, by microprocessors located at detectors (transceivers) 22 or at the surface, amplified (not shown in FIGS. 2A and 2B ), and rebroadcast, 26, into the bore hole and into the surrounding formation.
  • the signals follow their forward propagation directions in reverse, and focus at source 12 that also acts as a detector.
  • the phase relationships among the returning waves permit the constructive interference thereof resulting in space and time focusing.
  • the focused signal is large in amplitude and is effective for inducing an elastic nonlinear response in focal volume, 28. If a portion of focal volume 28 encompasses a high fluid or gas pressured region, the nonlinear response of generated harmonic frequencies (and potentially sum and difference frequencies), and time delays due to wave speed decreases, will be significantly greater than at an established baseline thereof. This nonlinear response is detected and interpreted at the surface or by microprocessors located behind the drilling string (not shown in FIGS. 2A and 2B ).
  • the diameter of the focal spot measured at the half maximum value is equal to one-half of the dominant wavelength. See, e.g., " Depth Profile Of A Time-Reversal Focus In An Elastic Solid," by Marcel C. Remillieux et al., Ultrasonics 58 (2015) 60-66 . Corrective action can then be taken by the drillers (placing blowout stops etc.).
  • FIGS. 3A - 3C A schematic representation of an embodiment of transceiver mount 24 is shown in FIGS. 3A - 3C.
  • FIGURE 3A shows a possible arrangement for low-frequency transceivers, 22a
  • FIGS. 3B and 3C show an arrangement for mid-frequency, 22b, and high-frequency transceivers, 22c, respectively.
  • the transducer mount includes a long metallic portion, 30, adapted to fit in a cased borehole having an inner diameter of 6 in.
  • Transducer mount 24 may also be constructed of sturdy plastics capable of withstanding high temperatures and caustic environments.
  • the long dimension of mount 24 is equal to or larger than five times the largest wavelength of the elastic waves propagating in the formation.
  • the actual shape of the mount can be optimized to improve the transfer of energy from the tool to the formation.
  • FIGURE 3D is a schematic representation of another embodiment of transceiver mount 24 where transceiver 12 is located approximately in the center of the mount. This design is suitable for measuring pore pressure in formations through well casings.
  • Transducers 22a mid-frequency transceivers 22b and high-frequency transceivers 22c controlled by digital synthesizers, 32, 34, and 36, respectively, which are directed by microcontroller and digital signal processor, 38, are affixed along mount 24 to provide the required excitation signals.
  • Transducers vary in size and relative spacing depending on the center frequency of excitation signal that is intended to be generated. For low frequency excitation, large transducers are distributed over the entire length of the tool. For high frequency excitation, smaller transducers are centered with a smaller span around the point where focus should be achieved (at transceiver 12 ).
  • FIGURE 4A illustrates the forward or calibration steps to achieve the classical time reversed process using the apparatus of FIGS. 2A and 2B , and 3
  • FIG. 4B illustrates that reverse or focusing steps for the time reversed process.
  • a predetermined calibration signal, P from memory storage is directed to an arbitrary waveform generator, which drives, Step 42, source transducer 12 at the focal position, S. Signals generated from source transducer 12 are received by all of the transceivers 22, Step 44, and stored in memory, Step 46.
  • Step 48 in FIG. 4B simultaneously directs the time reversed signals stored in memory in step 46 to an arbitrary waveform generator which drives transceivers 22, Step 50, the transmitted signals being focused onto transducer 12 in Step 52.
  • Nonlinear signals, F, generated, Step 54, in the focal area are stored in memory for later processing.
  • This calibration process would be undertaken every time measurements of ⁇ and ⁇ are to be made, since external conditions, such as increasing pressures of drilling mud, may change the calibration.
  • Noise from impulsive elastic waves generated from the action of drilling bit 14 on the materials in a formation can be used as a source for the classical time reversal measurements in place of acoustic source 12 in accordance with embodiments of the present invention.
  • the drilling bit would be stopped when the amplified time-reversed signals generated by transceivers 22 are employed to generate harmonics in front of drilling bit 14, the harmonic signals being correlated with the time-reversed signals from the drilling bit.
  • transceivers 22 located on transceiver mount 24 are caused to individually broadcast a signal for a 'training' step.
  • Transceiver 12 detects these signals that are broadcast one at a time.
  • the detected signals are time reversed, and amplified and reemitted from transceivers 22. They again focus on the detector 12. This process works due to the reciprocity of the wave equation, since the transfer function in one direction is the same as that in the other direction.
  • FIGURE 5 illustrates the 'training' process.
  • a predefined calibration signal, P, in memory storage is directed, Step 56, to an arbitrary waveform generator, which drives transceivers 22, Step 58, sequentially, the generated signals being focused onto transceiver 12, Step 60, and the received signals stored in memory for time reversal processing, Step 62, after each stored signal is time reversed.
  • an arbitrary waveform generator which drives transceivers 22, Step 58, sequentially, the generated signals being focused onto transceiver 12, Step 60, and the received signals stored in memory for time reversal processing, Step 62, after each stored signal is time reversed.
  • FIGURE 6 is a schematic representation of an apparatus for carrying out reciprocal time reversal measurements described above in the laboratory for a simulated down hole environment.
  • Pipe, 100 was embedded in block, 102, of Berea Sandstone.
  • Ten acoustic transceivers 22 were affixed to portion, 104, of pipe 100 emerging from block 102.
  • a reference signal for example, a pulsed (10-20 kHz) sinusoidal waveform having a 50 kHz bandwidth, is directed into at least one arbitrary waveform generator, 106, by computer 38. After amplification by at least one power amplifier, 108, each generator signal is directed to a single transceiver 22, one generator signal at a time.
  • the signal traveling through pipe 100 and block 102 is recorded by laser vibrometer, 110, after being received by fiber optic attachment, 112, disposed inside pipe 100.
  • the signal received for each emission is digitized by digitizer, 114, and directed to computer 38, which time reverses each of the received signals and programs the arbitrary waveform generators 106 with the time-reversed signals.
  • Reciprocal time reversal in its basic form, includes first sending the last digitized element, then the second to last digitized element and so forth until the entire waveform has been inverted. Focusing in pipe 100 may be improved (typically, to sharpen the focus), by filtering the time-reversed signal (not shown in FIG. 6 ). Signal averaging may also be performed on the detected signals.
  • the final step is for all of the arbitrary waveform generator to simultaneously direct the time-reversed signals to the transducers through the power amplifiers. During that phase, all the generators are synchronized via pulser, 116. While all the transceivers are emitting, laser vibrometer 110 records the signals generated in the focal volume that are digitized and analyzed for nonlinear components by computer 38.
  • the time-reversed signals may be broadcasted successively at different amplitudes to assist in the detection of the nonlinear signals.
  • the size of the region probed by focused waves in the formation depends of the wavelength used for the first reference signal.
  • FIGURE 7A is a graph of the detected pulse propagation down pipe 100 from a single transducer 22 without time reversal
  • FIG. 7B is a graph of the detected pulse propagation from all 10 transducers 22 without time reversal
  • FIGURE 7C is a graph of the detected focused pulse propagation down pipe 100 from a single transducer 22 using reciprocal time reversal
  • FIG. 7D is a graph of the detected focused pulse propagation down pipe 100 from all 10 transducers using classical time reversal.
  • the signal strength increases by a factor of 10 when using reciprocal time reversal over that resulting from the use of conventional sources.
  • FIGURE 8A is a graph of the delay in the arrival time of a detected pulse in the pipe of the apparatus of FIG. 6 as the amplitude of the signal pulse increases from (a) to (c). As mentioned above, this shift is related to ⁇ .
  • FIGURE 8B is a graph of the detected pulse as a function of frequency as the amplitude of the signal pulse increases from (a) to (c).
  • the fundamental as well as the second and third harmonics are readily observable. As discussed above, ⁇ may be obtained by monitoring the second harmonic. Monitoring the third harmonic is an alternative method for obtaining ⁇ . When the bandwidth of the fundamental is broader, however, the third harmonic can partially overlap with the second harmonic, which may make third harmonic measurements more difficult. Additionally, the third harmonic may be influenced by ⁇ (see Equations 5 above). Due to the relative sizes of ⁇ and ⁇ , this latter issue is generally not a serious problem, but having an alternative way of measuring ⁇ , such as by using the time delay, is advantageous.
  • FIGURE 9 is a graph of the pore pressure as a function of depth measured by conventional techniques for an actual formation, permitting the prediction of gas/water contacts (GWC) for locating hydrocarbon deposits.
  • GWC gas/water contacts
  • drilling can be redirected using apparatus 25 in FIG. 2A to accomplish the change in direction.
  • ⁇ can be quantified by monitoring the propagation speed of an elastic wave as a function of the strain amplitude, laboratory experiments were performed. Although the propagation of impulsive elastic waves remains the principal measurement, time reversal is not required to generate the strain since the measurements are restricted to a one-dimensional waveguide over a known propagation distance. The hysteretic nonlinearity parameter has never been measured in this manner, so the determination is validated using nonlinear resonant ultrasound spectroscopy.
  • FIGURE 10 is a schematic representation of the apparatus used to determine ⁇ in sample, 120, of Berea sandstone (Cleveland Quarries, Amherst, Ohio) having a length of 1794 mm (70.63 in) and a diameter of 39.6 mm (1.56 in).
  • the sample was supported by a foam pad, not shown in FIG. 10 , in order to simulate free (unconstrained) boundary conditions.
  • Elastic waves were generated using a piezoelectric transducer, 122, epoxied onto one flat end of sample 120.
  • Impulsive elastic waveforms generated in signal generator, 124, and amplified by voltage amplifier, 126, were propagated in sample 120 at different amplitudes.
  • the vibrational response of bar 120 was recorded on the surface of the sample using 3D Laser Doppler vibrometer, 128, and received by data acquisition apparatus, 130.
  • a vibrational motion sensor such as a piezoelectric transducer, 132, on wall, 134 of wellbore, 18, in FIGS. 2A and 2B , corresponding to pipe 100 in FIG. 6
  • a vibrational motion sensor such as a piezoelectric transducer, 132, on wall, 134 of wellbore, 18, in FIGS. 2A and 2B , corresponding to pipe 100 in FIG. 6
  • the peak strain EQU. 8
  • ⁇ can be evaluated. With ⁇ being available from measurement of the intensity of harmonics, the pore pressure can be determined. Note that since the metal walls or casings of wellbores exhibit linear response to vibrations, the particle velocity, which is the motion excited in the focal region is accurately measured through the wall.
  • FIGURES 11A and 11B show the time histories of the axial component of the particle velocity measured at 500 mm from the source 122 for 20 source amplitudes ranging from 1 to 20 Vpp in steps of 1Vpp (before amplification), with
  • FIG. 11B illustrating an expanded abscissa between 0.32 ms and 0.37 ms.
  • the original sinusoidal waveform observed at the lowest amplitudes progressively evolves into a triangular wave, as a result of hysteretic nonlinearity, with additional distortion caused by classical nonlinearity. It is observed that the waveforms experience a significant delay in arrival time as the source amplitude increases. This delay is observed not only at the extrema (peaks of amplitude) but also at the zero crossings.
  • Classical nonlinearity would not induce any delay at the zero crossings (where the strain is equal to zero), because it produces an instantaneous variation of the modulus without time constants involved. Therefore, the observed delay is a direct consequence of hysteretic nonlinearity.
  • the data shown in FIGS. 11A and 11B can be further reduced to quantify hysteretic nonlinearity.
  • the signal measured at the source amplitude i-1 is taken as a reference to compute a relative time delay, ⁇ t i / i- 1 / t 0 .
  • the relative time delay between two signals is estimated by cross-correlation.
  • the relative time delay between the signals is also equal to the relative change in speed of the longitudinal wave, ⁇ c i /0 / c 0 .
  • the relative changes in the elastic modulus over the propagation path of the waveform can be followed as a function of the maximum strain amplitude at the measurement point.
  • the strain component of interest is ⁇ xx , where x is axial direction.
  • FIGS. 12A and 12B Reduced data from FIGS. 11A and 11B are shown in FIGS. 12A and 12B .
  • Three distinct regimes may be identified in the waveform: (i) the initial portion of the waveform, where the dynamic strain at the measurement point transitions from zero to steady state, and labeled with a "star"; (ii) the center portion of the waveform, where the absolute value of the peak amplitude of the dynamic strain is quasi-constant (that is, steady state), and labeled with a "0"; and (iii) the final portion of the waveform, where the dynamic strain at the measurement point transitions from steady state back to zero, that is, a path similar to the earlier portion of the waveform but in reverse, and also labeled with a "+".
  • the elastic modulus decreases progressively from its undisturbed equilibrium value to a smaller value (ultimately up to 5.3% smaller at the largest source amplitude employed as the steady state is approached). This is the portion typically referred to as conditioning.
  • the evolution of the material softening with the strain amplitude is essentially independent of the wave cycle selected for the analysis (that is, similar results are obtained from the 9 th to 21 st peaks of the waveforms).
  • the slope of the curve is 10480 ⁇ 600 (dimensionless). This value quantifies the hysteretic nonlinearity, but it will be determined in EXAMPLE 3 below how this value compares to the parameter ⁇ .
  • the apparatus may be utilized for resonant ultrasound spectroscopy.
  • Sample, 120 of Berea sandstone (Cleveland Quarries, Amherst, Ohio) having a length of 1794 mm (70.63 in) and a diameter of 39.6 mm (1.56 in) is shown.
  • the sample was supported by a foam pad, not shown in FIG. 10 .
  • Piezoelectric transducer 122 was driven with a sequence of harmonic voltage signals. Each harmonic signal was applied for 55 ms, and the transient vibrational response was recorded during the last 40 ms of the source signal, to ensure that steady state conditions had been reached. Vibrational spectra were constructed from the harmonic responses.
  • Frequencies ranging between 0.3 kHz and 7 kHz in steps of 2.5 Hz were employed, and at 22 excitation amplitudes ranging between 0.25 to 10 Vpp (before amplification).
  • the axial component of the acceleration was measured on the flat end opposite to the source by an accelerometer using transducer, 132, the output of which was processed by signal conditioner, 135, and analyzed using data acquisition apparatus 136.
  • the resonance frequencies are not computed for the first three modes because of the poor signal-to-noise ratio.
  • any resonance mode can be selected to quantify hysteretic nonlinearity as long as the mode type is purely longitudinal.
  • the vibrational spectra for this experiment are shown in FIGS. 13A and 13B .
  • Material softening a drop in the effective elastic moduli resulting from the fact that elastic moduli are dependent on the strain amplitude and strain rate, is observed when the drive amplitude of the source becomes sufficiently large (Equ. 3).
  • such softening can be quantified by plotting the relative frequency shift as a function of the maximum strain in the sample, the slope of which is the nonlinear parameter ⁇ .
  • the maximum strain in the sample may be inferred analytically from the measured vibrational response.
  • the strain component of interest is ⁇ xx .
  • Eq. (8) can also be used in the context of a resonance experiment to relate the maximum amplitude of the axial component of the particle velocity at the free end of the sample (where data is acquired) to the maximum amplitude of the axial component of the strain in the sample.
  • the relative frequency shift varies almost linearly with the maximum strain beyond 4 microstrains.
  • the Young's modulus is approximately constant below 4 microstrains in the conditioning phase and varies linearly with strain above this value, with a sharp transition between the two regimes (see FIGS. 11A and 11B ).
  • the resonance frequency also varies linearly with strain beyond 4 microstrains but experiences a smooth transition to this regime.
  • classical nonlinearity plays a substantial role.
  • Classical nonlinearity can induce a frequency shift in the resonance experiment, but cannot induce a time delay in the pulse propagation experiment; hence, the sharp transitions in FIGS. 11A and 11B . In this sense, it is possible to decouple the contributions from classical nonlinearity and nonequilibrium dynamics in the pulse propagation experiment.
  • TR NEWS time reversal and elastic nonlinearity

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Force Measurement Appropriate To Specific Purposes (AREA)
  • Measuring Fluid Pressure (AREA)

Claims (11)

  1. Vorrichtung zum Messen des Porendrucks in einer Formation, umfassend:
    einen ersten Signalgenerator, der konfiguriert ist, um ein gepulstes Sinussignal zu erzeugen, wobei das gepulste Sinussignal eine Frequenz aufweist;
    einen ersten Transceiver, der in einem Bohrloch in der Formation angeordnet ist und konfiguriert ist, um das gepulste sinusförmige Signal zu empfangen und ein gepulstes sinusförmiges akustisches Signal gemäß dem empfangenen gepulsten sinusförmigen Signal zu senden;
    einen zweiten Transceiver, der im Bohrloch angeordnet ist und konfiguriert ist, um das gepulste sinusförmige akustische Signal zu empfangen und daraus ein erstes elektrisches Signal zu erzeugen;
    einen oder mehrere Mikroprozessoren, die konfiguriert sind, um das erste elektrische Signal zu empfangen und das erste elektrische Signal zeitlich umzukehren, um ein zeitumgekehrtes elektrisches Signal zu erzeugen;
    einen zweiten Signalgenerator, der dazu konfiguriert ist, das zeitlich umgekehrte elektrische Signal zu empfangen und das zeitlich umgekehrte elektrische Signal an den zweiten Transceiver zu leiten, wobei der zweite Transceiver ferner dazu konfiguriert ist, ein zweites akustisches Signal gemäß dem zeitlich umgekehrten elektrischen Signal zu senden, wobei das zweite akustische Signal ein Fokusvolumen bildet, das auf dem ersten Transceiver zentriert ist, wobei das zweite akustische Signal zweite und dritte Harmonische der Frequenz enthält, wobei der erste Transceiver ferner dazu konfiguriert ist, ein zweites elektrisches Signal basierend auf dem Empfang des zweiten akustischen Signals zu erzeugen, wobei das zweite elektrische Signal eine Amplitude aufweist; und
    eine Speicherablage, die dazu konfiguriert ist, die Amplitude des zweiten elektrischen Signals und/oder von der Amplitude des zweiten elektrischen Signals abgeleitete Informationen zu speichern,
    wobei ein Hysterese-Nichtlinearitätsparameter α in Bezug auf die Hysteresestärke bestimmt wird, basierend auf der Amplitude der dritten Harmonischen des zweiten elektrischen Signals,
    wobei ein nichtlinearer elastischer Parameter β bezüglich der Akustoelastizität bestimmt wird, basierend auf der Amplitude der zweiten Harmonischen des zweiten elektrischen Signals und der Porendruck in der Formation bestimmt wird, basierend auf dem hysteretischen Nichtlinearitätsparameter α und dem nichtlinearen elastischen Parameter β.
  2. Vorrichtung gemäß Anspruch 1, wobei der erste Transceiver im Bohrloch über einem Bohrmeißel angeordnet ist.
  3. Vorrichtung gemäß Anspruch 1, wobei die Frequenz variiert wird, wodurch die Größe des Fokus des zweiten akustischen Signals variiert wird.
  4. Verfahren zum Messen des Porendrucks in einer Formation, wobei das Verfahren von einem System durchgeführt wird, das umfasst einen in einem Bohrloch in der Formation angeordneten ersten Transceiver, einen in dem Bohrloch angeordneten zweiten Transceiver, einen oder mehrere Mikroprozessoren und eine Speicherablage, das Verfahren umfassend
    Erzeugen, mit dem ersten Transceiver (12), eines gepulsten sinusförmigen akustischen Signals, wobei das gepulste sinusförmige akustische Signal eine Frequenz hat;
    Empfangen, mit dem zweiten Transceiver (22), des gepulsten sinusförmigen akustischen Signals;
    Erzeugen, mit dem einen oder den mehreren Mikroprozessoren (38). eines zeitumgekehrten Signals durch Zeitumkehren des empfangenen Signals;
    Senden, mit dem zweiten Transceiver (22), eines zeitlich umgekehrten akustischen Signals entsprechend dem zeitlich umgekehrten Signal, wobei das zeitlich umgekehrte akustische Signal ein auf dem ersten Transceiver (12) zentriertes Fokusvolumen bildet;
    Empfangen, mit dem ersten Transceiver (12), von akustischen Signalen innerhalb des Fokusvolumens, wobei solche akustische Signale zweite und dritte harmonische Signale der Frequenz des gepulsten sinusförmigen akustischen Signals umfassen, wobei die zweiten und dritten harmonischen Signale im Fokusvolumen auf der ersten Transceiver (12) erzeugt werden, wobei die zweiten und dritten harmonischen Signale Amplituden aufweisen; und
    Speichern, in der Speicherablage, der Amplituden der zweiten und dritten harmonischen Signale,
    wobei ein Hysterese-Nichtlinearitätsparameter α in Bezug auf die Hysteresestärke bestimmt wird, basierend auf der Amplitude des dritten harmonischen Signals,
    wobei ein nichtlinearer elastischer Parameter β bezüglich der Akustoelastizität bestimmt wird, basierend auf der Amplitude des zweiten harmonischen Signals, und der Porendruck innerhalb der Formation bestimmt wird, basierend auf dem hysteretischen Nichtlinearitätsparameter α und dem nichtlinearen elastischen Parameter β.
  5. Verfahren gemäß Anspruch 4, ferner umfassend den Schritt des Variierens der Frequenz des gepulsten sinusförmigen akustischen Signals, wodurch die Größe des Fokus des zeitlich umgekehrten Signals variiert wird.
  6. Verfahren zum Messen des Porendrucks in einer Formation durch ein Bohrloch mit einer Metallverrohrung, wobei das Verfahren von einem System durchgeführt wird, das einen im Bohrloch angeordneten ersten Sendeempfänger (12) umfasst, einen im Bohrloch über dem ersten Sendeempfänger (12) angeordneten zweiten Sendeempfänger (22), einen Wandler (132), der auf dem Metallgehäuse angeordnet ist, einen ersten Prozessor und einen zweiten Prozessor, das Verfahren umfassend
    Erzeugen, mit dem ersten Transceiver (12), eines gepulsten sinusförmigen akustischen Signals mit einer Frequenz;
    Empfangen, mit dem zweiten Transceiver (22), des gepulsten sinusförmigen akustischen Signals;
    Erzeugen, mit dem ersten Prozessor, eines zeitumgekehrten Signals durch Zeitumkehren des vom zweiten Transceiver (22) empfangenen gepulsten sinusförmigen akustischen Signals;
    Senden, mit einer ausgewählten Intensität mit dem zweiten Transceiver (22), eines zeitlich umgekehrten akustischen Signals in Übereinstimmung mit dem zeitlich umgekehrten Signal, wobei das zeitlich umgekehrte akustische Signal ein Fokusvolumen bildet, das auf dem ersten Transceiver zentriert ist;
    Empfangen, mit dem ersten Transceiver (12), von Energie bei einer zweiten Harmonischen der Frequenz des gepulsten sinusförmigen akustischen Signals, wobei die Energie bei der zweiten Harmonischen eine Amplitude hat;
    Bestimmen eines nichtlinearen elastischen Parameters β bezüglich der Akustoelastizität mit dem zweiten Prozessor, basierend auf der Amplitude der vom ersten Transceiver (12) bei der zweiten Harmonischen empfangenen Energie;
    Empfangen, mit dem Wandler (132), eines zweiten akustischen Signals, das auf eine Schwingungsanregung im Fokusvolumen anspricht; und
    Variieren, mit dem zweiten Transceiver (22), der Intensität des gesendeten zeitumgekehrten akustischen Signals,
    wobei ein Hysterese-Nichtlinearitätsparameter α bezüglich der Hysteresestärke bestimmt wird durch Messen einer Zeitverzögerung des zweiten akustischen Signals relativ zum gesendeten zeitumgekehrten akustischen Signal als Funktion der Intensität des gesendeten zeitumgekehrten akustischen Signals, und der Porendruck in der Formation bestimmt wird, basierend auf dem hysteretischen Nichtlinearitätsparameter α und dem nichtlinearen elastischen Parameter β.
  7. Verfahren gemäß Anspruch 6, ferner umfassend den Schritt des Verstärkens des übertragenen zeitumgekehrten Signals.
  8. Verfahren gemäß Anspruch 6 oder 7, ferner umfassend den Schritt des Variierens der Frequenz des gepulsten sinusförmigen akustischen Signals, wodurch die Größe des Fokus des zeitlich umgekehrten Signals variiert wird.
  9. Vorrichtung zum Messen des Porendrucks in einer Formation durch ein Bohrloch mit einem Metallgehäuse, umfassend
    einen ersten Signalgenerator, der konfiguriert ist, um ein gepulstes Sinussignal mit einer Frequenz bereitzustellen;
    einen ersten im Bohrloch in der Formation angeordneten Transceiver (12), der konfiguriert ist, um das gepulste sinusförmige Signal vom ersten Signalgenerator zu empfangen und ein gepulstes sinusförmiges akustisches Signal gemäß dem gepulsten sinusförmigen Signal vom ersten Signalgenerator zu senden;
    einen zweiten im Bohrloch angeordneten Transceiver (22), der konfiguriert ist, um das gepulste sinusförmige akustische Signal zu empfangen und daraus ein erstes elektrisches Signal zu erzeugen;
    einen ersten Prozessor, der dazu konfiguriert ist, das erste elektrische Signal zu empfangen und ein zeitlich umgekehrtes elektrisches Signal durch zeitliches Umkehren des ersten elektrischen Signals zu erzeugen;
    einen zweiten Signalgenerator, der konfiguriert ist, um das zeitlich umgekehrte elektrische Signal zu empfangen, um daraus ein zweites elektrisches Signal mit einer Intensität zu erzeugen, die gesteuert wird, um zu variieren, und um das zweite elektrische Signal an den zweiten Transceiver (22) zu leiten, so dass der zweite Transceiver (22) ein zweites akustisches Signal sendet mit einer Intensität, die entsprechend der Intensität des zweiten elektrischen Signals variiert;
    wobei das zweite akustische Signal ein auf dem ersten Transceiver (12) zentriertes Fokalvolumen bildet, wobei der erste Transceiver (12) ein zweites harmonisches Signal bei der zweiten Harmonischen der Frequenz des gepulsten Sinussignals empfängt;
    wobei ein zweiter Prozessor konfiguriert ist, um einen nichtlinearen elastischen Parameter β bezüglich der Akustoelastizität zu bestimmen, basierend auf einer Amplitude des empfangenen zweiten harmonischen Signals; und
    einen Wandler (132), der im Bohrloch auf dem Metallgehäuse angeordnet ist, wobei der Wandler (132) dazu konfiguriert ist, ein drittes akustisches Signal als Reaktion auf eine Schwingungsanregung im Fokusvolumen zu empfangen,
    wobei ein Hysterese-Nichtlinearitätsparameter α bezüglich der Hysteresestärke bestimmt wird, basierend auf der Zeitverzögerung des dritten akustischen Signals relativ zum zweiten akustischen Signal als Funktion der Intensität des zweiten akustischen Signals, und der Porendruck in der Formation bestimmt wird basierend auf dem hysteretischen Nichtlinearitätsparameter α und dem nichtlinearen elastischen Parameter β.
  10. Vorrichtung gemäß Anspruch 9, ferner umfassend mindestens einen Verstärker, der konfiguriert ist, um das zweite akustische Signal zu verstärken.
  11. Vorrichtung gemäß Anspruch 9 oder 10, wobei die Frequenz so variiert wird, dass die Größe des Fokusvolumens des zweiten akustischen Signals variiert wird.
EP17764114.9A 2016-03-09 2017-03-09 Zeitumgekehrte nichtlineare akustik für den bohrlochdruckmessungen Active EP3427087B1 (de)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201662306037P 2016-03-09 2016-03-09
US201662367337P 2016-07-27 2016-07-27
PCT/US2017/021606 WO2017156292A1 (en) 2016-03-09 2017-03-09 Time-reversed nonlinear acoustics for downhole pressure measurements

Publications (3)

Publication Number Publication Date
EP3427087A1 EP3427087A1 (de) 2019-01-16
EP3427087A4 EP3427087A4 (de) 2019-10-23
EP3427087B1 true EP3427087B1 (de) 2021-12-29

Family

ID=59789889

Family Applications (1)

Application Number Title Priority Date Filing Date
EP17764114.9A Active EP3427087B1 (de) 2016-03-09 2017-03-09 Zeitumgekehrte nichtlineare akustik für den bohrlochdruckmessungen

Country Status (5)

Country Link
US (1) US11079505B2 (de)
EP (1) EP3427087B1 (de)
BR (1) BR112018068254B8 (de)
CA (1) CA3017196C (de)
WO (1) WO2017156292A1 (de)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3427087B1 (de) 2016-03-09 2021-12-29 Triad National Security, LLC Zeitumgekehrte nichtlineare akustik für den bohrlochdruckmessungen
US11086040B2 (en) 2016-03-09 2021-08-10 Triad National Security, Llc Time-reversed nonlinear acoustics for wellbore integrity characterization
CA3079020C (en) * 2017-10-13 2022-10-25 Exxonmobil Upstream Research Company Method and system for performing communications using aliasing
EP3850188B1 (de) * 2018-09-14 2022-08-17 ENI S.p.A. Verfahren zum schätzen eines porendruckwertes in geologischen formationen, die mit einer bohrvorrichtung durchbohrt werden sollen
CN110927790A (zh) * 2018-09-20 2020-03-27 中国石油化工股份有限公司 基于叠前弹性参数的反演地层压力预测方法及系统
US11650346B2 (en) * 2019-08-15 2023-05-16 Halliburton Energy Services, Inc. Downhole acoustic measurement
US11366049B2 (en) * 2020-07-23 2022-06-21 Baker Hughes Oilfield Operations Llc Estimation of objective driven porous material mechanical properties
US11808136B2 (en) 2021-05-27 2023-11-07 Halliburton Energy Services, Inc. Through-tubing, cased-hole sealed material evaluation using acoustic measurements
CN114894903B (zh) * 2022-05-23 2023-09-29 中国石油大学(华东) 一种各向同性材料三阶弹性常数测定系统及方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120075951A1 (en) * 2010-09-29 2012-03-29 Chevron U.S.A., Inc. Device and method for imaging of non-linear and linear properties of formations surrounding a borehole

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6269198B1 (en) 1999-10-29 2001-07-31 Litton Systems, Inc. Acoustic sensing system for downhole seismic applications utilizing an array of fiber optic sensors
EP1410072A4 (de) 2000-10-10 2005-08-31 Exxonmobil Upstream Res Co Verfahren zur bohrlochmessung von formationseigenschaften
US20020112540A1 (en) 2000-12-20 2002-08-22 Schlumberger Technology Corporation Acoustic method for estimating mechanical properties of a material and apparatus therefor
WO2004095077A1 (en) * 2003-04-23 2004-11-04 Commonwealth Scientific And Industrial Research Organisation Method for predicting pore pressure
US7735579B2 (en) 2005-09-12 2010-06-15 Teledrift, Inc. Measurement while drilling apparatus and method of using the same
US7639563B2 (en) 2007-12-27 2009-12-29 Schlumberger Technology Corporation Method for sonic indication of voids in casing cement
US7753118B2 (en) 2008-04-04 2010-07-13 Schlumberger Technology Corporation Method and tool for evaluating fluid dynamic properties of a cement annulus surrounding a casing
CA2772256A1 (en) * 2009-12-15 2011-07-14 Spectraseis Ag Energy density and stress imaging conditions for source localization and characterization
US9678236B2 (en) * 2010-04-27 2017-06-13 Halliburton Energy Services, Inc. Fracture characterization by interferometric drillbit imaging, time reversal imaging of fractures using drill bit seismics, and monitoring of fracture generation via time reversed acoustics and electroseismics
US9103196B2 (en) 2010-08-03 2015-08-11 Baker Hughes Incorporated Pipelined pulse-echo scheme for an acoustic image tool for use downhole
US8553495B2 (en) 2010-09-29 2013-10-08 Chevron U.S.A. Inc. Device and method for imaging of non-linear and linear properties of formations surrounding a borehole
US9046620B2 (en) * 2010-11-12 2015-06-02 Los Alamos National Security Llc System and method for investigating sub-surface features and 3D imaging of non-linear property, compressional velocity VP, shear velocity VS and velocity ratio VP/VS of a rock formation
FR2977671B1 (fr) * 2011-07-08 2013-11-29 01Db Metravib Procede et dispositif pour controler des structures par retournement temporel
WO2013074112A1 (en) 2011-11-18 2013-05-23 Halliburton Energy Services, Inc. Methods and systems for analyzing formation properties when performing subterranean operations
EP3126620A4 (de) * 2014-04-04 2017-11-29 Services Pétroliers Schlumberger Bohrlochzementsimulator
US10012749B2 (en) 2014-10-20 2018-07-03 Schlumberger Technology Corporation Fast model based inversion of acoustic impedance of annulus behind casing
US10605954B2 (en) 2015-12-29 2020-03-31 Halliburton Energy Services, Inc. Simultaneous evaluation of the volume and the position of voids in downhole cement
EP3427087B1 (de) 2016-03-09 2021-12-29 Triad National Security, LLC Zeitumgekehrte nichtlineare akustik für den bohrlochdruckmessungen
US10408053B2 (en) 2016-05-19 2019-09-10 Baker Hughes, A Ge Company, Llc Encapsulated phased array segment for downhole applications

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120075951A1 (en) * 2010-09-29 2012-03-29 Chevron U.S.A., Inc. Device and method for imaging of non-linear and linear properties of formations surrounding a borehole

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
SUTIN ALEXANDER M ET AL: "Single-channel time reversal in elastic solids", THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA, AMERICAN INSTITUTE OF PHYSICS FOR THE ACOUSTICAL SOCIETY OF AMERICA, NEW YORK, NY, US, vol. 116, no. 5, 1 November 2004 (2004-11-01), pages 2779 - 2784, XP012072603, ISSN: 0001-4966, DOI: 10.1121/1.1802676 *

Also Published As

Publication number Publication date
BR112018068254B1 (pt) 2023-06-13
CA3017196A1 (en) 2017-09-14
EP3427087A1 (de) 2019-01-16
US11079505B2 (en) 2021-08-03
EP3427087A4 (de) 2019-10-23
CA3017196C (en) 2020-08-04
BR112018068254A2 (pt) 2019-01-15
WO2017156292A1 (en) 2017-09-14
US20190101660A1 (en) 2019-04-04
BR112018068254B8 (pt) 2023-10-17

Similar Documents

Publication Publication Date Title
EP3427087B1 (de) Zeitumgekehrte nichtlineare akustik für den bohrlochdruckmessungen
US11086040B2 (en) Time-reversed nonlinear acoustics for wellbore integrity characterization
US11029435B2 (en) Downhole nonlinear acoustics measurements in rock formations using dynamic acoustic elasticity and time reversal
US11067711B2 (en) Time-reversed nonlinear acoustic downhole pore pressure measurements
RU2684267C1 (ru) Геонавигация при бурении скважин с использованием распределенного акустического зондирования
US10094945B2 (en) Formation measurements using nonlinear guided waves
CN103140773A (zh) 用于对钻孔周围地层的非线性和线性特性进行成像的装置和方法
NO20161639A1 (en) Formation measurements using flexural modes of guided waves
Marsala et al. Sonic while drilling: Have you thought about cuttings?
EP1438573B1 (de) Verwendung von bohrlöchern für echtzeit-dämpfungsvorhersagen
Sirevaag A Study of the Ultrasonic Measurements for Logging Behind a Steel Pipe: Expanding the proessing to improve the differentiation between heavy fluids and light solids
Wang et al. Identification of the seismoelectric field induced by a slow compressional wave at an ultrasonic frequency
Yang et al. Extension of tube wave detection for quality evaluation of pile foundation
Potylitsyn et al. The Results of Observing the Seismoelectric Effect in a Gas Condensate Hydrocarbon Field Using a Source of Seismic Field
Burns Predicting relative and absolute variations of in-situ permeability from full-waveform acoustic logs
BR112019008243B1 (pt) Método e aparelho configurado para determinar parâmetros não lineares e pressão de poro em uma formação através de um furo de poço
Sakiyama et al. Classifying and removing unwanted arrivals propagating through a drill collar for semblance processing
Cheng et al. Acoustic-waveform Logging-advances Intheory And Application
AU2003254730A1 (en) Estimating formation stress using borehole monopole and cross-dipole acoustic measurements: theory and method

Legal Events

Date Code Title Description
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE

PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20181008

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

DAV Request for validation of the european patent (deleted)
DAX Request for extension of the european patent (deleted)
REG Reference to a national code

Ref country code: DE

Ref legal event code: R079

Ref document number: 602017051582

Country of ref document: DE

Free format text: PREVIOUS MAIN CLASS: G01V0001400000

Ipc: G01V0001000000

A4 Supplementary search report drawn up and despatched

Effective date: 20190920

RIC1 Information provided on ipc code assigned before grant

Ipc: G01V 1/30 20060101ALI20190916BHEP

Ipc: G01V 1/48 20060101ALI20190916BHEP

Ipc: G01V 1/00 20060101AFI20190916BHEP

RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: CHEVRON U.S.A. INC.

Owner name: TRIAD NATIONAL SECURITY, LLC

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

17Q First examination report despatched

Effective date: 20200828

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

INTG Intention to grant announced

Effective date: 20210809

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE PATENT HAS BEEN GRANTED

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: AT

Ref legal event code: REF

Ref document number: 1459098

Country of ref document: AT

Kind code of ref document: T

Effective date: 20220115

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602017051582

Country of ref document: DE

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG9D

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: RS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220329

REG Reference to a national code

Ref country code: NL

Ref legal event code: MP

Effective date: 20211229

REG Reference to a national code

Ref country code: AT

Ref legal event code: MK05

Ref document number: 1459098

Country of ref document: AT

Kind code of ref document: T

Effective date: 20211229

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: LV

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: HR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220330

REG Reference to a national code

Ref country code: NO

Ref legal event code: T2

Effective date: 20211229

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: NL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SM

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220429

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: PL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220429

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602017051582

Country of ref document: DE

REG Reference to a national code

Ref country code: DE

Ref legal event code: R119

Ref document number: 602017051582

Country of ref document: DE

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MC

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: DK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: AL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed

Effective date: 20220930

REG Reference to a national code

Ref country code: BE

Ref legal event code: MM

Effective date: 20220331

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20220309

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20220331

Ref country code: IE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20220309

Ref country code: DE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20221001

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20220331

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20220331

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO

Effective date: 20170309

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20240201

Year of fee payment: 8

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: NO

Payment date: 20240222

Year of fee payment: 8

Ref country code: FR

Payment date: 20240213

Year of fee payment: 8

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229