Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic sources”

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.

At least 55 records · Page 3

Seismic vibration source

A system for vibrating the earth in a location where seismic mapping is to take place is described. A relatively shallow hole formed in the earth, such as a hole 10 feet deep, placing a solid propellant in the hole, sealing a portion of the hole above the solid propellant with a device that can rapidly open and close to allow a repeatedly interrupted escape of gas. The propellant is ignited so that high pressure gas is created which escapes in pulses to vibrate the earth.

Dowler, W. L.↗

Quantifying the effect of 3D models on moment tensor results using synthetic data

Moment tensors provide vital information on seismic source properties for seismic events. Moment tensors require seismic wavespeed models to compute Green’s functions, which measure the impulse response between a given source and receiver. Traditionally, researchers have used one-dimensional velocity models to calculate Green’s functions since 1D Green’s functions are computationally cheap to compute. Local 1D velocity models can also accurately model waveforms at short distances (< 500 km). However, 1D velocity models do not account for lateral heterogeneity, which can cause significant misfit in tectonically complex regions such as the Middle East (Covellone and Savage, 2012). Green’s functions calculated using 3D seismic wavespeed models have been shown to perform better in tectonically complex regions (Covellone and Savage, 2012; Kintner and Modrak, 2022), so we are interested in quantifying the effect of considering 3D structure on moment tensor inversion results. The Middle East is an ideal study area for a synthetic moment tensor test for two reasons. Firstly, the Middle East is a tectonically complex region that has been heavily studied. Secondly, the Middle East has significant seismic activity throughout the region, but imperfect station coverage due to limited open data through large swaths of the domain. The tectonic complexity and uneven station coverage will test real-world performance even in a synthetic experiment.

58 GEOSCIENCES↗

Deep compressed seismic learning for fast location and moment tensor inferences with natural and induced seismicity

Fast detection and characterization of seismic sources is crucial for decision-making and warning systems that monitor natural and induced seismicity. However, besides the laying out of ever denser monitoring networks of seismic instruments, the incorporation of new sensor technologies such as Distributed Acoustic Sensing (DAS) further challenges our processing capabilities to deliver short turnaround answers from seismic monitoring. In response, this work describes a methodology for the learning of the seismological parameters: location and moment tensor from compressed seismic records. In this method, data dimensionality is reduced by applying a general encoding protocol derived from the principles of compressive sensing. The data in compressed form is then fed directly to a convolutional neural network that outputs fast predictions of the seismic source parameters. Thus, the proposed methodology can not only expedite data transmission from the field to the processing center, but also remove the decompression overhead that would be required for the application of traditional processing methods. An autoencoder is also explored as an equivalent alternative to perform the same job. We observe that the CS-based compression requires only a fraction of the computing power, time, data and expertise required to design and train an autoencoder to perform the same task. Implementation of the CS-method with a continuous flow of data together with generalization of the principles to other applications such as classification are also discussed.

54 ENVIRONMENTAL SCIENCES↗

Local-Distance Seismic Event Relocation and Relative Magnitude Estimation, Applications to Mining Related Seismicity in the Powder River Basin, Wyoming

Recent efforts to characterize small ( M w < 3 ) seismic events at local distances have become more important because of the increased observation of human-triggered and induced seismicity and the need to advance nuclear explosion monitoring capabilities. The signals generated by low-magnitude seismic sources necessitate the use of nearby short-period observations, which are sensitive to local geological heterogeneity. Local to near-regional distance ( < 300 km ) surface and shear waves can dominate short-period observations from small, shallow seismic sources. In this work, we utilize these observations to estimate precise, relative locations and magnitudes of ~ 700 industrial mining events in Wyoming, using nearly 360,000 observations. The precise, relative location estimates (with formal location uncertainty estimates of less than 1 km) collapse a diffuse collection of mining events into discrete clusters associated with individual blasting operations. We also invert the cross-correlation amplitudes to estimate precise, relative moment magnitude estimates, which help validate and identify disparities in the event sizes reported by regional network catalogs. Joint use of multiple phases allows for the inclusion of more seismic events due to the increase in the number of observations. In some cases, using a single phase allowed us to relocate only 50% of the original reported seismic events within a cluster. Combining shear- and surface-wave phases increased the number of events to above 90% of the original events, allowing us to characterize a broader range of event sizes, source to station distances, and event distributions. This analysis takes a step toward making a fuller characterization of small industrial seismic events observed at local distances.

58 GEOSCIENCES↗

Field Validation of MVA Technology for Offshore CCS: Novel Ultra-High-Resolution 3D Marine Seismic Technology (P-Cable) (Final Report)

The objectives of the proposed study were to deploy and validate a specific monitoring technology, high-resolution 3D marine seismic (HR3D), appropriate for large-demonstration and commercial-scale offshore CCS sites. The project accomplished successful acquisition two HR3D seismic surveys. The first HR3D dataset was over the offshore injection site of the Tomakomai, Japan integrated pilot CCS project, which at the time of survey acquisition was actively injecting CO 2 . The first survey also represented a successful international collaboration between the DOE NETL program and Japan’s national CCS program and was the first successful acquisition and use of HR3D over an active CO 2 injection site (Meckel, Feng et al. 2019). The Tomakomai HR3D survey successfully tested a novel 4-streamer HR3D system array in which, for the first time, no cross-cable (aka “P-Cable”) was utilized and only four GeoEel streamers were used instead of the standard 12-streamer configuration. Consequently, this was not, strictly speaking, a deployment of the “P-Cable” system of (Planke and Berndt 2004) but rather a modified version, thereof, and it is the first known demonstration of the modified system configuration. One very positive outcome from the Japanese collaboration earlier in the project was the ability to learn from the Japanese how they used tail buoys with GPS to determine the position of the seismic source and receivers in time and space. Based on that experience, GCCC designed and built six GPS receivers that could be used to position the streamer receivers and the seismic source via tail buoys. A fundamental advance that was made on the original design, was the ability to directly power the tail buoy GPS units and transfer data through the streamers (i.e., vs. the batteries used at Tomakomai). The bulkiness of the GPS batteries caused drag and episodic surging of the buoys, which affected data quality by lifting up the tail end of the streamers so the receivers were not at the same depth. The units were tested onshore for accuracy and functionality, and the design was subsequently and successfully tested in marine acquisition mode during the SLP survey acquisition. The marine acquisition test and survey satisfied Subtasks 2.2.2, Novel Positioning Technology Selection and Subtask 2.2.3, Novel Positioning Technology Deployment. Results of the novel positioning technology selection (Subtask 2.2.2) were considered successful and will be incorporated in future HR3D seismic acquisition projects to reduce costs, improve deployment safety at sea, and integrate both seismic and data recording via a single data transfer through the streamers to the recording system. The project also established a permitting process through NETL NEPA compliance, which included an Environmental Assessment in a marine setting and is required for conducting these types of surveys using Federal funding. The permitting process charted a “boilerplate,” which can allow future surveys related to other funded projects to move forward more expeditiously. Future improvements that could be considered are more robust seals on the GPS module and stronger materials (especially joints) on tail buoy fabrication. These would increase fixed costs, but would be advisable and probably more economic long-term if multiple HR3D surveys are planned. Project Accomplishments include: • Pre-survey Sensitivity Study • Marine geochemistry methods and data analysis • Successful HR3D seismic dataset acquired @ Tomakomai active CO 2 injection marine site • Developed advanced seismic processing techniques • No NRMS anomalies detected in overburden; Demonstration of containment • Repeatability study • Second survey collected @ San Luis Pass, TX • 4D application using positioning techniques developed in the project for monitoring were successful

3D seismic GPS positioning↗

Uncertainty Quantification of Geophysical Inversion Using Stochastic Partial Differential Equations (LDRD #218329)

This report summarizes work completed under the Laboratory Directed Research and Development (LDRD) project "Uncertainty Quantification of Geophysical Inversion Using Stochastic Differential Equations." Geophysical inversions often require computationally expensive algorithms to find even one solution, let alone propagating uncertainties through to the solution domain. The primary purpose of this project was to find more computationally efficient means to approximate solution uncertainty in geophysical inversions. We found multiple computationally efficient methods of propagating Earth model uncertainty into uncertainties in solutions of full waveform seismic moment tensor inversions. However, the optimum method of approximating the uncertainty in these seismic source solutions was to use the Karhunen-Love theorem with data misfit residuals. This method was orders of magnitude more computationally efficient than traditional Monte Carlo methods and yielded estimates of uncertainty that closely approximated those of Monte Carlo. We will summarize the various methods we evaluated for estimating uncertainty in seismic source inversions as well as work toward this goal in the realm of 3-D seismic tomographic inversion uncertainty.

58 GEOSCIENCES↗

Identifying geological structures through microseismic cluster and burst analyses complementing active seismic interpretation

At the Decatur carbon capture and storage site (IL, USA) CO 2 has been injected from 2011–2014 and from 2017 to present near the base of the Lower Mt. Simon Sandstone saline reservoir, resulting in microseismicity. Microseismicity is mainly located in the basement and distributed in distinct spatial clusters. The lack of significant impedance contrasts within the basement makes the interpretation of active-source seismic reflection data challenging, however, recent reprocessing allowed to resolve faults above and at the top of the basement. These faults generally do not coincide with the location of microseismic events and their continuation to the general depth of the seismic events cannot be assumed. This paper shows how the interpretation of the microseismicity can complement structural interpretations of active-source seismic reflection data. In particular, we analyze clusters and bursts (abrupt increases) of microseismicity, identify unresolved, smaller-scale weaknesses and extract statistical parameters. These parameters allow comparisons with the interpreted faults, and with fracture sets intercepted by boreholes. During injection at the Decatur site, the injection pressure was kept far below fracture pressure, nevertheless, seismic events were induced and spread far beyond the expected extent of the CO 2 plume. We argue that local stress transfers related to the CO 2 injection reactivated pre-existing fractures within the critically stressed basement. Finally, we conducted a slip tendency analysis for faults interpreted from active seismic, selected cluster, bursts and nodal planes from focal mechanisms to determine if the interpreted structures are optimally oriented with respect to the stress regime. Our results suggest that the orientation of fractures close to the injection well, generally shows slight deviations from the optimal orientation for slip. This might indicate either slight local deviations of the maximum horizontal stress azimuth from the average direction used in the analysis, or the lack of optimally oriented fractures at this location.

58 GEOSCIENCES↗

Estimating Explosion Yields Using Moment Tensor Solutions and Seismic Moment

We report seismic moment, a measurable and well-understood quantity of seismic sources, is used to estimate the yield of explosions. Application of such a method in the past, as in the manner of m b -derived yields, has been complicated by the effect of variations in the explosion working point, depth, and secondary source effects (such as spalling and tectonic release) on the observed moment. We start using the full (six-element) moment tensor solution, which can capture the relevant source physics and, at least in theory, better isolate the primary explosion source. The moment-to-yield ratio is then estimated using an explosion source model which, provided with emplacement conditions, can relate the two parameters. We discuss the major sources of uncertainty associated with the method, and calibrate it with chemical and nuclear explosions at the Nevada National Security Site. We then apply the method to published moment tensor solutions for the six declared North Korean nuclear explosions that occurred between 2006 and 2017. The results are mostly consistent with other yield estimates made using a variety of high-frequency methods. This technique is a new approach to estimating explosive yield and simple to implement, as much of the complexity is captured by the source models.

58 GEOSCIENCES↗

MTUQ: a framework for estimating moment tensors, point forces, and their uncertainties

SUMMARY We introduce MTUQ, an open-source Python package for seismic source estimation and uncertainty quantification, emphasizing flexibility and operational scalability. MTUQ provides MPI-parallelized grid search and global optimization capabilities, compatibility with 1-D and 3-D Green’s function database formats, customizable data processing, C-accelerated waveform and first-motion polarity misfit functions, and utilities for plotting seismic waveforms and visualizing misfit and likelihood surfaces. Applicability to a range of full- and constrained-moment tensor, point force, and centroid inversion problems is possible via a documented application programming interface, accompanied by example scripts and integration tests. We demonstrate the software using three different types of seismic events: (1) a 2009 intraslab earthquake near Anchorage, Alaska; (2) an episode of the 2021 Barry Arm landslide in Alaska; and (3) the 2017 Democratic People’s Republic of Korea underground nuclear test. With these events, we illustrate the well-known complementary character of body waves, surface waves, and polarities for constraining source parameters. We also convey the distinct misfit patterns that arise from each individual data type, the importance of uncertainty quantification for detecting multimodal or otherwise poorly constrained solutions, and the software’s flexible, modular design.

58 GEOSCIENCES↗

Laboratory investigation of hydraulic fracturing in granitic rocks using active and passive seismic monitoring

SUMMARY Knowledge of the fracturing processes can be important for the optimization of pressurized fluid injection operations in the deep underground rock mass. Active and passive seismic monitoring techniques have been used in the field for tracking or mapping the propagating hydraulic fracture. Although both these monitoring techniques provide valuable information about the generated fracture network, it is difficult for either technique to comprehensibly identify the different processes associated with hydraulic fracturing. The combined active and passive monitoring has the potential for better characterization of the complex hydraulic fracturing phenomena. In this study, laboratory hydraulic fracturing experiments with combined active and passive seismic monitoring were conducted on true triaxially loaded Barre granite cubes with different fluid injection rates. The seismic inelastic fracturing was detected by 16 passive acoustic emission sensors, where 3678 and 2370 seismic source events were detected for the high and low injection rate experiments, respectively. For active monitoring, strong variations in the attributes of signals were observed which were transmitted through four source–receiver pairs, placed both perpendicular and parallel to the generated hydraulic fracture. Positive velocity changes were observed for active sensor pairs with ray paths passing through the generated hydraulic fracture indicating fluid permeation, whereas isolated dry deformation was characterized by a slight but permanent velocity decrease. Compared to velocity, the energy of the active signals was 1–2 orders of magnitude more sensitive to different hydraulic fracturing processes. However, the sensitivity and signatures of the active signal attributes were found to be dependent on the frequency range and direction of ray path with respect to the location of the generated fracture network. Using the coupled evaluation of the active and passive signals we were able to systematically identify various hydraulic fracturing processes including: (1) aseismic deformation, (2) fracture initiation and fluid permeation, (3) pressure build-up, (4) fracture propagation and (5) pressure release and leak-off. The results of this study showed that combining the respective advantages of active and passive seismic techniques and using both of them to monitor the failure processes can facilitate a more comprehensive understanding and better control of the hydraulic stimulations in subsurface operations.

Geochemistry & Geophysics↗

Aerospace technology can be applied to exploration 'back on earth'

Applications of aerospace technology to petroleum exploration are described. Attention is given to seismic reflection techniques, sea-floor mapping, remote geochemical sensing, improved drilling methods and down-hole acoustic concepts, such as down-hole seismic tomography. The seismic reflection techniques include monitoring of swept-frequency explosive or solid-propellant seismic sources, as well as aerial seismic surveys. Telemetry and processing of seismic data may also be performed through use of aerospace technology. Sea-floor sonor imaging and a computer-aided system of geologic analogies for petroleum exploration are also considered.

Jaffe, L. D.↗

Real-Time Seismic System for Monitoring, Imaging, and Characterization (RT-SEISMIC)

The goal of this Laboratory Directed Research and Development (LDRD) project was to develop a borehole seismic source and sensor array to enable real-time seismic imaging at scales and conditions relevant to the energy industry including both fossil-energy and geothermal. In FY21 and FY22, we designed, built, and tested both a prototype impulse source module for generating seismic energy and a sensing module for recording ground motions generated by the source module array. A pneumatically driven vibratory source was also designed. The source modules were fabricated with all high temperature components and the team has worked to incorporate the current RT-SEISMIC electronics design into a commercially available, high temperature silicon-on-insulator chip integrated circuit. Several issues were identified during fabrication and lab testing that led to redesign of several system components and subsequent retesting. The final round of testing showed that while metal/graphite-based seals worked quite well for static seals, they were unable to provide an adequate gas seal for dynamic, reciprocating part movements which necessitated a final redesign using Kalrez. This change will result in a continuous temperature rating of approximately 275 degrees C for the system. While a field test of the RT-SEISMIC system was targeted in FY22, due to the extended lab testing and redesign efforts, field testing was not achieved. As a result of this LDRD investment, several sponsors have expressed interest in RT-SEISMIC and we expect to continue towards a field demonstration of the full system in the future.

58 GEOSCIENCES↗

Modeling the near-field ionospheric disturbances during earthquakes

We present the latest development of the Wave Perturbation-Global Ionosphere-Thermosphere Model (WP-GITM), a threedimensional physics-based numerical model for seismic/tsunami-ionosphere coupling via atmospheric acoustic-gravity waves. WP-GITM was previously applied to simulate the ionospheric perturbations resulted from the epicentral crustal movement by assuming spherical acoustic-gravity waves originated from a point source, which was specified by the seismic measurement at a single location nearby the epicenter. In this work, we extend WP-GITM to include the effects of Rayleigh surface waves and adapt WP-GITM to utilize seismic measurements from more than one location in the attempt to capture the radiation pattern of the seismic source. We apply the new WP-GITM to model the near-field co-seismic ionospheric perturbations during the 16 September 2015 Illapel earthquake. The comparison between the simulated ionospheric total electron content perturbations and the GPS observations shows promising results.

Ravanelli, Michela↗

Monitoring Fracture Hydromechanical Evolution in the Lab and Field Using Unsupervised Metric Learning

Fractures evolve in time through thermal‐hydraulic‐mechanical‐chemical (THMC) processes that alter their long‐range hydraulic transport properties and modify subsurface behavior and activities. The location of subsurface fractures makes it necessary to use remote sensing techniques such as passive or active seismic monitoring for fracture characterization. In this paper, we develop a machine learning approach to monitor the evolution of fracture properties using passive seismic sources in a laboratory setting and using active seismic monitoring from the Sanford Underground Research Facility in Lead, South Dakota, at a depth of 1.25 km in amphibolite rock during stimulation of natural fractures as well as during induced fracturing. The unsupervised metric learning technique applies tandem neural networks (twin (Siamese) or triplet) with contrastive loss and adaptive margins to track slowly varying systems for which class or similarity labels are not available. The approach adopts locality‐sensitive hashing to divide time‐ordered contiguous data into an arbitrary number of pseudo‐classes. Contrastive‐loss training with many hash bins generates an evolving latent‐space trajectory. This approach enables unsupervised metric learning for seismic data stacks under the condition of contiguous state sampling and slowly varying fracture properties. The displacement discontinuity theory provides a mechanistic foundation for the fracture‐dependent trajectories that are related to relaxation of fractures with time‐dependent specific stiffness responding to changes in stress or fluid saturation.

02 PETROLEUM↗

Relative Source Time Functions, Spectral Ratios, and Near‐Source Spallation in the Source Physics Experiment Phase I Chemical Explosions

The Source Physics Experiment (SPE) Phase I was composed of six chemical explosions at the Nevada National Security Sites (NNSS) between 2011 and 2016. The experiment provided a robust set of dense, local to regional distance observations suitable for characterizing shallow chemical explosions located within the same borehole. We investigate the time-dependent source processes of each chemical explosion using Relative Source Time Function (RSTF) estimation, derived from five linear seismometer profiles located between 100 and 2000 m of the shot location. The RSTF estimate provide a detailed characterization of each chemical explosion's time-history. Subsequent modeling efforts suggest these measurements may be useful for precise characterization of explosion processes and spallation. RSTF estimation allows us to expand our understanding of the SPE Phase I chemical explosion series. Application of these techniques contributes to our understanding of explosion source physics and facilitates future applications to additional seismic source experiments and characterization of explosion phenomenology.

58 GEOSCIENCES↗

Machine learning for seismic low-frequency extrapolation

The cycle-skipping problem that plagues full waveform inversion (FWI) can be at least partially mitigated if low frequencies (which encode the kinematics of wave propagation in seismic data) are recorded. However, seismic sources and receivers are band-limited, so seismic data does not generally include signals down to 0 Hz. To improve our ability to solve the seismic inverse problem, one can synthesize this missing low-frequency (LF) content from the recorded high-frequency (HF) data using machine learning (ML) models. Deep learning models such as convolutional neural networks (CNNs) demonstrate impressive ability to perform low frequency extrapolation. However, such models require powerful hardware (GPU machines) and careful training. We assess the extrapolation capabilities of three different ML models that do not require GPU machines, namely, random forest, Gaussian process regression and gradient boosting, on both synthetic and real data. Experimental results on two synthetic data sets (generated from a low velocity lens embedded in a homogeneous medium, and the Marmousi model) demonstrate that FWI applied to the extrapolated data consistently improves inversion accuracy relative to FWI applied to the original data sets that do not contain low frequencies. Application of low-frequency extrapolation to real data from the Northwest Shelf of Australia demonstrates that tree-based ML models such as gradient boosting can outperform CNNs in terms of both accuracy and computational cost on non-GPU architectures.

58 GEOSCIENCES↗

Energy, frequency, and distance of moonquakes at the Apollo 17 site

'Thermal' moonquakes have been detected at the Apollo 17 site. A typical event releases about 1 to 10 million ergs of energy. The annual seismic energy release for the events observed in the Taurus-Littrow Valley is estimated at 100 billion ergs. Such small events would not cause impossibly rapid erosion on either the North Massif or within the craters of the Central Cluster. Seismic events become increasingly frequent after sunrise and reach a maximum at sunset. The largest events, however, occur most commonly near lunar noon. Rise times of the seismic signals, after being calibrated by the well-located Lunar Seismic Profiling Experiment (LSPE) seismic sources, were used to estimate distances to the seismic events. Most seismic events (about 90%) appear to occur within 2.5 km of the seismometer array.

Cooper, M. R.↗