USGS ScienceSearch

USGS · 70121924

Experimental investigation of changes in methane adsorption of bitumen-free Woodford Shale with thermal maturation induced by hydrous pyrolysis

Abstract

This study quantifies the effects of organic-matter (OM) thermal maturity on methane (CH 4 ) sorption, on the basis of five samples that were artificially matured through hydrous pyrolysis achieved by heating samples of immature Woodford Shale under five different time–temperature conditions. CH 4 -sorption isotherms at 35 °C, 50 °C, and 65 °C, and pressures up to 14 MPa on dry, solvent-extracted samples of the artificially matured Woodford Shale were measured. The results showed that CH 4 -sorption capacity, normalized to TOC, varied with thermal maturity, following the trend: maximum oil (367 °C) > oil cracking (400 °C) > maximum bitumen/early oil (333 °C) > early bitumen (300 °C) > immature stage (130 °C). The Langmuir constants for the samples at maximum-oil and oil-cracking stages are larger than the values for the bitumen-forming stages. The total pore volume, determined by N 2 physisorption at 77 K, increases with increased maturation: mesopores, 2–50 nm in width, were created during the thermal conversion of organic-matter and a dramatic increase in porosity appeared when maximum-bitumen and maximum-oil generation stages were reached. A linear relationship between thermal maturity and Brunauer–Emmett–Teller (BET) surface area suggests that the observed increase in CH 4 -sorption capacity may be the result of mesopores produced during OM conversion. No obvious difference is observed in pore-size distribution and pore volume for samples with pores <2.0 nm with the increase of thermal maturity based on CO 2 physisorption at 273 K. The isosteric heat of adsorption and the standard entropy for artificially matured samples ranged from 17.9 kJ mol −1 to 21.9 kJ mol −1 and from −85.4 J mol −1 K −1 to −101.8 J mol −1 K −1 , respectively. These values are similar to the values of immature Woodford kerogen concentrate previously observed, but are larger than naturally matured organic-rich shales. High-temperature hydrous pyrolysis might have induced Lewis acid sites on both organic and mineral surfaces, resulting to some extent, in chemical interactions between the adsorption site and the methane C–H bonds. The formation of abundant mesopores (2–50 nm) within organic matter during organic-matter thermal maturation makes a great contribution to the increase in both BET surface area and pore volume, and a significant increase in 2–6 nm pores occurs at maximum-oil-generation and oil-cracking to gas, ultimately controlling the methane-adsorption capacity. Therefore, consideration of pore-size effects and thermal maturity is very important for gas in place (GIP) prediction in organic-rich shales.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Haiyan Hu, Tongwei Zhang, Jaclyn D. Wiggins-Camacho, Geoffrey S. Ellis, Michael D. Lewan, Xiayong Zhang. 2014. Experimental investigation of changes in methane adsorption of bitumen-free Woodford Shale with thermal maturation induced by hydrous pyrolysis. https://doi.org/10.1016/j.marpetgeo.2014.07.029

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Terrebonne Basin, Gulf of Mexico gas hydrate resource evaluation and 3-D modeling of basin-scale sedimentation, salt tectonics, and hydrate system evolution since the early Miocene

In assessing methane hydrate as a potential transitional energy source, quantification of in-place volumes of gas is a first step—and yet, global, regional, and even local estimates of gas volumes contained within hydrate are highly variable, including within the extensively-studied Gulf of Mexico (GoM) gas hydrate province. Here, we construct the first 3-D basin and hydrocarbon system model of the Terrebonne Basin (Walker Ridge protraction area, northern GoM) to estimate in-place resources contained as gas within hydrate, as well as to provide the first estimate of in-place resources contained within free gas accumulations in the basin. The resultant model-predicted 3-D distributions (and saturations) of hydrate are consistent with independent seismic interpretation and borehole observations. This study reveals both sedimentation-driven and salt diapir-driven gas hydrate recycling produce elevated hydrate saturations in the basin (“recycling” occurs when relative upward shifts of the gas hydrate stability zone (GHSZ) cause hydrate dissociation, generating buoyant free gas that may reaccumulate as higher-saturation hydrate at the new base of GHSZ). An important finding of this study (particularly given that, globally, the source of methane in marine hydrate systems tends to be poorly understood) is that microbial sources can explain most or all hydrate occurrences in Terrebonne Basin. We calculate that ∼32 × 10 9 m 3 (∼1.1 10 12 ft 3 ) of methane gas is trapped within hydrate accumulations throughout Terrebonne Basin, while nearly 2 × 10 9 m 3 (∼67 × 10 9 ft 3 ) of methane occurs as free gas at high saturations within sandy reservoirs beneath the GHSZ and ∼35 × 10 9 m 3 (∼1.2 × 10 12 ft 3 ) occurs as free gas at low saturations within muddy units beneath the GHSZ. In total, our calibrated 3-D model predicts that nearly 70 × 10 9 m 3 (∼2.5 × 10 12 ft 3 ) of microbial methane is trapped as gas hydrate and free gas in the Terrebonne Basin.

Marine and Petroleum Geology

Mapping the resistivity structure of Walker Ridge 313 in the Gulf of Mexico using the marine CSEM method

A marine controlled source electromagnetic (CSEM) campaign was carried out in the Gulf of Mexico to further develop marine electromagnetic techniques in order to aid the detection and mapping of gas hydrate deposits. Marine CSEM methods are used to obtain an electrical resistivity structure of the subsurface which can indicate the type of substance filling the pore space, such as gas hydrates which are more resistive. Results from the Walker Ridge 313 study (WR 313) are presented in this paper and compared with the Gulf of Mexico Gas Hydrate Joint Industry Project II (JIP2) logging while drilling (LWD) results and available seismic data. The hydrate, known to exist within sheeted sand deposits, is mapped as a resistive region in the two dimensional (2D) CSEM inversion models. This is consistent with the JIP2 LWD resistivity results. CSEM inversions that use seismic horizons provide more realistic results compared to the unconstrained inversions by providing sharp boundaries and architectural control on the location of the resistive and conductive regions in the CSEM model. The seismic horizons include: 1) the base of the gas hydrate stability zone (BGHSZ), 2) the top of salt, and 3) the top and bottom of a fine grained marine mud interval with near vertical hydrate filled fractures, to constrain the CSEM inversion model. The top of salt provides improved location for brines, water saturated salt, and resistive salt. Inversions of the CSEM data map the occurrence of a ‘halo’ of conductive brines above salt. The use of the BGHSZ as a constraint on the inversion helps distinguish between free gas and gas hydrate as well as gas hydrate and water saturated sediments.

Louisiana

Comparison of measured versus modeled TOC in the Tuscaloosa marine shale of Southwestern Mississippi, U.S.A.

This study presents a comparison of measured versus modeled total organic carbon (TOC) in the Upper Cretaceous Tuscaloosa marine shale (TMS) of southwestern Mississippi as a case study to evaluate the effects of mineralogy on the TOC estimated from the ΔlogR method. The ΔlogR method is utilized to calculate TOC, which involves baselining sonic transit time and resistivity log curves in a non-source rock section of the formation. In our application, the well log curves were baselined in the upper TMS, which is described as a non-source rock section, and in the lower TMS above the high resistivity zone (HRZ), which is described as having a higher carbonate content. The ΔlogR calculated TOC values from these two baselining approaches show that the lower baseline results in improved agreement between measured TOC and calculated TOC. This improvement is likely to be due to the lower baseline accounting for the increase in resistivity caused by higher carbonate content, in addition to any presence of TOC. The upper baseline, which has a lower carbonate content, does not account for this resistivity increase. Additionally, sample type appears to affect the comparison of measured and ΔlogR calculated TOC. Most of the samples used in this study are legacy cuttings that were not preserved during storage, exposing the high surface area cuttings to increased rates of oxidation, whereas geophysical logs record the rock properties in situ. To account for this oxidation effect, the difference between the medians of the TMS HRZ TOC core and cuttings values was added to each TOC measurement in this study, resulting in a median measured TOC value that is similar to the median of the lower TMS-baselined ΔlogR calculated TOC value. Overall, this study demonstrates that carbonate content and sample type can affect how well measured and ΔlogR-modeled TOC values compare.

Mississippi