Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic data”

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 811 records · Page 45

Diagnostic assessment of reservoir response to fracturing: a case study from Hydraulic Fracturing Test Site (HFTS) in Midland Basin

Abstract This paper outlines a data collection and diagnostics case study involving multiple horizontal shale wells. We look at well production profiles using rate transient analysis, differences in near wellbore complexity, geologic variations within the area of interest, as well as compositional differences in the rocks based on cores obtained from within the stimulated reservoir. The Hydraulic Fracturing Test Site is a multi-well experiment involving stimulation of unconventional shale wells in the southeastern Midland portion of the Permian Basin. The targeted formations include both the upper as well as the middle Wolfcamp formations, also referred alternatively as Wolfcamp A and Wolfcamp B. Data integration and analysis shared in this paper help us understand the various geologic controls impacting well productivity, particularly the wide variance observed between the Wolfcamp A and Wolfcamp B formations. Rate transient analysis indicates similar system permeabilities for stimulated wells. However, we observe higher effective fracture half-lengths for upper Wolfcamp wells. Using observations from 3D seismic interpretations (such as pad scale faults) as well as petrophysical and image log data, we highlight the substantial differences in stimulation as we move along the well laterals from the heel toward the toe sections. These differences are further reconciled with observations from zones with high data density at the core locations through stimulated rock, as well as independent data such as microseismic emissions. At the test site, Wolfcamp A was found to be relatively quartz rich with significant heterogeneity whereas Wolfcamp B is richer in clay and organic content. This impacts the geomechanical characteristics of the rock mass with much higher natural fracture density in the shallower interval. Thus, the fracture growth is more uniform in the deeper interval and more heterogeneous with branching likely in upper interval. Increased complexity also leads to consistently better productivity from the wells in the shallower interval as demonstrated from RTA results. This case study is unique because it provides valuable insights from actual sampling of the stimulated zones in hydraulically fractured wells and helps understand impact of various factors that contribute toward variability in well production. The findings from this study provides insights into need for optimization of completion designs in the various Wolfcamp landing zones, such as optimization of cluster or fracture spacing in various Wolfcamp intervals. In addition, it provides a useful template for data collection and research direction in future field test sites of similar nature in unconventional reservoirs.

Energy & Fuels↗

An estimation of continental crust magnetization and susceptibility from Magsat data for the conterminous United States

The Magsat-derived magnetic anomaly field over the United States was inverted to magnetization in a seismically determined variable-thickness lower crust. The results are quite similar to those obtained from a previous inversion of POGO data. The average magnetization is 3.45 + or - 1.0 A/m, which gives a bulk magnetic susceptibility of 8.7 + or - 2.4 x 10 to the -2nd (SI units) for the lower crustal layer consistent with a mafic composition. These values are also consistent with previous estimates of lower crustal magnetization from aeromagnetic data and with laboratory measurements of mafic rock susceptibilities.

Schnetzler, C. C.↗

Time-dependent gravity in Southern California, May 1974 to April 1979

The Southern California gravity survey, begun in May 1974 to obtain high spatial and temporal density gravity measurements to be coordinated with long-baseline three dimensional geodetic measurements of the Astronomical Radio Interferometric Earth Surveying project, is presented. Gravity data was obtained from 28 stations located in and near the seismically active San Gabriel section of the Southern California Transverse Ranges and adjoining San Andreas Fault at intervals of one to two months using gravity meters relative to a base station standard meter. A single-reading standard deviation of 11 microGal is obtained which leads to a relative deviation of 16 microGal between stations, with data averaging reducing the standard error to 2 to 3 microGal. The largest gravity variations observed are found to correlate with nearby well water variations and smoothed rainfall levels, indicating the importance of ground water variations to gravity measurements. The largest earthquake to occur during the survey, which extended to April, 1979, is found to be accompanied in the station closest to the earthquake by the largest measured gravity changes that cannot be related to factors other than tectonic distortion.

Whitcomb, J. H.↗

NETL Plastic Pipes Project (Final Report)

Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.

42 ENGINEERING↗

Uncertainty Quantification of Geophysical Inversion Using Stochastic Partial Differential Equations (LDRD #218329)

This report summarizes work completed under the Laboratory Directed Research and Development (LDRD) project "Uncertainty Quantification of Geophysical Inversion Using Stochastic Differential Equations." Geophysical inversions often require computationally expensive algorithms to find even one solution, let alone propagating uncertainties through to the solution domain. The primary purpose of this project was to find more computationally efficient means to approximate solution uncertainty in geophysical inversions. We found multiple computationally efficient methods of propagating Earth model uncertainty into uncertainties in solutions of full waveform seismic moment tensor inversions. However, the optimum method of approximating the uncertainty in these seismic source solutions was to use the Karhunen-Love theorem with data misfit residuals. This method was orders of magnitude more computationally efficient than traditional Monte Carlo methods and yielded estimates of uncertainty that closely approximated those of Monte Carlo. We will summarize the various methods we evaluated for estimating uncertainty in seismic source inversions as well as work toward this goal in the realm of 3-D seismic tomographic inversion uncertainty.

58 GEOSCIENCES↗

Dynamic Rupture Simulations of the Mw7.2 1992 Landers,California, Earthquake

An important part of the simulation technique implemented in the Recipe is the earthquake rupture characterization based on established empirical scaling relationships between source parameters and seismic moment. The scaling relationships used in the rupture characterization are derived from recorded data and are tested against dynamic and kinematic rupture models of crustal earthquakes. Currently, a revised three-stage source scaling relationship between the source rupture area and seismic moment of inland crustal earthquakes, initially proposed by Irikura and Miyake, is being used. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, including the PI, Arben Pitarka, indicate that in order to reproduce the near-fault long-period ground motions and the permanent displacements, observed during surface-rupturing earthquakes, the recipe needs modifications to the kinematic rupture parametrization in the shallow crustal weak zone. In particular, improvements are required in characterizing the slip rate function at shallow depths, above the sesimogenic zone. The main objective of this research project was the development of a dynamic rupture model for the Mw7.2 1992 Landers, California earthquake. The dynamic rupture model was used to derive kinematic rupture parameters, such as slip, slip-rate, rupture time distributions, and slip rate functions in the weak (0-4km) and seismogenic crustal zones. In addition to the knowledge gained from the dynamic rupture modeling of Chi-chi, Ridgecrest, Kumamoto and other crustal earthquakes, studied by the PI, the rupture characterization of the Landers earthquake can be helpful in updating the conventional recipe for broadband ground motion simulations of crustal earthquakes with shallow rupture. The project emphasizes the need for performing physics-based earthquake rupture modeling which ultimately will inform further improvements of rupture characterization. In particular, the modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, used in adjustments of empirical source scaling relationships for large crustal earthquakes.

58 GEOSCIENCES↗

Widespread deep seismicity in the Delaware Basin, Texas, is mainly driven by shallow wastewater injection

Significance The Delaware Basin, western Texas, has experienced a surge in the number of earthquakes. We investigate basin-wide earthquake, hydrogeologic, industrial, and geodetic data spanning 1993–2020 to identify the reasons for the increase. We use these data and models to quantify the stress changes from wastewater injection and hydrocarbon extraction. We demonstrate that the basin-wide seismicity is dominated by the stresses caused by shallow injection that are transmitted vertically by deforming the underlying crust. We analyze both injection and production data to understand how subsurface properties control earthquake locations and timings. We suggest that injecting wastewater into porous formations will minimize seismic hazard. We establish a framework explaining why earthquakes sometimes occur and sometimes do not occur within human activity-impacted areas.

54 ENVIRONMENTAL SCIENCES↗

Big Sky Regional Carbon Sequestration Partnership (Phase III Final Scientific/Technical Report)

The Big Sky Carbon Sequestration Partnership (BSCSP) pursued a Phase III demonstration project at Kevin Dome in north central Montana. Kevin Dome covers approximately 700 square miles and is a naturally occurring CO 2 reservoir that is flanked by oil and gas fields. The carbon dioxide (CO 2 ) is in the upper Devonian Duperow (carbonate) formation and does not reach the spill point of the dome; therefore, the dome has potential as a CO 2 sequestration reservoir, a CO 2 supply, or as both if anthropogenic sources and enhanced oil recovery (EOR) operations are associated with the dome. Kevin Dome could potentially act as a buffer to continue accepting anthropogenic CO 2 when EOR flooding operations are interrupted or completed. The project objective was to produce one million tonnes of CO 2 from the gas cap of Kevin Dome, pipe it laterally, inject, and re-store it in the brine leg of the same formation to test the hub / buffer storage concept. This was to be accomplished by drilling up to five production wells, building a short pipeline and compression facilities, and drilling an injection well and several monitoring wells. BSCSP commenced outreach and site characterization activities including acquiring baseline data for near-surface insurance monitoring, acquiring 3-dimensional, 9-component surface seismic over the project area, drilling two test wells (one in the production area and one in the injection area), coring key intervals, and performing comprehensive logging. Well tests of those wells revealed two barriers to the project. The production (Danielson 33-17) well was perforated in multiple zones but failed to produce any significant CO 2 . This was despite being drilled in the near vicinity of a historic well that had produced 3700 MCF per day in a drill stem test. Modeling indicated that this was likely due to a phase change during production causing a temperature drop resulting in hydrate and/or water ice formation that clogged the formation. Tests of the injection zone (Wallewein 22-1) well indicated total dissolved solids (TDS) slightly below the EPA required 10,000 parts per million (ppm) for a Class VI underground injection control permit. While the project was initiated before Class VI rules were promulgated, and this was an experimental project (seemingly qualified for a Class V permit), the Environmental Protection Agency (EPA) indicated that injection would require a Class VI permit. The low salinity result was unexpected as contours plotted based on regional formation water quality data indicated an expected TDS above 20,000 ppm, and wells between the recharge zone and the Wallewein well tested above 10,000 ppm. Faced with the inability to obtain an injection permit, the demonstration project could not proceed. However, the project had generated valuable samples and data on a large natural analog including 32 sq. mi. of 3-D, 9-C seismic, 430 ft. of carbonate core covering seven different depositional environments taken from areas with, and without the presence of CO 2 , 30 ft. of core of two caprocks, a tight carbonate and an anhydrite, a full set of modern logs on both wells, and well tests. DOE decided to re-scope the project around completing studies utilizing this data. This report covers both the initial scope and the re-scope (Task / Section numbers preceded with an R). While the report covers a wide range of project activities, highlights of this work include: Development of a geostatic model using neural nets to match well logs to facies and using multi-waveform seismic to inform reservoir heterogeneity; Unique mechanical testing of permeability – stress relationship in two caprock materials; Development of full waveform inversion to generate a high resolution velocity model; Model development for dual permeability (fracture and matrix) systems to better account for matrix-matrix interactions; Joint seismic wave inversion (including the first quadr-joint inversion) exhibiting better imaging of a challenging reservoir zone in stiff rock; Core flow and core flood results on a reactive carbonate; and Innovative laboratory measurements of seismic response of fractured core as a function of fluid fill.

54 ENVIRONMENTAL SCIENCES↗

Lunar seismicity, structure, and tectonics

Natural seismic events have been detected by the long-period seismometers at Apollo stations 16, 14, 15, and 12 at annual rates of 3300, 1700, 800, and 700, respectively, with peak activity at 13- to 14-day intervals. The data are used to describe magnitudes, source characteristics, and periodic features of lunar seismicity. In a present model, the rigid lithosphere overlies an asthenosphere of reduced rigidity in which present-day partial melting is probable. Tidal deformation presumably leads to critical stress concentrations at the base of the lithosphere, where moonquakes are found to occur. The striking tidal periodicities in the pattern of moonquake occurrence and energy release suggest that tidal energy is the dominant source of energy released as moonquakes. Thus, tidal energy is dissipated by moonquakes in the lithosphere and probably by inelastic processes in the asthenosphere.

Lammlein, D. R.↗

Utah FORGE: Fluid Injection Induced Shearing Experiments on Fractured Granitoid at Elevated Temperatures

This repository contains experimental data from a series of fluid injection-induced shearing tests conducted on Utah FORGE granitoid. The experiments were performed using an aluminum triaxial pressure vessel (TEMCO) apparatus at Pennsylvania State University. The primary aim was to investigate the impact of temperature on fault seismicity and to explore how different pre-stress ratios influence shearing behavior. The data includes results from experiments conducted under varying conditions of temperature, pre-stress ratios, and pore pressure increments. The rock samples used in these experiments were 60-grit granitoid with a single inclined fracture oriented at 60 degrees with respect to the horizontal cross-section. The dataset also includes measurements of normal stiffness at different temperatures and images of the experimental apparatus. The included README file details each experiential setup and outlines which data files represent which experimental conditions.

15 GEOTHERMAL ENERGY↗

Graph-learning approach to combine multiresolution seismic velocity models

SUMMARY The resolution of velocity models obtained by tomography varies due to multiple factors and variables, such as the inversion approach, ray coverage, data quality, etc. Combining velocity models with different resolutions can enable more accurate ground motion simulations. Toward this goal, we present a novel methodology to fuse multiresolution seismic velocity maps with probabilistic graphical models (PGMs). The PGMs provide segmentation results, corresponding to various velocity intervals, in seismic velocity models with different resolutions. Further, by considering physical information (such as ray path density), we introduce physics-informed probabilistic graphical models (PIPGMs). These models provide data-driven relations between subdomains with low (LR) and high (HR) resolutions. Transferring (segmented) distribution information from the HR regions enhances the details in the LR regions by solving a maximum likelihood problem with prior knowledge from HR models. When updating areas bordering HR and LR regions, a patch-scanning policy is adopted to consider local patterns and avoid sharp boundaries. To evaluate the efficacy of the proposed PGM fusion method, we tested the fusion approach on both a synthetic checkerboard model and a fault zone structure imaged from the 2019 Ridgecrest, CA, earthquake sequence. The Ridgecrest fault zone image consists of a shallow (top 1 km) high-resolution shear-wave velocity model obtained from ambient noise tomography, which is embedded into the coarser Statewide California Earthquake Center Community Velocity Model version S4.26-M01. The model efficacy is underscored by the deviation between observed and calculated traveltimes along the boundaries between HR and LR regions, 38 per cent less than obtained by conventional Gaussian interpolation. The proposed PGM fusion method can merge any gridded multiresolution velocity model, a valuable tool for computational seismology and ground motion estimation.

Geochemistry & Geophysics↗

Evaluating multistation phase picking algorithm phase neural operator (PhaseNO) on local seismic networks

Reliable automatic phase picking is important for many seismic applications. With the development of machine learning approaches, many algorithms are proposed, evaluated and applied to different areas. Many of these algorithms are single station based, while recent proposed methods start to combine surrounding stations into consideration in the problem of phase picking. Among these algorithms, the phase neural operator (PhaseNO) shows promising results on regional data sets comparing to existing algorithms. But there are many use cases for the local seismic networks in our community, therefore in this paper we evaluate the performance of PhaseNO on four different local data sets and compare the results to PhaseNet and EQTransformer. We used both individual phase picking metrics as well as association metrics to illustrate the performance of PhaseNO. By manually reviewing the newly detected events, we find that the PhaseNO model outperforms the single station-based approaches in the local-scale use cases due to its consideration of coherent signals from multiple stations. We also explored PhaseNO’s behaviours when only using one station, as well as gradually increasing the number of stations in the seismic network to better understand its behaviour. Overall, using the off-the-shelf machine learning based phase pickers, PhaseNO demonstrated its good performance on local-scale seismic networks.

58 GEOSCIENCES↗

Compilation of a Comprehensive Earthquake Catalog and Relocations in the Caucasus Region

Instrumental seismic monitoring has a long history in the Caucasus and started in 1899 when the first seismograph was installed in Tbilisi, Georgia. Much of the analog paper records from this time period are preserved in the Tbilisi archives because Georgia served as the regional data center. In the 1990s, due to the collapse of the Soviet Union and the political turmoil in the region, the analog networks and the communication between the newly formed national networks deteriorated. In Georgia, for the next 13 yr, the seismic network coverage was poor until the 2002 Tbilisi earthquake. Following this earthquake, the first permanent digital seismic station in Georgia was established in Tbilisi in 2003. The digital era progressively improved the ability to collect and archive data and today more than a hundred broadband seismic stations (including temporary arrays) are operating in the southern Caucasus. Until recently, the region lacked a coordinated effort to catalog all analog and digital era data collected by different countries into a single repository. As a result of collaboration between Lawrence Livermore National Laboratory, the Ilia State University, and the Republican Seismic Survey Center of Azerbaijan, a comprehensive earthquake catalog was compiled for the Caucasus and neighboring areas as part of a broader probabilistic seismic hazard assessment project. Here this project digitized Soviet-era paper bulletins, compiled a unified earthquake catalog from regional bulletins, developed 1D reference velocity model, and used it to relocate the events. The final catalog contains 16,963 events with magnitudes 3.7 and above, bringing together all the available data sets in the Caucasus region from 1900 to 2015, significantly improving locations, and generating the most complete earthquake catalog in the region, temporally and geographically.

58 GEOSCIENCES↗

Appendices for Geothermal Exploration Artificial Intelligence Report

The Geothermal Exploration Artificial Intelligence looks to use machine learning to spot geothermal identifiers from land maps. This is done to remotely detect geothermal sites for the purpose of energy uses. Such uses include enhanced geothermal system (EGS) applications, especially regarding finding locations for viable EGS sites. This submission includes the appendices and reports formerly attached to the Geothermal Exploration Artificial Intelligence Quarterly and Final Reports. The appendices below include methodologies, results, and some data regarding what was used to train the Geothermal Exploration AI. The methodology reports explain how specific anomaly detection modes were selected for use with the Geo Exploration AI. This also includes how the detection mode is useful for finding geothermal sites. Some methodology reports also include small amounts of code. Results from these reports explain the accuracy of methods used for the selected sites (Brady Desert Peak and Salton Sea). Data from these detection modes can be found in some of the reports, such as the Mineral Markers Maps, but most of the raw data is included the DOE Database which includes Brady, Desert Peak, and Salton Sea Geothermal Sites.

15 GEOTHERMAL ENERGY↗

Salton Sea Geodatabase for Geothermal Exploration Artificial Intelligence

These files contain the geodatabases related to Salton Sea Geothermal Field. It includes all input and output files used with the Geothermal Exploration Artificial Intelligence. Input and output files are sorted into three categories: raw data, pre-processed data, and analysis (post-processed data). In each of these categories there are six additional types of raster catalogs which are titled Radar, SWIR, Thermal, Geophysics, Geology, and Wells. The files are used with the Geothermal Exploration Artificial Intelligence for the Salton Sea Geothermal Site to identify indicators of blind geothermal systems. The included zip file is a geodatabase to be used with ArcGIS and the tar file is an inclusive database that encompasses the inputs and outputs for the Salton Sea Geothermal Site.

15 GEOTHERMAL ENERGY↗

Internal friction in rocks and its relationship to volatiles on the moon

The physical properties of lunar rocks were measured using the vibrating bar technique in order to provide data for interpretation of geophysical results such as those from seismic measurements. The effect of volatiles on the mechanical Q in lunar rocks was studied in addition to the effect of exposing a sample to controlled amounts of those gases most likely to be present in the lunar environment or likely to have been outgassed from the lunar interior. The moderate temperatures to which the sample was exposed during the thermal treatment and the small drop in resonant frequency during the course of the outgassing suggests that there was little change in microfracture density. The frequency, composition and texture dependence of the damping were investigated, to study the loss mechanism.

Tittmann, B. R.↗

Satellite Relaying of Geophysical Data

Data Collection Platforms (DCPs) for transmitting surface data to an orbiting satellite for relaying to a central data distribution center are being used in a number of geophysical applications. "Off-the-shelf" DCP's, transmitting through Landsat or GOES satellites, are fully capable of relaying data from low-data-rate instruments, such as tiltmeters or tide gauges. In cooperation with the Lamont-Doherty Geological Observatory, Goddard has successfully installed DCP systems on a tide gauge and tiltmeter array on Anegada, British Virgin Islands. Because of the high-data-rate requirements, a practical relay system capable of handling seismic information is not yet available. Such a system could become the basis of an operational hazard prediction system for reducing losses due to major natural catastrophies such as earthquakes, volcanic eruptions, landslides or tsunamis.

Allenby, R. J.↗

Surface deformation associated with the November 23, 1977, Caucete, Argentina, earthquake sequence

The 1977 Caucete (San Juan) earthquake considered in the present paper occurred near the Sierra Pie de Palo in the Sierras Pampeanas tectonic province of western Argentina. In the study reported, coseismic surface deformation is combined with seismic observations (main shock and aftershocks, both teleseismic and local data) to place constraints on the geometry and slip of the main fault responsible for the 1977 earthquake. The implications of the 1977 event for long-term crustal shortening and earthquake recurrence rates in this region are also discussed. It is concluded that the 1977 Caucete earthquake was accompanied by more than 1 m of vertical uplift.

Kadinsky-Cade, K.↗