Search NASA⌕ Search

SEARCH · Search NASA

Results for “Velocity Model”

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 163 records · Page 9

Making Invisible Visible: Data-Driven Seismic Inversion With Spatio-Temporally Constrained Data Augmentation

Deep learning and data-driven approaches have shown great potential in scientific domains. The promise of data-driven techniques relies on the availability of a large volume of high-quality training datasets. Due to the high cost of obtaining data through expensive physical experiments, instruments, and simulations, data augmentation techniques for scientific applications have emerged as a new direction for obtaining scientific data recently. However, existing data augmentation techniques originating from computer vision yield physically unacceptable data samples that are not helpful for the domain problems that we are interested in. In this article, we develop new data augmentation techniques based on convolutional neural networks. Specifically, our generative models leverage different physics knowledge (such as governing equations, observable perception, and physics phenomena) to improve the quality of the synthetic data. To validate the effectiveness of our data augmentation techniques, we apply them to solve a subsurface seismic full-waveform inversion using simulated CO 2 leakage data. Our interest is to invert for subsurface velocity models associated with very small CO 2 leakage. We validate the performance of our methods using comprehensive numerical tests. Here via comparison and analysis, we show that data-driven seismic imaging can be significantly enhanced by using our data augmentation techniques. Particularly, the imaging quality has been improved by 15% in test scenarios of general-sized leakage and 17% in small-sized leakage when using an augmented training set obtained with our techniques.

58 GEOSCIENCES↗

High-resolution Bayesian spatial autocorrelation (SPAC) quasi-3-D Vs model of Utah FORGE site with a dense geophone array

SUMMARY We expand the application of spatial autocorrelation (SPAC) from typical 1-D Vs profiles to quasi-3-D imaging via Bayesian Monte Carlo inversion (BMCI) using a dense nodal array (49 nodes) located at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site. Combinations of 4 and 9 geophones in subarrays provide for 36 and 25 1-D Vs profiles, respectively. Profiles with error bars are determined by calculating coherency functions that fit observations in a frequency range of 0.2–5 Hz. Thus, a high-resolution quasi-3-D Vs model from the surface to 2.0 km depth is derived and shows that surface-parallel sedimentary strata deepen to the west, consistent with a 3-D seismic reflection survey. Moreover, the resulting Vs profile is consistent with a Vs profile derived from distributed acoustic sensing (DAS) data located in a borehole at the FORGE site. The quasi-3-D velocity model shows that the base of the basin dips ∼22° to the west and topography on the basement interface coincident with the Mag Lee Wash suggests that the bedrock interface is an unconformity.

Geochemistry & Geophysics↗

Analysis of Seismic Wave Amplification in Sedimentary Basins Using 3D Wavefield Simulations: Nenana Basin, Central Alaska

We investigate seismic wave amplification in Nenana basin, central Alaska, using 3D seismic wavefield simulations. We quantify frequency-dependent amplification by comparing synthetic seismograms among four different seismic velocity models: (a) a tomographic model, (b) a tomographic model with the uppermost 6.5 km replaced by a 1D basin profile, (c) a tomographic model with an embedded ellipsoidal basin, and (d) a tomographic model with an embedded realistic basin. For each model we perform wavefield simulations accurate up to 2 Hz for a set of 10 different earthquakes, which provide a range of frequencies and incident angles for waves interacting with the basin. Average amplification ratios are four on the horizontal components and seven on the vertical component. The dominant influence of the amplification is the 3D geometry of the basin, not the slow-velocity profile near the surface. Our synthetic amplification ratios are in general agreement with amplification estimates obtained from 14 stations that recorded the same set of earthquakes. Our approach offers a general strategy for documenting frequency-dependent basin amplification for a region with realistic basin structures and local earthquakes.

Geosciences↗

Modeling Induction Stirring and Particle Tracking in Molten Uranium

Two independent numerical models have been developed to simulate the behavior of carbon impurities in molten uranium metal. Informed by experimental parameters, one model was created using the commercial software Star-CCM+ and compared with another developed using open-source codes, including OpenFOAM, Finite Element Method Magnetics (FEMM) and a First Passage Kinetic Monte Carlo (FPKMC) approach. The target experimental system features a 404g uranium metal charge containing an average carbon concentration of 139 ppm which was melted in a vacuum induction furnace at 1400° C then resolidified. The microstructures of the uranium and its impurities before and after melting have been characterized and reported separately. Prior to simulating the uranium-carbon system described, the numerical models were validated using a previously published nonradioactive experimental system to ensure agreement with expected output values. Focus has been placed on modeling velocity fields under induction stirring and impurity particle trajectories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Correction of lunar seismograms for instrumental and near-surface effects and constraints on the velocity structure of the lunar interior

Long-period lunar seismograms were studied with the aim of identifying consistent sets of direct shear and secondary wave arrivals, thus constraining the velocities in the lunar mantle and the depths of the velocity discontinuities. Two velocity models were used to locate the natural impacts and the shallow moonquakes and to obtain the travel time residuals. Seismic sections were made of the radial, transverse, and vertical components of ground motion for impacts, shallow, and deep moonquakes in order to search for consistent sets of secondary wave arrivals. No conclusive set of secondary arrivals could be recognized on the seismic sections and thus the velocities and depths of the velocity discontinuities cannot be severely constrained by secondary arrivals. It is likely that the crust is thinner than 50 km and that a first-order discontinuity separates the upper and lower crust at a depth of between 20 and 30 km.

Horvath, P.↗

Modeling Subsurface Explosions Recorded on a Distributed Fiber Optic Sensor

Fiber optic distributed acoustic sensors (DAS) are becoming a widely used tool for seismic sensing. Here we examine recordings of two subsurface chemical explosions, DAG-1 and DAG-3, each of which was about one metric ton (TNT equivalent), that were recorded from a helical fiber installed in two boreholes 80 m away from the source location. Several clear phases including the initial P wave, a weak S wave, and a surface reflected P wave are observed on the helical DAS data. We estimate a velocity model using arrival times measured from the fiber. The DAS waveform data were compared with colocated accelerometers at specific depths in both frequency and time domains. Furthermore, the spectra of the DAS data matched spectra estimated from the accelerometer records. Comparisons of observed waveform shape between the accelerometer records and the fiber measurements (strain-rate) show reasonable agreement except for the data near the event depth. The DAS data and the accelerometer agreed in relative amplitudes but we had difficulties in matching absolute amplitudes, possibly due to errors in metadata. Synthetic strain-rate waveforms were calculated using a 2D wavenumber algorithm and matched the waveform shape and relative amplitudes. In general, DAS is effective at recording strong ground motions at high spatial density. Comparison of the synthetic seismograms with observed data indicate that the waveforms are not consistent with a pure isotropic explosion source and that the observed S waves originate from very near the source region.

58 GEOSCIENCES↗

Earthquake source mechanisms and transform fault tectonics in the Gulf of California

The source parameters of 19 large earthquakes in the Gulf of California were determined from inversions of long-period P and SH waveforms. The goal was to understand the recent slip history of this dominantly transform boundary between the Pacific and North American plates as well as the effect on earthquake characteristics of the transition from young oceanic to continental lithosphere. For the better recorded transform events, the fault strike is resolved to + or - 4 deg at 90 percent confidence. The slip vectors thus provide important constraints on the direction of relative plate motion. Most centroid depths are poorly resolved because of tradeoffs between depth and source time function. On the basis of waveform modeling, historical seismicity, and other factors, it is appropriate to divide the Gulf into three distinct zones. The difference in seismic character among the three zones is likely the result of differing levels of maturity of the processes of rifting, generation of oceanic crust, and formation of stable oceanic transform faults. The mechanism of an earthquake on the Tres Marias Escarpment is characterized by thrust faulting and likely indicates the direction of relative motion between the Rivera and North American plates. This mechanism requires revision in plate velocity models which predict strike slip motion at this location.

Goff, John A.↗

Spectral and Temporal Characteristics of X-Ray-Bright Stars in the Pleiades

We follow up our deep ROSAT imaging survey of the Pleiades (Stauffer et al. 1994) with an analysis of the spectral and temporal characteristics of the X-ray-bright stars in the Pleiades. Raymond & Smith (1977) one and two-temperature models have been used to fit the position-sensitive proportional counter (PSPC) pulse-height spectra of the dozen or so brightest sources associated with late-type Pleiades members. The best-fit temperatures suggest hot coronal temperatures for K, M, and rapidly rotating G stars, and cooler temperatures for F and slowly rotating G stars. In order to probe the many less X-ray-luminous stars, we have generated composite spectra by combining net counts from all Pleiades members according to spectral type and rotational velocity. Model fits to the composite spectra confirm the trend seen in the individual spectral fits. Particularly interesting is the apparent dependence of coronal temperature on L(sub x)/L(sub bol). A hardness-ratio analysis also confirms some of these trends. The PSPC data have also revealed a dozen or so strong X-ray flares with peak X-ray luminosities in excess of approx. 10(exp 30) ergs/sec. We have modeled the brightest of these flares with a simple quasi-static cooling loop model. The peak temperature and emission measure and the inferred electron density and plasma volume suggest a very large scale flaring event. The PSPC data were collected over a period of approx. 18 months, allowing us to search for source variability on timescales ranging from less than a day (in the case of flares) to more than a year between individual exposures. On approximately year-long timescales, roughly 25% of the late-type stars are variable. Since the Pleiades was also intensively monitored by the imaging instruments on the Einstein Observatory, we have examined X-ray luminosity variations on the 10 yr timescale between Einstein and ROSAT and find that up to 40% of the late-type stars are X-ray variable. Since there is only marginal evidence for increased variability on decade-long timescales, the variability observed on long and short timescales may have a common physical origin.

Gagne, Marc↗

Multi-Scale Seismic Measurements for Site Characterization and CO2 Monitoring in an Enhanced Oil Recovery/Carbon Capture, Utilization, and Sequestration Project, Farnsworth Field, Texas

To address the challenges of climate change, significantly more geologic carbon sequestration projects are beginning. The characterization of the subsurface and the migration of the plume of supercritical carbon dioxide are two elements of carbon sequestration that can be addressed through the use of the available seismic methods in the oil and gas industry. In an enhanced oil recovery site in Farnsworth, TX, we employed three separate seismic techniques. The three-dimensional (3D) surface seismic survey required significant planning, design, and processing, but produces both a better understanding of the subsurface structure and a three-dimensional velocity model, which is essential for the second technique, a timelapse vertical seismic profile, and the third technique, cross-well seismic tomography. The timelapse 3D Vertical Seismic Profile (3D VSP) revealed both significant changes in the reservoir between the second and third surveys and geo-bodies that may represent the extent of the underground carbon dioxide. The asymmetry of the primary geo-body may indicate the preferential migration of the carbon dioxide. The third technique, cross-well seismic tomography, suggested a strong correlation between the well logs and the tomographic velocities, but did not observe changes in the injection interval.

Energy & Fuels↗

Challenges and Potential of Waveform Modeling for Crustal Scale Predictions

Waveform modeling is crucial to improving our understanding of observed seismograms. Forward simulation of wavefields provides quantitative methods of testing interactions between complicated source functions and the propagation medium. Here, we discuss three experiments designed to improve under standing of high frequency seismic wave propagation. First, we compare observed and predicted travel times of crustal phases for a set of real observed earthquakes with calculations and synthetic seismograms. Second, we estimate the frequency content of a series of nearly co-located earthquakes of varying magnitude for which we have a relatively well- known 1D velocity model. Third, we apply stochastic perturbations on top of a 3D tomographic model and qualitatively assess how those variations map to differences in the seismograms. While different in scope and aim, these three vignettes illustrate the current state of crustal scale waveform modeling and the potential for future studies to better constrain the structure of the crust.

58 GEOSCIENCES↗

Regional Earthquake Ground Motion Simulations for Southern California With EQSIM: Insights From the 2008 Chino Hills, 2024 Highland Park, and 2021 Carson Earthquakes

This study presents physics-based, 3D simulations using the EQSIM framework for several earthquakes in the Los Angeles region. The primary objective was to assess the ability of deterministic physics-based ground motion simulations to reproduce the observed motions from historical events. The selected events included the mathematical equation M w 5.4 2008 Chino Hills, the mathematical equation M w 4.4 2024 Highland Park, and the mathematical equation M w 4.3 2021 Carson events. The simulated motions were evaluated by comparing the recorded and simulated seismograms, as well as the Fourier amplitude spectra, across multiple seismic stations. The SCEC 3D velocity model, CVM-S4.26.M01, was used to represent the regional geology, and ground motion simulations were carried out with a resolution of up to 5 Hz. The results indicate that the simulated motions captured the recorded motions up to approximately 4 Hz. While careful iterations regarding source parameters and corner frequencies were required, and, for the case of the Highland Park event, some of the near-source stations had relatively low accuracy, the present study established a positive step toward the utilization of physics-based simulations in practical applications. The computational efficiencies exhibited by EQSIM, especially on GPU clusters, further supported this assertion, as wall-clock times of simulations involving more than 10 billion grid points were as low as mathematical equation minutes. This permits ensemble simulations for a considered scenario event so that modeling uncertainties (e.g., source and geology) can be bracketed.

EQSIM↗

The Pantheon+ analysis: Improving the redshifts and peculiar velocities of Type Ia supernovae used in cosmological analyses

We examine the redshifts of a comprehensive set of published Type Ia supernovae, and provide a combined, improved catalogue with updated redshifts. We improve on the original catalogues by using the most up-to-date heliocentric redshift data available; ensuring all red shifts have uncertainty estimates; using the exact formulae to convert heliocentric redshifts into the Cosmic Microwave Background (CMB) frame; and utilising an improved peculiar velocity model that calculates local motions in redshift-space and more realistically accounts for the external bulk flow at high-redshifts. We review 2607 supernova redshifts; 2285 are from unique supernovae and 322 are from repeat observations of the same supernova. In total, we updated 990 unique heliocentric redshifts, and found 5 cases of missing or incorrect heliocentric corrections, 44 incorrect or missing supernova coordinates, 230 missing heliocentric or CMB frame redshifts, and 1200 missing redshift uncertainties. The absolute corrections range between 10 -8 ≤ Δz ≤ 0.038, and RMS(Δz) ~ 3 × 10 -3 . The sign of the correction was essentially random, so the mean and median corrections are small: 4 × 10 -4 and 4 × 10 -6 respectively. We examine the impact of these improvements for H 0 and the dark energy equation of state w and find that the cosmological results change by ΔH 0 = -0.12 km s -1 Mpc -1 and Δw = 0.003, both significantly smaller than previously reported uncertainties for H 0 of 1.0 km s -1 Mpc -1 and w of 0.04 respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

Machine Learning‐Assisted Microearthquake Location Workflow for Monitoring the Newberry Enhanced Geothermal System

Abstract Enhanced geothermal systems (EGS) offer a sustainable energy source but face challenges in accurately locating microearthquakes induced during reservoir stimulation. Locating these microearthquakes provides reliable feedback on the stimulation progress. Current deep learning methods for locating earthquakes require extensive data sets for training, which is problematic as detected microearthquakes are often limited. To address the scarcity of training data, we propose a practical workflow using probabilistic multilayer perceptron (PMLP) which predicts microearthquake locations from cross‐correlation time lags in waveforms. Utilizing a 3D velocity model of Newberry site derived from ambient noise interferometry, we generate numerous synthetic microearthquakes and 3D acoustic waveforms for PMLP training. Accurate synthetic tests prompt us to apply the trained network to the 2012 and 2014 stimulation field waveforms. To enhance the accuracy of source localization, we carefully handpick the P‐arrival times. Predictions on the 2012 stimulation data set show major microseismic activity at depths of 0.5–1.2 km, correlating with a known casing leakage scenario. In the 2014 data set, the majority of predictions concentrate at 2.0–2.9 km depths, consistent with results obtained from conventional physics‐based inversion, and align with the presence of natural fractures from 2.0 to 2.7 km. We validate our findings by comparing the synthetic and field picks, demonstrating a satisfactory match for the first arrivals. By combining the benefits of quick inference speeds and accurate location predictions, we demonstrate the feasibility of using realistic synthetic data set to locate microseismicity for EGS monitoring.

15 GEOTHERMAL ENERGY↗

Relocating microseismicity from downhole monitoring of the Decatur CCS site using a modified double-difference algorithm

SUMMARY The injection of CO2 at the Decatur carbon capture and storage site has generated significant microseismic activity, which occurs in distinct spatial clusters up to approximately 2.2 km from the primary injection well. Accurate and precise event locations are vital for the characterization of the microseismicity to help understand the reservoir response to the CO2 injection, whilst enabling the identification of minor faults and fractures below the resolution of conventional active seismic imaging. However, microseismic monitoring of fluid injection sites, such as Decatur, is often performed using a network of borehole sensors often from a single well. While these downhole sensors have excellent detection capabilities, their poor azimuthal coverage limits the ability to precisely determine event locations. We have developed a modified double-difference relocation algorithm suitable for both 1-D and 3-D velocity models, and which incorporates differential back azimuth observations to allow the benefits of the original double-difference algorithm to be applicable to a downhole microseismic monitoring setting. Applying the modified double-difference algorithm to the microseismicity at Decatur, we have successfully relocated 4293 events. The relocation included over 59 million observations for 757 285 event pairs, split across seven geographic regions. Despite the majority of observations being recorded in only two boreholes, with an almost identical azimuthal coverage, the results have shown to be reliable with significantly reduced residuals and low uncertainties associated with the final locations. We have analysed the residuals in terms of their association with each geographic region, data type, station and individual events, to fully appreciate their influence in the inversion and the fit of the data to the final set of event locations. For each region, the relocated seismicity has become less diffuse with improved clustering, and with newly visible linear features often orientated in a NE–SW direction. These results show the potential improvements that can be made to microseismic event locations recorded by a borehole network with a limited and variable azimuthal distribution.

58 GEOSCIENCES↗

Time-lapse seismic data inversion for estimating reservoir parameters using deep learning

Geologic carbon sequestration involves the injection of captured carbon dioxide ([Formula: see text]) into subsurface formations for long-term storage. The movement and fate of the injected [Formula: see text] plume is of great concern to regulators because monitoring helps to identify potential leakage zones and determines the possibility of safe long-term storage. To address this concern, we design a deep-learning framework for [Formula: see text] saturation monitoring to determine the geologic controls on the storage of the injected [Formula: see text]. We use different combinations of porosities and permeabilities for a given reservoir to generate saturation and velocity models. We train the deep-learning model with a few time-lapse seismic images and their corresponding changes in saturation values for a particular [Formula: see text] injection site. The deep-learning model learns the mapping from the change in the time-lapse seismic response to the change in [Formula: see text] saturation during the training phase. We then apply the trained model to data sets comprising different time-lapse seismic image slices (corresponding to different time instances) generated using different porosity and permeability distributions that are not part of the training to estimate the [Formula: see text] saturation values along with the plume extent. Our algorithm provides a deep-learning assisted framework for the direct estimation of [Formula: see text] saturation values and plume migration in heterogeneous formations using the time-lapse seismic data. Our method improves the efficiency of time-lapse inversion by streamlining the large number of intermediate steps in the conventional time-lapse inversion workflow. This method also helps to incorporate the geologic uncertainty for a given reservoir by accounting for the statistical distribution of porosity and permeability during the training phase. Tests on different examples verify the effectiveness of our approach.

Geochemistry & Geophysics↗

Utah FORGE Project 3-2417: DAS Microseismic Event Catalog from the 16A/16B Circulation Test, 2023

This preliminary data archive includes the relocated microseismic event catalog, 1D velocity model, and methods report from DAS acquisition conducted during the Well 16A and 16B circulation test (July 19th and 20th, 2023) at Utah FORGE. The methods report describes all processing steps, including real-time event detection, hierarchical clustering, joint velocity/hypocenter inversion, and relocation. The resulting work is accepted and will be presented at IMAGE 2024. This dataset was acquired by the FOGMORE R&D project (Fiber Optic MOnitoring for Reservoir Evolution), Utah FORGE R&D Project 3-2417.

15 GEOTHERMAL ENERGY↗

Utah FORGE: GES Well 16A(78)-32 and Well 16B(78)-32 Stimulation Seismic Event Catalogs

This dataset contains seismic event catalogs from the hydraulic stimulation of wells 16A(78)-32 and 16B(78)-32 at the Utah FORGE site in April 2024. The data was collected by Geo Energy Suisse (GES) using a variety of seismic monitoring technologies, including 3-component (3C) geophones and distributed acoustic sensing (DAS) systems. These technologies were deployed across several locations, including wells 16A, 16B, and Delano-1, with sensor arrays at multiple depths to capture microseismic activity during the stimulations. The catalogs provide both real-time and manually checked seismic event locations, with detailed parameters such as trigger conditions, velocity models, and data acquisition settings. The dataset includes information on the stimulation stages, event rates, and hydraulic injection conditions for each well, with a report detailing the data acquisition configuration and seismic event location methodologies. Users will need to reference the included report for a complete understanding of the sensor network, data processing techniques, and accuracy considerations.

15 GEOTHERMAL ENERGY↗

Relocation of the Seismicity of the Caucasus Region

Our objective is to improve the view of the seismicity in the Caucasus region using instrumental data between 1951 and 2019. To create a comprehensive catalog, we combine the bulletins of local agencies and the International Seismological Centre, and use an advanced single-event location algorithm, iLoc, to obtain better locations. We show that relocations with iLoc, using travel-time predictions from the 3D upper mantle velocity model, Regional Seismic Travel Time, improve the locations. Then, using the iLoc results as initial locations and the ground-truth events identified in the iLoc results as fix points, we apply Bayesloc, a multiple-event location algorithm, to simultaneously relocate the entire seismicity of the Caucasus region. Here we demonstrate that the simultaneous relocation of the seismicity with Bayesloc clarifies the location and geometry of major active structures accommodating ongoing convergence between the Arabian and Eurasian continents between the Black and Caspian Seas. Among our major findings is the confirmation of widespread seismicity in the mantle beneath the northern flank of the Greater Caucasus and central Caspian, resulting from north-dipping subduction of the Kura and South Caspian basins and the identification of a discrete band of crustal seismicity beneath the southern flank of the Greater Caucasus.

58 GEOSCIENCES↗