<p dir="ltr">The Alpine Fault, located in the South Island, New Zealand, is a subject of intense geological study due to its potential to trigger large earthquakes. It encompasses a complex system with the interplay of mechanics, thermodynamics, and fluid. Gaining insights into these systems not only enhances our understanding of the fault but also holds the potential to guide risk mitigation efforts.</p><p dir="ltr">The damage extent and fracture networks within the metamorphic rock mass adjacent to the fault can be effectively characterized by seismic anisotropy, an elastic property of rock, where seismic waves travel at different speeds with variation directions. This thesis presents a comprehensive exploration of seismic anisotropy in the hanging wall immediately adjacent to the principal slip zone of the Alpine Fault in New Zealand. Leveraging the borehole seismic data from a unique scientific drilling project and advanced numerical modeling techniques, the ultimate goal is to invert and parameterize the bulk seismic anisotropy.</p><p dir="ltr">Motivated by these challenges, the thesis undertakes several key initiatives: The first effort focuses on gaining a comprehensive understanding of an innovative method for seismic measurement: Distributed Acoustic Sensing (DAS) – examining its operational principles, factors influencing observed wavelets, and how it contrasts with traditional point sensors for accurate interpretation. Subsequently, the research introduces the implementation of an open-source seismic wave solver designed for modeling elastic wave propagation in complicated anisotropic media. This solver is further optimized for computational efficiency with its performance rigorously benchmarked.</p><p dir="ltr">With this preparedness, the inversion is further facilitated by high-performance computing (HPC) and a deep-learning algorithm specifically designed for automatically picking transit times. The inverted bulk elastic constants, compared to the intact rock, reveal 28% to 35% reductions in qP-wave velocity, characterizing the damage due to mesoscale fracture. Further analysis sheds light on the existence of orthogonal fracture sets and an intricate geometrical arrangement that agree with the previous borehole image log. This represents an advancement in our ability to characterize and understand the geologic processes with seismic anisotropy.</p>

<p>Crustal stress measurements play a crucial role in understanding how the subsurface deforms. As one of the most popular methods for stress characterization in deep wellbores, borehole breakout analysis examines the shape of drilling-induced compressive failures to determine stress directions and magnitudes, assuming that the rock formation is both isotropic and linearly elastic. To ensure accurate stress interpretations, the dissertation investigates the validity of underlying presumptions from two perspectives: (1) the effect of rock anisotropy (i.e., elastic anisotropy, and strength anisotropy) on wellbore failure patterns; and (2) the characterization of rock nonlinear elastic mechanical behaviors. </p> <p>The developed computer program, <em><strong>EASAfail</strong></em>, has broad applicability in calculating wellbore failure patterns for a wide range of scenarios. It takes into account factors such as elastic stiffness matrices of the rock, stress tensors in the surrounding environment, and the presence of weak planes. The program's generality allows it to handle various rock types with different degrees of symmetry in their elastic properties, as well as weak planes that are weaker than the intact rock matrix. By analyzing these factors, the program reveals that the patterns of wellbore failure in elastic and strength anisotropic rock formations are highly influenced by the sliding of weak planes. Complications from two modes of borehole failure, either in the intact rock matrix or in the weak planes, can cause the breakout azimuth to deviate from the direction of the minimum horizontal stress. </p> <p>In addition to hypothetical scenarios generated from numerical models, a case study from the field is presented to underscore the impact of foliations on the anomalous rotations of breakout azimuths. The wellbore was located in Northeastern Alberta, Canada, transecting both the sedimentary column and crystalline basement. Breakout rotations identified from caliper and image logs were highly likely caused by the slippage along foliations, supported by the close correlation between breakout azimuths and dip directions of foliations as well as polarization directions analyzed from dipole sonic logs. Stress magnitudes constrained from Monte Carlo simulations further reveal a lower stress field when rock anisotropy is taken into account, compared to what is inferred conventionally. </p> <p>The characterization of rock nonlinear elasticity involves the utilization of the third-order elastic (TOE) model. To measure the TOE moduli in a static manner, test-specific protocols were proposed based on the nonlinear stress-strain behaviors of the rock. By arranging the stress-strain responses obtained from hydrostatic, uniaxial, and triaxial compressive tests into a linear system of equations, it becomes possible to invert the equations for the TOE moduli. These analytical equations were validated through calculations from finite element models. </p> <p>By employing the established protocols, the TOE moduli were derived for four different rock types with varying pore structures when subjected to hydrostatic and uniaxial compressions. The TOE model successfully captured the nonlinear stress-strain responses exhibited by Indiana limestone, Vif-type Fontainebleau sandstone, and Snake River Plain basalt. However, it was found to be inadequate for Franc-type Fontainebleau sandstone, which displayed noticeable hysteresis and experienced significant strains. Future geomechanical applications will undoubtedly gain advantages from utilizing the inverted TOE moduli obtained through static measurements, as they allow for the examination of the impacts of nonlinear elasticity in rocks. </p>

Rock glaciers play a significant role in the alpine debris transport system. For practical and engineering considerations, identifying the internal structure and its relationship to surface characteristics is significant in terms of how a rock glacier settles during periods of melting, and the mode of deformation. A better understanding of these factors is important for engineers, engineering geologists and geomorphologists who must make prudent evaluations of rock glaciers as potential sites for human development and uses. It is equally important for evaluating potential stores for water on other planets such as Mars. Ground penetrating radar (GPR) shows that the internal structure of a lobate rock glacier located in the San Juan Mountains of southwest Colorado consists of continuous to semi-continuous horizontal layers of ice-supersaturated sediments and coarse blocky rockslide debris which likely formed through catastrophic episodes of rockfall from the cirque headwall. Folds in the uppermost layers correspond to the surface expression of ridges and furrows, indicating that compressive stresses originating in the steep accumulation zone are transmitted downslope through the rock glacier. The rock glacier is a composite feature that formed by a process involving the development and overlap of discrete flow lobes that have overridden older glacial moraine and protalus rampart materials. The latter materials have been incorporated into the present flow structure of the rock glacier. The discovery of rock glacier-like features on Mars suggests the presence of flowing, or once-flowing ice-rock mixtures. These landforms, which include lobate debris aprons, concentric crater fill and lineated valley fill, hold significant promise as reservoirs of stored water ice that could be used as fuel sources for human exploration of Mars and provide a frozen record of the climatic history of the planet. To this end, the rock glacier in this study was used as a surrogate for similar Martian landforms. Liquid water, found to be abundant in this rock glacier, occurs within a network of interconnected channels that permeate throughout the landform. In terms of water storage within Martian analogs, consideration must include the possibility that some water ice may be stored in relatively pure form within lenses and vein networks that are supplied by seasonal frost accumulation and/or water influx from below.

<p>Probing the subsurface for evidence related to the degradation of porous mediums and the evolution of damage mechanisms has been a long-standing challenge in geophysics. As such imaging and predicting fracture network development has remained a difficult area for subsurface science for decades despite the seminal and significant works put forward by many researchers. While this has provide great understanding about the behaviours and properties of natural porous media, there is still much that needs to be explored particularly in regard to the mineralogical composition and chemistry of clay-rich rocks. Despite the fact that argillaceous rocks which consist of different types of clays and varied mineral composition are ubiquitous in nature and are often the target of several technologies (e.g. geotechnical engineering, nuclear waste storage and disposal,hydrocarbon exploration and extraction, carbon capture and sequestration, etc.), many studies focus primarily on the bulk properties or the percentage of components in the matrix. For these reason and due to the problems that can be encountered with natural rocks that contain a swelling clay component whether randomly distirbuted or localized in consolidated globs in zones of the matrix, the influence of clay chemistry in relation to fracture development which is not well characterized, especially during desaturation is investigated with analogue rock samples which were systematically fabricated for this purpose.</p> <p><br></p> <p>The research performed in this dissertation investigated, the applicability of  the fabrication protocol for developing synthetic rocks with desirable rock like features and behavior, the impact and relationship between the rock properties, the microstructural composition, water loss, and the macroscopic behavior of the analogue rocks, focusing on the structure and chemistry of the constituent clay materials. Synthetic rocks were fashioned with the necessary geometries, properties, and material compositions. On the macroscopic scale the fracture and drying behavior of the synthetic rocks were examined with 3D X-ray microscopy and further evaluated through the utility of acoustic emission monitoring, water loss monitoring, and unconfined compressive testing. On the finer scale (nano-microscale), the chemical and mechanical properties, and behavior of select clays was explored by exploiting several methods of material characterization which also included cation exchange experiments coupled with inductively coupled plasma – optical emission spectrometry (ICP-OES). </p> <p><br></p> <p>For the finer scale, experiments verified that calcined kaolinite clay had a different mineral structure and negligible to non-existence shrinkage abilities. In contrast, the montmorillonite clays possessed higher and similar moisture contents but, owing to the different principal cations these clays interacted a bit differently in the highly akaline environment, experienced varying degrees of shrinkage, and had observedly minor structural dissimilarities. For the relatively larger scale, the emergence of damage, extent of the damage network, and the patterns of the crack network mainly depended on the microstructural composition of the analogue rocks, particularly it's clay chemistry and/ or distribution. The location of damage depended on the emplacement and percentage of swelling clay in the matrix, and numerical investigations with peridynamics revealed that the observed damage was a consequence of the action of the swelling and non-swelling components of the matrix. Furthermore, if the microstructure consisted of no clay or calcined kaolinite the AE activity was solely attributed to interfacial processes that occurred during fluid front movement. If the microstructure consisted of a particular montmorillonite, the cracks propagated in the direction of the drying front. Conversely, for montmorillonite clay predominated by a different principal cation, the crack network developed and propagated differently during water loss. Additionally, on the laboratory core scale, properties and behavior similar to natural rocks were confirmed and the rock strength, porosity, AE activity, and velocities were primarily affected by the microstructural composition of the analogue rocks. </p> <p><br></p> <p>An added challenge for investigating and monitoring evolving systems and processes, whether on the laboratory or field scale, is the problem of extracting useful information from the physical data that can be used to identify signatures of developing processes, and changes in the properties or the behavior of a system. Here, data driven machine learning modeling and clustering techniques were undertaken to build a mechanistic understanding of the AE activity generated during drying. The intent is for this work to add to the fundamental research aimed at developing methods that will robustly detect and extract signatures related to evolutionary processes or features in the AE signals, and group them according to some degree of similarity. Such research will support reliable interpretations of the physical data for predictions of the behavior of systems, development of engineering controls, and improvement of the understanding of intrinsic dynamics related to complex processes particularly those that occur in clay-rich systems.</p> <p><br></p> <p>Combined chemical and mechanical investigations have great potential for unraveling practical challenges in subsurface science, especially regarding damage processes in clay-rich rock systems, and identifying and interpreting the presence of discontinuities from geophysical data. The present findings are useful for establishing a link between the constituent clay and observed damage, and improving our understanding of the development of damage in clay-bearing systems. These results provide insight on the influence of swelling clay and the chemistry of such clays on the generation of cracks and crack networks in rock like materials which can be useful for the characterization of damage in both laboratory and the field. The work presented here can also be a basis for further experiments that aim to uncover methods and protocols that will help with the indirect characterization of evolutionary processes, damage mechanisms, and damage in clay rich porous media. Additionally support for the use of analogue rocks in experimental rock physics, architected with specific material compositions, pore structures, crack systems, or clay fractions, is provided here. </p>