SEARCH · Search NASA
Results for “SEISMIC”
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.
Contributions to Improving Detection Capabilities for Global Seismic Networks, Characterization of Socorro Region Seismicity, and Modeling of Seismic Wave Propagation
Explore the source record for details and available documents.
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.
Isotropic seismic radiation from rock damage and dilatancy
SUMMARY We present semi-analytical and numerical results on isotropic radiation expected to be generated by seismic sources with damaged dilatant materials that sustain during failure an overall shear deformation. The results are derived for inclusions with different stiffness components, including those responsible for shear-volumetric coupling in dilatant material, and different geometries, in a surrounding elastic solid. The obtained relations indicate that the isotropic seismic potency and seismic moment components are proportional to the shear potency components multiplied by the aspect ratio between the thickness and length of the source zone. For a typical rupture zone with a thin tabular geometry, the expected small aspect ratio explains the general seismological observations of small isotropic radiation. However, isotropic radiation induced by dynamic dilation for short and thick sources is more pronounced and may reach for spherical sources sustaining shear deformation 18 per cent of the shear components. Sources that include also tensile components and explosions can produce larger isotropic radiation. The results provide rigorous relations that can be used to estimate different source properties in numerical simulations, experiments and field studies.
The effects of earth model uncertainty on the inversion of seismic data for seismic source functions
SUMMARY We use Monte Carlo simulations to explore the effects of earth model uncertainty on the estimation of the seismic source time functions that correspond to the six independent components of the point source seismic moment tensor. Specifically, we invert synthetic data using Green’s functions estimated from a suite of earth models that contain stochastic density and seismic wave-speed heterogeneities. We find that the primary effect of earth model uncertainty on the data is that the amplitude of the first-arriving seismic energy is reduced, and that this amplitude reduction is proportional to the magnitude of the stochastic heterogeneities. Also, we find that the amplitude of the estimated seismic source functions can be under- or overestimated, depending on the stochastic earth model used to create the data. This effect is totally unpredictable, meaning that uncertainty in the earth model can lead to unpredictable biases in the amplitude of the estimated seismic source functions.
Detection of random noise and anatomy of continuous seismic waveforms in dense array data near Anza California
We develop a methodology to separate continuous seismic waveforms into random noise (RN), not random noise (NRN) produced by earthquakes, wind, traffic and other sources of ground motions, and an undetermined mixture of signals. The analysis is applied to continuous records from a dense seismic array on the San Jacinto fault zone. To detect RN signals, we cut hourly waveforms into non-overlapping 1 s time windows and apply cross-correlations to separate RN candidates from outliers. The cross-correlation coefficients between different RN candidates fall into a tight range (i.e. 0.09–0.35), while cross-correlation coefficients of RN candidates with NRN signals (e.g. seismic or air-traffic events) are lower. The amplitude spectra of RN candidates have a well-defined level, while the amplitude spectra of other signals deviate from that level. Using these properties, we examine the amplitude spectra of moving time windows and cross-correlation coefficients with RN templates in each hour. The hourly RN is quasi-stationary and the results cluster tightly in the parameter space of cross-correlation coefficients and L2 norm deviations from the mean spectra of RN candidates. Time windows with parameters in this tight cluster are identified as RN, windows that deviate significantly from the RN cluster are identified as NRN and windows with values in between are identified as mixed signals. Several iterations on each hourly data are used to update and stabilize the selection of RN templates and mean noise spectra. For the days examined, the relative fractions of RN, NRN and mixed signals in local day (night) times are about 26 (42 percent), 40 (33 percent) and 34 percent (25 percent), respectively.
Enhancing trap and fault seal analyses by integrating observations from HR3D seismic data with well logs and conventional 3D seismic data, Texas inner shelf
An understanding of trap and fault seal quality is critical for assessing hydrocarbon prospectivity. To achieve this, modern analytical techniques leverage well data and conventional industry-standard 3D seismic data to evaluate the trap, and any faults displacing the reservoir and top seal intervals. Above all, geological interpretation provides the framework of trap and fault seal analyses, but can be hindered by the data resolution, quality and acquisition style of the conventional seismic data. Furthermore, limiting the analysis to only the petroleum system at depth may lead to erroneous perceptions because interpreting overburden features, such as shallow faults or gas chimneys, can provide valuable observations with respect to container performance, and can to help validate trap and fault seal predictions. Here, a supplement to conventional 3D data are high-resolution 3D seismic (HR3D) data, which provide detailed images of the overburden geology. This study utilizes an HR3D seismic volume in the San Luis Pass area of the Texas inner shelf, where shallow fault tips and a sizeable gas chimney are interpreted over an unsuccessful hydrocarbon prospect. Static post-drill fault seal and trap analyses suggest that the primary fault displacing the structural closure could have withheld columns of gas c. 100 m high, but disagree with our HR3D seismic interpretations and dry-well analyses. From our results, we hypothesize that tertiary gas migration through fault conduits reduced the hydrocarbon column in the prospective Early Miocene reservoir, and may have resulted from continued movement along the intersecting faults. Overall, this study reinforces the importance of understanding the overburden geology and geohistory of faulted prospects, and demonstrates the utility of pre-drill HR3D acquisition when conducting trap and fault seal analyses.
Integrating distributed acoustic sensing, borehole 3C geophone array, and surface seismic array data to identify long-period long-duration seismic events during stimulation of a Marcellus Shale gas reservoir
Microseismic monitoring by downhole geophones, surface seismic, fiber-optic distributed acoustic sensing (DAS), and distributed temperature sensing (DTS) observations were made during the hydraulic fracture stimulation of the MIP-3H well in the Marcellus Shale in northern West Virginia. DAS and DTS data measure the fiber strain and temperature, respectively, along a fiber-optic cable cemented behind the casing of the well. The presence of long-period long-duration (LPLD) events is evaluated in the borehole geophones, DAS data, and surface seismic data of one of the MIP-3H stimulated stages. LPLD events are generally overlooked during the conventional processing of microseismic data, but they represent significant nonbrittle deformation produced during hydraulic fracture stimulation. In a single stage that was examined, 160 preexisting fractures and two faults of suboptimal orientation are noted in the image logs. We identified two low-frequency ([Formula: see text]) events of large temporal duration (tens of seconds) by comparing the surface seismic data, borehole geophone data, and DAS amplitude spectra of one of the MIP-3H stages. Spectrograms of DAS traces in time and depth reveal that the first low-frequency event might be an injection noise that has footprints on all DAS channels above the stimulated stage. However, the surface seismic array indicates an LPLD event concurrent with the first low-frequency event on DAS. The second LPLD event on DAS data and surface seismic data is related to a local deformation and does not have footprints on all DAS channels. The interpreted events have duration less than 100 s with frequencies concentrated below 10 Hz, and are accompanied by microseismic events.
Integration of seismic-pressure-petrophysics inversion of continuous active-seismic monitoring data for monitoring and quantifying CO 2 plume (Final Report)
The overall objective of this project is to develop and validate an integrated package of joint seismic-pressure-petrophysics inversion (jSPPI) of continuous active-source seismic monitoring dataset capable of providing real-time monitoring of CO 2 plume during geologic carbon sequestration (GCS). The three specific developments include: (a) the methodologies for fast seismic full waveform inversion of continuous active source seismic monitoring, (CASSM) datasets for simultaneously estimating velocity and attenuation, and with data assimilation; (b) joint Bayesian petrophysical inversion of seismic models and pressure data for providing and updating CO 2 saturation models; (c) the methods using multiple datasets including (Crainfield and Frio-II borehole) synthetic, laboratory, and field CASSM datasets. The outcomes of jSPPI include (a) a workflow for processing CASSM data, (b) Bayesian inversion algorithms using CASSM data and pressure response data, and (c) integration with data assimilation algorithms for continuously updating site-specific models used for prediction and reservoir management. The validation of joint FWI will be conducted using synthetic models based on the Cranfield and Frio experiments as well as field CASSM datasets collected as part of the Frio-II pilot injection. To quantify and map the mass and distribution of CO 2 (saturation), we will jointly invert velocity and attenuation measurements from the FWI with a Bayesian approach using a rock physics model for attenuation (e.g., White’s attenuation model with two selected patch sizes (White, 1976; Dutta and Seriff, 1979)). The Bayesian inversion will be applied to each time step in the CASSM survey in an updating scheme, which integrates with an ensemble of reservoir simulations at each step. A more complete experimental validation dataset will be collected as part of a mesoscale (2-3 m) gas-CO 2 injection experiment utilizing a higher frequency version of the CASSM system developed for laboratory studies; the integrated inversion will be demonstrated using this dataset which will provide both a dense geometry as well as more precise secondary confirmation measurements (e.g. saturation) typically not available in the field. The resulting real-time map of CO 2 saturation is able to provide a deeper scientific understanding of the complex, time-varying dynamics of subsurface fluid flow migration path as well as the rapid detection of CO 2 leakage hazards.
Apollo 14 and 16 Active Seismic Experiments, and Apollo 17 Lunar Seismic Profiling
Seismic refraction experiments were conducted on the moon by Apollo astronauts during missions 14, 16, and 17. Seismic velocities of 104, 108, 92, 114 and 100 m/sec were inferred for the lunar regolith at the Apollo 12, 14, 15, 16, and 17 landing sites, respectively. These data indicate that fragmentation and comminution caused by meteoroid impacts has produced a layer of remarkably uniform seismic properties moonwide. Brecciation and high porosity are the probable causes of the very low velocities observed in the lunar regolith. Apollo 17 seismic data revealed that the seismic velocity increases very rapidly with depth to 4.7 km/sec at a depth of 1.4 km. Such a large velocity change is suggestive of compositional and textural changes and is compatible with a model of fractured basaltic flows overlying anorthositic breccias. 'Thermal' moonquakes were also detected at the Apollo 17 site, becoming increasingly frequent after sunrise and reaching a maximum at sunset. The source of these quakes could possibly be landsliding.
Sequences of seismic and aseismic slip on bimaterial faults show dominant rupture asymmetry and potential for elevated seismic hazard
Here we perform numerical simulations of sequences of earthquake and aseismic slip on planar rate and state faults separating dissimilar material within the 2-D plane strain approximation. We resolve all stages of the earthquake cycle from aseismic slip to fast ruptures while incorporating full inertia effects during seismic event propagation. We show that bimaterial coupling results in favorable nucleation site and subsequent asymmetric rupture propagation. We demonstrate that increasing the material contrast enhances this asymmetry leading to higher slip rates and normal stress drops in the preferred rupture propagation direction. The normal stress drop, induced by the bimaterial effect, leads to strong dynamic weakening of the fault and may destabilize the creeping region on a heterogeneous rate and state fault, resulting in extended rupture propagation. Such rupture penetration into creeping patches may lead to more frequent opening of earthquake gates, causing increased seismic hazard. Furthermore, bimaterial coupling may lead to irregular seismicity pattern in terms of event length, peak slip rates, and hypocenter location, depending on the properties of the creeping patches bordering the seismogenically active part of the fault. Our results highlight robust characteristics of bimaterial interfaces that persist over long sequence of events and suggest the need for further exploration of the role of material contrast in earthquake physics and models of seismic hazard.
Analysis of surface and seismic sources in dense array data with match field processing and Markov chain Monte Carlo sampling
Not provided.
Seismicity Rate Surge on Faults after Shut–in: Poroelastic Response to Fluid Injection
Injection of large amounts of fluid into the subsurface alters the states of pore pressure and stress in the formation, potentially inducing earthquakes. Increase in the seismicity rate after shut–in is often observed at fluid–injection operation sites, but mechanistic study of the rate surge has not been investigated thoroughly. Considering full poroelastic coupling of pore pressure and stress, the earthquake occurrence after shut–in can be driven by two mechanisms: (1) post shut–in diffusion of pore pressure into distant faults and (2) poroelastic stressing caused by fluid injection. Interactions of these mechanisms can depend on fault geometry, hydraulic and mechanical properties of the formation, and injection operation. In this work, a 2D aerial view of the target reservoir intersected by strike–slip basement faults is used to evaluate the impact of injection–induced pressure buildup on seismicity rate surge. A series of sensitivity tests are performed by considering the variation in (1) permeability of the fault zone, (2) locations and the number of faults with respect to the injector, and (3) well operations with time–dependent injection rates. Lower permeability faults have higher seismicity rates than more permeable faults after shut–in due to delayed diffusion and poroelastic stressing. Hydraulic barriers, depending on their relative location to injection, can either stabilize or weaken a conductive fault via poroelastic stresses. In conclusion, gradual reduction of the injection rate minimizes the coulomb stress change and the least seismicity rates are predicted due to slower relaxation of coupling–induced compression as well as pore–pressure dissipation.
Pre-Sinkhole Seismicity at the Napoleonville Salt Dome: Implications for Local Seismic Monitoring of Underground Caverns
The formation of a sinkhole at the Napoleonville salt dome located near Bayou Corne in Assumption Parish, Louisiana, was preceded by an increase in seismicity, as observed 10 km away by a continuously recording broadband seismic sensor. The seismicity rate in the days prior to the sinkhole formation contained hundreds of long-period events, some of which could be detected using a power detector (short-term to long-term energy ratio) tuned to the duration and frequency content of the seismicity. A characteristic of some of these events is a resonance near 0.4 Hz, which could be due to either an acoustic resonance within the brine-filled cavity or possibly due to unsteady fluid flow through a crack. In this work, a more exhaustive catalog was obtained using some of the events detected by the power detector as templates in a correlation detector. Furthermore, the full pre-sinkhole sequence of detections could offer lessons in remote monitoring of underground caverns.
Recovery and Calibration of Legacy Underground Nuclear Test Seismic Data from the Leo Brady Seismic Network
The Leo Brady Seismic Network (LBSN, originally the Sandia Seismic Network) was established in 1960 by Sandia National Laboratories to monitor underground nuclear tests (UGTs) at the Nevada National Security Site (NNSS, formerly named the Nevada Test Site). The LBSN has been in various configurations throughout its existence, but it has generally been comprised of four to six stations at regional distances (~150–400 km) from the NNSS with approximately evenly spaced azimuthal coverage. Between 1962 and the end of nuclear testing in 1992, the LBSN—and a sister network operated by Lawrence Livermore National Laboratories—was the most comprehensive United States source of regional seismic data of UGTs. Approximately 75% of all UGTs performed by the United States occurred in the predigital era. At that time, LBSN data were transmitted as frequency-modulated (FM) audio over telephone lines to a central location and recorded as analog waveforms on high-fidelity magnetic audio tapes. These tapes have been in dry temperature-stable storage for decades and contain the sole record of this irreplaceable data; full waveforms of LBSN-recorded UGTs from this era were not routinely digitized or otherwise published. We have developed a process to recover and calibrate data from these tapes. First, we play back and digitize the tapes as audio. Next, we demodulate the FM “audio” into individual waveforms. We then estimate the various instrument constants through careful measurement of “weight-lift” tests performed prior to each UGT on each instrument. Finally, these coefficients allow us to scale and shape the derived instrument response of the seismographs and compute poles and zeros. Finally, the result of this process is a digital record of the recorded seismic ground motion in a modern data format, stored in a searchable database. To date, we have digitized tapes from 592 UGTs.
Investigations of acoustic-seismic effects at long range - Early-arriving seismic waves from Apollo 16
A reasonably comprehensive technical effort is described dealing with the investigations of acoustically generated seismic waves of Apollo 16 and Apollo 17 origin along the eastern seabord of the United States. This expanded effort is a continuation of earlier, rather successful detections of rocket-generated seismic disturbances on Skidaway Island, Georgia. The more recent effort has yielded few positive results other than a recording of an early-arriving seismic wave from Apollo 16 that was detected in Jacksonville. Evaluation of the negative results obtained in the Fort Monmouth area, with earlier studies of infrasound, local weather conditions, and geology, could be advantageous in the process of trying to gain a better insight into the acoustic-seismic resonance mechanism requiring phase-velocity matching at the atmosphere-ground interface.
Estimation of the Seismic Moment Release Rate of Mars from InSight Seismic Data
Seismicity models for Mars usually estimate the long-term average annual seismic moment rate, and also the average annual event rate. This holds for estimations based on geological evidence (Golombek et al., 1992, Golombek, 2002, Taylor et al., 2013) as well as for models based on thermal evolution and cooling of the Martian interior (Phillips, 1991, Knapmeyer et al., 2006, Plesa et al., 2018). All studies are compatible with the conclusion based on the non-observation of any unambiguous event by Viking (Anderson et al., 1977, Goins & Lazarewicz, 1979) that Martian seismicity lies somewhere between that of the Moon and that of the Earth. We developed tools to derive reasonable estimations of the annual seismic moment rate from a number of events as small as one, provided that the observed events are beyond the global completeness threshold for observable events. Numerical tests as well as evaluation of terrestrial data shows the feasibility of the approach.