Formation property determination apparatus, methods, and systems

Patent metadata

Field	Value
Publication number	US-9513396-B2
Application number	US-201213878727-A
Country	US
Kind code	B2
Filing date	Aug 28, 2012
Priority date	Sep 12, 2011
Publication date	Dec 6, 2016
Grant date	Dec 6, 2016

Abstract

In some embodiments, apparatus and systems, as well as methods, may operate to record a plurality of acoustic waveforms comprising cross-dipole waveforms and corresponding to acoustic waves received 5 at azimuthally orthogonal dipole receiver arrays surrounded by a geological formation, the waves being generated by azimuthally orthogonal dipole transmitter arrays. Further activity may include defining an objective function dependent on the cross-dipole waveforms, eigenvalues, an azimuth angle corresponding to an 10 orientation of the transmitter and receiver arrays; minimizing the objective function with respect to said angle and the set of auxiliary parameters; and determining at least one property of the geological formation based on the global minimum. Additional apparatus, systems, and methods are disclosed.

First claim

What is claimed is:

1. An apparatus, comprising:
azimuthally orthogonal dipole receiver arrays and azimuthally orthogonal transmitter arrays;
logic to record a plurality of acoustic waveforms corresponding to acoustic waves received at the azimuthally orthogonal dipole receiver arrays, the acoustic waves being generated by the azimuthally orthogonal transmitter arrays; and
a signal processor to estimate a global minimum of an objective function with respect to an azimuth angle and a set of auxiliary parameters, wherein the azimuth angle corresponds to an orientation of the transmitter arrays and the receiver arrays relative to fast and slow principal flexural wave axes, with the signal processor configured to estimate fast and slow principal flexural waves propagating along the axes and characterized by the waveforms as ordered complex eigenvalues of a waveform matrix diagonalized by a similarity transformation at each data point of the waveforms, and the objective function depending on the waveforms, the ordered eigenvalues, the azimuth angle, and the set of auxiliary parameters, wherein the signal processor is also configured to remove existing ambiguities associated with the fast and slow principal flexural wave axes.

2. The apparatus of claim 1 , further comprising:
a memory to receive and store values associated with the global minimum of the objective function, to include at least the azimuth angle.

3. The apparatus of claim 1 , further comprising:
a telemetry transmitter to communicate values associated with the global minimum of the objective function to a surface logging facility.

4. The apparatus of claim 1 , wherein the similarity transformation is a unitary similarity transformation.

5. A system, comprising:
a down hole tool;
azimuthally orthogonal dipole receiver arrays and azimuthally orthogonal transmitter arrays, each of the arrays attached to the down hole tool;
logic to record a plurality of acoustic waveforms corresponding to acoustic waves received at the azimuthally orthogonal dipole receiver arrays, the acoustic waves being generated by the azimuthally orthogonal transmitter arrays; and
a signal processor to estimate a global minimum of an objective function with respect to an azimuth angle and a set of auxiliary parameters, wherein the azimuth angle corresponds to an orientation of the transmitter arrays and the receiver arrays relative to fast and slow principal flexural wave axes, with the signal processor configured to estimate fast and slow principal flexural waves propagating along the axes and characterized by the waveforms as ordered complex eigenvalues of a waveform matrix diagonalized by a similarity transformation at each data point of the waveforms, and the objective function depending on the waveforms, the ordered eigenvalues, the azimuth angle, and the set of auxiliary parameters, wherein the signal processor is also configured to remove existing ambiguities associated with the fast and slow principal flexural wave axes.

6. The system of claim 5 , wherein the downhole tool comprises one of a wireline tool or a measurement while drilling tool.

7. The system of claim 5 , further comprising:
a surface computer comprising the signal processor.

8. The system of claim 5 , wherein the similarity transformation is a unitary similarity transformation.

9. A processor-implemented method of estimating an azimuth angle, to execute on one or more processors that perform the method, comprising:
recording a plurality of acoustic waveforms corresponding to acoustic waves received at azimuthally orthogonal dipole receiver arrays surrounded by a geological formation, the waves being generated by azimuthally orthogonal transmitter arrays, the azimuth angle corresponding to an orientation of the receiver arrays and the transmitter arrays relative to fast and slow principal flexural wave axes, wherein fast and slow principal flexural waves propagating along the axes and characterized by the waveforms being estimated as the ordered eigenvalues of a waveform matrix diagonalized by a similarity transformation;
defining an objective function dependent on the waveforms, the ordered eigenvalues, the azimuth angle, and a set of auxiliary parameters;
minimizing the objective function with respect to the azimuth angle; and
determining at least one anisotropy angle of the geological formation based on the global minimum.

10. The method of claim 9 , wherein the matrix represents a set of data points in a space defined by receiver ring index position and time.

11. The method of claim 10 , wherein determining at least one anisotropy angle of the geological formation comprises:
implementing the defining, the minimizing, and the determining for a data point space defined over multiple frequency bands having different starting frequencies or by multiple different time periods having different starting times.

12. The method of claim 11 , further comprising:
adding 90° to the anisotropy angle estimated for the one of the bands or time periods if its eigenvectors and associated ones of the eigenvalues have a different ordering than other ones of the bands or time periods across the bands or time periods.

13. The method of claim 9 , wherein the matrix represents a set of data points in a space defined by receiver ring index position and frequency.

14. The method of claim 9 , further comprising:
diagonalizing the waveform matrix by a similarity transformation at each data point of the waveforms, by computing a matrix having columns comprising complex eigenvectors of the waveform matrix at each of the data points of the waveforms; and computing the eigenvalues of the eigenvectors.

15. The method of claim 14 , wherein diagonalizing the waveform matrix is accomplished using an analytical method or a numerical method.

16. The method of claim 14 , wherein the similarity transformation is a unitary similarity transformation.

17. The method of claim 9 , further comprising:
ordering the eigenvalues by ordering eigenvectors and corresponding ones of the eigenvalues across data points of the waveforms such that a difference in angle between the eigenvectors across the data points is minimized, and ordering the eigenvectors and the corresponding ones of the eigenvalues across the data points using a parameter correlated in a known way with the difference in the angle.

18. The method of claim 9 , further comprising:
ordering eigenvectors and corresponding ones of the eigenvalues across receiver rings at each frequency in a space defined by data points in the waveforms;
estimating eigenvalue propagators across the receiver rings from the ordered eigenvalues at each frequency in the space; and
ordering the eigenvectors, the corresponding ones of the eigenvalues, and the propagators over frequency by comparing relative phases associated with the estimates of the eigenvalue propagators at one frequency to relative phases associated with the estimates of the eigenvalue propagators at another frequency.

19. The method of claim 18 , wherein ordering the eigenvectors and corresponding ones of the eigenvalues comprises:
switching the eigenvectors, the corresponding ones of the eigenvalues, and the propagators; and
adding 90° to an anisotropy angle estimated for a frequency band at a first start frequency if a relative phase of the propagators has opposite sign when compared to a relative phase of the propagators at a second start frequency.

20. The method of claim 9 , wherein the objective function is defined by uti

Assignees

• Collins Mark V
• Cheng Arthur
• Halliburton Energy Services Inc

Inventors

• Collins Mark V
• Cheng Arthur

Classifications

• G01V1/50 Primary

  Analysing data · CPC title

• G01V1/36

  Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy · CPC title

• G01V2210/67

  Wave propagation modeling · CPC title

• G01V2210/1299

  Subsurface, e.g. in borehole or below weathering layer or mud line · CPC title

• G01V2210/27

  Other pre-filtering · CPC title