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 55 records · Page 3

Utah FORGE: 2D and 3D Seismic Data

This set of data contains raw and Initially-processed 2D and 3D seismic data from the Utah FORGE study area near Roosevelt Hot Springs. Reprocessed versions of these data can be accessed at the linked submission in the Resources section titled "Reprocessed Seismic Reflection Data." The zipped archives numbered from 1-100 to 1001-1122 contain 3D seismic uncorrelated shot gatherers SEG-Y files. The zipped archives numbered from 1-100C to 1001-1122C contain 3D seismic correlated shot gatherers SEG-Y files. Other data have intuitive names.

15 GEOTHERMAL ENERGY↗

Automatic Calibration of a Geomechanical Model from Sparse Data for Estimating Stress in Deep Geological Formations

Summary In this study, we demonstrate geomechanical modeling with fully automatic parameter calibration to estimate the full geomechanical stress fields of a prospective US carbon dioxide (CO2) storage site, based on sparse measurement data. The goal is to compute full stress tensor field estimates (principal stresses and orientations) that are maximally compatible with observations within the constraints of the model assumptions, thereby extending pointwise, incomplete partial stress measurement to a simulated full formation stress field, as well as a rough assessment of the associated error. We use the Perch site, located in Otsego County, Michigan, USA, as our case study. The input data consist of partial stress tensor information inferred from in-situ borehole tests, geophysical well logs, and processing of seismic data. A static earth model (SEM) of the site was developed, and geomechanical simulation functionality of the open-source MATLAB Reservoir Simulation Toolbox (MRST) was used to model the stress field. Adjoint-based nonlinear optimization was used to adjust boundary conditions and material properties to calibrate simulated results of observations. Results were interpreted through a Bayesian framework. The focus of this paper is to demonstrate how the fully automatic calibration procedure works and discuss the results obtained; it does not attempt a detailed analysis of the stress field in the context of the proposed CO2 storage initiatives. Our work is part of a larger effort to noninvasively determine in-situ stresses in deep formations considered for CO2 storage. Guided by previously published research on geomechanical model calibration, our work presents a novel calibration approach supporting a potentially large number of linear or nonlinear calibration parameters to produce results optimally agreeing with available measurements and thus extend partial pointwise estimates to full tensor fields compatible with the physics of the site.

Engineering↗

Processed Lab Data for Neural Network-Based Shear Stress Level Prediction

Machine learning can be used to predict fault properties such as shear stress, friction, and time to failure using continuous records of fault zone acoustic emissions. The files are extracted features and labels from lab data (experiment p4679). The features are extracted with a non-overlapping window from the original acoustic data. The first column is the time of the window. The second and third columns are the mean and the variance of the acoustic data in this window, respectively. The 4th-11th column is the the power spectrum density ranging from low to high frequency. And the last column is the corresponding label (shear stress level). The name of the file means which driving velocity the sequence is generated from. Data were generated from laboratory friction experiments conducted with a biaxial shear apparatus. Experiments were conducted in the double direct shear configuration in which two fault zones are sheared between three rigid forcing blocks. Our samples consisted of two 5-mm-thick layers of simulated fault gouge with a nominal contact area of 10 by 10 cm^2. Gouge material consisted of soda-lime glass beads with initial particle size between 105 and 149 micrometers. Prior to shearing, we impose a constant fault normal stress of 2 MPa using a servo-controlled load-feedback mechanism and allow the sample to compact. Once the sample has reached a constant layer thickness, the central block is driven down at constant rate of 10 micrometers per second. In tandem, we collect an AE signal continuously at 4 MHz from a piezoceramic sensor embedded in a steel forcing block about 22 mm from the gouge layer The data from this experiment can be used with the deep learning algorithm to train it for future fault property prediction.

15 GEOTHERMAL ENERGY↗

Elastic Stochastic Full Waveform Inversion (eSFWI)

This collaboration between Lawrence Livermore National Security, LLC (LLNS) as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Chevron USA Inc., acting through its Chevron Technical Center division, aimed at developing next-generation computational methods for the Elastic Stochastic Full Waveform Inversion (eSFWI). Seismic imaging is heavily used in the oil and gas industry for identifying and operating subsurface reservoirs. Improved seismic imaging methods can improve productivity, lower costs, and improve operational and environmental safety. This CRADA demonstrated that new high-performance computing (HPC) architectures being rolled out over the next five years can enable unprecedented seismic imaging resolution when using eSFWI techniques to process active seismic data. An open-source computational mini-application was developed, capable of demonstrating near-peak performance for eSFWI algorithms on CPU and GPU enabled HPC platforms. Performance was demonstrated on LLNL HPC systems such as Lassen, as well as on Chevron systems. This project benefited Chevron USA Inc. by demonstrating the potential computational efficiency of their full waveform inversion capabilities used to characterize oil/gas reservoirs, which in turn benefits the public through potential increases in capabilities to perform analysis of leasing sites.

04 OIL SHALES AND TAR SANDS↗

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

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

Jaffe, L. D.↗

Feasibility of Geophysical Monitoring of Carbon-Sequestrated Deep Saline Aquifers

As carbon dioxide (CO{sub 2}) is sequestered from the bottom of a brine reservoir and allowed to migrate upward, the effects of the relative permeability hysteresis due to capillary trapping and buoyancy driven migration tend to make the reservoir patchy saturated with different fluid phases over time. Seismically, such a patchy saturated reservoir induces an effective anisotropic behavior whose properties are primarily dictated by the nature of the saturation of different fluid phases in the pores and the elastic properties of the rock matrix. By combining reservoir flow simulation and modeling with seismic modeling, it is possible to derive these effective anisotropic properties, which, in turn, could be related to the saturation of CO{sub 2} within the reservoir volume any time during the post-injection scenario. Therefore, if time-lapse seismic data are available and could be inverted for the effective anisotropic properties of the reservoir, they, in combination with reservoir simulation could potentially predict the CO{sub 2} saturation directly from the time-lapse seismic data. It is therefore concluded that the time-lapse seismic data could be used to monitor the carbon sequestrated saline reservoirs. But for its successful implementation, seismic modeling and inversion methods must be integrated with the reservoir simulations. In addition, because CO{sub 2} sequestration induces an effective anisotropy in the sequestered reservoir and anisotropy is best detected using multicomponent seismic data compared to single component (P-wave) data, acquisition, processing, and analysis is multicomponent seismic data is recommended for these time-lapse studies. Finally, a successful implementation of using time-lapse seismic data for monitoring the carbon sequestrated saline reservoirs will require development of a robust methodology for inverting multicomponent seismic data for subsurface anisotropic properties.

54 ENVIRONMENTAL SCIENCES↗

Identification of Natural Oscillation Modes for Purposes of Seismic Assessment and Monitoring of HPP Dams

The paper puts forward a method for processing data from detailed seismic assessments of HPP dams (dynamic tests). A detailed assessment (hundreds of observation points in dam galleries) is performed with consideration of operating dam equipment and the microseismic noise. It is shown that dynamic oscillation characteristics (natural oscillation frequencies and modes in the main dam axes, the velocities of propagation of elastic waves with given polarization, and so on.) can be determined with sufficient accuracy by using complex transfer functions and pulse characteristics. Monitoring data is processed using data from a detailed assessment, taking account of identified natural oscillation modes and determined ranges of natural frequencies. The spectra of characteristic frequencies thus obtained are used to choose substitution models and estimate the elastic characteristics of the “dam – rock bed” construction system, viz., the modulus of elasticity (the Young modulus), the Poisson ratio, the dam section stiffness with respect to shear, tension and compression and the elastic characteristics of the rock foundation.

42 ENGINEERING↗

sRecovery Act: Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO 2 Storage

This study focuses on evaluating the feasibility and suitability of using the Jurassic/Triassic (J/TR) sediments of the South Georgia Rift basin (SGR) for CO 2 storage in southern South Carolina and southern Georgia The SGR basin in South Carolina (SC), prior to this project, was one of the least understood rift basin along the east coast of the U.S. In the SC part of the basin there was only one well (Norris Lightsey #1) the penetrated into J/TR. Because of the scarcity of data, a scaled approach used to evaluate the feasibility of storing CO 2 in the SGR basin. In the SGR basin, 240 km (~149 mi) of 2-D seismic and 2.6 km2 3-D (1 mi2) seismic data was collected, process, and interpreted in SC. In southern Georgia 81.3 km (~50.5 mi) consisting of two 2-D seismic lines were acquired, process, and interpreted. Seismic analysis revealed that the SGR basin in SC has had a very complex structural history resulting the J/TR section being highly faulted. The seismic data is southern Georgia suggest SGR basin has not gone through a complex structural history as the study area in SC. The project drilled one characterization borehole (Rizer # 1) in SC. The Rizer #1 was drilled but due to geologic problems, the project team was only able to drill to 1890 meters (6200 feet) instead of the proposed final depth 2744 meters (9002 feet). The drilling goals outlined in the original scope of work were not met. The project was only able to obtain 18 meters (59 feet) of conventional core and 106 rotary sidewall cores. All the conventional core and sidewall cores were in sandstone. We were unable to core any potential igneous caprock. Petrographic analysis of the conventional core and sidewall cores determined that the average porosity of the sedimentary material was 3.4% and the average permeability was 0.065 millidarcy. Compaction and diagenetic studies of the samples determined there would not be any porosity or permeability at depth in SC. In Georgia there appears to be porosity in the J/TR section based on neutron log porosity values. The only zones in Rizer #1 that appear to be porous were fractured diabase units where saline formation water was flowing into the borehole. Two geocellular models were created for the SC and GA study area. Flow simulation modeling was performed on the SC data set. The injection simulation used the newly acquired basin data as well as the Petrel 3-D geologic model that included geologic structure. Due to the new basin findings as a result of the newly acquired data, during phase two of the modeling the diabase unit was used as reservoir and the sandstone units were used as caprock. Conclusion are: 1) the SGR basin is composed of numerous sub-basins, 2) this study only looked at portions of two sub-basins, 3) in SC, 30 million tonnes of CO 2 can be injected into the diabase units if the fracture network is continuous through the units, 4) due to the severity of the faulting there is no way of assuring the injected CO 2 will not migrate upward into the overlying Coastal Plain aquifers, 5) in Georgia there appears to porous zones in the J/TR sandstones, 6) as in SC there is faulting in the sub-basin and the seismic suggest the faulting extends upward into the Coastal Plain making that area not suitable for CO 2 sequestration, 7) the complex faulting observed at both study areas appear to be associated with transfer fault zones (Heffner 2013), if sub-basins in the Georgia portion of the SGR basin can be located that are far away from the transfer fault zones there is a strong possibility of sequestering CO 2 in these areas, and 9) the SGR basin covers area in three states and this project only studied two small areas so there is enormous potential for CO 2 sequestration in other portions the basin and further research needs to be done to find these areas.

58 GEOSCIENCES↗

The initial appraisal of buried DAS system in CO2CRC Otway Project: the comparison of buried standard fibre-optic and helically wound cables using 2D imaging

This study aims to assess the ability of shallow distributed acoustic sensing (DAS) to serve as a cost-effective seismic sensor array for permanent monitoring applications. To this end, as part of the CO2CRC seismic monitoring program, a fibre-optic DAS array was deployed alongside a permanently buried geophone array at the Otway Project site (Victoria, Australia). The DAS array consisted of a standard commercially available tactical fibre-optic cable, which was deployed in 0.8m deep trenches. A custom-designed helically wound (HW) cable was also deployed in one of the DAS trenches for comparison of the cable designs. Simultaneous acquisition of the seismic data was carried out using 3000 vibroseis source points and geophones, DAS standard and HW cables. For initial assessment of the seismic images acquired with DAS and to compare different cable designs, preliminary 2D seismic reflection processing is conducted on both DAS cables and geophone data along a single 2D line. The geophone data processing guided processing of the DAS data. Several shallow structures (100-450ms) and some important reflectors at 450-600ms are observed on the final DAS images. Comparison of the two different DAS cable types demonstrated that seismic imaging would benefit DAS technology. However, the benefit of utilising HW cable is modest compared with the standard cable. The workflows and results of this study pave the way for processing of the 3D seismic data set acquired with the DAS array, as well as further detailed analysis of the DAS cables and the system itself.

47 OTHER INSTRUMENTATION↗

West Flank Coso FORGE: Seismic Reflection

PDFs of seismic reflection profiles 101,110, 111 local to the West Flank FORGE site. 45 line kilometers of seismic reflection data are processed data collected in 2001 through the use of vibroseis trucks. The initial analysis and interpretation of these data was performed by Unruh et al. (2001). Optim processed these data by inverting the P-wave first arrivals to create a 2-D velocity structure. Kirchhoff images were then created for each line using velocity tomograms (Unruh et al., 2001).

15 GEOTHERMAL ENERGY↗

Magnetotelluric Data Collected in 2016 over the San Emidio Geothermal Field in Nevada

This data set includes the magnetotelluric (MT) data collected from October 21 to November 9, 2016 over the San Emidio geothermal field in Nevada by Quantec Geoscience USA Inc. on behalf of US Geothermal Inc. as part of a project entitled "A Novel Approach to Map Permeability Using Passive Seismic Emission Tomography". This data set includes descriptions of the instrumentation, data acquisition and processing procedures, as well as the final processed data and digital archive formats. A total of 81 MT locations were surveyed (52 profile sites, and 29 MT sites). Data were processed and inspected for quality assurance on site, and reviewed daily by the geophysicist in charge of the project.

15 GEOTHERMAL ENERGY↗

Estimation of anisotropy parameters in organic-rich shale: Rock physics forward modeling approach

Anisotropy analysis becomes an important step in processing and interpretation of seismic data. One of the most important things in anisotropy analysis is anisotropy parameter estimation which can be estimated using well data, core data or seismic data. In seismic data, anisotropy parameter calculation is generally based on velocity moveout analysis. However, the accuracy depends on data quality, available offset, and velocity moveout picking. Anisotropy estimation using seismic data is needed to obtain wide coverage of particular layer anisotropy. In anisotropic reservoir, analysis of anisotropy parameters also helps us to better understand the reservoir characteristics. Anisotropy parameters, especially ε, are related to rock property and lithology determination. Current research aims to estimate anisotropy parameter from seismic data and integrate well data with case study in potential shale gas reservoir. Due to complexity in organic-rich shale reservoir, extensive study from different disciplines is needed to understand the reservoir. Shale itself has intrinsic anisotropy caused by lamination of their formed minerals. In order to link rock physic with seismic response, it is necessary to build forward modeling in organic-rich shale. This paper focuses on studying relationship between reservoir properties such as clay content, porosity and total organic content with anisotropy. Organic content which defines prospectivity of shale gas can be considered as solid background or solid inclusion or both. From the forward modeling result, it is shown that organic matter presence increases anisotropy in shale. The relationships between total organic content and other seismic properties such as acoustic impedance and Vp/Vs are also presented.

58 GEOSCIENCES↗

GRAIL Refinements to Lunar Seismic Structure

The present ]day internal structure of the Moon provides insight not only into its own formation and evolution, but also that of all rocky planetary bodies. The most direct way to probe a planet fs interior structure is through seismology. As part of the Apollo lunar missions, four seismometers were deployed on the nearside surface of the Moon between the years 1969 and 1972. These instruments operated continuously until 1977, forming the only substantial extraterrestrial seismic data set in existence. These data have been used to constrain various aspects of the seismic velocity and density structure of the Moon. Typical 1-D models recognize a 30-60 km thick crust overlying a nearly constant ]velocity mantle, and extend to a depth of approximately 1000 km, below which the lack of penetrating moonquake ray ]paths precludes the seismic determination of deeper structure. Previously, the lack of observed moonquakes from the far side of the Moon has been used to infer the presence of a highly attenuating (possibly molten) core. Indirect geophysical measurements such as moment of inertia, magnetic induction, lunar laser ranging, and elemental abundances of mare basalts also place varying constraints on core size and state. In combination with seismic studies, these indirect measurements have been used to arrive at a commonly accepted model of the Moon's deepest interior that includes a solid inner and fluid outer core, overlain by a partial melt boundary layer. We recently applied modern array seismology techniques to the Apollo data and revealed detailed core structure, including the first direct confirmation of the presence of a solid inner core. Our study focused on the identification of core ]reflected phases in deep moonquake seismograms. The resulting model of the Moon fs innermost structure was found to be consistent with the commonly accepted model. However, the modeled layer radii may vary by tens of kilometers, as is expected when accounting for uncertainties such as moonquake location, timing errors, and potential seismic heterogeneities. In addition, the modeled velocities may vary with a 1-to-1 trade ]off with the modeled reflector depth. The GRAIL (Gravity Recovery and Interior Laboratory) mission, launched in Sept. 2011, placed two nearly identical spacecraft in lunar orbit. The two satellites make extremely high-resolution measurements of the lunar gravity field, which can be used to constrain the interior structure of the Moon using a "crust to core" approach. GRAIL fs constraints on crustal thickness, mantle structure, core radius and stratification, and core state (solid vs. molten) will complement seismic investigations in several ways. Here we present a progress report on our efforts to advance our knowledge of the Moon fs internal structure using joint gravity and seismic analyses. We will focus on methodology, including 1) refinements to the seismic core constraint accomplished through array processing of Apollo seismic data, made by applying a set of travel time corrections based on GRAIL structure estimates local to each Apollo seismic station; 2) modeling deep lunar structure through synthetic seismograms, to test whether the seismic core model can reproduce the core reflections observed in the Apollo seismograms; and 3) a joint seismic and gravity inversion in which we attempt to fit a family of seismic structure models with the gravity constraints from GRAIL, resulting in maps of seismic velocities and densities that vary from a nominal model both laterally and with depth.

Weber, Renee↗

2D Seismic Reflection Data across Central Illinois

In a continuing collaboration with the Midwest Geologic Sequestration Consortium (MGSC) on the Evaluation of the Carbon Sequestration Potential of the Cambro-Ordovician Strata of the Illinois and Michigan Basins project, Schlumberger Carbon Services and WesternGeco acquired two-dimensional (2D) seismic data in the Illinois Basin. This work included the design, acquisition and processing of approximately 125 miles of (2D) seismic reflection surveys running west to east in the central Illinois Basin. Schlumberger Carbon Services and WesternGeco oversaw the management of the field operations (including a pre-shoot planning, mobilization, acquisition and de-mobilization of the field personnel and equipment), procurement of the necessary permits to conduct the survey, post-shoot closure, processing of the raw data, and provided expert consultation as needed in the interpretation of the delivered product. Three 2D seismic lines were acquired across central Illinois during November and December 2010 and January 2011. Traversing the Illinois Basin, this 2D seismic survey was designed to image the stratigraphy of the Cambro-Ordovician sections and also to discern the basement topography. Prior to this survey, there were no regionally extensive 2D seismic data spanning this section of the Illinois Basin. Between the NW side of Morgan County and northwestern border of Douglas County, these seismic lines ran through very rural portions of the state. Starting in Morgan County, Line 101 was the longest at 93 miles in length and ended NE of Decatur, Illinois. Line 501 ran W-E from the Illinois Basin – Decatur Project (IBDP) site to northwestern Douglas County and was 25 miles in length. Line 601 was the shortest and ran N-S past the IBDP site and connected lines 101 and 501. All three lines are correlated to well logs at the IBDP site. Originally processed in 2011, the 2D seismic profiles exhibited a degradation of signal quality below ~400 millisecond (ms) which made interpretation of the Mt. Simon and Knox sections difficult. The data quality also gradually decreased moving westward across the state. To meet evolving project objectives, in 2012 the seismic data was re-processed using different techniques to enhance the signal quality thereby rendering a more coherent seismic profile for interpreters. It is believed that the seismic degradation could be caused by shallow natural gas deposits and Quaternary sediments (which include abandoned river and stream channels, former ponds, and swamps with peat deposits) that may have complicated or changed the seismic wavelet. Where previously limited by seismic coverage, the seismic profiles have provided valuable subsurface information across central Illinois. Some of the interpretations based on this survey included, but are not limited to: - Stratigraphy generally gently dips to the east from Morgan to Douglas County. - The Knox Supergroup roughly maintains its thickness. There is little evidence for faulting in the Knox. However, at least one resolvable fault penetrates the entire Knox section. - The Eau Claire Formation, the primary seal for the Mt. Simon Sandstone, appears to be continuous across the entire seismic profile. - The Mt. Simon Sandstone thins towards the western edge of the basin. As a result, the highly porous lowermost Mt. Simon section is absent in the western part of the state. - Overall basement dip is from west to east. - Basement topography shows evidence of basement highs with on-lapping patterns by Mt. Simon sediments. - There is evidence of faults within the lower Mt. Simon Sandstone and basement rock that are contemporaneous with Mt. Simon Sandstone deposition. These faults are not active and do not penetrate the Eau Claire Shale. It is believed that these faults are associated with a possible failed rifting event 750 to 560 million years ago during the breakup of the supercontinent Rodinia.

Smith, Valerie↗