Coupled poroelastic stressing and pore-pressure accumulation along pre-existing faults in deep basement contribute to recent occurrence of seismic events at subsurface energy exploration sites. Our coupled fluid-flow and geomechanical model describes the physical processes inducing seismicity corresponding to the sequential stimulation operations in Pohang, South Korea. Simulation results show that prolonged accumulation of poroelastic energy and pore pressure along a fault can nucleate seismic events larger than Mw3 even after terminating well operations. In particular the possibility of large seismic events can be increased by multiple-well operations with alternate injection and extraction that can enhance the degree of pore-pressure diffusion and subsequent stress transfer through a rigid and low-permeability rock to the fault. This study demonstrates that the proper mechanistic model and optimal well operations need to be accounted for to mitigate unexpected seismic hazards in the presence of the site-specific uncertainty such as hidden/undetected faults and stress regime.
Measuring the size and orientation of borehole breakouts is one of the primary methods for determining the orientation and magnitudes of the in situ stresses in the subsurface. To better understand the effects of anisotropy on borehole breakouts, experiments were conducted on Mancos Shale, a finely laminated mudrock. A novel testing configuration was developed to conduct borehole breakout experiments in a standard triaxial vessel and load frame. Samples were prepared at three different orientations and deformed under 6.9 to 20.7 MPa confining pressure. The results show a variation of peak strength and breakout geometry depending on the lamination orientation. Samples deformed parallel to laminations failed at a higher maximum compressive stress than samples deformed perpendicular to laminations, which were stronger than inclined samples. These relationships are quantified by a cosine-based failure envelope. Observed breakout shapes in perpendicular samples are V-shaped and symmetric around the borehole, which advance as a series of fractures of increasing size into the sidewalls. In inclined samples, fractures form along weaker laminations planes and grow in an en echelon pattern towards the axial stress direction. In parallel samples, long fractures grow from the wellbore towards the axial stress direction. The observed geometries highlight potential sources of error in calculating in situ stresses from borehole breakouts.
The Nevada National Security Site (NNSS) serves as the geologic setting for a Source Physics Experiment (SPE) program. The SPE provides ground truth data to create and improve strong ground motion and seismic S-wave generation and propagation models. The NNSS was chosen as the test bed because it provides a variety of geologic settings ranging from relatively simple to very complex. Each series of SPE testing will comprise the setting and firing of explosive charges (source) placed in a central borehole at varying depths and recording ground motions in instrumented boreholes located in two rings around the source, positioned at different radii. Modeling using advanced simulation codes will be performed both before and after each test to predict ground response and to improve models based on acquired field data, respectively. A key component in the predictive capability and ultimate validation of the models is the full understanding of the intervening geology between the source and the instrumented boreholes including the geomechanical behavior of the site's rock/structural features. This report summarizes unconfined compression testing (UCS) from coreholes U-15n#12 and U-15n#13 and compares those datasets to UCS results from coreholes U-15n and U-15n#10. U-15n#12 corehole was drilled at -60° to the horizontal and U-15n#13 was drilled vertically in granitic rock (quartz monzonite) after the third SPE shot. Figure 1 illustrates at the surface, U 15n#12 and U-15n#13 coreholes were approximately 30 meters and 10 meters from the central SPE borehole (U-15n) respectively. Corehole U-15n#12 intersects the central SPE borehole (U 15n) at a core depth of 174 feet (approximately 150 feet vertical depth). The location of U 15n#12 and U-15n#13 is the site of the first, second and third SPE's, in Area 15 of the NNSS.
Dynamic Brazilian tension (DBR) tests from core hole U-15n are part of a larger material characterization effort for the Source Physics Experiment (SPE) project. This larger effort encompasses characterizing Climax Stock granite rock from the Nevada National Security Site (NNSS) both before and after each SPE shot. The current test series includes DBR tests on dry intact granite and fault material at depths of -85 and -150 ft.
Triaxial compression tests from core hole U-15n are part of a larger material characterization effort for the Source Physics Experiment (SPE) project. This larger effort encompasses characterizing Climax Stock granite rock from the Nevada National Security Site (NNSS) both before and after each SPE shot. The current test series includes triaxial compression tests on dry and saturated intact granite and fault material at 100, 200, 300, and 400 MPa confining pressure.
The Nevada National Security Site (NNSS) serves as the geologic setting for a Source Physics Experiment (SPE) program. The SPE provides ground truth data to create and improve strong ground motion and seismic S-wave generation and propagation models. The NNSS was chosen as the test bed because it provides a variety of geologic settings ranging from relatively simple to very complex. Each series of SPE testing will comprise the setting and firing of explosive charges (source) placed in a central borehole at varying depths and recording ground motions in instrumented boreholes located in two rings around the source, positioned at different radii. Modeling using advanced simulation codes will be performed both before and after each test to predict ground response and to improve models based on acquired field data, respectively. A key component in the predictive capability and ultimate validation of the models is the full understanding of the intervening geology between the source and the instrumented boreholes including the geomechanical behavior of the site's rock/structural features. This memorandum reports on an initial phase of unconfined compression testing from corehole U-15n#10. Specimens tested came from the U-15n#10 core hole, which was drilled at -60° to the horizontal in granitic rock (quartz monzonite) after the second SPE shot (SPE-2). Figure 1 illustrates at the surface, the core hole was approximately 90 feet from the central SPE borehole. Corehole U 15n#10 intersects the central SPE borehole (U-15n) at a core depth of 170 feet (approximately 150 feet vertical depth) which is within the highly damaged zone of SPE-2. The U-15n#10 location is the site of the first, second and third SPE's, in Area 15 of the NNSS.
The state of stress in the earth is complicated and it is difficult to determine all three components and directions of the stress. However, the state of stress affects all activities which take place in the earth, from causing earthquakes on critically stressed faults, to affecting production from hydraulically fractured shale reservoirs, to determining closure rates around a subterranean nuclear waste repository. Current state of the art methods commonly have errors in magnitude and direction of up to 40%. This is especially true for the intermediate principal stress. This project seeks to better understand the means which are used to determine the state of stress in the earth and improve upon current methods to decrease the uncertainty in the measurement. This is achieved by a multipronged experimental investigation which is closely coupled with advanced constitutive and numeric modeling.
In recent years, seismicity rates in the US have dramatically risen due to increased activity in onshore oil and gas production. This project attempts to tie observations about induced seismicity to dehydration reactions in laumontite, a common mineral found in fault gouge in crystalline basement formations. It is the hypothesis of this study that in addition to pressurerelated changes in the in situ stress state, the injection of wastewater pushes new fluids into crystalline fault fracture networks that are not in chemical equilibrium with the mineral assemblages, particularly laumontite in fault gouge. Experiments were conducted under hydrothermal conditions where samples of laumontite were exposed to NaC1 brines at different pH values. After exposure to different fluid chemistries for 8 weeks at 90° C, we did not observe substantial alteration of laumontite. In hydrostatic compaction experiments, all samples deformed similarly in the presence of different fluids. Pore pressure decreases were observed at the start of a 1 week hold at 85° C in a 1M NaC1 pH 3 solution, suggesting that acidic fluids might stabilize pore pressures in basement fault networks. Friction experiments on laumontite and kaolinite powders showed both materials have similar coefficients of friction. Mixtures with partial kaolinite content showed a slight decrease in the coefficient of friction, which could be sufficient to trigger slip on critically stressed basement faults.
The study described in this report involves heated and unheated pressurized slot testing to determine thermo-mechanical properties of the Tptpll (Tertiary, Paintbrush, Topopah Spring Tuff Formation, crystal poor, lower lithophysal) and Tptpul (upper lithophysal) lithostratigraphic units at Yucca Mountain, Nevada. A large volume fraction of the proposed repository at Yucca Mountain may reside in the Tptpll lithostratigraphic unit. This unit is characterized by voids, or lithophysae, which range in size from centimeters to meters, making a field program an effective method of measuring bulk thermal-mechanical rock properties (thermal expansion, rock mass modulus, compressive strength, time-dependent deformation) over a range of temperature and rock conditions. The field tests outlined in this report provide data for the determination of thermo-mechanical properties of this unit. Rock-mass response data collected during this field test will reduce the uncertainty in key thermal-mechanical modeling parameters (rock-mass modulus, strength and thermal expansion) for the Tptpll lithostratigraphic unit, and provide a basis for understanding thermal-mechanical behavior of this unit. The measurements will be used to evaluate numerical models of the thermal-mechanical response of the repository. These numerical models are then used to predict pre- and post-closure repository response. ACKNOWLEDGEMENTS The authors would like to thank David Bronowski, Ronnie Taylor, Ray E. Finley, Cliff Howard, Michael Schuhen (all SNL) and Fred Homuth (LANL) for their work in the planning and implementation of the tests described in this report. This is a reprint of SAND2004-2703, which was originally printed in July 2004. At that time, it was printed for a restricted audience. It has now been approved for unlimited release.
As a part of the U.S. Department of Energy (DOE) SubTER (Subsurface Technology and Engineering Research, Development and Demonstration) initiative, University of Wisconsin- Madison, Sandia National Laboratories, and Lawrence Berkeley National Laboratory conducted the Permeability (k) and Induced Seismicity Management for Energy Technologies (kISMET) project. The objectives of the project are to define the in situ status of stress in the Sanford Underground Research Facility (SURF) in Lead, South Dakota and to establish the relations between in situ stress and induced fracture through hydraulically stimulating the fracture. (SURF) in Lead, South Dakota. In situ tests are conducted in a 7.6 cm diameter and 100 long vertical borehole located in the 4850 Level West Access Drift near Davies Campus of SURF (Figure 1). The borehole is located in the zone of Precambrian Metamorphic Schist.
A helium leakage detection system was modified to measure gas permeability on extracted cores of nearly impermeable rock. Here we use a Helium - Mass - Spectrometry - Permeameter (HMSP) to conduct a constant pressure, steady state flow test through a sample using helium gas. Under triaxial stress conditions, the HMSP can measure flow and estimate permeability of rocks and geomaterials down to the nanodarcy scale (10-21 m2). In this study, measurements of flow through eight shale samples under hydrostatic conditions were in the range of 10-7 to 10-9 Darcy. We extend this flow measurement technology by dynamically monitoring the release of helium from a helium saturated shale sample during a triaxial deformation experiment. The helium flow, initially extremely low, consistent with the low permeability of shale, is observed to increase in advance of volume strain increase during deformation of the shale. This is perhaps the result of microfracture development and flow path linkage through the microfractures within the shale. Once microfracturing coalescence initiates, there is a large increase in helium release and flow. This flow rate increase is likely the result of development of a macrofracture in the sample, a flow conduit, later confirmed by post-test observations of the deformed sample. The release rate (flow) peaks and then diminishes slightly during subsequent deformation; however the post deformation flow rate is considerably greater than that of undeformed shale.
A helium leakage detection system was modified to measure gas permeability on extracted cores of nearly impermeable rock. Here we use a Helium - Mass - Spectrometry - Permeameter (HMSP) to conduct a constant pressure, steady state flow test through a sample using helium gas. Under triaxial stress conditions, the HMSP can measure flow and estimate permeability of rocks and geomaterials down to the nanodarcy scale (10-21 m2). In this study, measurements of flow through eight shale samples under hydrostatic conditions were in the range of 10-7 to 10-9 Darcy. We extend this flow measurement technology by dynamically monitoring the release of helium from a helium saturated shale sample during a triaxial deformation experiment. The helium flow, initially extremely low, consistent with the low permeability of shale, is observed to increase in advance of volume strain increase during deformation of the shale. This is perhaps the result of microfracture development and flow path linkage through the microfractures within the shale. Once microfracturing coalescence initiates, there is a large increase in helium release and flow. This flow rate increase is likely the result of development of a macrofracture in the sample, a flow conduit, later confirmed by post-test observations of the deformed sample. The release rate (flow) peaks and then diminishes slightly during subsequent deformation; however the post deformation flow rate is considerably greater than that of undeformed shale.
The essential characteristics of the issue of radioactive waste management can be conceptualized as complex, with a variety of facets and uncertainty. These characteristics tend to cause people to perceive the issue of radioactive waste management as a 'risk'. This study was initiated in response to a desire to understand the perceptions of risk that the Korean public holds towards radioactive waste and the relevant policies and policy-making processes. The study further attempts to identify the factors influencing risk perceptions and the relationships between risk perception and social acceptance.
There is a long history of testing crushed salt as backfill for the Waste Isolation Pilot Plant program, but testing was typically done at 100 C or less. Future applications may involve backfilling crushed salt around heat-generating waste packages, where near-field temperatures could reach 250 C or hotter. A series of experiments were conducted to investigate the effects of hydrostatic stress on run-of-mine salt at temperatures up to 250 C and pressures to 20 MPa. The results of these tests were compared with analogous modeling results. By comparing the modeling results at elevated temperatures to the experimental results, the adequacy of the current crushed salt reconsolidation model was evaluated. The model and experimental results both show an increase in the reconsolidation rate with temperature. The current crushed salt model predicts the experimental results well at a temperature of 100 C and matches the overall trends, but over-predicts the temperature dependence of the reconsolidation. Further development of the deformation mechanism activation energies would lead to a better prediction of the temperature dependence by the crushed salt reconsolidation model.
Three dimensional finite element analyses were performed to evaluate the structural integrity of the caverns located at the Bayou Choctaw (BC) site which is considered a candidate for expansion. Fifteen active and nine abandoned caverns exist at BC, with a total cavern volume of some 164 MMB. A 3D model allowing control of each cavern individually was constructed because the location and depth of caverns and the date of excavation are irregular. The total cavern volume has practical interest, as this void space affects total creep closure in the BC salt mass. Operations including both cavern workover, where wellhead pressures are temporarily reduced to atmospheric, and cavern enlargement due to leaching during oil drawdowns that use water to displace the oil from the caverns, were modeled to account for as many as the five future oil drawdowns in the six SPR caverns. The impacts on cavern stability, underground creep closure, surface subsidence, infrastructure, and well integrity were quantified.
This report represents the final product of a background literature review conducted for the Nuclear Waste Management Organization of Japan (NUMO) by Sandia National Laboratories, Albuquerque, New Mexico, USA. Internationally, research of hydrological and transport processes in the context of high level waste (HLW) repository performance, has been extensive. However, most of these studies have been conducted for sites that are within tectonically stable regions. Therefore, in support of NUMO's goal of selecting a site for a HLW repository, this literature review has been conducted to assess the applicability of the output from some of these studies to the geological environment in Japan. Specifically, this review consists of two main tasks. The first was to review the major documents of the main HLW repository programs around the world to identify the most important hydrologic and transport parameters and processes relevant in each of these programs. The review was to assess the relative importance of processes and measured parameters to site characterization by interpretation of existing sensitivity analyses and expert judgment in these documents. The second task was to convene a workshop to discuss the findings of Task 1 and to prioritize hydrologic and transport parameters in the context of the geology of Japan. This report details the results and conclusions of both of these Tasks.
Detailed statistical analysis of the experimental data from testing of alumina-loaded epoxy (ALOX) composites was conducted to better understand influences of the selected compositional properties on the compressive strength of these ALOX composites. Analysis of variance (ANOVA) for different models with different sets of parameters identified the optimal statistical model as, y{sub l} = -150.71 + 29.72T{sub l} + 204.71D{sub l} + 160.93S{sub 1l} + 90.41S{sub 2l}-20.366T{sub l}S{sub 2l}-137.85D{sub l}S{sub 1l}-90.08D{sub l}S{sub 2l} where y{sub l} is the predicted compressive strength, T{sub l} is the powder type, D{sub l} is the density as the covariate for powder volume concentration, and S{sub il}(i=1,2) is the strain rate. Based on the optimal statistical model, we conclude that the compressive strength of the ALOX composite is significantly influenced by the three main factors examined: powder type, density, and strain rate. We also found that the compressive strength of the ALOX composite is significantly influenced by interactions between the powder type and the strain rate and between the powder volume concentration and the strain rate. However, the interaction between the powder type and the powder volume concentration may not significantly influence the compressive strength of the ALOX composite.
3-D finite element analyses were performed to evaluate the structural integrity of caverns located at the Strategic Petroleum Reserve's Big Hill site. State-of-art analyses simulated the current site configuration and considered additional caverns. The addition of 5 caverns to account for a full site and a full dome containing 31 caverns were modeled. Operations including both normal and cavern workover pressures and cavern enlargement due to leaching were modeled to account for as many as 5 future oil drawdowns. Under the modeled conditions, caverns were placed very close to the edge of the salt dome. The web of salt separating the caverns and the web of salt between the caverns and edge of the salt dome were reduced due to leaching. The impacts on cavern stability, underground creep closure, surface subsidence and infrastructure, and well integrity were quantified. The analyses included recently derived damage criterion obtained from testing of Big Hill salt cores. The results show that from a structural view point, many additional caverns can be safely added to Big Hill.
Specimens of poled 'chem-prep' PNZT ceramic from batch HF803 were tested under hydrostatic, uniaxial, and constant stress difference loading conditions at three temperatures of -55, 25, and 75 C and pressures up to 500 MPa. The objective of this experimental study was to obtain the electro-mechanical properties of the ceramic and the criteria of FE (Ferroelectric) to AFE (Antiferroelectric) phase transformations so that grain-scale modeling efforts can develop and test models and codes using realistic parameters. The poled ceramic undergoes anisotropic deformation during the transition from a FE to an AFE structure. The lateral strain measured parallel to the poling direction was typically 35 % greater than the strain measured perpendicular to the poling direction. The rates of increase in the phase transformation pressures per temperature changes were practically identical for both unpoled and poled PNZT HF803 specimens. We observed that the retarding effect of temperature on the kinetics of phase transformation appears to be analogous to the effect of shear stress. We also observed that the FE-to-AFE phase transformation occurs in poled ceramic when the normal compressive stress, acting perpendicular to a crystallographic plane about the polar axis, equals the hydrostatic pressure at which the transformation otherwise takes place.
This document is the final report for the Compressed Air Energy Storage Monitoring to Support Refrigerated-Mined Rock Cavern Technology (CAES Monitoring to Support RMRCT) (DE-FC26-01NT40868) project to have been conducted by CAES Development Co., along with Sandia National Laboratories. This document provides a final report covering tasks 1.0 and subtasks 2.1, 2.2, and 2.5 of task 2.0 of the Statement of Project Objectives and constitutes the final project deliverable. The proposed work was to have provided physical measurements and analyses of large-scale rock mass response to pressure cycling. The goal was to develop proof-of-concept data for a previously developed and DOE sponsored technology (RMRCT or Refrigerated-Mined Rock Cavern Technology). In the RMRCT concept, a room and pillar mine developed in rock serves as a pressure vessel. That vessel will need to contain pressure of about 1370 psi (and cycle down to 300 psi). The measurements gathered in this study would have provided a means to determine directly rock mass response during cyclic loading on the same scale, under similar pressure conditions. The CAES project has been delayed due to national economic unrest in the energy sector.
To establish mechanical material properties of cellular concrete mixes, a series of quasi-static, compression and tension tests have been completed. This report summarizes the test methods, set-up, relevant observations, and results from the constitutive experimental efforts. Results from the uniaxial and triaxial compression tests established failure criteria for the cellular concrete in terms of stress invariants I{sub 1} and J{sub 2}. {radical}J{sub 2} (MPa) = 297.2 - 278.7 exp{sup -0.000455 I}{sub 1}{sup (MPa)} for the 90-pcf concrete {radical}J{sub 2} (MPa) = 211.4 - 204.2 exp {sup -0.000628 I}{sub 1}{sup (MPa)} for the 60-pcf concrete
Sandia is currently developing a lead-zirconate-titanate ceramic 95/5-2Nb (or PNZT) from chemically prepared ('chem-prep') precursor powders. Previous PNZT ceramic was fabricated from the powders prepared using a 'mixed-oxide' process. The specimens of unpoled PNZT ceramic from batch HF803 were tested under hydrostatic, uniaxial, and constant stress difference loading conditions within the temperature range of -55 to 75 C and pressures to 500 MPa. The objective of this experimental study was to obtain mechanical properties and phase relationships so that the grain-scale modeling effort can develop and test its models and codes using realistic parameters. The stress-strain behavior of 'chem-prep' PNZT under different loading paths was found to be similar to that of 'mixed-oxide' PNZT. The phase transformation from ferroelectric to antiferroelectric occurs in unpoled ceramic with abrupt increase in volumetric strain of about 0.7 % when the maximum compressive stress, regardless of loading paths, equals the hydrostatic pressure at which the transformation otherwise takes place. The stress-volumetric strain relationship of the ceramic undergoing a phase transformation was analyzed quantitatively using a linear regression analysis. The pressure (P{sub T1}{sup H}) required for the onset of phase transformation with respect to temperature is represented by the best-fit line, P{sub T1}{sup H} (MPa) = 227 + 0.76 T (C). We also confirmed that increasing shear stress lowers the mean stress and the volumetric strain required to trigger phase transformation. At the lower bound (-55 C) of the tested temperature range, the phase transformation is permanent and irreversible. However, at the upper bound (75 C), the phase transformation is completely reversible as the stress causing phase transformation is removed.
We have conducted five hydraulic fracturing stress measurement campaigns in Korea, involving 13 test holes ranging in depth from 30 to 250 m, at locations from North Seoul to the southern coast of the peninsula. The measurements reveal consistent crustal stress magnitudes and directions that suggest persistence throughout western and southern Korea. The maximum horizontal stress {sigma}{sub H} is oriented between ENE-WSW and E-W, in accord with plate movement and deformation, and with directions indicated by both focal mechanism solutions from earthquakes inland and offshore as well as borehole breakouts in mainland China close to its eastern coast. With respect to magnitudes, the vertical stress is the overall minimum stress at all tested locations, suggesting a thrust faulting regime within the relatively shallow depths reached by our tests. Typically, such a stress regime becomes one favoring strike-slip at greater depths, as is also indicated by the focal mechanism solutions around Korea.
A suite of laboratory triaxial compression and triaxial steady-state creep tests provide quasi-static elastic constants and damage criteria for bedded rock salt and dolomite extracted from Cavern Well No.1 of the Tioga field in northern Pennsylvania. The elastic constants, quasi-static damage criteria, and creep parameters of host rocks provides information for evaluating a proposed cavern field for gas storage near Tioga, Pennsylvania. The Young's modulus of the dolomite was determined to be 6.4 ({+-}1.0) x 10{sup 6} psi, with a Poisson's ratio of 0.26 ({+-}0.04). The elastic Young's modulus was obtained from the slope of the unloading-reloading portion of the stress-strain plots as 7.8 ({+-}0.9) x 10{sup 6} psi. The damage criterion of the dolomite based on the peak load was determined to be J{sub 2}{sup 0.5} (psi) = 3113 + 0.34 I{sub 1} (psi) where I{sub 1} and J{sub 2} are first and second invariants respectively. Using the dilation limit as a threshold level for damage, the damage criterion was conservatively estimated as J{sub 2}{sup 0.5} (psi) = 2614 + 0.30 I{sub 1} (psi). The Young's modulus of the rock salt, which will host the storage cavern, was determined to be 2.4 ({+-}0.65) x 10{sup 6} psi, with a Poisson's ratio of 0.24 ({+-}0.07). The elastic Young's modulus was determined to be 5.0 ({+-}0.46) x 10{sup 6} psi. Unlike the dolomite specimens under triaxial compression, rock salt specimens did not show shear failure with peak axial load. Instead, most specimens showed distinct dilatancy as an indication of internal damage. Based on dilation limit, the damage criterion for the rock salt was estimated as J{sub 2}{sup 0.5} (psi) = 704 + 0.17 I{sub 1} (psi). In order to determine the time dependent deformation of the rock salt, we conducted five triaxial creep tests. The creep deformation of the Tioga rock salt was modeled based on the following three-parameter power law as {var_epsilon}{sub s} = 1.2 x 10{sup -17} {sigma}{sup 4.75} exp(-6161/T), where {var_epsilon}{sub s} is the steady state strain rate in s{sup -1}, {sigma} is the applied axial stress difference in psi, and T is the temperature in Kelvin.