US20140366621A1 - Gas Sorption Analysis Of Unconventional Rock Samples - Google Patents

Gas Sorption Analysis Of Unconventional Rock Samples Download PDF

Info

Publication number
US20140366621A1
US20140366621A1 US14/475,436 US201414475436A US2014366621A1 US 20140366621 A1 US20140366621 A1 US 20140366621A1 US 201414475436 A US201414475436 A US 201414475436A US 2014366621 A1 US2014366621 A1 US 2014366621A1
Authority
US
United States
Prior art keywords
sample
gas
shale
pore volume
pore
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Abandoned
Application number
US14/475,436
Inventor
John J. Valenza
Nicholas J. Drenzek
Flora Marques
Hendrik Grotheer
Dean M. Willberg
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Schlumberger Technology Corp
Original Assignee
Schlumberger Technology Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Schlumberger Technology Corp filed Critical Schlumberger Technology Corp
Priority to US14/475,436 priority Critical patent/US20140366621A1/en
Publication of US20140366621A1 publication Critical patent/US20140366621A1/en
Abandoned legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/082—Investigating permeability by forcing a fluid through a sample
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
    • G01N2015/0866—Sorption

Definitions

  • This disclosed subject matter is generally related to gas shale formations, and more particularly to methods for accessing the maturity, free gas content, and transport characteristics or pore connectivity of gas shale formations.
  • OM organic matter
  • Rock-Eval pyrolysis is one of these techniques and consists of subjecting a small shale aliquot to a temperature cycle and monitoring the amount of gas-phase products evolving as a function of temperature. Typically, three generation events occur.
  • the first event, 51 occurs as the sample is exposed to a temperature of 300° C. and consists of free hydrocarbons (gas and oil) that are volatilized.
  • the second event, S2 occurs as the temperature increases to 550° C. at a rate of 25° C. per minute and consists of cracking non-volatile hydrocarbons, e.g. kerogen.
  • T max The temperature at which this generation reaches a maximum
  • T max may also be influenced by other OM characteristics, for example, initial composition and catalytic association with various minerals. Typical ranges of T max for immature are 400-430° C., mature 435° C.-450° C. and over-mature 450° C. A large proportion of gas shales are over-mature, but T max is relevant for assessing the degree to which the shale is over-mature.
  • S3 is a measure of the amount of carbon dioxide that is associated with hydrocarbon cracking and S3 is useful for estimating the oxygen content of OM.
  • a second technique for assessing OM maturity consists of observing the percentage of light reflected from a shale sample, Ro %, where the percentage is calibrated against a standard that reflects 100% of the light. OM macerals become more glass-like and therefore reflect more light with increasing maturity. Typical ranges of Ro % are as follows: depositional (immature)—0.2-0.7; oil producing (immature-mature)—0.7-1.2 and gas producing (mature to over-mature) 1.2 ⁇ Ro % ⁇ 5.0.
  • This technique was originally devised for determining the rank of coal and may be tedious, subjective, and imprecise. Moreover, due to the absence of vitrinite-producing land plants before ca. 360 million years ago, this technique cannot be applied to gas shales of older provenance.
  • a technique for measuring the free gas content consists of measuring the high pressure, high temperature methane capacity of a shale by fitting successive measurements of gas uptake as a function of pressure with a Langmuir isotherm to quantify total gas and adsorbed gas, where the difference is the free gas.
  • employing the Langmuir isotherm implies that analysis gas chemisorbs on the solid.
  • methane does not chemisorb on the material constituents of shale.
  • information obtained by applying the Langmuir analysis may not be representative.
  • Measuring the high pressure, high temperature methane capacity is time consuming, cumbersome and necessitates high pressure, for example, reservoir pressure, rendering it unattractive for applications at the wellsite.
  • systems and methods of analyzing rock samples from an unconventional hydrocarbon reservoir include performing gas sorption on a sample of rock from the unconventional hydrocarbon reservoir; and determining at least one characteristic associated with the sample based at least in part on the gas sorption.
  • a method for determining a location of a fracture zone for stimulation in a shale reservoir includes obtaining a shale sample from the shale reservoir and grinding and sieving the shale sample. The method further includes degassing the shale sample and performing gas sorption on the shale sample from the shale reservoir to determine at least one characteristic associated with the shale sample based at least in part on the gas sorption. Finally, the method includes determining the fracture zone from the one or more characteristics.
  • a method of analyzing a rock sample includes obtaining the rock sample and measuring a pore volume or porosity on a plurality of particles of varying sizes from the rock sample. The method further includes, comparing the measured pore volume or porosity to determine a correlation length or pore connectivity in the rock sample.
  • unconventional reservoir includes reservoirs having an unconventional microstructure, such as having submicron pore size, and/or substantial amounts of primary organic matter such as kerogen.
  • unconventional reservoirs include hydrocarbon-bearing shales such as gas shales and oily shales.
  • FIG. 1 shows a sampling tool being deployed in a wellbore and an analysis facility, according to some embodiments
  • FIG. 2 is a plot of A s, N2 normalized by total organic content (TOC) versus vitrinite reflectance (VR 0 ) or the bitumen equivalent (VR E );
  • FIG. 3 is a plot of the change in surface area after bitumen extraction and combustion as a function of VR 0 or VR E ;
  • FIG. 4 is a plot of isotherms for shale of varying maturity
  • FIG. 5 is a plot of the natural log of V A , normalized by V M from BET analysis, versus the double natural log of P/P 0 ;
  • FIG. 6 shows sorption results using water vapor as the analysis gas, A S.H2O , normalized by the mass of illite+smectite (I+S) in the sample;
  • FIG. 7 is a plot of permeability, k, versus the difference between the pore volume measured on two different particle sizes.
  • FIG. 8 is a flow chart of an embodiment of the subject disclosure.
  • methods are disclosed for accessing the quality of a shale reservoir.
  • the information derived from the methods disclosed may be utilized for completion design.
  • Stimulation is the most costly process (in general at least 50% of a well cost) utilized for gas shale plays so gaining insight into the location of the best zones for stimulation is important.
  • zones with the most mature organic matter as mature organic matter contains oil-wet pores which may represent unimpeded pathways to the fractures formed during stimulation.
  • the gas in shale is in one of two states which are (a) tightly absorbed on the rock surface (liquid-like) or (b) free in the body of the pores (gas-like). Laboratory measurements indicate that initially the free-gas is produced and the adsorbed gas is produced when the reservoir pressure is reduced. It is therefore an advantage to identify mature zones with the greatest proportion of free gas.
  • the subject disclosure utilizes conventional gas sorption for characterizing the microstructure of gas shale.
  • Gas sorption is applicable for characterizing nanoporosity ( ⁇ 200 nm) and therefore is an ideal technique for characterizing the microstructure of gas shale.
  • a gas sorption study is performed at the boiling point of the analysis gas ( ⁇ 196° C. for nitrogen). Under these conditions the saturation vapor pressure, P 0 , of the analysis gas is equivalent to ambient pressure ( ⁇ 0.1 MPa).
  • Other analysis gases which may be used include argon, xenon, neon or any molecular or atomic probe known to be useful in adsorption studies.
  • the disclosure implicitly covers the normal thermodynamic conditions with respect to range of temperature and pressure, regularly imposed for the various analysis gases.
  • FIG. 1 shows a sampling tool being deployed in a wellbore and an analysis facility, according to some embodiments.
  • the sampling tool is a core sampling tool.
  • Wireline truck 110 is deploying wireline cable 112 into well 130 via well head 120 .
  • Wireline tool 114 is disposed on the end of the cable 112 in an unconventional subterranean formation 100 .
  • formation 100 is an unconventional reservoir, such as a hydrocarbon bearing shale reservoir.
  • Tool 114 includes a sampling tool as shown, in a non-limiting example a core sampling tool.
  • a wireline sampling tool is shown, according to other embodiments, other types of sampling tools are used such as while drilling and/or coiled tubing conveyed tools.
  • Samples 132 are from an unconventional rock formation 100 and are retrieved at the surface from the tool 114 and transported to an analysis facility 160 .
  • the analysis facility 160 can be located at the wellsite (which can be onshore or offshore) or it can be located remotely from the wellsite.
  • Facility 160 includes a gas sorption apparatus ( 152 ), one or more central processing units 140 , storage system 144 , communications and input/output modules 140 , a user display 146 and a user input system 148 .
  • Input/output modules 140 include modules to communicate with and control the gas sorption machine ( 152 ).
  • Facility 160 may also include a helium pycnometry (He-pyc), mercury intrusion porosimetry (MICP) or other apparatuses as known to those skilled in the art for characterizing rock samples.
  • He-pyc helium pycnometry
  • MIMP mercury intrusion porosimetry
  • microstructural characteristics of shale are dependent on the shale maturity.
  • the correlation between the microstructure and maturity is due to pores which form in the organic matter as hydrocarbons are generated and ultimately expelled. These pores are characterized by length scales that range from 10 ⁇ 9 -10 ⁇ 7 m.
  • the high surface to volume ratio of these nanopores is responsible for an increase in specific area, A s , with maturity.
  • An increase in maturity leads to a concomitant increase in porosity and surface area due to nanometric pores in the total organic content.
  • These pores increase pore connectivity and are formed as a result of hydrocarbon generation and expulsion. This increase in specific area, A s , with maturity is depicted in FIG. 2 .
  • FIG. 2 is a plot of A s, N2 normalized by TOC versus virtinite reflectance (VR 0 ) or the bitumen equivalent (VR E ). Thermal maturity is indicated by vitrinite reflectance, VR 0 , or the equivalent determined on bitumen, VR E (for higher maturity), which are acquired by standard organic petrology. A similar related metric commonly used is Tmax.
  • the plot shows that the surface area increases with maturity. The open circles are the raw data, while the filled surfaces show the surface area attributed to organic matter alone. Since A s, N2 represents that due to organic and clay (open symbols), there is considerable scatter in the data, especially in the samples that are moderately mature.
  • a s, N2 /TOC is reduced by an amount attributable to clay (filled symbols).
  • the proportion of A s, N2 attributable to clay is given by the XRD mineralogy, the relationship displayed in Equation 1 (free exponent) below for the proportion of clay accessible to the analysis gas, and G the ratio of the surface area of pure clay as measured by nitrogen to that measured by water.
  • a S , N ⁇ ⁇ 2 ⁇ ( no ⁇ ⁇ clay ) A S , N ⁇ ⁇ 2 ⁇ ( with ⁇ ⁇ clay ) - B * ( I + S ) * G * ( TOC 100 ) - 0.998 Equation ⁇ ⁇ 1
  • FIG. 3 is a plot showing the change in surface area after bitumen extraction, and combustion as a function of VR 0 or VR E . Also shown is the volume of micropores V ⁇ determined using, in a non-limiting example, the t-plot analysis on gas sorption data from the as received material. V ⁇ could also be measured directly by performing gas sorption with carbon dioxide. The high specific surface area of the kerogen is due to the increase in nanoporosity, including V ⁇ , with maturity.
  • FIG. 4 is a graph of isotherms for shale of varying maturity. The presence of micropores elevates V A (P/P 0 ⁇ 0) in the mature shale. This observation provides a rapid maturity index.
  • the horizontal dashed line in FIG. 4 is an example of a lower bound on V A (P/P 0 ⁇ 0) above which a sample would be considered mature.
  • V A adsorbed gas
  • P/P 0 ⁇ 0 e.g. maturity
  • a S surface area
  • V P pore volume
  • these quantities can be utilized to estimate the total gas in place or hydrogen index given the reservoir pressure and the associated gas density profile in the pore (e.g. higher density in the adsorbed phase near the pore wall, decaying to bulk density in the center of the pore.).
  • the Brunauer, Emmet and Teller analysis (BET) is utilized to determine, A s given V A at several pressures below P/P 0 ⁇ 0.3.
  • BET is a statistical analysis of the buildup of gas multilayer's on a solid surface. The primary result is an expression for the number of gas molecules in a monolayer.
  • the t-plot technique is used to determine the volume of micropores, V ⁇ , given V A over the entire range of pressure.
  • This approach utilizes a transformation of variables between P and the thickness of adsorbed gas, t A , to deduce V ⁇ , where the transformation of variables is determined on a chemically similar non-porous sample.
  • V ⁇ 0
  • a plot of V A versus t A passes through the origin.
  • FIG. 5A is a plot of the natural log of V A , normalized by V M from BET analysis, versus the double natural log of P/P 0 .
  • the plot is linear over nearly two orders of magnitude in pressure, which corresponds to nearly two orders of magnitude in pore size.
  • Shales are distinguished from conventional reservoir rock by their more abundant clay and OM content.
  • a s, H2O normalized by the mass of illite and smectite (I+S) exhibits a power law dependence on TOC.
  • FIG. 6 plots A s, H2O normalized by the mass of illite and smectite (I+S) versus TOC. The data corresponds to a number of samples across the spectrum of maturity.
  • FIG. 6 there are two power law fits to the A s, H2O /(I+S) data, one with the exponent as a free parameter (gray) and the other with the exponent fixed (black). Water vapor is used as it will not condense on organic surfaces.
  • a s, H2O inferred from these measurements is indicative of that due to exposed inorganic surfaces.
  • inorganic surface area is dominated by clay, where the specific surface area of clay as measured by water is A s, H2O ⁇ 300-400 m 2 /g.
  • the correlation length can also be determined by running the gas sorption measurement on a single sample by varying the particle size used in the analysis. This will give the discrete range of particle sizes in which ⁇ lies (e.g., bounded below by the largest particle size which yields a consistent A S , and above by the first particle size which yields lower A S ).
  • Embodiments of the subject disclosure comprise techniques for characterizing shale using gas sorption.
  • the permeability of shale is very low (10 ⁇ 19 -10 ⁇ 20 m 2 ) so the particle size is reduced (0.001-10 mm) so the gas molecules may probe the pore network rapidly.
  • samples are degassed before running experiments.
  • a shale sample is retrieved and the sample is grinded and then sieved ensuring all material passes through a sieve with 150 ⁇ m openings.
  • the sample size is on the order of 1 g.
  • the sample is then loaded into an analysis tube and the mass of the sample is noted.
  • the analysis tube is then loaded onto an analyzer to degas. Gas sorption analysis is then performed.
  • V A is noted for a plurality of P/P 0 in the pressure range 0.001 ⁇ P/P 0 ⁇ 0.995. This includes a point very near the upper bound.
  • V A is also noted in the limit, P/P 0 ⁇ 0.
  • the pore volume (V p ) may be noted from P/P 0 and the pore volume available to free gas (V F ), and the gas in place or hydrogen index may also be determined as indicated above.
  • the particle size may be varied to quantify the associated effect on access to the porosity. If the particle size is varied over a broad enough range, the data will provide insight into the correlation length of the shale and the total pore volume or porosity.
  • the particle size may be varied and helium pycnometry (He-pyc), mercury intrusion porosimetry (MICP), or other techniques as known to those skilled in the art for measuring pore volume or porosity are used to yield information and also to determine the correlation length of shale.
  • He-pyc determines the volume of solids in a sample of known mass, by allowing a known amount of Helium gas to expand into a container of known volume that holds the sample of interest.
  • the volume of solids is given by the equilibrium pressure using the ideal gas law (corrected for non-ideality).
  • the skeletal density is taken as the ratio between the sample mass and the volume of solids that make up the porous sample.
  • MICP forces mercury into a porous body by incrementally increasing the pressure applied to the mercury.
  • the size of pores the mercury can access is inversely proportional to the pressure.
  • the sample volume is the difference between the volume of the sample holder and the volume of the mercury required to surround but not invade the sample.
  • the pore size distribution can be determined from the volume of mercury injected at each incremental increase in pressure, and the well known relationship between the pore size and pressure as given by the Washburn equation.
  • the pore volume is determined from the total volume of mercury injected at the highest pressure, and the skeletal density is given by the mass of the sample and the volume of solids that make up the porous medium. The latter quantity is taken as the difference between the volume of the sample holder, and the total amount of mercury required to surround the sample, and fill the accessible pores.
  • Bulk density, ⁇ B is measured by weighing in air and the volume is determined by measuring the displacement of a liquid.
  • Gas permeability is determined using a modified pycnometer technique by analyzing the kinetics of pressure equilibration at a plurality of pressures. With this information, the liquid permeability, k, is inferred from the y-intercept of a linear fit to the correlation with the experimental pressure (See Klinkenberg, “The permeability of porous media to liquids and gases,” Drill, Proc. API (1941), 200-213).
  • the skeletal density ⁇ S is determined from the equilibrium pressure during the gas permeability measurements utilizing the ideal gas law.
  • FIG. 7 is a plot of permeability k versus the difference between the measured and calculated pore volume.
  • the measured pore volume is measured by nitrogen gas sorption, V P,meas , and the calculated pore volume is determined from ⁇ B and ⁇ S , V P,calc :
  • V P , calc 1 ⁇ B - 1 ⁇ S Equation ⁇ ⁇ 2
  • the measurement of ⁇ S is on a material that has a characteristic dimension that is at least one order of magnitude greater than that used in the gas sorption experiment, in a non-limiting example 10 ⁇ 3 m vs. 10 ⁇ 4 m.
  • the difference in size of the particles can be visually seen in FIG. 7 with the particles ( 701 ) used in the measured gas sorption being smaller than the particles ( 703 ) used where ⁇ S is measured.
  • any combination of the results from gas sorption, and similar results from other characterization techniques e.g., He-pyc or MIP using a plurality of varying particle sizes (0.001-50 mm) will yield information which may be used to determine the correlation length of shale.
  • the microstructure of mature shale has nanometer sized pores in the kerogen network. Therefore, since the surface area is expected to drop upon removal of organic matter from the shale matrix, a comparison of A s measured on samples in the native and organic free (via high temperature combustion) states is an alternative means for identifying mature zones (See FIG. 2 combustion data).
  • the samples are analyzed before and after combustion or bitumen extraction.
  • the shale samples are homogenized and split into two samples, a first and a second sample.
  • the two samples are grinded and then sieved ensuring all material passes through a sieve with 150 ⁇ m openings.
  • the sample sizes are on the order of 1 g.
  • a first sample is then combusted at 450° C. in air for at least 16 hours to remove organic matter or subjected to solvent reflux for several days to remove bitumen.
  • the first sample is dried.
  • a second sample is placed in a vacuum oven overnight.
  • the total organic content (TOC) is estimated in the first sample from mass loss due to combustion or bitumen removal.
  • the samples are then loaded into a clean sample tube and the mass of each of the samples is noted. Gas sorption analysis is then performed on both samples noting the V A of P/P 0 in the pressure range 0.01 ⁇ P/P 0 ⁇ 0.995, including a point at the upper bound.
  • the surface area (A s ) is determined by the BET method on both samples and the results are compared. If the surface area (A s ) is reduced by 200 m 2 /g/TOC after combustion the sample may be considered mature (see FIG. 3 —filled symbols). If the surface area (A s ) increases markedly upon bitumen extraction, the sample may be considered moderately mature (see FIG. 3 —open symbols).
  • the gas sorption data from the native-state sample is analyzed using the t-plot method to check for microporosity, if the sample contains micropores, the sample is mature.
  • the pore volume is determined from the volume of adsorbed gas (V A ) at the highest pressure and the pore volume available to free gas (V F ). (V F ) may also be determined as indicated above.
  • FIG. 8 depicts a flow chart of embodiments of the subject disclosure.
  • the shale sample ( 801 ) is grinded ( 803 ) and then sieved ( 805 ) ensuring all material passes through a sieve with a predefined opening, in a non-limiting example the opening is 150 ⁇ m.
  • the shale sample ( 801 ) is then degassed ( 807 ) in preparation for gas sorption analysis ( 809 ). From gas sorption analysis ( 809 ) a plurality of microstructural characteristics ( 823 ) are determined which include pore volume, surface area, pore size distribution and fractal dimension and correlation length or relative pore connectivity.
  • the shale sample ( 801 ) is grinded ( 803 ) and then sieved ( 805 ) and the sample is then analyzed using helium pycnometry (He-pyc) and bulk density.
  • He-pyc helium pycnometry
  • the pore volume ( 813 ) can be determined with this approach.
  • the shale sample ( 801 ) is grinded ( 803 ) and then sieved ( 805 ) and the sample is then analyzed using MICP ( 817 ). From MICP ( 817 ) a plurality of physical and microstructural characteristics ( 823 ) are determined which include bulk density, skeletal density, pore volume, surface area and pore size distribution.
  • the pore volume or porosity may also be measured on varying sized particles in the rock sample using any of the techniques disclosed or any techniques known to those skilled in the art for measuring the pore volume and the porosity in a rock sample and the results are compared to determine the correlation length or relative pore connectivity in the system.
  • Comparing any output 813 , 815 , 819 , 821 or 823 yields information regarding pore connectivity, more specifically the dimension over which the porosity is connected in 3D.
  • MIDP mercury intrusion porosimetry

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Dispersion Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Geology (AREA)
  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Solid-Sorbent Or Filter-Aiding Compositions (AREA)
  • Sampling And Sample Adjustment (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigating Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

Systems and methods for gas sorption analysis, or analogous practices, of samples from unconventional reservoirs are described. The described analysis of samples is used to determine various properties of unconventional reservoirs, which are used in evaluating their worth and producibility.

Description

    CROSS-REFERENCE TO RELATED APPLICATIONS
  • This application is a continuation of U.S. patent application Ser. No. 13/359,121 filed on Jan. 26, 2012.
    • FIELD
  • This disclosed subject matter is generally related to gas shale formations, and more particularly to methods for accessing the maturity, free gas content, and transport characteristics or pore connectivity of gas shale formations.
    • BACKGROUND
  • A number of techniques exist for accessing organic matter (OM) maturity in shales. Rock-Eval pyrolysis is one of these techniques and consists of subjecting a small shale aliquot to a temperature cycle and monitoring the amount of gas-phase products evolving as a function of temperature. Typically, three generation events occur. The first event, 51, occurs as the sample is exposed to a temperature of 300° C. and consists of free hydrocarbons (gas and oil) that are volatilized. The second event, S2, occurs as the temperature increases to 550° C. at a rate of 25° C. per minute and consists of cracking non-volatile hydrocarbons, e.g. kerogen. The temperature at which this generation reaches a maximum, is called Tmax, and may be interpreted as a measure of the shale maturity. Tmax may also be influenced by other OM characteristics, for example, initial composition and catalytic association with various minerals. Typical ranges of Tmax for immature are 400-430° C., mature 435° C.-450° C. and over-mature 450° C. A large proportion of gas shales are over-mature, but Tmax is relevant for assessing the degree to which the shale is over-mature. The third event, S3, is a measure of the amount of carbon dioxide that is associated with hydrocarbon cracking and S3 is useful for estimating the oxygen content of OM.
  • A second technique for assessing OM maturity consists of observing the percentage of light reflected from a shale sample, Ro %, where the percentage is calibrated against a standard that reflects 100% of the light. OM macerals become more glass-like and therefore reflect more light with increasing maturity. Typical ranges of Ro % are as follows: depositional (immature)—0.2-0.7; oil producing (immature-mature)—0.7-1.2 and gas producing (mature to over-mature) 1.2<Ro %<5.0. This technique was originally devised for determining the rank of coal and may be tedious, subjective, and imprecise. Moreover, due to the absence of vitrinite-producing land plants before ca. 360 million years ago, this technique cannot be applied to gas shales of older provenance.
  • A technique for measuring the free gas content consists of measuring the high pressure, high temperature methane capacity of a shale by fitting successive measurements of gas uptake as a function of pressure with a Langmuir isotherm to quantify total gas and adsorbed gas, where the difference is the free gas. In general, employing the Langmuir isotherm implies that analysis gas chemisorbs on the solid. However, methane does not chemisorb on the material constituents of shale. In this regard, information obtained by applying the Langmuir analysis may not be representative. Measuring the high pressure, high temperature methane capacity is time consuming, cumbersome and necessitates high pressure, for example, reservoir pressure, rendering it unattractive for applications at the wellsite.
    • Summary
  • According to some embodiments, systems and methods of analyzing rock samples from an unconventional hydrocarbon reservoir are described. The methods include performing gas sorption on a sample of rock from the unconventional hydrocarbon reservoir; and determining at least one characteristic associated with the sample based at least in part on the gas sorption.
  • According to some other embodiments, a method for determining a location of a fracture zone for stimulation in a shale reservoir is described. The method includes obtaining a shale sample from the shale reservoir and grinding and sieving the shale sample. The method further includes degassing the shale sample and performing gas sorption on the shale sample from the shale reservoir to determine at least one characteristic associated with the shale sample based at least in part on the gas sorption. Finally, the method includes determining the fracture zone from the one or more characteristics.
  • According to some embodiments, a method of analyzing a rock sample is described. The method includes obtaining the rock sample and measuring a pore volume or porosity on a plurality of particles of varying sizes from the rock sample. The method further includes, comparing the measured pore volume or porosity to determine a correlation length or pore connectivity in the rock sample.
  • As used herein the term “unconventional” reservoir includes reservoirs having an unconventional microstructure, such as having submicron pore size, and/or substantial amounts of primary organic matter such as kerogen. Examples of unconventional reservoirs include hydrocarbon-bearing shales such as gas shales and oily shales.
  • Further features and advantages of the invention will become more readily apparent from the following detailed description when taken in conjunction with the accompanying drawings.
    • BRIEF DESCRIPTION OF THE DRAWINGS
  • The present invention is further described in the detailed description which follows, in reference to the noted plurality of drawings by way of non-limiting examples of exemplary embodiments of the present invention, in which like reference numerals represent similar parts throughout the several views of the drawings, and wherein:
  • FIG. 1 shows a sampling tool being deployed in a wellbore and an analysis facility, according to some embodiments;
  • FIG. 2 is a plot of As, N2 normalized by total organic content (TOC) versus vitrinite reflectance (VR0) or the bitumen equivalent (VRE);
  • FIG. 3 is a plot of the change in surface area after bitumen extraction and combustion as a function of VR0 or VRE;
  • FIG. 4 is a plot of isotherms for shale of varying maturity;
  • FIG. 5 is a plot of the natural log of VA, normalized by VM from BET analysis, versus the double natural log of P/P0;
  • FIG. 6 shows sorption results using water vapor as the analysis gas, AS.H2O, normalized by the mass of illite+smectite (I+S) in the sample;
  • FIG. 7 is a plot of permeability, k, versus the difference between the pore volume measured on two different particle sizes; and
  • FIG. 8 is a flow chart of an embodiment of the subject disclosure.
    •  DETAILED DESCRIPTION
  • The particulars shown herein are by way of example and for purposes of illustrative discussion of the embodiments of the present invention only and are presented in the cause of providing what is believed to be the most useful and readily understood description of the principles and conceptual aspects of the present invention. In this regard, no attempt is made to show structural details of the present invention in more detail than is necessary for the fundamental understanding of the present invention, the description taken with the drawings making apparent to those skilled in the art how the several forms of the present invention may be embodied in practice. Further, like reference numbers and designations in the various drawings indicate like elements.
  • In embodiments, methods are disclosed for accessing the quality of a shale reservoir. The information derived from the methods disclosed may be utilized for completion design. Stimulation is the most costly process (in general at least 50% of a well cost) utilized for gas shale plays so gaining insight into the location of the best zones for stimulation is important. These include zones with the most mature organic matter as mature organic matter contains oil-wet pores which may represent unimpeded pathways to the fractures formed during stimulation. Generally, the gas in shale is in one of two states which are (a) tightly absorbed on the rock surface (liquid-like) or (b) free in the body of the pores (gas-like). Laboratory measurements indicate that initially the free-gas is produced and the adsorbed gas is produced when the reservoir pressure is reduced. It is therefore an advantage to identify mature zones with the greatest proportion of free gas.
  • The subject disclosure utilizes conventional gas sorption for characterizing the microstructure of gas shale. Gas sorption is applicable for characterizing nanoporosity (<200 nm) and therefore is an ideal technique for characterizing the microstructure of gas shale. In a non-limiting example a gas sorption study is performed at the boiling point of the analysis gas (−196° C. for nitrogen). Under these conditions the saturation vapor pressure, P0, of the analysis gas is equivalent to ambient pressure (˜0.1 MPa). Other analysis gases which may be used include argon, xenon, neon or any molecular or atomic probe known to be useful in adsorption studies. The disclosure implicitly covers the normal thermodynamic conditions with respect to range of temperature and pressure, regularly imposed for the various analysis gases. Experiments consist of measuring the volume of adsorbed gas, VA, from low (˜1 kPa) to ambient pressure at constant temperature (commonly referred to as an isotherm). At the low end of the pressure range, gas molecules stack in multiple layers on the internal surfaces of the shale, permitting determination of surface area AS. At higher pressure, liquid nitrogen condenses in the fine pores of the shale, permitting determination of the pore size distribution for pores smaller than 200 nm, and the pore volume VP.
  • FIG. 1 shows a sampling tool being deployed in a wellbore and an analysis facility, according to some embodiments. In a non-limiting example the sampling tool is a core sampling tool. Wireline truck 110 is deploying wireline cable 112 into well 130 via well head 120. Wireline tool 114 is disposed on the end of the cable 112 in an unconventional subterranean formation 100. According to some embodiments, formation 100 is an unconventional reservoir, such as a hydrocarbon bearing shale reservoir. Tool 114 includes a sampling tool as shown, in a non-limiting example a core sampling tool. Although a wireline sampling tool is shown, according to other embodiments, other types of sampling tools are used such as while drilling and/or coiled tubing conveyed tools. Samples 132 are from an unconventional rock formation 100 and are retrieved at the surface from the tool 114 and transported to an analysis facility 160. Note that the analysis facility 160 can be located at the wellsite (which can be onshore or offshore) or it can be located remotely from the wellsite. Facility 160 includes a gas sorption apparatus (152), one or more central processing units 140, storage system 144, communications and input/output modules 140, a user display 146 and a user input system 148. Input/output modules 140 include modules to communicate with and control the gas sorption machine (152). Facility 160 may also include a helium pycnometry (He-pyc), mercury intrusion porosimetry (MICP) or other apparatuses as known to those skilled in the art for characterizing rock samples.
  • By using conventional gas sorption it can be shown that the microstructural characteristics of shale are dependent on the shale maturity. The correlation between the microstructure and maturity is due to pores which form in the organic matter as hydrocarbons are generated and ultimately expelled. These pores are characterized by length scales that range from 10−9-10−7 m. The high surface to volume ratio of these nanopores is responsible for an increase in specific area, As, with maturity. An increase in maturity leads to a concomitant increase in porosity and surface area due to nanometric pores in the total organic content. These pores increase pore connectivity and are formed as a result of hydrocarbon generation and expulsion. This increase in specific area, As, with maturity is depicted in FIG. 2.
  • FIG. 2 is a plot of As, N2 normalized by TOC versus virtinite reflectance (VR0) or the bitumen equivalent (VRE). Thermal maturity is indicated by vitrinite reflectance, VR0, or the equivalent determined on bitumen, VRE (for higher maturity), which are acquired by standard organic petrology. A similar related metric commonly used is Tmax. The plot shows that the surface area increases with maturity. The open circles are the raw data, while the filled surfaces show the surface area attributed to organic matter alone. Since As, N2 represents that due to organic and clay (open symbols), there is considerable scatter in the data, especially in the samples that are moderately mature. Therefore, As, N2/TOC is reduced by an amount attributable to clay (filled symbols). The proportion of As, N2 attributable to clay is given by the XRD mineralogy, the relationship displayed in Equation 1 (free exponent) below for the proportion of clay accessible to the analysis gas, and G the ratio of the surface area of pure clay as measured by nitrogen to that measured by water.
  • A S , N 2 ( no clay ) = A S , N 2 ( with clay ) - B * ( I + S ) * G * ( TOC 100 ) - 0.998 Equation 1
  • As can be seen in FIG. 2 the specific surface area of the organic matter (filled symbols) increases markedly over the spectrum of maturity investigated. Mature (VR0>2) TOC has a specific surface area of 400-600 m2/g which is consistent with previously published values from artificially matured shale, and is similar to that of activated carbon. The plot also shows the average As, N2 from samples with 1<VR0 or E (immature), 1 <VR0 or E<1.5 (moderate) and 1.5<VR0 or E (mature). The average value for mature shale is more than an order of magnitude greater than that for immature shale, and nearly twice that of moderately mature shale.
  • Further, it has been found that micropores (pore size <2 nm as defined by the International Union of Pure and Applied Chemistry (IUPAC)) are found in moderate to highly mature shale as depicted in FIG. 3. FIG. 3 is a plot showing the change in surface area after bitumen extraction, and combustion as a function of VR0 or VRE. Also shown is the volume of micropores Vμ determined using, in a non-limiting example, the t-plot analysis on gas sorption data from the as received material. Vμ could also be measured directly by performing gas sorption with carbon dioxide. The high specific surface area of the kerogen is due to the increase in nanoporosity, including Vμ, with maturity. These extremely small micropores (<2 nm) absorb a high volume of gas at low pressure and this may be used for identifying mature zones by simply measuring the amount of gas adsorbed at very low pressure (P/P0<0.05), as depicted in FIG. 4. FIG. 4 is a graph of isotherms for shale of varying maturity. The presence of micropores elevates VA(P/P0→0) in the mature shale. This observation provides a rapid maturity index. The horizontal dashed line in FIG. 4 is an example of a lower bound on VA(P/P0→0) above which a sample would be considered mature.
  • Further, adsorbed gas (VA) is measured for a plurality of increments, in a non-limiting example, seven increments of pressure up to P0, where the last data is taken at P/P0˜0.995. VA(P/P0→0) (e.g. maturity), AS (surface area), and VP (pore volume) can then be determined from this information. The latter two quantities may be combined to provide an estimate of the pore volume available to free gas, VF, given the thickness of the layer of natural gas adsorbed in practice, t˜1 nm; VF=VP−ASt. Moreover, these quantities can be utilized to estimate the total gas in place or hydrogen index given the reservoir pressure and the associated gas density profile in the pore (e.g. higher density in the adsorbed phase near the pore wall, decaying to bulk density in the center of the pore.). In a non-limiting example the Brunauer, Emmet and Teller analysis (BET) is utilized to determine, As given VA at several pressures below P/P0<0.3. BET is a statistical analysis of the buildup of gas multilayer's on a solid surface. The primary result is an expression for the number of gas molecules in a monolayer. In a non-limiting example, the t-plot technique is used to determine the volume of micropores, Vμ, given VA over the entire range of pressure. This approach utilizes a transformation of variables between P and the thickness of adsorbed gas, tA, to deduce Vμ, where the transformation of variables is determined on a chemically similar non-porous sample. The idea is that if Vμ=0, a plot of VA versus tA passes through the origin.
  • Shale is a fractal material, and as such access to the microstructure is described by percolation concepts. FIG. 5A is a plot of the natural log of VA, normalized by VM from BET analysis, versus the double natural log of P/P0. The plot is linear over nearly two orders of magnitude in pressure, which corresponds to nearly two orders of magnitude in pore size. The slope of the line, m, gives the fractal dimension DF=m+3. DF increases with maturity, indicating that at a given characteristic length, a greater pore volume is accessible with increasing maturity.
  • Shales are distinguished from conventional reservoir rock by their more abundant clay and OM content. As, H2O normalized by the mass of illite and smectite (I+S) exhibits a power law dependence on TOC. FIG. 6 plots As, H2O normalized by the mass of illite and smectite (I+S) versus TOC. The data corresponds to a number of samples across the spectrum of maturity. As can be seen in FIG. 6 there are two power law fits to the As, H2O/(I+S) data, one with the exponent as a free parameter (gray) and the other with the exponent fixed (black). Water vapor is used as it will not condense on organic surfaces. Therefore, As, H2O inferred from these measurements is indicative of that due to exposed inorganic surfaces. In shale, inorganic surface area is dominated by clay, where the specific surface area of clay as measured by water is As, H2O≈300-400 m2/g. Normalizing As, H2O by the total clay content of the sample (I+S=Illite+Smectite) as determined by x-ray diffraction, the resulting value is strongly correlated to TOC, and that in the limit of TOC→0, the value As, H2O/(I+S) is consistent with the specific surface area of clay. These results are very useful because they yield several key insights. First, the power law fit to the data in FIG. 6 has an exponent that is similar to that for 3D conductivity in the system. This is important because it is indicative that percolation concepts apply for gauging access to the porosity. Moreover, the TOC in the mature systems contributes numerous pores which provide access to the porosity which yields access to the inorganic surfaces. Applying percolation concepts, we expect the correlation length (χ) to vary with a similar power law, χ=χ0(p−pc)−υ, where υ=0.88 for 3D percolation, χ0 is a fundamental length, like the characteristic dimension of the organic matter, and p is the volume fraction of a component in the sample, in this case TOC. Since we use a consistent particle size in our experiments, we conclude that the relative correlation length decreases with TOC content, or a reduction in maturity. For the oilfield this information is useful because smaller x necessitates greater stimulation to access the natural resource trapped inside the rock. Of course, the correlation length can also be determined by running the gas sorption measurement on a single sample by varying the particle size used in the analysis. This will give the discrete range of particle sizes in which χ lies (e.g., bounded below by the largest particle size which yields a consistent AS, and above by the first particle size which yields lower AS).
  • Embodiments of the subject disclosure comprise techniques for characterizing shale using gas sorption. The permeability of shale is very low (10−19-10−20 m2) so the particle size is reduced (0.001-10 mm) so the gas molecules may probe the pore network rapidly. To achieve low pressure, samples are degassed before running experiments.
  • According to some embodiments a shale sample is retrieved and the sample is grinded and then sieved ensuring all material passes through a sieve with 150 μm openings. According to some embodiments the sample size is on the order of 1 g. The sample is then loaded into an analysis tube and the mass of the sample is noted. The analysis tube is then loaded onto an analyzer to degas. Gas sorption analysis is then performed. VA is noted for a plurality of P/P0 in the pressure range 0.001<P/P0<0.995. This includes a point very near the upper bound. VA is also noted in the limit, P/P0→0. The sample is determined to be mature or not by comparing to lower bound for mature samples: VA(P/P0→0)=3-4 cc/g, see also FIG. 4. If the sample is mature, a t-plot analysis may be used to determine As and the volume of micropores, otherwise a BET analysis may be used to determine As. The pore volume (Vp) may be noted from P/P0 and the pore volume available to free gas (VF), and the gas in place or hydrogen index may also be determined as indicated above.
  • According to some embodiments the particle size may be varied to quantify the associated effect on access to the porosity. If the particle size is varied over a broad enough range, the data will provide insight into the correlation length of the shale and the total pore volume or porosity.
  • According to some embodiments the particle size may be varied and helium pycnometry (He-pyc), mercury intrusion porosimetry (MICP), or other techniques as known to those skilled in the art for measuring pore volume or porosity are used to yield information and also to determine the correlation length of shale. He-pyc determines the volume of solids in a sample of known mass, by allowing a known amount of Helium gas to expand into a container of known volume that holds the sample of interest. The volume of solids is given by the equilibrium pressure using the ideal gas law (corrected for non-ideality). The skeletal density is taken as the ratio between the sample mass and the volume of solids that make up the porous sample. MICP forces mercury into a porous body by incrementally increasing the pressure applied to the mercury. The size of pores the mercury can access is inversely proportional to the pressure. By monitoring the volume of mercury injected with each increase in pressure, the bulk density can be determined from the sample mass, and the sample volume. The sample volume is the difference between the volume of the sample holder and the volume of the mercury required to surround but not invade the sample. Also the pore size distribution can be determined from the volume of mercury injected at each incremental increase in pressure, and the well known relationship between the pore size and pressure as given by the Washburn equation. The pore volume is determined from the total volume of mercury injected at the highest pressure, and the skeletal density is given by the mass of the sample and the volume of solids that make up the porous medium. The latter quantity is taken as the difference between the volume of the sample holder, and the total amount of mercury required to surround the sample, and fill the accessible pores.
  • In a non-limiting example, Bulk density, ρB is measured by weighing in air and the volume is determined by measuring the displacement of a liquid. Gas permeability is determined using a modified pycnometer technique by analyzing the kinetics of pressure equilibration at a plurality of pressures. With this information, the liquid permeability, k, is inferred from the y-intercept of a linear fit to the correlation with the experimental pressure (See Klinkenberg, “The permeability of porous media to liquids and gases,” Drill, Proc. API (1941), 200-213). The skeletal density ρS is determined from the equilibrium pressure during the gas permeability measurements utilizing the ideal gas law.
  • FIG. 7 is a plot of permeability k versus the difference between the measured and calculated pore volume. The measured pore volume is measured by nitrogen gas sorption, VP,meas, and the calculated pore volume is determined from ρB and ρS, VP,calc:
  • V P , calc = 1 ρ B - 1 ρ S Equation 2
  • The measurement of ρS is on a material that has a characteristic dimension that is at least one order of magnitude greater than that used in the gas sorption experiment, in a non-limiting example 10−3 m vs. 10−4 m. The difference in size of the particles can be visually seen in FIG. 7 with the particles (701) used in the measured gas sorption being smaller than the particles (703) used where ρS is measured. As a result of this difference, we measure the same or more VP,meas in the gas sorption measurement than we determine from the densities using equation (2). This difference may be attributed to the fact that the gas used in the ρS measurement is not able to access all of the porosity. For those samples where VP,meas>VP,calc, correlation length (χ) is smaller than the particle size used in the ρS measurement (10−3). On the other hand, when VP,meas≈VP,calc, correlation length (χ) is greater than or equivalent to the particle size used in the ρS measurement (10−4). As can be seen in FIG. 7 the mature samples have higher k than the moderate or lower maturity systems over the range in VP,meas-VP,calc. The exception occurs when VP,meas≈VP,calc, in this case the moderate maturity systems have k that is similar to the highest k of the mature systems. This observation is consistent with the expectation of percolation theory where the highest permeabilities are correlated to the fact that the measured porosity/pore volume does not increase when the particle size is reduced. In any case, k decreases as VP,meas becomes greater than VP,calc. A simple permeability model, like the Carmen-Kozeny equation (3) may be used to better understand the relationship between the microstructure and k further:
  • k = φ r H 2 κ Equation 3
  • Where Ø=VPØB is the porosity, rH is the hydraulic radius, and k is the kozeny constant, which is roughly equivalent to the tortuosity. It has been shown that rH decreases by a factor of 2 while VP increases by the same factor over the range of maturity investigated here. According to equation 3 these microstructural variations suggests that if k is constant, k should decrease by a factor of 2 as maturity increases. FIG. 7 indicates that k increases by a factor of 2-4 with increasing maturity. So κ drops by as much as factor of 4-8 with an increase in maturity. This observation is also consistent with the predictions of percolation theory. This decrease in κ is due to the opening of pores in the organic matter (OM). In terms of reservoir quality, the ideal systems are those with large correlation length (χ) as these systems require less stimulation to fully tap the natural hydrocarbon reserve. At the same times these systems also have more attractive transport properties.
  • The data shows that the gas sorption measurement sees more pore volume than that determined from the densities. Also, the permeability of moderate and highly mature samples increases as the pore volume determined by the independent technique comes into agreement.
  • According to further embodiments any combination of the results from gas sorption, and similar results from other characterization techniques e.g., He-pyc or MIP using a plurality of varying particle sizes (0.001-50 mm) will yield information which may be used to determine the correlation length of shale.
  • As discussed above, the microstructure of mature shale has nanometer sized pores in the kerogen network. Therefore, since the surface area is expected to drop upon removal of organic matter from the shale matrix, a comparison of As measured on samples in the native and organic free (via high temperature combustion) states is an alternative means for identifying mature zones (See FIG. 2 combustion data).
  • In further embodiments the samples are analyzed before and after combustion or bitumen extraction. The shale samples are homogenized and split into two samples, a first and a second sample. The two samples are grinded and then sieved ensuring all material passes through a sieve with 150 μm openings. According to some embodiments the sample sizes are on the order of 1 g. A first sample is then combusted at 450° C. in air for at least 16 hours to remove organic matter or subjected to solvent reflux for several days to remove bitumen. In order to run gas sorption analysis the first sample is dried. A second sample is placed in a vacuum oven overnight. The total organic content (TOC) is estimated in the first sample from mass loss due to combustion or bitumen removal. The samples are then loaded into a clean sample tube and the mass of each of the samples is noted. Gas sorption analysis is then performed on both samples noting the VA of P/P0 in the pressure range 0.01<P/P0<0.995, including a point at the upper bound. The surface area (As) is determined by the BET method on both samples and the results are compared. If the surface area (As) is reduced by 200 m2/g/TOC after combustion the sample may be considered mature (see FIG. 3—filled symbols). If the surface area (As) increases markedly upon bitumen extraction, the sample may be considered moderately mature (see FIG. 3—open symbols). The gas sorption data from the native-state sample is analyzed using the t-plot method to check for microporosity, if the sample contains micropores, the sample is mature. The pore volume is determined from the volume of adsorbed gas (VA) at the highest pressure and the pore volume available to free gas (VF). (VF) may also be determined as indicated above.
  • FIG. 8 depicts a flow chart of embodiments of the subject disclosure. In an embodiment the shale sample (801) is grinded (803) and then sieved (805) ensuring all material passes through a sieve with a predefined opening, in a non-limiting example the opening is 150 μm. The shale sample (801) is then degassed (807) in preparation for gas sorption analysis (809). From gas sorption analysis (809) a plurality of microstructural characteristics (823) are determined which include pore volume, surface area, pore size distribution and fractal dimension and correlation length or relative pore connectivity.
  • In a further embodiment the shale sample (801) is grinded (803) and then sieved (805) and the sample is then analyzed using helium pycnometry (He-pyc) and bulk density. The pore volume (813) can be determined with this approach.
  • In a further embodiment the shale sample (801) is grinded (803) and then sieved (805) and the sample is then analyzed using MICP (817). From MICP (817) a plurality of physical and microstructural characteristics (823) are determined which include bulk density, skeletal density, pore volume, surface area and pore size distribution.
  • The pore volume or porosity may also be measured on varying sized particles in the rock sample using any of the techniques disclosed or any techniques known to those skilled in the art for measuring the pore volume and the porosity in a rock sample and the results are compared to determine the correlation length or relative pore connectivity in the system.
  • Comparing any output 813, 815, 819, 821 or 823 yields information regarding pore connectivity, more specifically the dimension over which the porosity is connected in 3D. When mercury intrusion porosimetry (MICP) is used the size of pores controlling transport and the pore morphology may also be determined.
  • Whereas many alterations and modifications of the present invention will no doubt become apparent to a person of ordinary skill in the art after having read the foregoing description, it is to be understood that the particular embodiments shown and described by way of illustration are in no way intended to be considered limiting. Further, the invention has been described with reference to particular preferred embodiments, but variations within the spirit and scope of the invention will occur to those skilled in the art. It is noted that the foregoing examples have been provided merely for the purpose of explanation and are in no way to be construed as limiting of the present invention. While the present invention has been described with reference to exemplary embodiments, it is understood that the words, which have been used herein, are words of description and illustration, rather than words of limitation. Changes may be made, within the purview of the appended claims, as presently stated and as amended, without departing from the scope and spirit of the present invention in its aspects. Although the present invention has been described herein with reference to particular means, materials and embodiments, the present invention is not intended to be limited to the particulars disclosed herein; rather, the present invention extends to all functionally equivalent structures, methods and uses, such as are within the scope of the appended claims.

Claims (14)

What is claimed is:
1. A method of analyzing rock samples from an unconventional hydrocarbon reservoir comprising:
grinding and sieving a sample of rock from the unconventional hydrocarbon reservoir;
degassing the sample of rock;
performing gas sorption on the sample of rock; and
determining at least one characteristic associated with the sample based at least in part on the gas sorption.
2. The method of claim 1, wherein the sample is at least one of a core sample or a cutting sample from the unconventional hydrocarbon reservoir.
3. The method of claim 1, wherein the at least one characteristic includes a measurement of one or more types selected from a group consisting of: total organic content, surface area, thermal maturity, pore size distribution, pore volume, volume of micropores, free gas capacity, permeability, fractal dimension, and correlation length.
4. The method of claim 1, wherein the sample includes a plurality of particles.
5. The method of claim 4, wherein a size of the plurality of particles is selected from a range of sizes.
6. The method of claim 1, further comprising:
calculating a bulk density, a skeletal density, a pore volume, a surface area or a pore size distribution of the sample using mercury intrusion porosimetry.
7. The method of claim 1, further comprising:
calculating a skeletal density of the sample;
calculating a bulk density of the sample; and
determining a pore volume of the sample from the calculated skeletal density and the bulk density.
8. The method of claim 7, further comprising:
using pycnometry to calculate the skeletal density.
9. The method of claim 1, further comprising:
extracting bitumen from the sample using a solvent;
extracting kerogen from the sample using combustion; and
calculating a pore volume, surface area and pore size distribution of the sample.
10. The method of claim 3, wherein the at least one characteristic is used to determine a location for initiating a fracture in the unconventional hydrocarbon reservoir.
11. The method of claim 1, wherein the unconventional hydrocarbon reservoir includes hydrocarbon-bearing shales.
12. The method of claim 3, wherein a size of the micropores is smaller than 2 nm
13. The method of claim 1, wherein performing gas sorption on the sample uses one of nitrogen, water, carbon dioxide, argon, xenon, neon, or combinations thereof.
14. The method of claim 1, further comprising:
extracting the sample from an unconventional hydrocarbon reservoir at a plurality of depths in the reservoir.
US14/475,436 2011-01-27 2014-09-02 Gas Sorption Analysis Of Unconventional Rock Samples Abandoned US20140366621A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US14/475,436 US20140366621A1 (en) 2011-01-27 2014-09-02 Gas Sorption Analysis Of Unconventional Rock Samples

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201161436819P 2011-01-27 2011-01-27
US13/359,121 US8881587B2 (en) 2011-01-27 2012-01-26 Gas sorption analysis of unconventional rock samples
US14/475,436 US20140366621A1 (en) 2011-01-27 2014-09-02 Gas Sorption Analysis Of Unconventional Rock Samples

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
US13/359,121 Continuation US8881587B2 (en) 2011-01-27 2012-01-26 Gas sorption analysis of unconventional rock samples

Publications (1)

Publication Number Publication Date
US20140366621A1 true US20140366621A1 (en) 2014-12-18

Family

ID=46576210

Family Applications (2)

Application Number Title Priority Date Filing Date
US13/359,121 Active 2033-03-18 US8881587B2 (en) 2011-01-27 2012-01-26 Gas sorption analysis of unconventional rock samples
US14/475,436 Abandoned US20140366621A1 (en) 2011-01-27 2014-09-02 Gas Sorption Analysis Of Unconventional Rock Samples

Family Applications Before (1)

Application Number Title Priority Date Filing Date
US13/359,121 Active 2033-03-18 US8881587B2 (en) 2011-01-27 2012-01-26 Gas sorption analysis of unconventional rock samples

Country Status (5)

Country Link
US (2) US8881587B2 (en)
CN (1) CN103339488B (en)
AU (1) AU2012211109B2 (en)
CA (1) CA2822696A1 (en)
WO (1) WO2012103479A2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016130314A1 (en) 2015-02-11 2016-08-18 Schlumberger Technology Corporation Porosimetry transition region adjustment
WO2016160076A1 (en) * 2015-04-01 2016-10-06 Schlumberger Canada Limited Determination of free volume of a rock sample using high pressure adsorption data
KR101826132B1 (en) 2016-11-01 2018-03-22 한국지질자원연구원 Prediction apparatus and the method of desorbed gas volume of shale using geophysical logging data

Families Citing this family (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013066549A1 (en) * 2011-10-31 2013-05-10 Baker Hughes Incorporated Hydrocarbon determination in unconventional shale
US20130262069A1 (en) * 2012-03-29 2013-10-03 Platte River Associates, Inc. Targeted site selection within shale gas basins
US8857243B2 (en) 2012-04-13 2014-10-14 Schlumberger Technology Corporation Methods of measuring porosity on unconventional rock samples
US8967249B2 (en) 2012-04-13 2015-03-03 Schlumberger Technology Corporation Reservoir and completion quality assessment in unconventional (shale gas) wells without logs or core
US20140088878A1 (en) * 2012-09-27 2014-03-27 Jinhong Chen Isotherm and gas-in-place estimation considering capillary condensation in shale gas reservoir
BR112015027750A2 (en) * 2013-05-02 2017-07-25 Schlumberger Technology Bv land sample analysis method, land sample analysis system, and land sample analysis system
BR112015032469A2 (en) 2013-06-25 2017-07-25 Geoservices Equipements Sas method of quantifying bitumen quantities in a soil sample, and system for quantifying bitumen quantities in a soil sample.
US10385259B2 (en) 2013-08-07 2019-08-20 Schlumberger Technology Corporation Method for removing bitumen to enhance formation permeability
CN103822849B (en) * 2013-11-22 2016-06-29 中国地质科学院地质力学研究所 The scene of a kind of non-conventional gas resolves instrument system automatically
CN103900928B (en) * 2014-04-22 2014-11-19 中国地质大学（武汉） Full-automatic shale gas content testing analyzer
CN104198228B (en) * 2014-08-18 2017-01-11 中国石油天然气股份有限公司 Nanoscale unconventional reservoir rock adsorption gas extracting and hydrocarbon analysis method
CN105628574B (en) * 2014-10-29 2019-04-09 中国石油天然气股份有限公司 The measurement method of rock sample porosity
CN104964993A (en) * 2015-06-29 2015-10-07 中国科学院青海盐湖研究所 X ray diffraction testing method for lake sediment lamina thickness single lamina samples
CN105136643B (en) * 2015-08-25 2018-07-17 成都理工大学 The implementation method of one-way flow dynamic pressure mercury system under effective overburden pressure
CN105136646A (en) * 2015-09-15 2015-12-09 东北大学 Permeability evolution experimental device and method considering shale adsorption and fracturing process
CN105445159B (en) * 2015-11-11 2021-08-13 成都理工大学 Method for obtaining pore size distribution curve and sample specific surface area
CN105486621B (en) * 2015-12-29 2018-05-04 中国石油天然气股份有限公司 A kind of shale pore-size distribution test method
CN107543775B (en) * 2017-05-12 2020-02-11 河海大学 Method for determining rockfill material filling index and detecting field filling quality based on fractal theory
CN107228816B (en) * 2017-06-01 2018-11-13 中国石油大学(华东) Different type pore evolution evaluation method in a kind of mud shale
CN108982287B (en) * 2017-06-02 2020-09-08 中国石油化工股份有限公司 Method and device for determining free adsorption ratio of shale gas reservoir
CN108829950B (en) * 2018-05-31 2021-04-27 中国科学院力学研究所 Unconventional reservoir permeability evaluation method based on core image
CN108843310B (en) * 2018-06-02 2021-11-12 东北石油大学 Method for determining law of residual oil in water drive of matrix-high permeability strip by applying fractal theory
CN108843311B (en) * 2018-06-02 2021-08-03 东北石油大学 Method for determining distribution and migration rule of water-drive residual oil by applying fractal theory
CN108678738B (en) * 2018-06-02 2021-11-12 东北石油大学 Method for determining distribution and migration rule of matrix-hypertonic strip functional polymer flooding residual oil by applying fractal theory
CN108825222B (en) * 2018-06-02 2021-08-03 东北石油大学 Method for determining distribution and migration rule of functional polymer flooding residual oil by applying fractal theory
CN108952696B (en) * 2018-06-02 2021-08-03 东北石油大学 Method for determining chemical flooding residual oil distribution and migration rule by applying fractal theory
CN108444895B (en) * 2018-06-14 2020-08-25 长安大学 Efficient method for obtaining unsaturated permeability parameters of loess body
CN108846218B (en) * 2018-06-22 2020-07-28 中国石油大学(北京) Method and system for predicting internal specific surface areas of different types of pores of shale reservoir
CN109635410B (en) * 2018-12-06 2021-08-17 中国华能集团清洁能源技术研究院有限公司 Permeability discrete element simulation method based on pore system
CN110364227A (en) * 2019-07-18 2019-10-22 中国石油大学(华东) A kind of mud shale supercritical methane isotherm adsorption model based on variable density
CN111007233B (en) * 2019-12-25 2022-03-11 西南石油大学 Method for analyzing movement behavior of methane-carbon dioxide in micro pores of shale
CN111287738B (en) * 2020-02-10 2023-06-30 中国石油天然气股份有限公司 Shale gas data determination method and device for shale gas well
CN111537416B (en) * 2020-04-10 2021-10-22 中国地质大学（北京） Method for judging change of nanoscale pore structure and application
CN113324889B (en) * 2021-06-01 2022-05-06 西南石油大学 Device for evaluating shale oil in-situ pyrolysis exploitation displacement efficiency and testing method
CN115598023B (en) * 2022-11-08 2023-04-07 昆明理工大学 Method for determining three-dimensional tortuosity of loose broken rock-soil mass
CN117388113B (en) * 2023-08-18 2024-03-19 中国地质大学（北京） Shale gas associated helium content acquisition method, terminal and medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8522599B2 (en) * 2008-04-17 2013-09-03 Isotech Laboratories, Inc. Low dead-volume core-degassing apparatus
US20130270011A1 (en) * 2012-04-13 2013-10-17 Schlumberger Technology Corporation Reservoir and completion quality assessment in unconventional (shale gas) wells without logs or core
US8714004B2 (en) * 2010-01-13 2014-05-06 Santos Ltd. Measuring gas content of unconventional reservoir rocks
US8857243B2 (en) * 2012-04-13 2014-10-14 Schlumberger Technology Corporation Methods of measuring porosity on unconventional rock samples

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3349625A (en) * 1964-01-10 1967-10-31 Gulf Research Development Co Adsorption measuring apparatus and method
US4489593A (en) * 1982-09-09 1984-12-25 Omicron Technology Corporation Method and apparatus for determining the amount of gas adsorbed or desorbed from a solid
US4856320A (en) * 1988-09-08 1989-08-15 Universite Du Quebec A Trois-Rivieres, Societe Quebecoise D'initiatives Petrolieres And Gaz Metropolitain Inc. Method and apparatus for measuring physical adsorption of gases based on dielectric measurements
JP2881242B2 (en) * 1990-02-27 1999-04-12 富士シリシア化学株式会社 Device for measuring adsorption and desorption and method for measuring adsorption and desorption
FR2661990B1 (en) * 1990-05-09 1992-07-24 Inst Francais Du Petrole METHOD AND DEVICE FOR MEASURING THE ADSORPTION AND DESORPTION OF A GAS ABSORBED BY A SOLID SAMPLE AS WELL AS ITS USE.
US5342580A (en) * 1990-04-17 1994-08-30 Alan Brenner Apparatus and method for measuring the amount of gas adsorbed on or desorbed from a solid and reactions of a gas with a solid
US5408864A (en) * 1993-06-18 1995-04-25 Coulter Corporation Method of determining the amount of gas adsorbed or desorbed from a solid
JP3243079B2 (en) * 1993-09-16 2002-01-07 株式会社日立製作所 Atmospheric pressure heated gas desorption equipment
AU2001251019A1 (en) * 2000-03-27 2001-10-08 Peter J. Ortoleva Method for simulation of enhanced fracture detection in sedimentary basins
WO2002034931A2 (en) 2000-10-26 2002-05-02 Guyer Joe E Method of generating and recovering gas from subsurface formations of coal, carbonaceous shale and organic-rich shales
CN2469437Y (en) * 2001-03-30 2002-01-02 西南石油学院 Gas adsorption testing device for core with high temp. and high pressure
AU2002353888B1 (en) 2001-10-24 2008-03-13 Shell Internationale Research Maatschappij B.V. In situ thermal processing of a hydrocarbon containing formation using a natural distributed combustor
US6595036B1 (en) * 2002-02-27 2003-07-22 Bel Japan, Inc. Method and apparatus for measuring amount of gas adsorption
US7429358B1 (en) * 2002-05-20 2008-09-30 Hy-Energy, Llc Method and apparatus for measuring gas sorption and desorption properties of materials
CN1200262C (en) * 2003-01-10 2005-05-04 清华大学 Method and equipment for measuring amount of gas adsorbed by powder
FR2853071B1 (en) * 2003-03-26 2005-05-06 Inst Francais Du Petrole METHOD AND DEVICE FOR EVALUATING PHYSICAL PARAMETERS OF A UNDERGROUND DEPOSIT FROM ROCK DEBRIS WHICH ARE TAKEN THEREFROM
JP3846800B2 (en) * 2003-08-18 2006-11-15 株式会社リガク Evolved gas analysis method and apparatus
JP4162138B2 (en) * 2003-10-27 2008-10-08 株式会社リガク Thermal desorption gas analyzer
US20080162056A1 (en) * 2006-12-29 2008-07-03 Keith Howarth Greaves Method and apparatus for determination of gas in place
CN100575921C (en) * 2007-09-28 2009-12-30 北京科技大学 A kind of method for detecting iron ore air hole characteristic
US8573298B2 (en) 2008-04-07 2013-11-05 Baker Hughes Incorporated Method for petrophysical evaluation of shale gas reservoirs
EP2116690A1 (en) * 2008-04-09 2009-11-11 Bp Exploration Operating Company Limited Geochemical surveillance of gas production from tight gas fields
WO2011060494A1 (en) * 2009-11-19 2011-05-26 Ian Gray System for analysing gas from strata being drilled under high mud flows

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8522599B2 (en) * 2008-04-17 2013-09-03 Isotech Laboratories, Inc. Low dead-volume core-degassing apparatus
US8714004B2 (en) * 2010-01-13 2014-05-06 Santos Ltd. Measuring gas content of unconventional reservoir rocks
US20130270011A1 (en) * 2012-04-13 2013-10-17 Schlumberger Technology Corporation Reservoir and completion quality assessment in unconventional (shale gas) wells without logs or core
US8857243B2 (en) * 2012-04-13 2014-10-14 Schlumberger Technology Corporation Methods of measuring porosity on unconventional rock samples

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016130314A1 (en) 2015-02-11 2016-08-18 Schlumberger Technology Corporation Porosimetry transition region adjustment
US10371619B2 (en) 2015-02-11 2019-08-06 Schlumberger Technology Corporation Porosimetry transition region adjustment
WO2016160076A1 (en) * 2015-04-01 2016-10-06 Schlumberger Canada Limited Determination of free volume of a rock sample using high pressure adsorption data
US10156507B2 (en) 2015-04-01 2018-12-18 Schlumberger Technology Corporation Determination of free volume of a rock sample using high pressure adsorption data
KR101826132B1 (en) 2016-11-01 2018-03-22 한국지질자원연구원 Prediction apparatus and the method of desorbed gas volume of shale using geophysical logging data
WO2018084394A1 (en) * 2016-11-01 2018-05-11 한국지질자원연구원 Apparatus for predicting amount of desorption gas in shale gas layer by using geophysical well logging data analysis and method therefor

Also Published As

Publication number Publication date
CN103339488A (en) 2013-10-02
AU2012211109A1 (en) 2013-07-04
CN103339488B (en) 2015-09-02
AU2012211109B2 (en) 2014-12-18
WO2012103479A3 (en) 2012-11-01
US8881587B2 (en) 2014-11-11
US20120192639A1 (en) 2012-08-02
CA2822696A1 (en) 2012-08-02
WO2012103479A2 (en) 2012-08-02

Similar Documents

Publication Publication Date Title
US8881587B2 (en) Gas sorption analysis of unconventional rock samples
Yang et al. Measurement of the surface diffusion coefficient for adsorbed gas in the fine mesopores and micropores of shale organic matter
Kang et al. Impact of fracturing fluid on multi-scale mass transport in coalbed methane reservoirs
Aljamaan et al. Experimental investigation and Grand Canonical Monte Carlo simulation of gas shale adsorption from the macro to the nano scale
Pang et al. Experimental measurement and analytical estimation of methane absorption in shale kerogen
Wang et al. The effect of moisture on the methane adsorption capacity of shales: A study case in the eastern Qaidam Basin in China
AU2006223089B2 (en) Determination of coal bed natural gas production factors and a system to determine same
Sondergeld et al. Petrophysical considerations in evaluating and producing shale gas resources
Peng et al. Permeability measurements in mudrocks using gas-expansion methods on plug and crushed-rock samples
Yang et al. Review of shale gas sorption and its models
Wang et al. Experimental and numerical investigation of dynamic gas adsorption/desorption–diffusion process in shale
Achang et al. The influence of particle size, microfractures, and pressure decay on measuring the permeability of crushed shale samples
Yang et al. Laboratory study of cryogenic treatment induced pore-scale structural alterations of Illinois coal and their implications on gas sorption and diffusion behaviors
CN105866002A (en) Method for accurately measuring nuclear magnetic resonance porosity of oil-containing shale
Heller Multiscale investigation of fluid transport in gas shales
US10436727B2 (en) Prediction of gas production rates from time-dependent NMR measurements
Gao et al. Water saturation-driven evolution of helium permeability in Carboniferous shale from Qaidam Basin, China: An experimental study
Hatami et al. Poroelastic effects on gas transport mechanisms and influence on apparent permeability in shale
Mahmoud Effect of gas adsorption on the estimation of gas in place (GIP) in conventional and unconventional reservoirs
Yang et al. Modeling water imbibition and penetration in shales: New insights into the retention of fracturing fluids
Zhao et al. Experimental investigation on effect of water film thickness in unsaturated sandstone cores on CO2 transport during geologic storage
Guo et al. Characterizing shale oil occurrence in the Yanchang Formation of the Ordos Basin, assisted by petrophysical and geochemical approaches
Ghanizadeh et al. Controls on liquid hydrocarbon permeability of tight oil and liquid-rich gas reservoirs: Examples from Bakken and Montney formations (Canada)
Zan et al. Experimental investigation of spontaneous water imbibition into methane-saturated shales under different methane pressures
Zamirian New Experimental Approach to Measure Petrophysical Properties of Organic-Rich Shales

Legal Events

Date Code Title Description
STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION