Deep Learning Denoising Applied to the University of Utah Seismic Station Network Data.
Abstract not provided.
SEARCH · Engineering Papers
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.
Abstract not provided.
Seismic refraction data, obtained at the Apollo 14 and 16 sites, when combined with other lunar seismic data, allow a compressional wave velocity profile of the lunar near-surface and crust to be derived. The regolith, although variable in thickness over the lunar surface, possesses surprisingly similar seismic properties. Underlying the regolith at both the Apollo 14 Fra Mauro site and the Apollo 16 Descartes site is low-velocity brecciated material or impact derived debris. Key features of the lunar seismic velocity profile are: (1) velocity increases from 100 to 300 m/sec in the upper 100 m to about 4 km/sec at 5 km depth, (2) a more gradual increase from about 4 km/sec to about 6 km/sec at 25 km depth,(3) a discontinuity at a depth of 25 km, and (4) a constant value of about 7 km/sec at depths from 25 km to about 60 km.
To probe a planet's interior, seismology provides the most direct constraints on the variables that govern the dynamic properties of the body. However, the GRAIL (Gravity Recovery and Interior Laboratory) mission's high-resolution measurements of the lunar gravity field provide constraints on crustal thickness, mantle structure, core radius and stratification, and core state (solid vs. molten). These data complement seismic investigations, and joint interpretation permits improved constraints on the Moon's internal structure. Joint interpretation of disparate geophysical datasets helps reduce drawbacks that can result from analyzing them individually. The Apollo seismic network was situated on the lunar nearside surface in a roughly equilateral triangle having sides approximately 1000 km long, with stations 12/14 nearly co-located at one corner. Due to this limited geographical extent, near-surface ray coverage from moonquakes is low, but increase with depth. In comparison, gravity surveys and their resulting gravity anomaly maps have traditionally offered optimal resolution at crustal depths. Gravimetric maps and seismic data sets are therefor well suited to joint inversion, since the complementary information reduces inherent model ambiguity. We will perform a joint inversion of Apollo seismic delay times and gravity data collected by GRAIL lunar gravity mission, in order to recover seismic velocity and density as a function of latitude, longitude and depth within the Moon. We will relate density (rho) to seismic velocity (v) using a linear relationship that is allowed to be depth-dependent. The corresponding coefficient (B) can reflect a variety of material properties that vary with depth, including temperature and composition. The inversion seeks to recover the set of rho, v, and B perturbations that minimize (in a least-squares sense) the difference between the observed and calculated data.
Abstract Faults play an important role in recharging many geothermal reservoirs, and seismic information can image the locations of these faults. The value of information (VOI) metric is used to objectively quantify and compare the value of two types of seismic receiver data via a machine learning approach. The demonstrated VOI methodology is novel by including spatial models from seismic data and obtaining the information statistics from machine learning. Our two-dimensional numerical experiments compare images created from sparsely spaced (80 m), two-component geophone sampling to high spatial resolution (1 m), single-component DAS. We used a three-fold cross validation of a U-Net convolutional neural networks to achieve average classification statistics. The results suggest that when horizontal sources are utilized, geophones and DAS identify reflectors and non-reflectors at roughly the same rate. The average F1 score for horizontal DAS is 0.939 and 0.931 for geophones. For images created from a vertical source, DAS performed marginally better (F1 = 0.919) than geophones (F1 = 0.877). Our transferrable methodology can provide guidance on which acquisition scenarios can improve images of important structures in the subsurface and present an efficient method for obtaining reliability statistics from high-dimensional, spatial data.
Data relevant to the shallow structure of the moon obtained at the Apollo seismic stations are compared with previously published results of the active seismic experiments. It is concluded that the lunar surface is covered by a layer of low seismic velocity which appears to be equivalent to the lunar regolith defined previously by geological observations. This layer is underlain by a zone of distinctly higher seismic velocity at all of the Apollo landing sites. The regolith thicknesses at the Apollo 11, 12, and 15 sites are estimated from the shear-wave resonance to be 4.4, 3.7, and 4.4 m, respectively. These thicknesses and those determined at the other Apollo sites by the active seismic experiments appear to be correlated with the age determinations and the abundances of extralunar components at the sites.
The Silverpit crater is a recently discovered, 60-65 Myr old complex crater, which lies buried beneath the North Sea, about 150 km east of Britain. High-resolution images of Silverpit's subsurface structure, provided by three-dimensional seismic reflection data, reveal an inner-crater morphology similar to that expected for a 5-8 km diameter terrestrial crater. The crater walls show evidence of terracestyle slumping and there is a distinct central uplift, which may have produced a central peak in the pristine crater morphology. However, Silverpit is not a typical 5-km diameter terrestrial crater, because it exhibits multiple, concentric rings outside the main cavity. External concentric rings are normally associated with much larger impact structures, for example Chicxulub on Earth, or Orientale on the Moon. Furthermore, external rings associated with large impacts on the terrestrial planets and moons are widely-spaced, predominantly inwardly-facing, asymmetric scarps. However, the seismic data show that the external rings at Silverpit represent closely-spaced, concentric fault-bound graben, with both inwardly and outwardly facing faults-carps. This type of multi-ring structure is directly analogous to the Valhalla-type multi-ring basins found on the icy satellites. Thus, the presence and style of the multiple rings at Silverpit is surprising given both the size of the crater and its planetary setting.
The Silverpit crater is a recently discovered, 60-65 Myr old complex crater, which lies buried beneath the North Sea, about 150 km east of Britain. High-resolution images of Silverpit's subsurface structure, provided by three-dimensional seismic reflection data, reveal an inner-crater morphology similar to that expected for a 5-8 km diameter terrestrial crater. The crater walls show evidence of terrace-style slumping and there is a distinct central uplift, which may have produced a central peak in the pristine crater morphology. However, Silverpit is not a typical 5-km diameter terrestrial crater, because it exhibits multiple, concentric rings outside the main cavity. External concentric rings are normally associated with much larger impact structures, for example Chicxulub on Earth, or Orientale on the Moon. Furthermore, external rings associated with large impacts on the terrestrial planets and moons are widely-spaced, predominantly inwardly-facing, asymmetric scarps. However, the seismic data show that the external rings at Silverpit represent closely-spaced, concentric faultbound graben, with both inwardly and outwardly facing fault-scarps. This type of multi-ring structure directly analogous to the Valhalla-type multi-ring basins found on the icy satellites. Thus, the presence and style of the multiple rings at Silverpit is surprising given both the size of the crater and its planetary setting. A further curiosity of the Silverpit structure is that the external concentric rings appear to be extensional features on the West side of the crater and compressional features on the East side. The crater also lies in a local depression, thought to be created by postimpact movement of a salt layer buried beneath the crater.
On 24 September 2023, the Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer sample return capsule (SRC) entered Earth’s atmosphere after successfully collecting samples from an asteroid. The known trajectory and timing of this return provided a rare opportunity to strategically instrument sites to record geophysical signals produced by the capsule because it travelled at hypersonic speeds through the atmosphere. For this work, we deployed two distributed acoustic sensing (DAS) interrogators to sample over 12 km of surface-draped fiber-optic cables, as well as six collocated seismometer-infrasound sensor pairs, spread across two sites near Eureka, Nevada. This campaign-style rapid deployment is the first reported recording of an SRC re-entry with any distributed fiber-optic sensing technology. The DAS interrogators recorded an impulsive arrival with an extended coda which had features that were similar to recordings from both the seismometers and infrasound sensors. While the signal-to-noise ratio of the DAS data was lower than the seismic-infrasound data, the extremely dense spacing of fiber-optic sensors allowed for more phases to be clearly distinguished and the visualization of the continuous transformation of the wavefront as it impacted the ground. Unexpectedly, the DAS recordings contain less low-frequency content than is present in both the seismic and infrasound data. The deployment conditions strongly affected the recorded DAS data; in particular, we observed that fiber selection and placement exert strong controls on data quality.
MTtime is a software tool developed for time domain inversion of complete seismic waveform data to obtain the seismic moment tensor. It supports both deviatoric and full moment tensor inversions as well as 1-D and 3-D basis Green's functions.
This document presents final technical findings for the project Using Natural Gas Liquids to Recover Unconventional Oil and Gas Resources (FE0031782). The project is part of the U.S. Department of Energy Oil and Gas Program to develop and advance technologies that can significantly improve the recovery efficiencies of unconventional oil and gas resources. The overall objective of this project is to improve the ultimate recovery from unconventional oil and gas (UOG) resources in the United States by developing a method for using unrefined natural gas liquids (NGLs) as treatment fluids to improve hydrocarbon production. Horizontal extended-lateral drilling coupled with high volume hydraulic fracturing has significantly increased production from UOG resources in the U.S. However, the recovery efficiency is low compared to the estimated oil and gas in place. Recent data indicate that less than 10% of the oil in the liquid-rich UOG reservoirs is produced. Alternative completion methods using NGLs could increase production (Battelle, 2016; Wan et al., 2013; Wan, 2013; Downey et al, 2021); however, field validation tests are needed to develop an approach that is economical, efficient, and compatible in the UOG setting to advance towards commercial deployment. This project aims to develop and field test a method to improve recovery of oil resources in UOG shale plays by using Y-Grade NGLs, or a similar combination of NGLs, as treatment fluids. Refined NGLs have been used as a hydraulic treatment fluid in UOG plays for decades and are shown to be particularly effective because their miscibility with oil allows oil to flow more freely; however, the use of Y-Grade (unrefined) NGLs has not been studied. The use of Y-Grade NGLs would be advantageous over refined NGLs because Y-Grade NGLs do not require infrastructure or investment in refining and are already being produced from many UOG reservoirs. The concept was tested and monitored in the field at a commercial well site owned by project partner Hopco, Ltd. The project team, which consists of multiple oil and gas operators, Linde Gas North America LLC (Linde) and the Ohio Division of Geological Survey (ODGS), has extensive experience with oil and gas production in the Appalachian basin and the ability to work together quickly to solve technical issues and research needs. A key part of the proposed work was the use of existing wells for field testing and monitoring. A total of four wells (three vertical and one horizontal) were available for this project. One of the vertical wells was utilized as the test well for the NGL treatment test. A nearby vertical well was used for microseismic monitoring. The remaining vertical well and the horizontal well provided a baseline for typical UOG production in the oil window of the Utica/Point Pleasant (UPP). Major technical tasks of the project include characterization of the geotechnical properties of the UPP with an emphasis on the field site; design and planning for the NGLs testing; field testing and monitoring; analysis and integration of field data; and economic and resource/reserve assessment. The Shoman monitoring well and Doughty NGL treatment well were successfully plugged back during September-October 2020 in preparation for treatment and monitoring. A nitrogen diagnostic fracture injection test (DFIT) was completed on the Doughty well on July 22, 2021, consisting of 133,000 scf (91 Bbl.) of nitrogen. A nitrogen foam frac was completed in the Utica-Point Pleasant interval on August 17, 2021 with funding from outside sources. A microseismic monitoring array was installed in the Shoman well and monitored microseismic activity during the Doughty well frac job. Y-Grade NGL injection commenced on August 26, 2021. A total of 215 Bbl. was injected but the job was shut down due to a small leak on the suction hose on the pump truck. The Y-Grade treatment resumed on August 27, 2021 and an additional 726 Bbl. of Y-Grade NGL was injected at a well head pressure of 3850 psi. Total volume of injected Y-Grade over the two days of injection was 941 Bbl. The well was shut-in for 17 days following injection to allow the Y-Grade NGLs to soak on the formation. Y-Grade treatment flow back commenced on 9/13/2021 on a weekly basis. Production data, including surface pressures, oil, nitrogen, natural gas, and flow times was measured and recorded. Periodic gas samples were collected and analyzed to determine composition of flowback gas. As of July 2022, the treatment well had produced 726 Bbl of oil and 2,888 mcf gas. In August 2022, tubing and packer in the well were removed and a pump was installed to enhance oil recovery. Currently, the operator is producing the well about 2 days a week for a few hours. The performance of the NGL treatment test was evaluated based on reservoir simulations of the treatment process, processing of well testing data, analysis of micro seismic monitoring data, and production data analysis. This analysis suggested that oil production in the small test would continue through 2025. Upscaling the treatment to a horizontal Utica Point-Pleasant well would allow more oil recovery, but the process would involve more investment, services, and operational support. An economic analysis was conducted for scenarios aimed at upscaling the NGL treatment process for more typical horizontal UPP wells in the Appalachian Basin.
A newly discovered feature in the midcontinent - a gravity low that begins at a break in the midcontinent gravity high in SE Nebraska, extends across Missouri in a NW-SE direction, and intersects the Mississippi Valley graben to form the Pascola arch - is discussed. The anomaly varies from 120 to 160 km in width, extends approximately 700 km, and is best expressed in southern Missouri, where it has a Bouguer amplitude of about -34 mGal. It is noted that the magnitude of the anomaly cannot be explained on the basis of a thickened section of Paleozoic sedimentary rock. The gravity data and the sparse seismic refraction data for the region are found to be consistent with an increased crustal thickness beneath the gravity low. It is thought that the gravity anomaly is probably the present expression of a failed arm of a rifting event, perhaps one associated with the spreading that led to or preceded formation of the granite and rhyolite terrain of southern Missouri.
Results from various age dating techniques, seismic reflection profiling hydrocarbon maturation studies, and structural analysis were used to evaluate the Cenozoic deformation in the Ruby Mountains and adjoining ranges (pinyon Range and Cortez Range) in Elko and Eureka Counties, Nevada. Age dating techniques used include potassium-argon ages of biotites from granites published by Kistler et al. (1981) and fission track ages from apatite and zircon. Fission track ages from apatite reflect a closing temperature of 100 plus or minus 20 deg C. Zircon fission track ages reflect a closing temperature of 175 plus or minus 25 deg C and potassium-argon ages from brotite reflect a closing temperature of 250 plus or minus 30 deg C. Thus these results allow a reasonably precise tracking of the evolution of the ranges during the Cenozoic. Seismic reflection data are available from Huntington Valley. Access to seismic reflection data directly to the west of the Harrison Pass Pluton in the central Ruby Mountains was obtained. In addition results are available from several deep exploration holes in Huntington Valley.
The SOV/DAS software program focuses on seismic monitoring and data processing of permanent seismic sources Surface Orbital Vibrators (SOVs) and fiber-optics sensing Distributed Acoustic Sensing (DAS). Its key features include automated data processing of the continuous seismic monitoring data acquired with DAS, and output of processed shot gathers and QC plots. The software supports timelapse seismic for long-term reservoir monitoring, and is used for monitoring of geological carbon storage sites, geothermal reservoirs, and oil and gas, and in general subsurface resource management.
Understanding hydraulic fracturing is crucial to improving the stimulation of unconventional reservoirs and increasing fluid production. This study develops a novel seismic monitoring technology using distributed acoustic sensing (DAS) and surface orbital vibrators (SOV) to capture fracture seismic response and mechanical properties at high temporal intervals. We analyze continuous time-lapse vertical seismic profiling (VSP) data acquired every hour during the first nine days of treatment of an unconventional reservoir in the Austin Chalk/Eagle Field Laboratory. The VSP data contain clear seismic signals scattered from the activated fractures. The spatiotemporal changes of the fracture reflectivity revealed by the SOV/DAS data correlate well with the observations of fracture locations inferred from low-frequency DAS data. These results capture the fracture opening and closure processes, as well as highlight potential prestage activations of the fractures due to hydraulic connectivity with preexisting fracture systems. Therefore, analysis of the presented data set provides a unique opportunity to understand fracture initiation and subsequent evolution, not only in the context of unconventional resources but also in enhanced geothermal systems.
Deep moonquakes, occurring at depths about halfway to the center of the Moon, were discovered during the Apollo lunar landing missions, 1969-1972. Their near-monthly occurrence with nearly identical waveforms at any given seismic station suggests that they are strongly influenced by tides, caused by the Earth and the Sun, acting on certain limited regions of the deep lunar interior. However, much about them is unknown, why they are restricted to certain depths and to limited source regions (nests), and what they tell us about the nature of the material and dynamics of the interior of the Moon. A piece of information helpful to decipher their true nature is the spatial extent and distribution of their hypocenters. The occurrence of nearly identical waveforms suggests groups of hypocenters appear in close proximity to one another, but can we tell how closely they are located and how they are distributed? The Apollo PSE (Passive Seismic Experiment) data from all stations of the seismic network were received in real time at a common receiving station on Earth in digital form. This provided extremely high inter-station timing accuracy not achievable for most earthquake data on Earth at that time. We took advantage of this to compute relative locations of waveform-matched events, and concluded that deep moonquake foci in the A1 moonquake nest were concentrated on a nearly horizontal plane of less than 1 km in diameter. This earlier study was based on deep moonquake events visually identified on seismograms. A recent reanalysis of earlier unidentified seismic events fully utilizing the high capability of present-day computers expanded the list of positively identified deep moonquakes by more than a factor of five. The new list contains many events that are not visually matched in waveforms, yet correlated at a significant level when cross-correlated with a computer. Thus it became imperative to reexamine the spatial distribution of deep moonquake hypocenters including the newly identified events. Additional information is included in the original extended abstract.
Included are experimental data recorded from shear experiments that specifically explore the link between fluid-injection rate and seismic moment resulting from shear reactivation of laboratory faults. Raw mechanical data from three experiments are included alongside corresponding MATLAB scripts that import and plot the data, as well as use it to calculate shear and normal stress. Experiments are performed on 2.5-3 inch long granitoid cores from the Utah FORGE EGS demonstration site, containing a single inclined fracture with small-scale roughness added to the fracture surface. The raw data included here were recorded from an aluminum triaxial pressure vessel (TEMCO) configured with three independent servo-controlled pumps, with distilled water used as the working fluid. The pumps control confining pressure, upstream pore pressure, and axial pressure, with each pump connected to a LabView interface to record applied pressures, cumulative injected water volumes, and pump flow rates. The downstream outlet from the fracture is closed to allow pressurization, which is measured by an external pressure transducer. A linear variable differential transformer (LVDT) attached to the axial piston measures axial displacement, from which we calculate shear displacement along the fracture. Additionally, P-wave transducers are used to record acoustic signatures, where acoustic emission events and maximum amplitudes are compared against seismic moment and shear slip velocity. Fluid injection rates range between 0.05 mL/min, 0.25 mL/min, and 0.75 mL/min for each experiment. Along-fault pressure distributions are progressively less uniform as injection rates increase, representing a switch from steady-state to transient conditions. Triggered shear displacement is used as a proxy for seismic moment and is indexed against cumulative injection volume and rate. Each experiment is performed under constant shear stress conditions, and the sample is fully saturated with DI water. Axial and confining stresses are applied to 3 MPa through pressure-stepping in 500 kPa increments. The pore pressure is held constant at 200 kPa prior to initiating the experiment, and initial axial displacement is recorded. The axial stress is then increased to initiate shear mobilization during the loading phase (run-in) until a peak steady state is achieved. The initial shear stress is reduced to approximately 80% of the peak shear stress by decreasing the axial stress, then held constant for the duration of each experiment.
A new 1.5 m diameter impact crater was discovered on Mars only ~40 km from the InSight lander. Context camera images constrained its formation between 21 February and 6 April 2019; follow-up High Resolution Imaging Science Experiment images resolved the crater. During this time period, three seismic events were identified in InSight data. We derive expected seismic signal characteristics and use them to evaluate each of the seismic events. However, none of them can definitively be associated with this source. Atmospheric perturbations are generally expected to be generated during impacts; however, in this case, no signal could be identified as related to the known impact. Using scaling relationships based on the terrestrial and lunar analogs and numerical modeling, we predict the amplitude, peak frequency, and duration of the seismic signal that would have emanated from this impact. The predicted amplitude falls near the lowest levels of the measured seismometer noise for the predicted frequency. Hence, it is not surprising this impact event was not positively identified in the seismic data. Finding this crater was a lucky event as its formation this close to InSight has a probability of only ~0.2, and the odds of capturing it in before and after images are extremely low. We revisit impact-seismic discriminators in light of real experience with a seismometer on the Martian surface. Using measured noise of the instrument, we revise our previous prediction of seismic impact detections downward, from ~a few to tens, to just ~2 per Earth year, still with an order of magnitude uncertainty.
The results of monitoring of carbon dioxide (CO 2 ) injection at the Illinois Basin—Decatur Project (IBDP) and the companion Illinois Industrial Carbon Capture and Sequestration Sources (IL-ICCS) project—have shown that reservoir response to fluid pressure changes can vary significantly at different injection locations within the same reservoir. Predrill reservoir characterization is important to identify potentially seismogenic faults. However, interpretations of newly reprocessed 3D seismic reflection data illustrate the challenges related to their identification in a region dominated by faulting with small vertical offsets. Faults interpreted in the 3D seismic volume range from ~300 to 1200 m wide and are in the same size range as faults that could have been the source of historical events up to Mw 2.7 in central Illinois. The array of monitoring sensors that was installed for the IBDP continues to collect data, as injection operates in IL-ICCS, the second injection well. CO 2 injection rates for the IL-ICCS well are on average 1.7 times the rates injected in the IBDP well, but a significantly reduced rate of induced seismicity is observed. This article presents results of passive seismic monitoring for the duration of the project to date, integrating active and passive seismic data to develop a new interpretation of the subsurface structure at the Decatur site that explicitly identifies pathways for fluid flow into the basement leading to induced seismicity, and provides a geological explanation for the sharp reduction of induced seismicity during injection at higher rates into the second well. The use of seismic moment to estimate the length of seismogenic slip planes in the local subsurface suggests that faults large enough to produce felt seismicity are unlikely to be present at or near the Decatur site.