Method for determining elastic-wave attenuation of rock formations using computer tomograpic images thereof

1. A method for estimating at least one elastic-wave-related property of a porous material from a sample thereof, comprising:

making a three dimensional tomographic image of the sample of the material;

segmenting the image into pixels each representing pore space or rock grain;

determining bulk modulus and shear modulus of the porous material from the segmented image at a frequency corresponding to mobile fluid;

determining bulk modulus and shear modulus of the porous material from the segmented image at a frequency corresponding to immobile fluid;

determining the at least one elastic-wave-related property from the mobile fluid and immobile fluid moduli, wherein the at least one elastic-wave-related property comprises maximum inverse quality factor for shear elastic waves or maximum inverse quality factor for compressional elastic waves or both; and

at least one of storing and displaying the at least one elastic-wave-related property.

2. The method of claim 1 wherein the tomographic image comprises a value of gray scale allocated to each of a plurality of pixels in the image.

3. The method of claim 2 wherein the segmenting comprises:

(a) determining an initial gray scale threshold for each of pore space and rock grain;

(b) allocating each pixel in the image to pore space or rock grain for each pixel meeting threshold criteria for each of the pore space and rock grain thresholds, respectively;

(c) interrogating pixels adjacent to each seed;

(d) allocating the interrogated adjacent pixels previously not allocated to a seed to the pore space seed or the rock grain seed based on threshold criteria; and

(e) repeating (c) and (d) until substantially all pixels in the image are allocated to the rock grain or the pore space.

4. The method of claim 3 wherein the determining initial gray scale thresholds comprises histogram analysis of the tomographic image.

5. The method of claim 3 wherein the allocating interrogated adjacent pixels comprises determining a difference between a gray scale value and a gray scale value of the adjacent pixel in the seed, and allocating the interrogated pixel to the seed if the difference falls below a selected threshold.

6. The method of claim 5 further comprising:

determining whether unallocated pixels exist in the image;

increasing the selected difference threshold; and

repeating the interrogating adjacent and allocating the interrogate pixels having gray scale values wherein the difference is below the increased difference threshold.

7. The method of claim 1 wherein the material comprises a rock formation.

8. The method of claim 1 further comprising determining bulk and shear moduli for immobile fluid and mobile fluid for the rock formation in the presence of at least partial hydrocarbon saturation in pore spaces thereof, and determining a maximum inverse quality factor for the rock formation in the presence of partial hydrocarbon saturation of the pore spaces thereof.

9. The method of claim 1 wherein the determining the bulk and shear moduli for mobile fluid comprises determining bulk and shear moduli for the rock formation wherein no material is disposed in pore spaces thereof, and using Gassmann's fluid substitution to estimate the bulk and shear moduli in the presence of mobile fluid in the pore spaces.

10. The method of claim 1 wherein the determining the bulk and shear moduli for immobile fluid comprises using a bulk modulus for water or brine solutions thereof in pore spaces of the rock formation.

11. The method of claim 1 wherein the determining the mobile fluid moduli and the immobile fluid moduli comprises performing finite element analysis of the segmented image.