Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic attenuation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Seismic attenuation system for a nuclear reactor

A system for attenuating seismic forces includes a reactor pressure vessel containing nuclear fuel and a containment vessel that houses the reactor pressure vessel. Both the reactor pressure vessel and the containment vessel include a bottom head. Additionally, the system includes a base support to contact a support surface on which the containment vessel is positioned in a substantially vertical orientation. An attenuation device is located between the bottom head of the reactor pressure vessel and the bottom head of the containment vessel. Seismic forces that travel from the base support to the reactor pressure vessel via the containment vessel are attenuated by the attenuation device in a direction that is substantially lateral to the vertical orientation of the containment vessel.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Seismic attenuation system for a nuclear reactor

A system for attenuating seismic forces includes a reactor pressure vessel containing nuclear fuel and a containment vessel that houses the reactor pressure vessel. Both the reactor pressure vessel and the containment vessel include a bottom head. Additionally, the system includes a base support to contact a support surface on which the containment vessel is positioned in a substantially vertical orientation. An attenuation device is located between the bottom head of the reactor pressure vessel and the bottom head of the containment vessel. Seismic forces that travel from the base support to the reactor pressure vessel via the containment vessel are attenuated by the attenuation device in a direction that is substantially lateral to the vertical orientation of the containment vessel.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

A viscoelastic model for seismic attenuation using fractal mechanical networks

SUMMARY Seismic attenuation (quantified by the quality factor Q) has a significant impact on the seismic waveforms, especially in the fluid-saturated rocks. This dissipative process can be phenomenologically represented by viscoelastic models. Previous seismological studies show that the Q value of Earth media exhibits a nearly frequency-independent behaviour (often referred to as constant-Q in literature) in the seismic frequency range. Such attenuation can be described by the mathematical Kjartansson constant-Q model, which lacks of a physical representation in the viscoelastic sense. Inspired by the fractal nature of the pore fluid distribution in patchy-saturated rocks, here we propose two fractal mechanical network (FMN) models, that is, a fractal tree model and a quasi-fractal ladder model, to phenomenologically represent the frequency-independent Q behaviour. As with the classic viscoelastic models, the FMN models are composed of mechanical elements (spring and dashpots) arranged in different hierarchical patterns. A particular parametrization of each model can produce the same complex modulus as in the Kjartansson model, which leads to the constant-Q. Applying the theory to several typical rock samples, we find that the seismic attenuation signature of these rocks can be accurately represented by either one of the FMN models. Besides, we demonstrate that the ladder model in particular exhibits the realistic multiscale fractal structure of the saturated rocks. Therefore, the FMN models as a proxy could provide a new way to estimate the microscopic rock structure property from macroscopic seismic attenuation observation.

Geochemistry & Geophysics↗
DOE OSTI2017

Spatiotemporal changes of seismic attenuation caused by injected CO 2 at the Frio-II pilot site, Dayton, TX, USA

A continuous active source seismic monitoring data set was collected with crosswell geometry during CO 2 injection at the Frio-II brine pilot, near Liberty, TX. Previous studies have shown that spatiotemporal changes in the P wave first arrival time reveal the movement of the injected CO 2 plume in the storage zone. To further constrain the CO 2 saturation, particularly at higher saturation levels, we investigate spatial-temporal changes in the seismic attenuation of the first arrivals. The attenuation changes over the injection period are estimated by the amount of the centroid frequency shift computed by local time-frequency analysis. We observe that (1) at receivers above the injection zone seismic attenuation does not change in a physical trend; (2) at receivers in the injection zone attenuation sharply increases following injection and peaks at specific points varying with distributed receivers, which is consistent with observations from time delays of first arrivals; then, (3) attenuation decreases over the injection time. The attenuation change exhibits a bell-shaped pattern during CO 2 injection. Under Frio-II field reservoir conditions, White's patchy saturation model can quantitatively explain both the P wave velocity and attenuation response observed. We have combined the velocity and attenuation change data in a crossplot format that is useful for model-data comparison and determining patch size. Our analysis suggests that spatial-temporal attenuation change is not only an indicator of the movement and saturation of CO 2 plumes, even at large saturations, but also can quantitatively constrain CO 2 plume saturation when used jointly with seismic velocity.

58 GEOSCIENCES↗
DOE OSTI2017

Global seismic attenuation imaging using full-waveform inversion: a comparative assessment of different choices of misfit functionals

We present the results of synthetic tests that aim at evaluating the relative performance of three different definitions of misfit functionals in the context of 3-D imaging of shear wave attenuation in the earth’s upper mantle at the global scale, using long-period full-waveform data. The synthetic tests are conducted with simple hypothetical upper-mantle models that contain Q μ anomalies centred at different depths and locations, with or without additional seismic velocity anomalies. To build synthetic waveform data sets, we performed simulations of 50 events in the hypothetical (target) models, using the spectral element method, filtered in the period range 60–400 s. The selected events are chosen among 273 events used in the development of radially anisotropic model SEMUCB-WM1 and recorded at 495 stations worldwide. The synthetic Z-component waveforms correspond to paths and time intervals (fundamental mode and overtone Rayleigh waves) that exist in the real waveform data set. The inversions for shear attenuation structure are carried out using a Gauss–Newton optimization scheme in which the gradient and Hessian are computed using normal mode perturbation theory. The three different misfit functionals considered are based on time domain waveform (WF) and waveform envelope (E-WF) differences, as well as spectral amplitude ratios (SA), between observed and predicted waveforms. We evaluate the performance of the three misfit functional definitions in the presence of seismic noise and unresolved S-wave velocity heterogeneity and discuss the relative importance of physical dispersion effects due to 3-D Q μ structure. We observed that the performance of WF is poorer than the other two misfit functionals in recovering attenuation structure, unless anelastic dispersion effects are taken into account in the calculation of partial derivatives. WF also turns out to be more sensitive to seismic noise than E-WF and SA. Overall, SA performs best for attenuation imaging. Our tests show that it is important to account for 3-D elastic effects (focusing) before inverting for Q μ . Additionally, we show that including high signal-to-noise ratio overtone wave packets is necessary to resolve Q μ structure at depths greater than 250 km.

58 GEOSCIENCES↗
NASA NTRS1980

Seismic velocity, attenuation and rheology of the upper mantle

Seismic and rheological properties of the upper mantle in the vicinity of the low-velocity zone are expressed in terms of relaxation by dislocation glide. Dislocation bowing in the glide plane explains seismic velocities and attenuation. Climbing at higher stresses for longer periods of time give the observed viscosity, and explain the low velocity and high temperature attenuation found at seismic frequencies. Due to differing parameters, separate terms for thermal, seismic and rheological lithospheres are proposed. All three lithospheres, however, are related and are functions of temperature, and must be specified by parameters such as period, stress, and stress duration.

Anderson, D. L.↗
DOE OSTI2024

Experimental Observation of a New Attenuation Mechanism in hcp ‐Metals That May Operate in the Earth's Inner Core

Abstract Seismic observations show the Earth's inner core has significant and unexplained variation in seismic attenuation with position, depth and direction. Interpreting these observations is difficult without knowledge of the visco‐ or anelastic dissipation processes active in iron under inner core conditions. Here, a previously unconsidered attenuation mechanism is observed in zinc, a low pressure analog ofhcp‐iron, during small strain sinusoidal deformation experiments. The experiments were performed in a deformation‐DIA combined with X‐radiography, at seismic frequencies (∼0.003–0.1 Hz), high pressure and temperatures up to ∼80% of melting temperature. Significant dissipation (0.077 ≤ Q −1 (ω) ≤ 0.488) is observed along with frequency dependent softening of zinc's Young's modulus and an extremely small activation energy for creep (⩽7 kJ mol −1 ). In addition, during sinusoidal deformation the original microstructure is replaced by one with a reduced dislocation density and small, uniform, grain size. This combination of behavior collectively reflects a mode of deformation called “internal stress superplasticity”; this deformation mechanism is unique to anisotropic materials and activated by cyclic loading generating large internal stresses. Here we observe a new form of internal stress superplasticity, which we name as “elastic strain mismatch superplasticity.” In it the large stresses are caused by the compressional anisotropy. If this mechanism is also active inhcp‐iron and the Earth's inner‐core it will be a contributor to inner‐core observed seismic attenuation and constrain the maximum inner‐core grain‐size to ≲10 km.

Geochemistry & Geophysics↗
DOE OSTI2024

A Regional Phase Amplitude Model of 2-D Attenuation for North America

We analyzed seismic attenuation patterns across the North America using over 70,000,000 Lg wave amplitudes recorded by the various networks at frequencies of 0.1-32 Hz. Our inversion solved for laterally varying attenuation, site terms, moments, and apparent stress following Phillips et al., (2016). The inversion was anchored by independently constrained: corner frequencies (via coda spectral ratios) to control the attenuation-stress tradeoff and moment measurements of teleseismic (GCMT, USGS) and regional (St. Louis University and UC Berkeley) earthquakes to provided absolute scaling. The quality factor (Q) shows clear regional patterns: low values in coastal, volcanic, and tectonically active regions, and high values in stable areas like the Great Plains and major plateaus throughout North America. These 2-D Q models enable improved regional source characterization, magnitude estimation, and yield determination.

58 GEOSCIENCES↗
NASA NTRS1977

The distribution of seismic velocities and attenuation in the earth

Estimates of the radial distribution of seismic velocities and density and of seismic attenuation within the earth are obtained through inversion of body wave, surface wave, and normal mode data. The effect of attenuation related dispersion on gross earth structure, and on the reliability of eigenperiod identifications is discussed. The travel time baseline discrepancies between body waves and free oscillation models are examined and largely resolved.

Hart, R. S.↗
DOE OSTI2018

Subsurface Structure of Long Valley Caldera Imaged With Seismic Scattering and Intrinsic Attenuation

Here we image seismic intrinsic (Q i -1 ) and scattering (Q s -1 ) attenuations in Long Valley Caldera, California, by analyzing more than 1,700 vertical component waveforms from local earthquakes. Observed energy envelopes are fit to the diffusion model and seismic attenuation images are produced using two-dimensional space weighting functions. Low intrinsic and low scattering attenuation anomalies in the center south of the caldera correspond to the location of an earthquake swarm in 2014. We identify high intrinsic and high scattering attenuation anomalies in the fluid-rich western and eastern areas of the caldera. From a comparison with other geophysical images (magnetotellurics and seismic tomography) we attribute these anomalies to a hydrothermal system (high attenuation). Average to high attenuation values are also observed at the adjacent Mammoth Mountain (southwest of the caldera) and may also have a hydrothermal origin. High intrinsic attenuation at low frequencies to the west of the Hartley Springs Fault may be produced by the magmatic system that produced the Inyo Craters. Finally, seismic attenuation imaging provides insights into subsurface structures that are complementary to velocity and conductivity images.

58 GEOSCIENCES↗
DOE OSTI2024

Advancing attenuation estimation through integration of the Hessian in multiparameter viscoacoustic full-waveform inversion

Accurate seismic attenuation models of subsurface structures not only enhance subsequent migration processes by improving fidelity, resolution, and facilitating amplitude-compliant angle gather generation but also provide valuable constraints on subsurface physical properties. Leveraging full-wavefield information, multiparameter viscoacoustic full-waveform inversion ( Q-FWI) simultaneously estimates seismic velocity and attenuation ( Q) models. However, a major challenge in Q-FWI is the contamination of crosstalk artifacts, where inaccuracies in the velocity model are mistakenly mapped to the inverted attenuation model. While incorporating the Hessian is expected to mitigate these artifacts, the explicit implementation is prohibitively expensive due to its formidable computational cost. In this study, we formulate and develop a Q-FWI algorithm via the Newton-conjugate gradient (CG) framework, where the search direction at each iteration is determined through an internal CG loop. In particular, the Hessian is integrated into each CG step in a matrix-free fashion using the second-order adjoint-state method. We find through synthetic experiments that our Newton-CG Q-FWI significantly mitigates crosstalk artifacts compared with the limited-memory Broyden-Fletcher-Goldfarb-Shanno method and the CG method, albeit with a notable computational cost. In the discussion of several key implementation details, we also determine the significance of the approximate Gauss-Newton Hessian, the second-order adjoint-state method, and the two-stage inversion strategy.

Geochemistry & Geophysics↗
DOE OSTI2023

Compensating attenuation effects in full-waveform inversion with dissipation-dispersion decoupling

Seismic attenuation poses challenges to velocity model building from acoustic/elastic full-waveform inversion (FWI). In particular, when constructing the FWI gradient for velocity inversion, recent studies have indicated that the high-attenuation structure could distort the gradient by damping the amplitudes and shifting the kinematic phases, resulting in an imbalanced update and thus unreliable velocity model. These Q effects are particularly significant in the reflection acquisition geometry due to the “double-damping” issue. Here, we develop a Q-compensated FWI algorithm for constructing a Q-free FWI gradient. By using a recently developed viscoacoustic wave propagator, this compensation can be done conveniently by keeping the dispersion term and flipping the dissipation term in the wave equation when we simulate the forward and time-reversed adjoint wavefields. The resultant gradient obtained by interacting these two wavefields has correct kinematics and Q-free amplitudes. This Q-compensated FWI can balance the update between Q- and no- Q-areas, which we determine using synthetic examples. In addition, we illustrate how to take advantage of the dissipation-dispersion decoupling to determine the anomalous Q value in the Q model building workflow via Q-compensated reverse time migration.

Geochemistry & Geophysics↗
DOE OSTI2022

Ambient Noise Attenuation Tomography Reveals an Asymmetric Damage Zone Across San Jacinto Fault Near Anza, California

Here, we perform seismic attenuation tomography of the shallow structure for the San Jacinto Fault in the Ramona Reservation of southern California. The study uses ambient seismic noise recorded by a linear array of 65 3-C sensors across the fault. We extract amplitude decay information, with uncertainty, from noise interferometry functions. To account for strong heterogeneities in the complex shallow fault zone structure, we apply a frequency-dependent amplitude correction for focusing/defocusing effects using inverted phase velocity maps obtained by solving the transport equation. We then estimate the attenuation structure based on the linear station-triplet method for both Love and Rayleigh waves. The attenuation tomography indicates strong attenuation correlated with known San Jacinto Fault surface traces. The Love wave attenuation tomography reveals an asymmetric damage zone that exists primarily on the fault side with faster seismic velocity, consistent with earthquake ruptures on a fault bimaterial interface with preferred propagation direction to the northwest.

58 GEOSCIENCES↗
DOE OSTI2020

Modulation of seismic noise near the San Jacinto fault in southern California: origin and observations of the cyclical time dependence and associated crustal properties

SUMMARY We examine the cyclic amplitude variation of seismic noise recorded by continuous three-component broad-band seismic data with durations spanning 91–713 d (2008–2011) from three different networks: Anza seismic network, IDA network and the transportable seismic array. These stations surround the San Jacinto Fault Zone (SJFZ) in southern California. We find the seismic noise amplitudes exhibit a cyclical variation between 0.3 and 7.2 Hz. The high-frequency (≥0.9 Hz) noise variations can be linked to human activity and are not a concern. Our primary interest is signals in the low frequencies (0.3–0.9 Hz), where the seismic noise is modulated by semi-diurnal tidal mode M2. These long-period (low-frequency) variations of seismic noise can be attributed to a temporal change of the ocean waves breaking at the shoreline, driven by ocean tidal loading. We focus on the M2 variation of seismic noise at f = 0.6 Hz, travelling distances of ∼92 km through the crust from offshore California to the inland Anza, California, region. Relative to the shoreline station, data from the inland stations show a phase lag of ∼ –12°, which we attribute to the cyclic change in M2 that can alter crustal seismic attenuation. We also find that for mode M2 at 0.6 Hz, the amplitude variations of the seismic quality factor (Q) depend on azimuth and varies from 0.22 per cent (southeast to northwest) to 1.28 per cent (northeast to southwest) with Q = 25 for Rayleigh waves. We propose the direction dependence of the Q variation at 0.6 Hz reflects the preferred orientation of subfaults parallel to the main faulting defined by the primarily N45°W strike of the SJFZ.

Geochemistry & Geophysics↗