Estimating SHmax azimuth with P sources and vertical geophones: Use P-P reflection amplitudes or use SV-P reflection times?

We compared two methods for extracting the azimuth of maximum horizontal stress (SHmax) from 3D land-based seismic data generated by a P source and recorded with vertical geophones. In the first method, we used the direct-SV mode that is produced by all land-based P sources. P sources generate SV illumination that radiates in all azimuth directions from a source station and creates SV-P reflections that are recorded by vertical geophones. Unless stratigraphy has steep dip, SV-P raypaths recorded by vertical geophones are the reverse of P-SV raypaths recorded by horizontal geophones. Thus, SV-P data provide the same S-wave sensitivity to stress fields as popular P-SV data do. In the second method, we retrieved P-P reflections and then performed an amplitude-versus-incident-angle (AVA) analysis of the amplitude-gradient behavior of P-P reflection wavelets. We did this analysis in narrow azimuth corridors to determine the gradient of reflection-wavelet amplitudes as a function of azimuth. This P-P AVA amplitude-gradient method has been of great interest in the reflection seismology community since it was introduced in the late 1990s. Each of these methods, AVA analysis of the gradient of P-P reflection amplitudes and azimuth-dependent arrival times of SV-P reflections can be used to determine the azimuth of SHmax stress. We compare the results of the two methods with ground truth measurements of SHmax azimuth at a CO2 sequestration site in the Michigan Basin. SHmax azimuths were determined from P-P and SV-P data at three major boundaries at depths of approximately 3500 ft (1067 m), 5500 ft (1676 m), and 7500 ft (2286 m). Two estimates of SHmax azimuth (one using SV-P data and one using P-P data) were made at each stacking bin inside a 24 mi2 (62 km2) image space. The result was approximately 98,000 estimates of SHmax azimuth across each of these three boundaries for each of these two prediction strategies. Histogram displays of PP AVA gradient estimates had peaks at correct azimuths of SHmax at all three depths, but the spread of the distributions widened with depth and split into two peaks at the deepest boundary. In contrast, each histogram of SHmax azimuth predicted by azimuth-dependent SV-P traveltimes had a single, definitive peak that was positioned at the correct SHmax azimuth at all three boundary depths.

Fluid Histories of Middle Ordovician fault–fracture hydrothermal dolomite oil fields in the southern Michigan Basin, U.S.A.

ABSTRACT Dolomitized fault–fracture structures in the Trenton and Black River formations (TBR) are the type example for “hydrothermal” petroleum reservoirs world-wide. However, fluid histories of these structures are only partially understood. Trenton and Black River reservoirs in the southern Michigan Basin are composed of fault-associated, vertical dolomite bodies that are highly fractured and brecciated. Open spaces are partially to completely filled by saddle dolomite and less frequently by calcite cement. Cathodoluminescence microstratigraphies of void-filling carbonate cements are not correlatable between oil fields. Fluid inclusion homogenization temperatures (Th) measured in carbonate cements indicate two fluid endmembers: a warm fluid (∼ 80° to 180° C) and a hot fluid (180° to ∼ 260° C). Increasing Th proximal to the underlying Proterozoic Mid-Michigan Rift (MMR) suggest that the hot fluids emanated from the rift area. Included fluids are saline (16.1–49.4 wt. % NaCl equivalent), and salinity likely is sourced from overlying Silurian Salina Group evaporites. First melting temperatures (Tfm), interpreted as eutectic temperatures (Te), of fluids range from –112° C to –50° C, indicating a complex Na–Ca–KCl brine; the expected composition of dissolved Salina salts. Lower Te proximal to the MMR suggest the rift as a source of additional complexing ions. C and O isotope values for carbonate cements are depleted with respect to δ18O (–6.59 to –12.46‰ VPDB) relative to Ordovician seawaters, and somewhat depleted with respect to δ13C (–1.22 to +1.18‰ VPDB). Equilibrium calculations from δ18O and Th values indicate that cement precipitating waters were highly evolved (+1.3 to +14.4‰ δ18O‰ VSMOW) compared to Ordovician and Silurian seawaters (–5.5‰ δ18O‰ VSMOW). Strontium isotope values indicate two fluid sources: Proterozoic basement and Late Silurian evaporites. Values of 87Sr/86Sr for cements in the Freedom, Napoleon, Reading, and Scipio fields (0.7086–0.7088) are influenced by warm water sourced from Silurian strata, and values for cements in the Albion, Branch County, and Northville fields (0.7091–0.7110) record continental basement signatures. Cement precipitating fluids in TBR oil fields likely have similar sources and timing. However, water–rock interactions along fault pathways modified source waters, giving each oil field a unique petrographic and geochemical signature. Fluid movement in TBR oil fields likely were initiated by reactivation of basement faulting during Silurian–Devonian tectonism.

Enhancing Paleoreef Reservoir Characterization through Machine Learning and Multi-Attribute Seismic Analysis: Silurian Reef Examples from the Michigan Basin

Historically, Silurian pinnacle reef complexes in the Michigan Basin have been largely identified using 2D seismic with very little research on the reservoir characterization of these reefs using 3D seismic data. By incorporating a high-resolution 3D dataset constrained by a well-studied and data-rich paleoreef reservoir, the Puttygut reef, seismic attributes were correlated to petrophysical properties through machine learning and self-organizing maps (SOMs). A suite of structural and frequency-based attributes was calculated from pre-stack time migrated (PSTM) seismic data, with only a subset of them selected as SOM inputs. Structural attributes enhanced details in the reef but frequency attributes were overall more useful for correlating with reservoir quality. A strong relationship between certain combination percentages of attributes and certain sections of the reef with porosity and permeability was found after the SOM results were compared to wireline log and core analysis data. Areas with high permeability and porosity correlated with the average frequency and spectral decomposition at 29 and 81 Hz. Areas with high porosity and varying permeability correlated with the average frequency and spectral decomposition at 29, 57, and 81 Hz. Areas with intermediate porosity correlated with the average frequency and spectral decomposition at 29 and 57 Hz. The efficacy of the procedure was then demonstrated on two nearby reefs with very similar results.

Variation in Assemblages of Freshwater Mussels Downstream of Dams and Dam Removals in the Lake Michigan Basin, Michigan, USA

Impacts of barriers, including dam removals, on aquatic diversity are poorly understood. We assessed freshwater mussel assemblages and microhabitat downstream of uncontrolled and controlled low-head dams and low-head dam removals in Michigan, USA. The objectives of this study were to quantify whether downstream mussel assemblages and microhabitat parameters differ by anthropogenic barrier and along a downstream gradient, and to determine parameters that were predictors of mussel diversity and density. Sampling consisted of standardized timed-searches and quadrat excavations. Results suggest that areas downstream of dams had higher mussel diversity than dam removals, and mussel assemblages differed along a downstream gradient for uncontrolled and controlled dams. Indicator Species Analyses determined mussel species representative of downstream river reaches from uncontrolled low-head dams and removals. Predictor variables for mussel assemblages included substrate classes and total suspended solids. Controlled dams contained the least fine substrates (%) and highest coarse substrates (%) in downstream reaches. This study suggests that rivers with uncontrolled low-head dams and removals provide downstream habitat that support viable mussel assemblages. Results from this study also suggest that evidence of mussel assemblage recovery following dam removal may take many years. Quantification of barrier-related impacts, as shown in this study, are imperative to guide conservation of aquatic fauna.

Bentonite swelling characteristics with a hypersaline multi-component pore fluid

Swelling characteristics of compacted bentonite when hydrated with a hypersaline pore fluid (332 g/L total dissolved solids; 6.6 mol/L ionic strength) are reported. The pore fluid mimics the multiple constituents and their concentrations found for the Cobourg limestone of the Michigan Basin and is dominated by sodium (25% mole fraction) with some potassium, calcium, and magnesium (10%, 5%, and 4% mole fractions). Measurements of swell pressure for two sodium bentonites when hydrated under conditions of zero volume increase are reported. Swell pressure reached a peak within 10–30 h from the onset of hydration, followed by a continual decrease over 1 year of testing from chemical interaction between the bentonite and pore fluid. After 1 year, the swell pressure of the MX-80 bentonite tested decreased by a factor of nine relative to the peak swell pressure with deionized water when the dry density was 1.6 Mg/m3. Swell pressures increased as dry density increased. However, chemical interactions appear to have more influence on swell pressure than density for the pore fluid examined as a swell pressure of just under 1200 kPa was measured for MX-80 after 1.8 years of hydration when compacted to the highest dry density of 1.8 Mg/m3 examined.