Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Seismic Data Compression”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

64 records · Page 4

Evidence of shock-compressed stishovite above 300 GPa

SiO 2 is one of the most fundamental constituents in planetary bodies, being an essential building block of major mineral phases in the crust and mantle of terrestrial planets (1–10 M E ). Silica at depths greater than 300 km may be present in the form of the rutile-type, high pressure polymorph stishovite (P4 2 /mnm) and its thermodynamic stability is of great interest for understanding the seismic and dynamic structure of planetary interiors. Previous studies on stishovite via static and dynamic (shock) compression techniques are contradictory and the observed differences in the lattice-level response is still not clearly understood. Here, laser-induced shock compression experiments at the LCLS- and SACLA XFEL light-sources elucidate the high-pressure behavior of stishovite on the lattice-level under in situ conditions on the Hugoniot to pressures above 300 GPa. We find stishovite is still (meta-)stable at these conditions, and does not undergo any phase transitions. This contradicts static experiments showing structural transformations to the CaCl 2 , α-PbO 2 and pyrite-type structures. However, rate-limited kinetic hindrance may explain our observations. These results are important to our understanding into the validity of EOS data from nanosecond experiments for geophysical applications.

58 GEOSCIENCES↗

Source characteristics and geological implications of the January 2016 induced earthquake swarm near Crooked Lake, Alberta

On 2016 January 12, an intraplate earthquake with an initial reported local magnitude (M L ) of 4.8 shook the town of Fox Creek, Alberta. While there were no reported damages, this earthquake was widely felt by the local residents and suspected to be induced by the nearby hydraulic-fracturing (HF) operations. We determine the earthquake source parameters using moment tensor inversions, and then detect and locate the associated swarm using a waveform cross-correlation based method. The broad-band seismic recordings from regional arrays suggest a moment magnitude (M) 4.1 for this event, which is the largest in Alberta in the past decade. Similar to other recent M ~ 3 earthquakes near Fox Creek, the 2016 January 12 earthquake exhibits a dominant strike-slip (strike = 184°) mechanism with limited non-double-couple components (~22 per cent). This resolved focal mechanism, which is also supported by forward modelling and P-wave first motion analysis, indicates an NE–SW oriented compressional axis consistent with the maximum compressive horizontal stress orientations delineated from borehole breakouts. Further detection analysis on industry-contributed recordings unveils 1108 smaller events within 3 km radius of the epicentre of the main event, showing a close spatial-temporal relation to a nearby HF well. The majority of the detected events are located above the basement, comparable to the injection depth (3.5 km) on the Duvernay shale Formation. The spatial distribution of this earthquake cluster further suggests that (1) the source of the sequence is an N–S-striking fault system and (2) these earthquakes were induced by an HF well close to but different from the well that triggered a previous (January 2015) earthquake swarm. Reactivation of pre-existing, N-S oriented faults analogous to the Pine Creek fault zone, which was reported by earlier studies of active source seismic and aeromagnetic data, are likely responsible for the occurrence of the January 2016 earthquake swarm and other recent events in the Crooked Lake area.

58 GEOSCIENCES↗

Radiation-induced aging of PDMS Elastomer TR-55: a summary of constitutive, mesoscale, and population-based models

Filled and cross-linked elastomeric rubbers are versatile network materials with a multitude of applications ranging from artificial organs and biomedical devices to cushions, coatings, adhesives, interconnects, and seismic-isolation-, thermal-, and electrical barriers. External factors like mechanical stress, temperature fluctuations, or radiation are known to create chemical changes in such materials that can directly affect the molecular weight distribution (MWD) of the polymer between cross-links and alter the structural and mechanical properties. From a Materials Science point of view it is highly desirable to understand, effect, and manipulate such property changes in a controlled manner. In this report we summarize our modeling efforts on a polysiloxane elastomer TR-55, which is an important component in several of our systems, and representative of a wide class of filled rubber materials. The primary aging driver in this work has been γ-radiation, and a variety of modeling approaches have been employed, including constitutive, mesoscale, and population-based models. The work utilizes diverse experimental data, including mechanical stress-strain and compression set measurements, as well as MWD measurements using multiquantum NMR.

36 MATERIALS SCIENCE↗

Subtask 1.4 - Techno-Economic Assessment of Regional Carbon Utilization Scenarios and Attendant Monitoring Technology

The U.S. Department of Energy (DOE) Carbon Storage Program seeks to develop technologies to significantly improve the effectiveness of carbon capture and storage (CCS), reduce the cost of implementation, and support the readiness of widespread commercial deployment by 2035. To accomplish this goal, technical and economic barriers must be addressed and the results used to inform regulators and industry on the safety, permanence, and viability of carbon storage. To help overcome these barriers, the Energy & Environmental Research Center (EERC) conducted three activities: 1) investigate the regional techno-economic life cycle of CO2 enhanced oil recovery (EOR) in North Dakota, 2) assist DOE in determining the prospective CO2 storage resource of residual oil zones (ROZs), and 3) develop a seismic technique for monitoring geologic storage of CO2. These three activities address two strategic goals enumerated in the DOE–EERC Cooperative Agreement, namely, 1) to advance foundational science and 2) to inform data-driven policies that enhance U.S. economic growth and job creation, energy security, and environmental quality. Activity 1 In North Dakota, there is synergistic incentive to develop and implement CO2 EOR technologies. Regional industry partners are actively seeking options that can cost-effectively improve the efficiency of power production, while also coproducing CO2 for use in EOR applications. A full-value-chain techno-economic assessment that links lignite resources, power generation, and EOR is needed to determine the implications for a state or region to continue to use these resources or look elsewhere for other energy sources and revenue streams. The integration of North Dakota’s lignite power industry with the state’s oil industry has the potential to revitalize legacy oil fields and produce up to 1 billion barrels of additional oil using CO2 EOR while providing a market for a CO2 commodity captured at coal-fired power plants. Revitalization of North Dakota’s legacy oil fields will bring new industry investment to the state, resulting in the growth of oil-related jobs and revenue through oil sales, severance tax, and associated oilfield development activities (drilling, pipelines, etc.). CO2 demand for EOR and incremental oil production were calculated for 201 conventional oil fields in North Dakota. Infrastructure development in the fields was estimated using 1) a traditional vertical well arrangement of CO2 injectors and oil producers and 2) an advanced case where horizontal wells were used exclusively. Pipeline, compression, and recycling infrastructure requirements were also estimated. When viewed from the CO2 capture and CO2 EOR sectors, the advanced case showed a positive investment scenario when including potential benefits of the 45Q tax credit policy. An economic evaluation of a full implementation of CO2 EOR showed an average $60 billion potential state economic impact estimated over the 37-year construction and production time frame could be realized with an assumed oil price of $65/bbl. Oil price ranges of $35–$90/bbl cause an estimated variation of ±$5 billion over the potential effort. The total net state government revenues were estimated over the 37-year construction and production time frame. With an assumed oil price of $65/bbl, total potential fiscal impact could average $7.4 billion. This result ranges about ±$2.5 billion for a potential EOR implementation effort for oil price ranges of $35–$90/bbl. This broad-brush assessment was focused on “magnitude-of-effect” perspectives of what CO2 capture and use for EOR could mean to the state. Whether CO2 EOR ever becomes widely adopted is not the issue at this stage of evaluation. The question being answered is, How valuable could CO2 capture and CO2 EOR be for the state and local jurisdictions in western North Dakota? The results of the current economic assessment remove any ambiguity to that question. Full implementation of CO2 EOR as modeled in this study provides substantial economic benefits to the state, even at low oil prices (i.e., $35/bbl). Job creation and government revenues could be substantial, depending on the scale of implementation. With a better understanding of what CO2 EOR might mean for the state’s economy, the focus can now shift to the policies and financial mechanisms needed to make private investment in CO2 capture and use as EOR competitive with other energy-related pursuits. Activity 2 ROZs contain remnants of oil not swept away by natural waterflooding. Although nearly all of the identified ROZs are found in West Texas, there is evidence that this category of oil reservoir can be found in other major oil-producing basins of the United States. Because the formation of these oil reservoirs is distinct from other conventional reservoirs, the method to estimate the prospective CO2 storage resource in ROZs may be equally distinct from the method used for conventional fields. As such, the EERC collaborated with DOE’s National Energy Technology Laboratory (NETL) to investigate the feasibility of CO2 storage in ROZs and determine what specific modification (if any) needs to be made to the current DOE volumetric method to estimate prospective CO2 storage resource in oil reservoirs. This activity built on a previous EERC investigation to identify ROZs in the Williston and Powder River Basins and upon collaborations between the EERC and NETL to develop storage resource estimation methods for different types of geologic CO2 storage. Review of the volumetric approach raised questions about the applicability of such an approach to hydrocarbon-bearing zones without the inclusion of potential hydrocarbon production. A production-based approach was presented by the EERC to inform the working group of additional methods of increasing the CO2 storage resource potential. Because the production-based method focuses more heavily on oil recovery than CO2 storage, it was not included in the working group manuscript. Activity 3 The EERC has used sparse and irregular seismic arrays to detect CO2 injected into a reservoir. The variability of seismic noise recorded by the sensors directly impacts the seismic interpretation. However, designing noise attenuation workflows effectively can be a more challenging task using sparse and irregular arrays for data acquisition. Based on the hypothesis that noncoherent noise is stationary during weekly data acquisition, a new machine-learning algorithm was developed to automatically estimate noise conditions using passive records, noise from active records, and ancillary information. The new noise attenuation processing workflow can improve the time-lapse detection of CO2 from sparse and irregular array acquisitions. Results show spatial and temporal noise patterns can be related to the location of the sensors and the weather conditions during the data acquisition and that our attenuation processing workflow can improve the time-lapse detection of CO2 from sparse and irregular array acquisitions. This effort was funded in part through the EERC–DOE Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Department of Commerce.

Peck, Wesly↗

Mechanical Testing Results on Core from Borehole U-15n, NNSS, in support of SPE

The Nevada National Security Site (NNSS) will serve as the geologic setting for a Source Physics Experiment (SPE) program. The SPE will provide 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 bore hole at varying depths and recording ground motions in instrumented bore holes located in two rings around the source positioned at different radii. Modeling using advanced simulation codes will be performed both a priori 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 bore holes including the geomechanical behavior of the site rock/structural features. This report presents a limited scope of work for an initial phase of primarily unconfined compression testing. Samples tested came from the U-15n core hole, which was drilled in granitic rock (quartz monzonite). The core hole was drilled at the location of the central SPE borehole, and thus represents material in which the explosive charges will be detonated. The U-15n location is the site of the first SPE, in Area 15 of the NNSS.

58 GEOSCIENCES↗

Seismic Analysis of Spatio-Temporal Fracture Generation During EGS Resource Development

This DOE GTO funded project was a collaborative effort between Array Information Technology and the University of California at Berkeley. The project goal is to improve technology to assess in-situ stress magnitude and orientation, kinematic fracture parameters, rupture size, rupture orientation as well as temporal changes and volumetric distribution of the injected fluid during Enhanced Geothermal System (EGS) resource development. The project leveraged high-frequency seismic data recorded by the LBNL 34-station permanent geophone network and seismic broadband data recorded by a temporary 33-station seismometer network that operated in The Geysers during the injection phase of the DOE GTO funded EGS demonstration project at the well Prati-32. The operation of the broadband network with high station density and good azimuthal coverage to monitor the injection phase produced an unprecedented dataset that is typically unavailable for EGS operations. The two datasets offer the advantages of studying the parameters under investigation over a broader frequency band and the higher sensitivity of the broadband seismometers resulted in a richer dataset with a lower magnitude of completeness. Moment tensor (MT) inversion have been employed to develop a catalog of earthquake source solutions including moment tensor and finite source rupture area. This resulted in a 168-event waveform-based seismic moment tensor catalog for events 0.7≤Mw≤3.9 based on a semi-automated moment tensor analysis technique. The MT solutions were used for in-situ stress estimation during the EGS injection phases. We found an approximate 15-degree counterclockwise rotation of the least compressive stress sigma-3, and a rotation of the maximum compressive stress sigma-1 toward the vertical as the injected volume of water increased. We developed a rupture-area magnitude scaling relationship of The Geysers earthquakes obtained from finite-source inversion for fault slip that reveals a very high correlation to published results. Based on the rupture area-magnitude relationship we mapped Mw from the LBNL earthquake catalog to estimate the activated fracture area, its location and orientation within the injection volume of Prati-32. The double-difference Wadati (DDW) technique was applied to high-accuracy P- and S-wave differential travel times derived from waveform cross correlation to image fluid saturation in the subsurface based on Vp/Vs-ratio. The recently developed technique has been applied to estimate temporal changes of fluid distribution in fault zones and in volcanic regions but not yet in engineered geo-reservoirs. The improvement of the technology includes high-resolution waveform cross correlation applied to the high-quality datasets with lower magnitude of completeness resulting in a higher number of recorded events. The improvement allowed to estimate both, the spatial and temporal distribution of injected fluids and the state of the injected fluid during reservoir stimulation. Results of the DDW analysis revealed the development of an injection-derived steam plume during the early injection phase, particularly in the northern part of the Prati-32 study area.

15 GEOTHERMAL ENERGY↗

Stratigraphically Controlled Stress Variations at the Hydraulic Fracture Test Site-1 in the Midland Basin, TX

We investigated the relationship between stratigraphy, stress, and microseismicity at the Hydraulic Fracture Test Site-1. The site comprises two sets of horizontal wells in the Wolfcamp shale and a deviated well drilled after hydraulic fracturing. Regional stresses indicate normal/strike-slip faulting with E-W compression. Stress measurements in vertical and horizontal wells show that the minimum principal stress varies with depth. Strata with high clay and organic content show high values of the least compressive stress, consistent with the theory of viscous stress relaxation. By integrating data from core, logs, and the hydraulic fracturing stages, we constructed a stress profile for the Wolfcamp sequence, which predicts how much pressure is required for hydraulic fracture growth. We applied the results to fracture orientation data from image logs to determine the population of pre-existing faults that are expected to slip during stimulation. We also determined microseismic focal plane mechanisms and found slip on steeply dipping planes striking NW, consistent with the orientations of potentially active faults predicted by the stress model. This case study represents a general approach for integrating stress measurements and rock properties to predict hydraulic fracture growth and the characteristics of injection-induced microseismicity.

microseismicity↗

High pressure thermoelasticity and sound velocities of Fe-Ni-Si alloys

The Earth's iron-dominant core is known to contain nickel from cosmochemical analysis and some amount of light elements from geophysical constraints on density and seismic wave velocities. Although there have been several studies to constrain thermoelastic properties of iron-alloys, there has been no systematic study on the effects of nickel and light elements on properties of iron using the same experimental methods and data analysis approach. We conducted nuclear resonant inelastic X-ray scattering and X-ray diffraction experiments on body-centered cubic and hexagonal close-packed (hcp) Fe 0.91 Ni 0.09 and Fe 0.8 Ni 0.1 Si 0.1 up to 104 GPa and 86 GPa, respectively, and compare to similar measurements conducted on hcp-Fe up to 171 GPa. Specifically, we determine the Debye sound velocity from the low-energy transfer region of the (partial) phonon density of states (DOS) using the equation of state determined for each material and a new approach which utilizes information criteria and probability distributions. Nickel decreases the shear velocity of iron, while 10 at% Si has little to no effect on the shear velocity of Fe 0.91 Ni 0.09 . We observe that the shape of the phonon DOS of these alloys remains similar with increasing pressure. In the measured compression range, we therefore apply a generalized scaling law to describe the volume dependence of the phonon DOS and find that the vibrational Grüneisen parameters of hcp-Fe 0.91 Ni 0.09 are nearly indistinguishable from those hcp-Fe and those for Fe 0.8 Ni 0.1 Si 0.1 trend lower. From the vibrational free energy, we constrain the harmonic vibrational component of thermal pressure, which shows a significant positive deviation from theoretical calculations of hcp-Fe at pressures and temperatures of Earth's core. Collectively, our results demonstrate that the effects of nickel should be considered when modeling iron-rich planetary cores.

36 MATERIALS SCIENCE↗

Phase I Project: Fiber Optic Distributed Acoustic Sensing for Periodic Hydraulic Tests

The extraction of heat from hot rock requires circulation of fluid through fracture networks. Because the geometry and connectivity of these fractures determines the efficiency of fluid circulation, many tools are used to characterize fractures before and after development of the reservoir. Under this project, a new tool was developed that allows hydraulic connectivity between geothermal boreholes to be identified. Nanostrain in rock fractures is measured using fiber optic distributed acoustic sensing (DAS). This strain is measured in one borehole in response to periodic pressure pulses induced in another borehole. The strain in the fractures represents hydraulic connectivity between wells. DAS is typically used at frequencies of Hz to kHz, but strain at mHz frequencies were measured for this project. The tool was demonstrated in the laboratory and in the field. In the laboratory, strain in fiber optic cables was measured in response to compression due to oscillating fluid pressure. DAS recorded strains as small as 10 picometer/m in response to 1 cm of water level change. At a fractured crystalline rock field site, strain was measured in boreholes. Fiber-optic cable was mechanically coupled borehole walls using pressured flexible liners. In one borehole 30 m from the oscillating pumping source, pressure and strain were measured simultaneously. The DAS system measured fracture displacement at frequencies of less than 1 mHz (18 min periods) and amplitudes of less than 1 nm, in response to fluid pressure changes of less 20 Pa (2 mm of water). The attenuation and phase shift of the monitored strain signal is indicative of the permeability and storage (compliance) of the fracture network that connects the two wells. The strain response as a function of oscillation frequency is characteristic of the hydraulic structure of the formation. This is the first application of DAS to the measurement of low frequency strain in boreholes. It has enormous potential for monitoring geothermal reservoirs for purposes of understanding reservoir compliance and for assuring security of injection fluids. Periodic pressure pulses can be induced by oscillating injection or pumping during operation so the system could provide real time reservoir data. DAS cable may already be installed at a site as it is becoming increasingly used for seismic observation. Simulations conducted for this project indicate that strain should be propagated through borehole cements, so observations can be made outside of well casing. In uncased holes, the cable would need to be mechanically coupled to the borehole wall to provide measurements. One option would be to install fiber into cemented and abandoned boreholes to extend their utility.

15 GEOTHERMAL ENERGY↗

Nonlinear Time Domain Seismic Soil-Structure Interaction (SSI) Deep Soil Site Methodology Development

Currently the Department of Energy (DOE) and the nuclear industry perform seismic soil-structure interaction (SSI) analysis using equivalent linear numerical analysis tools. For lower levels of ground motion, these tools should produce reasonable in-structure response values for evaluation of existing and new facilities. For larger levels of ground motion these tools likely overestimate the in-structure response (and therefore structural demand) since they do not consider geometric nonlinearities (such as gaping and sliding between the soil and structure) and are limited in the ability to model nonlinear soil behavior. The current equivalent linear SSI (SASSI) analysis approach either joins the soil and structure together in both tension and compression or releases the soil from the structure for both tension and compression. It also makes linear approximations for material nonlinearities and generalizes energy absorption with viscous damping. This produces the potential for inaccurately establishing where the structural concerns exist and/or inaccurately establishing the amplitude of the in-structure responses. Seismic hazard curves at nuclear facilities have continued to increase over the years as more information has been developed on seismic sources (i.e. faults), additional information gathered on seismic events, and additional research performed to determine local site effects. Seismic hazard curves are used to develop design basis earthquakes (DBE) that are used to evaluate nuclear facility response. As the seismic hazard curves increase, the input ground motions (DBE’s) used to numerically evaluation nuclear facility response increase causing larger in-structure response. As ground motions increase so does the importance of including nonlinear effects in numerical SSI models. To include material nonlinearity in the soil and geometric nonlinearity using contact (gaping and sliding) it is necessary to develop a nonlinear time domain methodology. This methodology will be known as, NonLinear Soil-Structure Interaction (NLSSI). In general NLSSI analysis should provide a more accurate representation of the seismic demands on nuclear facilities their systems and components. INL, in collaboration with a Nuclear Power Plant Vender (NPP-V), will develop a generic Nuclear Power Plant (NPP) structural design to be used in development of the methodology and for comparison with SASSI. This generic NPP design has been evaluated for the INL soil site because of the ease of access and quality of the site specific data. It is now being evaluated for a second site at Vogtle which is located approximately 15 miles East-Northeast of Waynesboro, Georgia and adjacent to Savanna River. The Vogtle site consists of many soil layers spanning down to a depth of 1058 feet. The reason that two soil sites are chosen is to demonstrate the methodology across multiple soil sites. The project will drive the models (soil and structure) using successively increasing acceleration time histories with amplitudes. The models will be run in time domain codes such as ABAQUS, LS-DYNA, and/or ESSI and compared with the same models run in SASSI. The project is focused on developing and documenting a method for performing time domain, non-linear seismic soil structure interaction (SSI) analysis. Development of this method will provide the Department of Energy (DOE) and industry with another tool to perform seismic SSI analysis.

58 GEOSCIENCES↗