Method for estimating pore structure of porous materials and its application to determining physical properties of the materials
A method is disclosed for modeling a pore structure of a porous material. The method includes generating an initial model of the porous material. The initial model includes an interconnected network of passages including a first geometric shape and a second geometric shape. At least one petrophysical property is calculated from the initial model. At least one geometric parameter in at least one passage of the first geometric shape is adjusted. The at least one petrophysical property is recalculated, and the adjusting and recalculating are repeated until a difference between the calculated petrophysical property and a measured petrophysical property reaches a minimum. In some embodiments, the adjusting and recalculating are repeated until an indeterminacy condition is satisfied.
1. A method for modeling a pore structure of porous materials, comprising:
generating an initial model of the porous material, the initial model comprising an interconnected network of passages including a first geometric shape and a second geometric shape;
calculating at least one petrophysical property from the initial model;
adjusting at least one geometric parameter of the first geometric shape of at least one passage;
recalculating the at least one petrophysical property;
repeating the adjusting and recalculating until a difference between the calculated petrophysical property and a measured petrophysical property reaches a minimum and an indeterminacy condition is satisfied; and
at least one of recording and tabulating results of the repeated adjusting and calculating.
2. The method of claim 1 wherein the indeterminacy condition comprises a minimum mean size of the passages in the interconnected network.
3. The method of claim 1 wherein the indeterminacy condition comprises a maximum mean aspect ratio of the first geometric shape.
4. The method of claim 1 further comprising determining an entropy of the initial model based on the first geometric shape, and repeating the adjusting and recalculating until the entropy reaches a maximum.
5. The method of claim 1 wherein the initial model comprises average length.
6. The method of claim 5 wherein the average length is related to a geometric parameter of the second geometric shape.
7. The method of claim 6 wherein the second geometric shape comprises a sphere, and the geometric parameter of the second geometric shape comprises average radius.
8. The method of claim 1 wherein the at least one geometric parameter comprises aspect ratio.
9. The method of claim 1 wherein the at least one geometric parameter comprises a cross-sectional dimension.
10. The method of claim 1 further comprising setting the initial model to include a fractional volume of total pore space occupied by features of the first geometric shape and a fractional volume of features occupied by the second geometric shape.
11. The method of claim 1 wherein the first geometric shape comprises tubes.
12. The method of claim 1 wherein the second geometric shape comprises spheres.
13. The method of claim 1 wherein the second geometric shape is adapted to have substantially no effect on transport properties of the Earth formation.
14. The method of claim 13 wherein the transport properties comprise at least one of fluid permeability and formation resistivity factor.
15. The method of claim 1 wherein the at least one petrophysical property comprises at least one of nuclear magnetic resonance relaxation distribution, formation resistivity factor, fluid permeability, mercury injection capillary pressure, compressional wave acoustic velocity and shear wave acoustic velocity.
16. The method of claim 1 further comprising adding pore space to the initial model, the added pore space comprising voids which are modeled as substantially separate from the network of interconnected passages, the voids substantially not in hydraulic communication with each other.
17. The method of claim 16 wherein the adjusting comprises changing an amount of the additional porosity.
18. The method of claim 1 further comprising adding pore space to the initial model, added pore space comprising a second network of interconnected passages, the second network substantially separate from the network of interconnected passages including the first and second geometric shapes.
19. The method of claim 18 wherein the second network of passages comprises the first geometric shape, and the adjusting comprises at least one of changing amount of the additional porosity and changing at least one geometric parameter of the first geometric shape in the second network.
20. The method of claim 1 further comprising:
varying the at least one petrophysical property;
repeating the generating an initial model of the Earth formation, repeating the calculating at least one petrophysical property from the initial model, repeating the adjusting at least one geometric parameter of the first geometric shape of at least one passage, repeating the recalculating the at least one petrophysical property, repeating the adjusting and recalculating until a difference between the calculated petrophysical property and a measured petrophysical property reaches a minimum; and
determining a degree of correspondence between the varied petrophysical property, and the calculated petrophysical property.
21. The method of claim 20 further comprising removing the at least one petrophysical property as a constraint on the initial model and the adjusted model if the degree of correspondence falls below a selected threshold.
22. The method of claim 1 further comprising:
varying at least one model parameter;
repeating the generating an initial model of the Earth formation, repeating the calculating at least one petrophysical property from the initial model, repeating the adjusting at least one geometric parameter of the first geometric shape of at least on passage, repeating the recalculating the at least one petrophysical property, repeating the adjusting and recalculating until a difference between the calculated petrophysical property and a measured petrophysical property reaches a minimum; and
determining a degree of correspondence between the at least one petrophysical property, and the calculated petrophysical property.
23. A computer program stored in a computer readable medium, the program having logic operable to cause a programmable computer to perform steps comprising:
generating an initial model of a porous material, the initial model comprising an interconnected network of passages including a first geometric shape and a second geometric shape;
calculating at least one petrophysical property from the initial model;
adjusting at least one geometric parameter of the first geometric shape of at least one passage;
recalculating the at least one petrophysical property;
repeating the adjusting and recalculating until a difference between the calculated petrophysical property and a measured petrophysical property reaches a minimum and an indeterminacy condition is satisfied; and
at least one of recording and tabulating results of the repeated adjusting and calculating.
24. The computer program of claim 23 wherein the indeterminacy condition comprises a minimum mean size of the passages in the interconnected network.
25. The computer program of claim 23 wherein the indeterminacy condition comprises a maximum mean aspect ratio of the first geometric shape.
26. The computer program of claim 23 further comprising logic operable to cause the computer to perform determining an entropy of the model based on the first geometric shape, and repeating the adjusting and recalculating until the entropy reaches a maximum.
27. The computer program of claim 23 wherein the initial model comprises average length.
28. The computer program of claim 27 wherein the average length is related to a geometric parameter of the second geometric shape.
29. The computer program of claim 28 wherein the second geometric shape comprises a sphere, and the geometric parameter of the second geometric shape comprises average radius.
30. The computer program of claim 23 wherein the at least one geometric parameter comprises aspect ratio.
31. The computer program of claim 23 wherein the at least one geometric parameter comprises and cross-sectional dimension.
32. The computer program of claim 23 further comprising logic operable to cause the computer to perform setting the initial model to include a fractional volume of total pore space occupied by features of the first geometric shape and a fractional volume of features occupied by the second geometric shape.
33. The computer program of claim 23 wherein the first geometric shape comprises tubes.
34. The computer program of claim 23 wherein the second geometric shape comprises spheres.
35. The computer program of claim 23 wherein the second geometric shape is adapted to have substantially no effect on transport properties of the porous material.
36. The computer program of claim 35 wherein the transport properties comprise at least one of fluid permeability and formation resistivity factor.
37. The computer program of claim 23 wherein the at least one petrophysical property comprises at least one of nuclear magnetic resonance relaxation distribution, formation resistivity factor, fluid permeability, mercury injection capillary pressure, compressional wave acoustic velocity and shear wave acoustic velocity.
38. The computer program of claim 23 further comprising logic operable to cause the computer to perform adding pore space to the initial model, the added pore space comprising voids which are substantially separate from the network of interconnected passages, the voids substantially not in hydraulic communication with each other.
39. The computer program of claim 38 wherein the adjusting comprises changing an amount of the additional porosity.
40. The computer program of claim 23 further comprising logic operable to cause the computer to perform adding pore space to the initial model, added pore space comprising a second network of interconnected passages, the second network substantially separate from the network of interconnected passages including the first and second geometric shapes.
41. The computer program of claim 40 wherein the second network of passages comprises the first geometric shape, and the adjusting comprises at least one of changing amount of the additional porosity and changing at least one geometric parameter of the first geometric shape in the second network.
42. The computer program of claim 23 further comprising logic operable to cause the computer to perform:
varying the at least one petrophysical property;
repeating the generating an initial model of the porous material, the initial model comprising an interconnected network of passages including a first geometric shape and a second geometric shape, repeating the calculating at least one petrophysical property from the initial model, repeating the adjusting at least one geometric parameter of the first geometric shape of at least one passage, repeating the recalculating the at least one petrophysical property, repeating the adjusting and recalculating until a difference between the calculated petrophysical property and a measured petrophysical property reaches a minimum; and
determining a degree of correspondence between the varied petrophysical property and the calculated petrophysical property.
43. The computer program of claim 42 further comprising logic operable to cause the computer to remove the at least one petrophysical property as a constraint on the initial model and the adjusted model if the degree of correspondence falls below a selected threshold.
44. The computer program of claim 23 further comprising logic operable to cause the computer to perform:
varying at least one model parameter;
repeating the generating an initial model of the porous material, repeating the calculating at least one petrophysical property from the initial model, repeating the adjusting at least one geometric parameter of the first geometric shape of at least one passage, repeating the recalculating the at least one petrophysical property, repeating the adjusting and recalculating until a difference between the calculated petrophysical property and a measured petrophysical property reaches a minimum; and
determining a degree of correspondence between the at least one petrophysical property and the calculated petrophysical property.