System and method of determining rock properties using terahertz-band dielectric measurements

Patent metadata

Field	Value
Publication number	US-9513213-B2
Application number	US-201313962816-A
Country	US
Kind code	B2
Filing date	Aug 8, 2013
Priority date	Aug 8, 2013
Publication date	Dec 6, 2016
Grant date	Dec 6, 2016

Abstract

A method of analyzing physical properties of a sample includes obtaining the sample and obtaining an electromagnetic spectrum of the sample using terahertz spectroscopy. A sample complex permittivity is computed from the electromagnetic spectrum of the sample. The method further includes estimating the constituents and the constituent fractions and computing an estimated effective complex permittivity based upon a model and the constituent fractions. The method further includes comparing the computed sample complex permittivity with the estimated effective complex permittivity in order to determine the physical properties the sample.

First claim

What is claimed is:

1. A method of analyzing a physical property of a sample, the method comprising:
operating a terahertz spectrometer to obtain an electromagnetic spectrum of the sample, wherein the electromagnetic spectrum comprises a terahertz spectrum;
computing a sample complex permittivity from the electromagnetic spectrum of the sample;
estimating types of constituents and constituent fractions;
computing an estimated effective complex permittivity based upon a model and the constituent fractions; and
comparing the computed sample complex permittivity with the estimated effective complex permittivity to determine the physical property of the sample.

2. The method of claim 1 , wherein the model is determined using a dielectric mixing law that combines known reference complex permitivities for each of the constituents using the constituent fractions.

3. The method of claim 2 , wherein the comparing further comprises:
minimizing a difference between the estimated effective complex permittivity and the sample complex permittivity by varying the estimated constituent fractions; and
determining the physical property of the sample using the estimated constituent fractions that minimize the difference between the estimated effective complex permittivity and the sample complex permittivity.

4. The method of claim 3 , wherein minimizing the difference comprises using a nonlinear least squares fit of the sample complex permittivity to the estimated effective complex permittivity using at least one estimated constituent fraction as a fitting parameter.

5. The method of claim 2 , wherein one of the estimated constituent fractions is a volume fraction of a rock matrix of the sample.

6. The method of claim 5 , wherein one of the estimated constituent fractions is selected from a group consisting of water volume fraction, oil volume fraction, and gas volume fraction.

7. The method of claim 2 , wherein one of the estimated fractions is a volume fraction of a mineral constituent of the sample.

8. The method of claim 7 , wherein the mineral constituent of the sample is selected from a group consisting of smectite and illite.

9. The method of claim 2 , wherein one of the estimated constituent fractions is a volume fraction of smectite and another of the estimated fractions is a volume fraction of illite.

10. The method of claim 2 , wherein an initial estimate of at least one of the estimated fractions is obtained by Fourier transform spectroscopy of the sample.

11. The method of claim 2 , wherein the dielectric mixing law is selected from a group consisting of a Maxwell Garnett (MG) model, Polder and van Santen (PvS) model, Bruggeman model, and a Landau, Lifshitz, Looyenga (LLG) model.

12. The method of claim 1 , further comprising:
creating the sample, wherein the creating comprises:
obtaining a raw earth sample;
crushing the raw earth sample to create a crushed earth sample;
mixing the crushed earth sample with a matrix material to create a crushed sample-matrix mixture; and
pressing the crushed sample-matrix mixture to form the sample.

13. The method of claim 1 , comprising obtaining the electromagnetic spectrum by terahertz time domain spectrometry.

14. A system for terahertz (“THz”) band dielectric measurements, the system comprising:
a THz spectrometer configured to obtain a THz spectrum of a sample; and
an analysis module configured to:
receive the THz spectrum of the sample;
compute a sample complex permittivity from the THz spectrum of the sample;
estimate types of constituents and constituent fractions;
compute an estimated effective complex permittivity based upon a model and the constituent fractions; and
compare the computed sample complex permittivity with the estimated effective complex permittivity to determine a physical property of the sample.

15. The system of claim 14 , wherein the analysis module further comprises:
a sample modeling engine configured determine the model of the sample using a dielectric mixing law that combines a plurality of known reference complex permitivities for each of the constituents using the estimated constituent fractions.

16. The system of claim 15 , wherein the analysis module further comprises:
an error estimation engine configured to:
minimize a difference between the estimated effective complex permittivity and the sample complex permittivity by varying the estimated constituent fractions; and
determine the physical properties of the sample using the estimated constituent fractions that minimize the difference between the estimated effective complex permittivity and the sample complex permittivity.

17. The system of claim 15 , wherein one of the estimated constituent fractions is a volume fraction of a rock matrix of the sample.

18. The system of claim 17 , wherein one of the estimated constituent fractions is selected from a group consisting of water volume fraction, oil volume fraction, and gas volume fraction.

19. The system of claim 14 , wherein one of the estimated constituent fractions is a volume fraction of smectite and another of the estimated fractions is a volume fraction of illite.

20. The system of claim 14 , wherein one of the estimated constituent fractions is a volume fraction of a mineral constituent of the sample.

21. The system of claim 20 , wherein the mineral constituent of the sample is selected from a group consisting of smectite and illite.

22. The system of claim 14 , wherein the THz spectrometer is configured as a THz time domain spectrometer.

23. A non-transitory computer readable medium comprising computer readable program code embodied therein, that, when executed on a processor, causes the processor to:
receive an electromagnetic spectrum of a sample that comprises a terahertz spectrum;
compute a sample complex permittivity from the electromagnetic spectrum of the sample;
estimate types of constituents and constituent fractions;
compute an estimated effective complex permittivity based upon a model and the constituent fractions; and
compare the computed sample complex permittivity with the estimated effective complex permittivity to determine a physical property the sample.

24. The non-transitory computer readable medium of claim 23 , wherein the model is determined using a dielectric mixing law that combines a plurality of known reference complex permitivities for each of the constituents using the constituent fractions.

25. The non-transitory computer readable medium of claim 24 , further comprising computer readable program code embodied therein, that, when executed on a processor, causes the processor to:
compare the computed sample complex permittivity with the estimated effective complex permittivity by minimizing a difference between the estimated effective complex permittivity and the sample complex permittivity by varying the estimated constituent fractions; and
determine the physical property of the sample using the estimated constituent fractions that minimize the difference between the estimated effective complex permittivity and the sample complex permittivity.

26. The non-transitory computer readable medium of claim 24 , wherein one of the estimated fractions is a volume fraction of a rock matrix of the sample.

27. The non-transitory computer readable medium of claim 24 , further comprising computer readable program code embodied therein, that, when

Assignees

• Chieffo Logan
• Averitt Richard
• Speck Andrew
• Herron Michael M
• Andrews A Ballard
• Van Hal Ronald E G
• Schlumberger Technology Corp

Inventors

• Chieffo Logan
• Averitt Richard
• Speck Andrew
• Herron Michael M
• Andrews A Ballard
• Van Hal Ronald E G

Classifications

• G01N21/3586 Primary

  by Terahertz time domain spectroscopy [THz-TDS] · CPC title

• G01N33/241

  for hydrocarbon content · CPC title