Connecting microstructural attributes and permeability from 3D tomographic images of in situ shear-enhanced compaction bands using multiscale computations

Tomographic images taken inside and outside a compaction band in a field specimen of Aztec sandstone are analyzed by using numerical methods such as graph theory, level sets, and hybrid lattice Boltzmann/finite element techniques. The results reveal approximately an order of magnitude permeability reduction within the compaction band. This is less than the several orders of magnitude reduction measured from hydraulic experiments on compaction bands formed in laboratory experiments and about one order of magnitude less than inferences from two-dimensional images of Aztec sandstone. Geometrical analysis concludes that the elimination of connected pore space and increased tortuosities due to the porosity decrease are the major factors contributing to the permeability reduction. In addition, the multiscale flow simulations also indicate that permeability is fairly isotropic inside and outside the compaction band

Submarine landslides in the Santa Barbara Channel as potential tsunami sources

International audienceRecent investigations using the Monterey Bay Aquarium Research Institutes (MBARI) Remotely Operated Vehicles (ROVs) "Ventana" and "Tiburon" and interpretation of MBARI's EM 300 30 kHz multibeam bathymetric data show that the northern flank of the Santa Barbara Basin has experienced massive slope failures. Of particular concern is the large (130 km2) Goleta landslide complex located off Coal Oil Point near the town of Goleta, that measures 14.6-km long extending from a depth of 90 m to nearly 574 m deep and is 10.5 km wide. We estimate that approximately 1.75 km3 has been displaced by this slide during the Holocene. This feature is a complex compound submarine landslide that contains both surfical slump blocks and mud flows in three distinct segments. Each segment is composed of a distinct head scarp, down-dropped head block and a slide debris lobe. The debris lobes exhibit hummocky topography in the central areas that appear to result from compression during down slope movement. The toes of the western and eastern lobes are well defined in the multibeam image, whereas the toe of the central lobe is less distinct. Continuous seismic reflection profiles show that many buried slide debris lobes exist and comparison of the deformed reflectors with ODP Drill Site 149, Hole 893 suggest that at least 200 000 years of failure have occurred in the area (Fisher et al., 2005a). Based on our interpretation of the multibeam bathymetry and seismic reflection profiles we modeled the potential tsunami that may have been produced from one of the three surfical lobes of the Goleta slide. This model shows that a 10 m high wave could have run ashore along the cliffs of the Goleta shoreline. Several other smaller (2 km2 and 4 km2) slides are located on the northern flank of the Santa Barbara Basin, both to the west and east of Goleta slide and on the Conception fan along the western flank of the basin. One slide, named the Gaviota slide, is 3.8 km2, 2.6 km long and 1.7 km wide. A distinct narrow scar extends from near the eastern head wall of this slide for over 2km eastward toward the Goleta slide and may represent either an incipient failure or a remnant of a previous failure. Push cores collected within the main head scar of this slide consisted of hydrogen sulfide bearing mud, possibly suggesting active fluid seepage and a vibra-core penetrated ~50 cm of recent sediment overlying colluvium or landslide debris confirming the age of ~300 years as proposed by Lee et al. (2004). However, no seeps or indications of recent movement were observed during our ROV investigation within this narrow head scar indicating that seafloor in the scar is draped with mud

Detecting Compaction Disequilibrium with Anisotropy of Magnetic Susceptibility

In clay-rich sediment, microstructures and macrostructures influence how sediments deform when under stress. When lithology is fairly constant, anisotropy of magnetic susceptibility (AMS) can be a simple technique for measuring the relative consolidation state of sediment, which reflects the sediment burial history. AMS can reveal areas of high water content and apparent overconsolidation associated with unconformities where sediment overburden has been removed. Many other methods for testing consolidation and water content are destructive and invasive, whereas AMS provides a nondestructive means to focus on areas for additional geotechnical study. In zones where the magnetic minerals are undergoing diagenesis, AMS should not be used for detecting compaction state. By utilizing AMS in the Santa Barbara Basin, we were able to identify one clear unconformity and eight zones of high water content in three cores. With the addition of susceptibility, anhysteretic remanent magnetization, and isothermal remanent magnetization rock magnetic techniques, we excluded 3 out of 11 zones from being compaction disequilibria. The AMS signals for these three zones are the result of diagenesis, coring deformation, and burrows. In addition, using AMS eigenvectors, we are able to accurately show the direction of maximum compression for the accumulation zone of the Gaviota Slide

Increased fluid flow activity in shallow sediments at the 3 km Long Hugin Fracture in the central North Sea

The North Sea hosts a wide variety of seafloor seeps that may be important for transfer of chemical species, such as methane, from the Earth's interior to its exterior. Here we provide geochemical and geophysical evidence for fluid flow within shallow sediments at the recently discovered, 3-km long Hugin Fracture in the Central North Sea. Although venting of gas bubbles was not observed, concentrations of dissolved methane were significantly elevated (up to six-times background values) in the water column at various locations above the fracture, and microbial mats that form in the presence of methane were observed at the seafloor. Seismic amplitude anomalies revealed a bright spot at a fault bend that may be the source of the water column methane. Sediment porewaters recovered in close proximity to the Hugin Fracture indicate the presence of fluids from two different shallow (<500m) sources: (i) a reduced fluid characterized by elevated methane concentrations and/or high levels of dissolved sulfide (up to 6 mmol L−1), and (ii) a low-chlorinity fluid (Cl ∼305 mmol L−1) that has low levels of dissolved methane and/or sulfide. The area of the seafloor affected by the presence of methane-enriched fluids is similar to the footprint of seepage from other morphological features in the North Sea

Fracture-pattern growth in the deep, chemically reactive subsurface

Arrays of natural opening-mode fractures show systematic patterns in size and spatial arrangement. The controls on these factors are enigmatic, but in many cases the depth of formation appears to be critical. Physical, potentially depth-dependent factors that could account for these variations include confining stress, fluid pressure, and strain rate; these factors are common inputs to existing fracture models. However, temperature-dependent chemical processes likely exert an equally important control on patterns, and such processes have not yet been rigorously incorporated into models of fracture formation. Here we present a spring-lattice model that simulates fracturing in extending sedimentary rock beds, while explicitly accounting for cementation during opening of fractures, and for rock failure via both elastic and time-dependent failure criteria. Results illustrate three distinct fracturing behaviors having documented natural analogs, which we here term fracture facies. “Exclusionary macrofracturing” occurs at shallow levels and produces large, widely spaced, uncemented fractures; “multi-scale fracturing” occurs at moderate depth and produces partially cemented fractures having a wide range of sizes and spacings; and “penetrative microfracturing” occurs at great depth and produces myriad narrow, sealed fractures that are closely and regularly spaced. The effect of depth is primarily to accelerate both dissolution and precipitation reactions via increased temperature and porewater salinity; the specific depth range of each fracture facies will vary by host-rock lithology, grain size, strain rate, and thermal history

Fracture growth during exhumation in low-permeability rock formations&amp;#8212;the role of fluid PVT properties

&amp;lt;p&amp;gt;Detailed fluid inclusion analyses of fracture cements in tightly cemented hydrocarbon-bearing sandstones and shales reveal that natural fractures tend to form under conditions approaching maximum burial, coinciding with hydrocarbon generation, and during incipient exhumation. Fluid inclusion analyses also reveal that these fractures form under abnormal (above-hydrostatic) pore fluid pressures. While compaction disequilibrium can account for elevated pore fluid pressures that promote fracture growth during early prograde burial, hydrocarbon maturation is likely the primary driver for fracture growth under peak burial conditions. Tectonic processes and thermal stresses provide secondary drivers. Thermal contraction with exhumation and cooling of the rock mass can promote fracture growth depending on the PVT properties of the fluid phase. The possible contribution of hydrocarbon generation after peak burial as a driver for fracture growth during incipient exhumation is discussed.&amp;lt;/p&amp;gt; </jats:p

Poroelastic modeling of basement fault reactivation caused by saltwater disposal near Venus, Johnson County, Texas

Presented at the 54th US Rock Mechanics/Geomechanics Symposium held online, 28 June-1 July 2020To assess the potential for fault reactivation in response to wastewater injection near Venus, Texas, we conducted fully coupled 3D poroelastic finite element simulations for a site-specific geomechanical analysis. We find that simulations using the best estimates of in situ stress azimuth and fault orientation, and the Byerlee’s friction coefficient of 0.6 do not predict fault reactivation, in contrast to observed seismicity patterns. Increasing the maximum horizontal stress azimuth by 10° and the basement fault dip by 5°, both within the uncertainty space of the input parameters, and lowering the friction coefficient of the fault in the basement to 0.35, leads to fault reactivation in the basement as observed. Using the same model geometry but a friction coefficient of 0.6 leads to fault reactivation within Ellenburger, which is inconsistent with observed hypocenter depths at Venus. Expanding the model domain from 5 disposal wellbores to 35 increases the excess pore pressure which favors basement fault reactivation. These simulations demonstrate the sensitivity of geomechanical models of fault reactivation in response to injection requiring high-quality field parameters and help refine the azimuths of the horizontal stresses, the basement fault dip, and the basement fault friction coefficient on the basis of earthquake hypocenter temporal distribution.Bureau of Economic Geolog