Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Seismic Data Compression”

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 37 records · Page 2

Alsep data processing: How we processed Apollo Lunar Seismic Data

The Apollo lunar seismic station network gathered data continuously at a rate of 3 x 10 to the 8th power bits per day for nearly eight years until the termination in September, 1977. The data were processed and analyzed using a PDP-15 minicomputer. On the average, 1500 long-period seismic events were detected yearly. Automatic event detection and identification schemes proved unsuccessful because of occasional high noise levels and, above all, the risk of overlooking unusual natural events. The processing procedures finally settled on consist of first plotting all the data on a compressed time scale, visually picking events from the plots, transferring event data to separate sets of tapes and performing detailed analyses using the latter. Many problems remain especially for automatically processing extraterrestrial seismic signals.

Latham, G. V.↗

How we processed Apollo lunar seismic data

The Apollo lunar seismic station network gathered data continuously at a rate of 3 x 10 to the 8th bits per day for nearly eight years until termination in September, 1977. The data were processed and analyzed using a PDP-15 minicomputer. On average, 1500 long-period seismic events were detected yearly. Automatic event detection and identification schemes proved unsuccessful because of occasional high noise levels and, above all, the risk of overlooking unusual natural events. The processing procedures which were finally chosen consist of plotting all the data on a compressed time scale, visually picking events from the plots, transferring event data to separate sets of tapes and performing detailed analyses using the latter. Many problems remain, especially in the automatic processing of extraterrestrial seismic signals.

Nakamura, Y.↗

Surface uplift and time-dependent seismic hazard due to fluid injection in eastern Texas

Observations that unequivocally link seismicity and wastewater injection are scarce. We show that wastewater injection in eastern Texas causes uplift, detectable in radar interferometric data up to >8 kilometers from the wells. Using measurements of uplift, reported injection data, and a poroelastic model, we computed the crustal strain and pore pressure. We infer that an increase of >1 megapascal in pore pressure in rocks with low compressibility triggers earthquakes, including the 4.8–moment magnitude event that occurred on 17 May 2012, the largest earthquake recorded in eastern Texas. Seismic activity increased even while injection rates declined, owing to diffusion of pore pressure from earlier periods with higher injection rates. Finally, induced seismicity potential is suppressed where tight confining formations prevent pore pressure from propagating into crystalline basement rocks.

58 GEOSCIENCES↗

Reducing Disk Storage of Full-3D Seismic Waveform Tomography (F3DT) Through Lossy Online Compression

Full-3D seismic waveform tomography (F3DT) is the latest seismic tomography technique that can assimilate broadband, multi-component seismic waveform observations into high-resolution 3D subsurface seismic structure models. The main drawback in the current F3DT implementation, in particular the scattering-integral implementation (F3DT-SI), is the high disk storage cost and the associated I/O overhead of archiving the 4D space-time wavefields of the receiver- or source-side strain tensors. The strain tensor fields are needed for computing the data sensitivity kernels, which are used for constructing the Jacobian matrix in the Gauss-Newton optimization algorithm. In this study, we have successfully integrated a lossy compression algorithm into our F3DT SI workflow to significantly reduce the disk space for storing the strain tensor fields. The compressor supports a user-specified tolerance for bounding the error, and can be integrated into our finite-difference wave-propagation simulation code used for computing the strain fields. The decompressor can be integrated into the kernel calculation code that reads the strain fields from the disk and compute the data sensitivity kernels. During the wave-propagation simulations, we compress the strain fields before writing them to the disk. To compute the data sensitivity kernels, we read the compressed strain fields from the disk and decompress them before using them in kernel calculations. Experiments using a realistic dataset in our California statewide F3DT project have shown that we can reduce the strain-field disk storage by at least an order of magnitude with acceptable loss, and also improve the overall I/O performance of the entire F3DT-SI workflow significantly. The integration of the lossy online compressor may potentially open up the possibilities of the wide adoption of F3DT-SI in routine seismic tomography practices in the near future.

97 MATHEMATICS AND COMPUTING↗

The Complicated Link Between Material Properties and Microfracture Density for an Underground Explosion in Granite

Geophysical techniques are often implemented as quick and inexpensive ways to locate and characterize fractures in the subsurface, which is important for a number of geoscience fields. Seismic velocities are the most widely used proxies for identification of fractures, but the correlation is not always well-defined. In this study we present material property data: unconfined compressive strength (UCS), bulk density (ρ), Young's modulus (E), Poisson's ratio (ν), P wave velocity (V p ), and S wave velocity (V s ), in conjunction with microfracture densities measured on samples of granite collected before and after underground chemical explosions. Results indicate the relationship between fractures and material properties is complex, even in this single-lithology environment. We interpret that this complexity arises from varying fracture mechanisms (e.g. dilation-inducing fractures vs compression-inducing fractures) in different parts of the core, due to differences in stress conditions. Additional complexity may result from chemical interactions between the fresh fractures and the fluids in the area. Water content appears to have a significant, if not dominant, role in the unconfined compressive strength (UCS) of the samples. We suggest caution when using elastic property measurements as a proxy for fracturing in areas of explosion-induced damage, or in other areas where a variety of mechanisms induce fracturing.

58 GEOSCIENCES↗

Lateral variations in lower mantle seismic velocity

To obtain a theoretical model which provides a rationale for the observed high values of velocity variations, the effect of a 0.1 to 0.2 percent partially molten volatile-rich material in various geometries which are heterogeneously dispersed in the lower mantle is examined. Data obtained indicate that, depending on aspect ratio and geometry, 0.1-0.2 percent partial melting in conjunction with about 100 K thermal anomalies can explain the seismic variations provided the compressibility of the melt differs by less than about 20 percent from the surrounding solid.

Duffy, Thomas S.↗

Unconfined Compression Results on Core from Borehole U-15n#10, NNSS, in support of SPE

The Nevada National Security Site (NNSS) serves as the geologic setting for a Source Physics Experiment (SPE) program. The SPE provides ground truth data to create and improve strong ground motion and seismic S-wave generation and propagation models. The NNSS was chosen as the test bed because it provides a variety of geologic settings ranging from relatively simple to very complex. Each series of SPE testing will comprise the setting and firing of explosive charges (source) placed in a central borehole at varying depths and recording ground motions in instrumented boreholes located in two rings around the source, positioned at different radii. Modeling using advanced simulation codes will be performed both before and after each test to predict ground response and to improve models based on acquired field data, respectively. A key component in the predictive capability and ultimate validation of the models is the full understanding of the intervening geology between the source and the instrumented boreholes including the geomechanical behavior of the site's rock/structural features. This memorandum reports on an initial phase of unconfined compression testing from corehole U-15n#10. Specimens tested came from the U-15n#10 core hole, which was drilled at -60° to the horizontal in granitic rock (quartz monzonite) after the second SPE shot (SPE-2). Figure 1 illustrates at the surface, the core hole was approximately 90 feet from the central SPE borehole. Corehole U 15n#10 intersects the central SPE borehole (U-15n) at a core depth of 170 feet (approximately 150 feet vertical depth) which is within the highly damaged zone of SPE-2. The U-15n#10 location is the site of the first, second and third SPE's, in Area 15 of the NNSS.

36 MATERIALS SCIENCE↗

Unconfined Compression Results on Core from Boreholes U-15n#12 and U-15n#13, NNSS in support of SPE

The Nevada National Security Site (NNSS) serves as the geologic setting for a Source Physics Experiment (SPE) program. The SPE provides ground truth data to create and improve strong ground motion and seismic S-wave generation and propagation models. The NNSS was chosen as the test bed because it provides a variety of geologic settings ranging from relatively simple to very complex. Each series of SPE testing will comprise the setting and firing of explosive charges (source) placed in a central borehole at varying depths and recording ground motions in instrumented boreholes located in two rings around the source, positioned at different radii. Modeling using advanced simulation codes will be performed both before and after each test to predict ground response and to improve models based on acquired field data, respectively. A key component in the predictive capability and ultimate validation of the models is the full understanding of the intervening geology between the source and the instrumented boreholes including the geomechanical behavior of the site's rock/structural features. This report summarizes unconfined compression testing (UCS) from coreholes U-15n#12 and U-15n#13 and compares those datasets to UCS results from coreholes U-15n and U-15n#10. U-15n#12 corehole was drilled at -60° to the horizontal and U-15n#13 was drilled vertically in granitic rock (quartz monzonite) after the third SPE shot. Figure 1 illustrates at the surface, U 15n#12 and U-15n#13 coreholes were approximately 30 meters and 10 meters from the central SPE borehole (U-15n) respectively. Corehole U-15n#12 intersects the central SPE borehole (U 15n) at a core depth of 174 feet (approximately 150 feet vertical depth). The location of U 15n#12 and U-15n#13 is the site of the first, second and third SPE's, in Area 15 of the NNSS.

58 GEOSCIENCES↗

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↗

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↗

Sound velocity of Fe-S liquids at high pressure: Implications for the Moon's molten outer core

Sound velocities of Fe and Fe–S liquids were determined by combining the ultrasonic measurements and synchrotron X-ray techniques under high pressure–temperature conditions from 1 to 8 GPa and 1573 K to 1973 K. Four different liquid compositions were studied including Fe, Fe–10 wt% S, Fe–20 wt% S, and Fe–27 wt% S. Our data show that the velocity of Fe-rich liquids increases upon compression and decreases with increasing sulfur content, whereas temperature has negligible effect on the velocity of Fe–S liquids. The sound velocity data were combined with ambient-pressure densities to fit the Murnaghan equation of state (EOS). Compared to the lunar seismic model, our velocity data constrain the sulfur content at 4±3 wt%, indicating a significantly denser (6.5±0.5 g/cm 3 ) and hotter (1870-70+100 K) outer core than previously estimated. A new lunar structure model incorporating available geophysical observations points to a smaller core radius. Our model suggests a top–down solidification scenario for the evolution of the lunar core. Such “iron snow” process may have been an important mechanism for the growth of the inner core.

Jing, Zhicheng↗

An autoencoder compression approach for accelerating large-scale inverse problems

Abstract Partial differential equation (PDE)-constrained inverse problems are some of the most challenging and computationally demanding problems in computational science today. Fine meshes required to accurately compute the PDE solution introduce an enormous number of parameters and require large-scale computing resources such as more processors and more memory to solve such systems in a reasonable time. For inverse problems constrained by time-dependent PDEs, the adjoint method often employed to compute gradients and higher order derivatives efficiently requires solving a time-reversed, so-called adjoint PDE that depends on the forward PDE solution at each timestep. This necessitates the storage of a high-dimensional forward solution vector at every timestep. Such a procedure quickly exhausts the available memory resources. Several approaches that trade additional computation for reduced memory footprint have been proposed to mitigate the memory bottleneck, including checkpointing and compression strategies. In this work, we propose a close-to-ideal scalable compression approach using autoencoders to eliminate the need for checkpointing and substantial memory storage, thereby reducing the time-to-solution and memory requirements. We compare our approach with checkpointing and an off-the-shelf compression approach on an earth-scale ill-posed seismic inverse problem. The results verify the expected close-to-ideal speedup for the gradient and Hessian-vector product using the proposed autoencoder compression approach. To highlight the usefulness of the proposed approach, we combine the autoencoder compression with the data-informed active subspace (DIAS) prior showing how the DIAS method can be affordably extended to large-scale problems without the need for checkpointing and large memory.

Mathematics↗

Missile impact craters (White Sands Missile Range, New Mexico) and applications to lunar research: Contributions to astrogeology

Craters in natural materials at White Sands Missile Range, N. Mex., were produced by the impact of high-velocity to hypervelocity missiles traveling along oblique trajectories with kinetic energies between 2.1 and 81 × 1014 ergs. The oblique impacts produce craters 2 to 10 m across with morphologies and ejecta that are bilaterally symmetrical with respect to the plane of the missile trajectory. Rims are high and the amount of ejecta large in down-trajectory and lateral directions, whereas rims are low to nonexistent and ejecta thin to absent up-trajectory. Symmetry development and modifications of the symmetry are a function of target material, local topography, and angle of impact. Seven mappable units can be recognized in and around the craters. Three of these are ejecta: thick ejecta near the crater, thin to discontinuous ejecta at greater distances, and scattered ejecta at the greatest distances to the limit of throwout. These ejecta units may be absent on the up-trajectory side; if present, they are rarely as thick or continuous as on other sides of the crater. Three units are target materials: undeformed target material exposed in local patches through thin to discontinuous ejecta and everywhere between the fragments of scattered ejecta, tilted and broken target material exposed in upper crater walls, and shattered and fractured target material exposed on the up-trajectory crater wall. The seventh unit is slope material composed of talus and fallback within the crater. Development, character, and exposure of these units varies chiefly with the target material. Ejecta from the craters is chiefly broken but relatively undeformed target material that may range in size from very fine grained debris to large blocks. Where the target is porous, significant amounts of the ejecta are composed of sheared and compressed fragments, some coated with dark layers of mixed projectile pieces, powder, and fused metal mixed with crushed target material. For layered targets, the original stratigraphic sequence is crudely preserved and in inverted order in thick ejecta. Secondary impact craters are produced by the impact of ejected fragments when the surrounding surface materials are sufficiently weak. A wide variety of secondary impact crater relations may result. Secondary craters nearest the primary crater have blocks in them that are larger than or the same size as the crater they produced. Farther from the primary crater, the fragments are generally smaller than the secondary crater and are ejected from it. Excavation of four craters revealed a mixed breccia beneath the crater floor composed of missile pieces, sheared and compressed target material, and crushed debris. Banded, disaggregated target material and nonmixed breccia surrounded the mixed breccia, and these breccias were surrounded by a zone of conjugate fractures. Beneath the ejecta on the lateral and down-trajectory crater flanks, the target materials were tilted upward and broken. Up-trajectory, open fractures and downward displacement occurred in two of the craters. No displacement was observed for the other two. Beneath the down-trajectory rims of craters with distinct layering, overturned synclines were observed. Missile breakup and behavior during cratering are a function of target and missile properties. Missile breakup depends on missile velocity and is most extensive at high velocities, where the missile is fragmented, powdered, and partly fused. Burial of missile or its fragmented, powdered, and fused remains is greatest for porous targets and least for dense cohesive targets. For very porous targets, camouflet structures containing the fragmented missile may form. Least squares fit to the data on craters in dry to moist targets indicate V(a) = 10^(-11.433)E(p)^(1.205) where V(a) is the volume of the apparent crater and E(p) is the kinetic energy of the missile. This equation is consistent with expectations of the equations relating apparent depth and radius to kinetic energy. Extrapolation of displaced masses and kinetic energies for laboratory impacts with sand and rock converge near 10^(15) to 10^(16) ergs, where the extrapolations are near the data on missile impact craters, corrected for impact angle. Displaced masses of craters produced by missile impacts and by chemical explosives with small scaled depths of burial are about the same when the kinetic energies of the missiles (corrected for angle of impact) are equal to the TNT equivalent energy of the explosive. The problem of equivalent scaled depth of burst for an impact crater is complicated and not entirely resolved, however. Both missile impact craters and chemical explosive craters in water-saturated targets are larger than their counterparts in dry to moist materials. Data collected during the study of missile impact craters have helped resolve a number of problems in lunar research: (1) the soillike nature of lunar surface materials was predicted, (2) sizes of craters produced by artificial impacts were correctly predicted, (3) certain features imaged by Surveyor were found to be analogous to features associated with missile impact craters, {4) missile impacts were used in support of the Apollo passive seismic experiment, (5) craters seen in Apollo orbital photographs were found to be similar to some missile impact craters, (6) missile impact craters supplied data on sample collection and crater phenomenology used in training astronauts, and (7) some returned lunar samples are similar to coated, sheared, and compressed fragments ejected from missile impact craters.

H. J. Moore↗

Strength and compressibility of returned lunar soil.

Two oedometer and three direct shear tests have been performed in vacuum on a 200 g sample of lunar soil from Apollo 12 (12001, 119). The compressibility data have been used to calculate bulk density and shear wave velocity versus depth on the lunar surface. The shear wave velocity was found to increase approximately with the one-fourth power of the depth, and the results suggest that the Apollo 14 Active Seismic Experiment may not have detected the Fra Mauro formation at a depth of 8.5 m, but only naturally consolidated lunar soil. The shear data indicate that the strength of the lunar soil sample is about 65% that of a ground basalt simulant at the same void ratio.

Carrier, W. D., III↗

Data User's Note: Apollo seismological investigations

Seismological objectives and equipment used in the passive seismic, active seismic, lunar seismic profiling, and the lunar gravimeter experiments conducted during Apollo 11, 12, 14, 15, 16, and 17 missions are described. The various formats in which the data form these investigations can be obtained are listed an an index showing the NSSDC identification number is provided. Tables show manned lunar landing missions, lunar seismic network statistics, lunar impact coordinate statistics, detonation masses and times of EP's, the ALSEP (Apollo 14) operational history; compressed scale playout tape availability, LSPE coverage for one lunation, and experimenter interpreted events types.

Vostreys, R. W.↗

Geodynamics of the Baikal-Stanovoy seismic belt

The convection generated tensional stress field in the earth, as inferred from satellite gravity data, reveals an anomalous lens of upwelling mantle rocks under the Baikal rift zone. The point of no strain at 56 deg N 116 deg E forms a seismic gap along the Baikal-Stanovoy seismic belt. East of this point, the stress field changes from extension to compression. Therefore, the position of no strain at the eastern termination of the rift accounts for the dying-out of the rift zone and for the appearance of a compressive structure in the Stanovoy Range.

Liu, H.-S.↗

Anomalous density and elastic properties of basalt at high pressure: Reevaluating of the effect of melt fraction on seismic velocity in the Earth's crust and upper mantle

Independent measurements of the volumetric and elastic properties of Columbia River basalt glass were made up to 5.5 GPa by high-pressure X-ray microtomography and GHz-ultrasonic interferometry, respectively. The Columbia River basalt displays P and S wave velocity minima at 4.5 and 5 GPa, respectively, violating Birch’s law. These data constrain the pressure dependence of the density and elastic moduli at high pressure, which cannot be modeled through usual equations of state nor determined by stepwise integrating the bulk sound velocity as is common practice. We propose a systematic variation in compression behavior of silicate glasses that is dependent on the degree of polymerization and arises from the flexibility of the aluminosilicate network. Likewise, this behavior likely persists into the liquid state for basaltic melts resulting in weak pressure dependence for P wave velocities perhaps to depths of the transition zone. By modeling the effect of partial melt on P wave velocity reductions it is suggested that melt fraction determined by seismic velocity variations may be significantly overestimated in the crust and upper mantle.

58 GEOSCIENCES↗