Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic data processing”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Produced Water Treatment Using the Switchable Polarity Solvent Forward Osmosis (SPS FO) Desalination Process: Preliminary Engineering Design Basis

Switchable Polarity Solvent Forward Osmosis (SPS FO) is a semi-permeable membrane-based water treatment technology. INL is currently advancing SPS FO technology such that a prototype unit can be designed and demonstrated for the purification of produced water from oil and gas production operations. The SPS FO prototype unit will used the thermal energy in the produced water as a source of process heat, thereby reducing the external process energy demands. Treatment of the produced water stream will reduce the volume of saline wastewater requiring disposal via injection, an activity that is correlated with undesirable seismic events, as well as generate a purified product water stream with potential beneficial uses. This paper summarizes experimental data that has been collected in support of the SPS FO scale-up effort, and describes how this data will be used in the sizing of SPS FO process equipment. An estimate of produced water treatment costs using the SPS FO process is also provided.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Causal mechanism of injection-induced earthquakes through the M w 5.5 Pohang earthquake case study

Causal mechanisms for fluid injection-induced earthquakes remain a challenge to identify. Past studies largely established spatiotemporal correlations. Here, we propose a multi-process causal mechanism for injection-induced earthquakes through a case study of the 2017 M w 5.5 induced earthquake near Pohang Enhanced Geothermal System, Korea, where detailed hydraulic stimulation and on-site seismicity monitoring data provide an unprecedented opportunity. Pore pressure modeling reveals that pore pressure changes initiate seismicity on critically stressed faults and Coulomb static stress transfer modeling reveals that earthquake interactions promote continued seismicity, leading to larger events. On the basis of these results, we propose the following causal mechanism for induced seismicity: pore pressure increase and earthquake interactions lead to fault weakening and ultimately triggering larger earthquakes later in the process. We suggest that it is prudent that pore pressure change, initial seismicity locations, and Coulomb static stress transfer from seismicity earlier in the sequence are assessed in real-time.

58 GEOSCIENCES↗
DOE OSTI2016

Seismic Fracture Characterization Methodologies for Enhanced Geothermal Systems

Executive Summary The overall objective of this work was the development of surface and borehole seismic methodologies using both compressional and shear waves for characterizing faults and fractures in Enhanced Geothermal Systems. We used both surface seismic and vertical seismic profile (VSP) methods. We adapted these methods to the unique conditions encountered in Enhanced Geothermal Systems (EGS) creation. These conditions include geological environments with volcanic cover, highly altered rocks, severe structure, extreme near surface velocity contrasts and lack of distinct velocity contrasts at depth. One of the objectives was the development of methods for identifying more appropriate seismic acquisition parameters for overcoming problems associated with these geological factors. Because temperatures up to 300º C are often encountered in these systems, another objective was the testing of VSP borehole tools capable of operating at depths in excess of 1,000 m and at temperatures in excess of 200º C. A final objective was the development of new processing and interpretation techniques based on scattering and time-frequency analysis, as well as the application of modern seismic migration imaging algorithms to seismic data acquired over geothermal areas. The use of surface seismic reflection data at Brady's Hot Springs was found useful in building a geological model, but only when combined with other extensive geological and geophysical data. The use of fine source and geophone spacing was critical in producing useful images. The surface seismic reflection data gave no information about the internal structure (extent, thickness and filling) of faults and fractures, and modeling suggests that they are unlikely to do so. Time-frequency analysis was applied to these data, but was not found to be significantly useful in their interpretation. Modeling does indicate that VSP and other seismic methods with sensors located at depth in wells will be the most effective seismic tools for getting information on the internal structure of faults and fractures in support of fluid flow pathway management and EGS treatment. Scattered events similar to those expected from faults and fractures are seen in the VSP reported here. Unfortunately, the source offset and well depth coverage do not allow for detailed analysis of these events. This limited coverage also precluded the use of advanced migration and imaging algorithms. More extensive acquisition is needed to support fault and fracture characterization in the geothermal reservoir at Brady's Hot Springs. The VSP was effective in generating interval velocity estimates over the depths covered by the array. Upgoing reflection events are also visible in the VSP results at locations corresponding to reflection events in the surface seismic. Overall, the high temperature rated fiber optic sensors used in the VSP produced useful results. Modeling has been found useful in the interpretation of both surface reflection seismic and VSP data. It has helped identify possible near surface scattering in the surface seismic data. It has highlighted potential scattering events from deeper faults in the VSP data. Inclusion of more detailed fault and fracture specific stiffness parameters are needed to fully interpret fault and fracture scattered events for flow properties (Pyrak-Nolte and Morris, 2000, Zhu and Snieder, 2002). Shear wave methods were applied in both the surface seismic reflection and VSP work. They were not found to be effective in the Brady's Hot Springs area. This was due to the extreme attenuation of shear waves in the near surface at Brady's. This does not imply that they will be ineffective in general. In geothermal areas where good shear waves can be recorded, modeling suggests they should be very useful for characterizing faults and fractures.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Probabilistic Assessment and Uncertainty Analysis of CO2 Storage Capacity of the Morrow B Sandstone—Farnsworth Field Unit

This paper presents probabilistic methods to estimate the quantity of carbon dioxide (CO2) that can be stored in a mature oil reservoir and analyzes the uncertainties associated with the estimation. This work uses data from the Farnsworth Field Unit (FWU), Ochiltree County, Texas, which is currently undergoing a tertiary recovery process. The input parameters are determined from seismic, core, and fluid analyses. The results of the estimation of the CO2 storage capacity of the reservoir are presented with both expectation curve and log probability plot. The expectation curve provides a range of possible outcomes such as the P90, P50, and P10. The deterministic value is calculated as the statistical mean of the storage capacity. The coefficient of variation and the uncertainty index, P10/P90, is used to analyze the overall uncertainty of the estimations. A relative impact plot is developed to analyze the sensitivity of the input parameters towards the total uncertainty and compared with Monte Carlo. In comparison to the Monte Carlo method, the results are practically the same. The probabilistic technique presented in this paper can be applied in different geological settings as well as other engineering applications.

03 NATURAL GAS↗
DOE OSTI2020

Integrated Modeling Report for CO 2 Storage with Enhanced Oil Recovery in Northern Michigan

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. MRCSP has been investigating various reservoir characterization and modeling technologies related to Carbon Capture, Utilization, and Storage (CCUS) in conjunction with carbon dioxide enhanced oil recovery (CO 2 -EOR) operations in multiple Silurian-age (Niagaran), oil-bearing, carbonate pinnacle reefs in northern Michigan, USA. This report provides a comprehensive discussion of reservoir modeling studies that were conducted for tracking oil production, forecasting CO 2 plume migration, and estimating associated storage in a number of reefs that were at different stages of their CO 2 -EOR life cycle. The modeling process for simulating oil production, CO 2 injection, and associated storage in these reefs entails two phases. The first phase, geologic framework modeling, integrates all pertinent geological and geophysical data (from logs, cores and seismic surveys) about reservoir structure, geometry, rock types, and property distributions (porosity, permeability, water saturation) into a 3-D distributed grid-based static earth model (SEM). The second phase, dynamic reservoir modeling, uses the SEM as a platform to simulate the movement of oil, gas, water, and CO 2 within the reservoir during primary hydrocarbon production, as well as during subsequent phases such as CO 2 -injection assisted EOR, plume migration, and associated storage. In addition, an assessment of coupled process effects is also carried out, where the impacts of geochemical and geo-mechanical processed induced by CO 2 injection are studied.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Vector elastic deconvolution migration with dual wavefield decomposition

Here, we report elastic-wave imaging using multicomponent data can provide more useful subsurface information than acoustic-wave imaging, but it is usually algorithmically challenging. We have developed a vector elastic deconvolution migration method for high-resolution imaging of subsurface structures in isotropic and anisotropic elastic media. Our new method uses a vector deconvolution imaging condition based on dual wavefield decomposition, including an explicit directional wavefield separation using the Hilbert transform and a P/S vector wavefield decomposition using the low-rank decomposition method. Using three elastic models, we numerically determine that our new method produces notably higher resolution and more amplitude-balanced elastic images compared with a crosscorrelation-based vector elastic reverse time migration method.

58 GEOSCIENCES↗
DOE OSTI2013

Natural time analysis of critical phenomena: The case of pre-fracture electromagnetic emissions

Criticality of complex systems reveals itself in various ways. One way to monitor a system at critical state is to analyze its observable manifestations using the recently introduced method of natural time. Pre-fracture electromagnetic (EM) emissions, in agreement to laboratory experiments, have been consistently detected in the MHz band prior to significant earthquakes. It has been proposed that these emissions stem from the fracture of the heterogeneous materials surrounding the strong entities (asperities) distributed along the fault, preventing the relative slipping. It has also been proposed that the fracture of heterogeneous material could be described in analogy to the critical phase transitions in statistical physics. In this work, the natural time analysis is for the first time applied to the pre-fracture MHz EM signals revealing their critical nature. Seismicity and pre-fracture EM emissions should be two sides of the same coin concerning the earthquake generation process. Therefore, we also examine the corresponding foreshock seismic activity, as another manifestation of the same complex system at critical state. We conclude that the foreshock seismicity data present criticality features as well.

58 GEOSCIENCES↗
DOE OSTI2013

3d Wave Propagation and Site Effects in the Humboldt Bay Area Using Strong Ground Motion Records from the M6.5 2010 Ferndale Earthquake

In this study we simulate and analyze strong ground motion waveform data recorded in the Humboldt Bay and Eureka areas during the M6.5 Ferndale earthquake of January 2010. The scope of the presented work was two-fold. First, we investigated the main aspects of seismic wave generation and propagation, including kinematic rupture process and 3D wave propagation scattering. Our goal is to analyze their potential effects on seismic motion recorded at free field stations across Humboldt Bay and Eureka, California, and test the performance of a standard broadband strong ground motion simulation technique. Second, we evaluated the implication of the non-linear response of shallow sediments on ground motion recorded at the Humboldt Bay geotechnical array. Using synthetic and recorded data we compared the performance of numerical and empirical techniques in simulating local non-linear site response. Our study provides insight into the composition of the wave field during the earthquake and an improved understanding of how the strong ground motion is affected by the 3D structure and non-linear response of shallow sedimentary layers in the Humboldt Bay.

58 GEOSCIENCES↗
DOE OSTI2024

Evaluating Physics-Informed Neural Network Performance for Seismic Discrimination between Earthquakes and Explosions

In this article, we evaluate adding a weak physics constraint, that is, a physics‐based empirical relationship, to the loss function with a physics‐informed manner in local distance explosion discrimination in the hope of improving the generalization capability of the machine learning (ML) model. We compare the proposed model with the two‐branch model we previously developed, as well as with a pure data‐driven model. Unexpectedly, the proposed model did not consistently outperform the pure data‐driven model. By varying the level of inconsistency in the training data, we find this approach is modulated by the strength of the physics relationship. In conclusion, this result has important implications for how to best incorporate physical constraints in ML models.

58 GEOSCIENCES↗
DOE OSTI2019

Distributed Acoustic Sensing Using Dark Fiber for Near-Surface Characterization and Broadband Seismic Event Detection

We present one of the first case studies demonstrating the use of distributed acoustic sensing deployed on regional unlit fiber-optic telecommunication infrastructure (dark fiber) for broadband seismic monitoring of both near-surface soil properties and earthquake seismology. We recorded 7 months of passive seismic data on a 27 km section of dark fiber stretching from West Sacramento, CA to Woodland, CA, densely sampled at 2 m spacing. This dataset was processed to extract surface wave velocity information using ambient noise interferometry techniques; the resulting VS profiles were used to map both shallow structural profiles and groundwater depth, thus demonstrating that basin-scale variations in hydrological state could be resolved using this technique. The same array was utilized for detection of regional and teleseismic earthquakes and evaluated for long period response using records from the M8.1 Chiapas, Mexico 2017, Sep 8th event. The combination of these two sets of observations conclusively demonstrates that regionally extensive fiber-optic networks can effectively be utilized for a host of geoscience observation tasks at a combination of scale and resolution previously inaccessible.

58 GEOSCIENCES↗
DOE OSTI2015

Vertical deformation at western part of Sumatra

This research tries to make advancement in GPS signal processing to estimate the interseismic vertical deformation field at western part of Sumatra Island. The data derived by Continuous Global Positioning System (CGPS) from Badan Informasi Geospasial (BIG) between 2010 and 2012. GPS Analyze at Massachusetts Institute of Technology (GAMIT) software and Global Kalman Filter (GLOBK) software are used to process the GPS signal to estimate the vertical velocities of the CGPS station. In order to minimize noise due to atmospheric delay, Vienna Mapping Function 1 (VMF1) is used as atmospheric parameter model and include daily IONEX file provided by the Center for Orbit Determination in Europe (CODE) as well. It improves GAMIT daily position accuracy up to 0.8 mm. In a second step of processing, the GLOBK is used in order to estimate site positions and velocities in the ITRF08 reference frame. The result shows that the uncertainties of estimated displacement velocity at all CGPS stations are smaller than 1.5 mm/yr. The subsided deformation patterns are seen at the northern and southern part of west Sumatra. The vertical deformation at northern part of west Sumatra indicates postseismic phase associated with the 2010 and 2012 Northern Sumatra earthquakes and also the long-term postseismic associated with the 2004 and 2005 Northern Sumatra earthquakes. The uplifted deformation patterns are seen from Bukit Tinggi to Seblat which indicate a long-term interseismic phase after the 2007 Bengkulu earthquake and 2010 Mentawai earthquake. GANO station shows a subsidence at rate 12.25 mm/yr, indicating the overriding Indo-Australia Plate which is dragged down by the subducting Southeast Asian Plate.

58 GEOSCIENCES↗
DOE OSTI2015

Bradys Hot Springs Ambient Noise Correlation Functions (Initial Waveforms)

These files are ambient noise correlation (ANC) functions calculated for 11 days of continuous seismic data recorded by the Lawrence Berkeley network in the Brady geothermal field. These are SAC formatted seismic waveforms. The stations included are BPB04, BPB05, BPB07, BPB08, BPRT1, BPRT2, BPRT3, BPRT5, BRB10, BRP01, BRP02, BRP03, BRP04, BPR06, and BRP09 The original data were cut into hour long traces and processed by differentiating, removing the mean, removing the trend, applying a 1% taper, whitening, removing the mean and trend again, and converting to single bit traces. The data were then correlated with those from other stations and stacked. The resulting files were then named according to the convention: STA1.STA2.CHAN1_CHAN2.NHOURS.stacked.sac The days included are from records during 2013 (julian days, 200,220-229). The ANC correlations were calculated on the raw data traces (without instrument corrections applied) to assess the quality of the signal as a function of frequency throughout the network. The data were recorded at 500 Hz. We observe high quality signals 30 Hz on all traces, and measurable signal up to 80 Hz on a subset of the traces.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

The Role of Dislocations in the Anelasticity of the Upper Mantle

Dislocation‐based dissipation mechanisms potentially control the viscoelastic response of Earth's upper mantle across a variety of geodynamic contexts, including glacial isostatic adjustment, postseismic creep, and seismic‐wave attenuation. However, there is no consensus on which dislocation‐based, microphysical process controls the viscoelastic behavior of the upper mantle. Although both intergranular (plastic anisotropy) and intragranular (backstress) mechanisms have been proposed, there is currently insufficient laboratory data to discriminate between those mechanisms. Here, we present the results of forced‐oscillation experiments in a deformation‐DIA apparatus at confining pressures of 3–7 GPa and temperatures of 298–1370 K. Our experiments tested the viscoelastic response of polycrystalline olivine—the main constituent of the upper mantle—at stress amplitudes from 70 to 2,800 MPa. Mechanical data are complemented by microstructural analyses of grain size, crystallographic preferred orientation, and dislocation density. We observe amplitude‐ and frequency‐dependent attenuation and modulus relaxation and find that numerical solutions of the backstress model match our results well. Therefore, we argue that interactions among dislocations, rather than intergranular processes (e.g., plastic anisotropy or grain boundary sliding), control the viscoelastic behavior of polycrystalline olivine in our experiments. In addition, we present a linearized version of the constitutive equations of the backstress model and extrapolate it to conditions typical of seismic‐wave propagation in the upper mantle. Our extrapolation demonstrates that the backstress model can explain the magnitude of seismic‐wave attenuation in the upper mantle, although some modification is required to explain the weak frequency dependence of attenuation observed in nature and in previous experimental work.

Hein, Diede [Univ. of Minnesota, Minneapolis, MN (↗
DOE OSTI2015

Earthquake hypocenter relocation using double difference method in East Java and surrounding areas

Determination of precise hypocenter location is very important in order to provide information about subsurface fault plane and for seismic hazard analysis. In this study, we have relocated hypocenter earthquakes in Eastern part of Java and surrounding areas from local earthquake data catalog compiled by Meteorological, Climatological, and Geophysical Agency of Indonesia (MCGA) in time period 2009-2012 by using the double-difference method. The results show that after relocation processes, there are significantly changes in position and orientation of earthquake hypocenter which is correlated with the geological setting in this region. We observed indication of double seismic zone at depths of 70-120 km within the subducting slab in south of eastern part of Java region. Our results will provide useful information for advance seismological studies and seismic hazard analysis in this study.

58 GEOSCIENCES↗
DOE OSTI2025

Reducing the Parameter Dependency of Phase-Picking Neural Networks with Dice Loss

Training a neural network for picking seismic phase arrivals has been commonly posed as a segmentation problem. It is a highly imbalanced segmentation problem in the sense that the background vastly dominates the foreground because we are trying to pick the optimal single sample point that represents the arrival of a seismic phase in a many seconds long time window. Here, we test the Dice loss, which is a preferred loss function for highly imbalanced image segmentation problems. We show that phase-picking neural networks trained on the Dice loss behave in a binary fashion for which the prediction output is almost always either nearly 1 or nearly 0. This feature removes the strong dependence of data processing workflows on the prediction score threshold, which is an otherwise critical parameter to determine when using neural networks trained on the cross-entropy loss. When strategically used, models trained on the Dice loss can reduce the parameter dependency of machine learning-based seismic monitoring.

58 GEOSCIENCES↗
DOE OSTI2019

Induced seismicity during the 2012 Newberry EGS stimulation: Assessment of two advanced earthquake detection techniques at an EGS site

To improve the microearthquake detection capability during the 2012 Newberry enhanced geothermal system (EGS) demonstration, we applied two advanced earthquake detection algorithms to the continuous seismic data. These techniques have the power to detect microearthquakes within the seismic noise that would otherwise be missed using traditional techniques. We compare the efficacy of a matched field processing algorithm with a multi-channel waveform correlation algorithm. Both methods identified between 60%–80% more events, with each slightly outperforming the other in different seismic noise contexts. Incorporating either template-matching technique into a seismic processing routine could therefore effectively increase the sensitivity of the seismic network without changes to the existing hardware. Additionally, lowering the magnitude of completeness would significantly add to the number of events that could subsequently inform future decision making at the reservoir level.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Central Asia Seismic Hazard Assessment (Casha) "Complete" Catalog of Earthquakes

Probabilistic Seismic Hazard Assessments (PSHA) greatly benefit from reliable earthquake catalogs. Central Asia, while rich in seismic monitoring history, reports a fraction of their earthquake data to international data collection centres. Furthermore, since the collapse of Soviet Union until now, monitoring and data collection efforts in Central Asia have not been centralized, making monitoring and reporting uneven and disconnected in this time period. Lawrence Livermore National Laboratory (LLNL) has undertaken a multi-year effort to systematically work with national seismic monitoring and data centres and seismic networks in Kazakhstan, Kyrgyzstan, and Tajikistan to compile a unified seismological bulletin for Central Asia. In tandem with this effort, an earthquake catalog that is processed to be PSHA ready was also generated. The processing for this effort included further duplicate identification and removal, magnitude harmonization and catalog completeness analyses. This report describes the “complete” PSHA-ready earthquake catalog, including the completeness intervals that were applied.

58 GEOSCIENCES↗