Geology Reports⌕ Search
Home/Archive/USGS

USGS · 70015494

Subsurface temperatures and surface heat flow in the Michigan Basin and their relationships to regional subsurface fluid movement

R. Vugrinovich
View original source ↗Export metadata ↓

Abstract

Linear regression of 405 bottomhole temperature (BHT) measurements vs. associated depths from Michigan's Lower Peninsula results in the following equation relating BHT and depth: BHT(??C) = 14.5 + 0.0192 ?? depth(m) Temperature residuals, defined as (BHT measured)-(BHT calculated), were determined for each of the 405 BHT's. Areas of positive temperature residuals correspond to areas of regional groundwater discharge (determined from maps of equipotential surface) while areas of negative temperature residuals correspond to areas of regional groundwater recharge. These relationships are observed in the principal aquifers in rocks of Devonian and Ordovician age and in a portion of the principal aquifer in rocks of Silurian age. There is a similar correspondence between high surface heat flow (determined using the silica geothermometer) and regional groundwater discharge areas and low surface heat flow and regional groundwater recharge areas. Post-Jurassic depositional and tectonic histories suggest that the observed coupling of subsurface temperature and groundwater flow systems may have persisted since Jurassic time. Thus the higher subsurface palaeotemperatures (and palaeogeothermal gradients) indicated by recent studies most likely pre-date the Jurassic. ?? 1989.

Explore related subjects

Marine and Petroleum GeologyArticle

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

R. Vugrinovich. 1989. Subsurface temperatures and surface heat flow in the Michigan Basin and their relationships to regional subsurface fluid movement. https://doi.org/10.1016/0264-8172(89)90076-7

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

KEEP EXPLORING

Related USGS reports

USGS2025

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↗
USGS2024

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↗
USGS2024

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↗