Method of estimating elastic properties of kerogen using multi-scale data integration
The present disclosure is directed to numerically estimating the shear modulus of Kerogen by using a combination of mineralogy from digital image analysis and sonic log analysis, when measured data on only one elastic constant (Bulk, Young's or P-wave modulus) is available. In some instances, elastic properties predicted from the digital images are compared with sonic, shear, and density logs, to estimate the shear modulus of kerogen. As a one-to-one correspondence is not expected between the core sub-samples and the rock unit sampled by the well logs, cross-property relations can be used to identify the suitability of the effective medium models and to iteratively determine the shear modulus of kerogen.
1. A method comprising:
developing, by a computing device, a first numerical model of a rock sample, the first numerical model comprising modeled elastic properties of the rock sample including a first estimated elastic property of the rock sample;
calculating, by the computing device, elastic properties of the rock sample;
obtaining, by the computing device, a first result by comparing a first one of the calculated elastic properties of the rock sample with a second one of the calculated elastic properties of the rock sample;
obtaining, by the computing device, a second result by comparing a first one of the modeled elastic properties of the rock sample with a second one of the modeled elastic properties of the rock sample;
determining, by the computing device, that the first result and the second result are within a first acceptable threshold range, and in response to determining that the first result and the second result are not within the first acceptable threshold range, repeating calculations of elastic properties using a revised estimate for shear modulus of kerogen, obtaining the first result, and obtaining the second result until a determination is made that the first result and the second result are within the first acceptable threshold range;
comparing, by the computing device, based upon the determining that the first result and the second result are within a first acceptable threshold range, at least one of a second numerical model with the calculated elastic properties of the rock sample, wherein the at least one of a second numerical model comprises elastic properties of another rock sample;
determining, by the computing device, that the elastic properties of the second numerical model and the calculated elastic properties of the rock sample are within a second acceptable threshold range, and in response to determining that the elastic properties of the second numerical model and the calculated elastic properties of the rock sample are not within the second acceptable threshold range, repeating comparing the at least one of the second numerical model with the calculated elastic properties of the rock sample until a determination is made that the elastic properties of the second numerical model and the calculated elastic properties of the rock sample are within the second acceptable threshold range; and
identifying, by the computing device, a particular one of the elastic properties of the at least one of a second numerical model,
wherein a type of the particular one of the elastic properties of the at least one of a second numerical model is substantially the same as a type of the first estimated elastic property of the rock sample.
2. The method according to claim 1 , wherein the developing a first numerical model of a rock sample includes determining a volume fraction of minerals comprising the rock sample and a volume fraction of organic matter comprising the rock sample.
3. The method according to claim 2 , wherein the organic matter is kerogen.
4. The method of claim 3 , where the type of the particular one of the elastic properties of the at least one of a second numerical model is a shear modulus of the kerogen.
5. The method of claim 1 , wherein the calculating elastic properties of the rock sample includes calculating bulk modulus, shear modulus, compressional velocity, shear velocity, and density.
6. The method of claim 1 , wherein the first set of calculated elastic properties of the rock sample includes a density of the rock sample, and the first portion of the first numerical model of the rock sample includes a volume fraction of kerogen.
7. The method of claim 6 , wherein the second set of the calculated elastic properties of the rock sample includes a calculated compressional velocity and a calculated shear velocity of the rock sample, and the second portion of the first numerical model of the rock sample includes a numerical modeled compressional velocity and a numerical modeled shear velocity of the rock sample.
8. The method of claim 1 , wherein the at least one of a second numerical model includes a plurality of effective-medium models.
9. The method of claim 8 , wherein each of the plurality of effective-medium models includes a third set of calculated elastic properties of another rock sample.
10. The method of claim 9 , wherein a type of the second set of calculated elastic properties of the rock sample is substantially the same as a type of the third set of calculated elastic properties of another rock sample.
11. The method of claim 1 , further comprising:
determining that no second numerical models have been identified as having the calculated elastic properties of the rock sample; and
widening the acceptable threshold range until a minimum number of the second numerical models are identified.
12. A non-transitory computer readable medium storing instructions to cause a processor to perform operations comprising:
developing, by a computing device, a first numerical model of a rock sample, the first numerical model comprising modeled elastic properties of the rock sample including a first estimated elastic property of the rock sample;
calculating, by the computing device, elastic properties of the rock sample;
obtaining, by the computing device, a first result by comparing a first one of the calculated elastic properties of the rock sample with a second one of the calculated elastic properties of the rock sample;
obtaining, by the computing device, a second result by comparing a first one of the modeled elastic properties of the rock sample with a second one of the modeled elastic properties of the rock sample;
determining, by the computing device, that the first result and the second result are within a first acceptable threshold range, and in response to determining that the first result and the second result are not within the first acceptable threshold range, repeating calculations of elastic properties using a revised estimate for shear modulus of kerogen, obtaining the first result, and obtaining the second result until a determination is made that the first result and the second result are within the first acceptable threshold range;
comparing, by the computing device, based upon the determining that the first result and the second result are within a first acceptable threshold range, at least one of a second numerical model with the calculated elastic properties of the rock sample, wherein the at least one of a second numerical model comprises elastic properties of another rock sample;
determining, by the computing device, that the elastic properties of the second numerical model and the calculated elastic properties of the rock sample are within a second acceptable threshold range, and in response to determining that the elastic properties of the second numerical model and the calculated elastic properties of the rock sample are not within the second acceptable threshold range, repeating comparing the at least one of the second numerical model with the calculated elastic properties of the rock sample until a determination is made that the elastic properties of the second numerical model and the calculated elastic properties of the rock sample are within the second acceptable threshold range; and
identifying, by the computing device, a particular one of the elastic properties of the at least one of a second numerical model,
wherein a type of the particular one of the elastic properties of the at least one of a second numerical model is substantially the same as a type of the first estimated elastic property of the rock sample.
13. The computer readable medium according to claim 12 , wherein the developing a first numerical model of a rock sample includes determining a volume fraction of minerals comprising the rock sample and a volume fraction of organic matter comprising the rock sample.
14. The computer readable medium according to claim 13 , wherein the organic matter is kerogen.
15. The computer readable medium of claim 14 , where the type of the particular one of the elastic properties of the at least one of a second numerical model is a shear modulus of the kerogen.
16. The computer readable medium of claim 12 , wherein the calculating elastic properties of the rock sample includes calculating bulk modulus, shear modulus, compressional velocity, shear velocity, and density.
17. The computer readable medium of claim 12 , wherein the first set of calculated elastic properties of the rock sample includes a density of the rock sample, and the first portion of the first numerical model of the rock sample includes a volume fraction of kerogen.
18. The computer readable medium of claim 17 , wherein the second set of the calculated elastic properties of the rock sample includes a calculated compressional velocity and a calculated shear velocity of the rock sample, and the second portion of the first numerical model of the rock sample includes a numerical modeled compressional velocity and a numerical modeled shear velocity of the rock sample.
19. The computer readable medium of claim 12 , wherein the at least one of a second numerical model includes a plurality of effective-medium models.
20. The computer readable medium of claim 19 , wherein each of the plurality of effective-medium models includes a third set of calculated elastic properties of another rock sample.
21. The computer readable medium of claim 20 , wherein a type of the second set of calculated elastic properties of the rock sample is substantially the same as a type of the third set of calculated elastic properties of another rock sample.