Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Algorithms”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57Linked to original sources
USGS1989

A terracing operator for physical property mapping with potential field data

The terracing operator works iteratively on gravity or magnetic data, using the sense of the measured field's local curvature, to produce a field comprised of uniform domains separated by abrupt domain boundaries. The result is crudely proportional to a physical-property function defined in one (profile case) or two (map case) horizontal dimensions. This result can be extended to a physical-property model if its behavior in the third (vertical) dimension is defined, either arbitrarily or on the basis of the local geologic situation. The terracing algorithm is computationally fast and appropriate to use with very large digital data sets. Where gravity and magnetic data are both available, terracing provides an effective means by which the two data sets can be compared directly. Results of the terracing operation somewhat resemble those of conventional susceptibility (or density) mapping. In contrast with conventional susceptibility mapping, however, the terraced function is a true step function, which cannot be depicted by means of contour lines. Magnetic or gravity fields calculated from the physical-property model do not, in general, produce an exact fit to the observed data. By intent, the terraced map is more closely analogous to a geologic map in that domains are separated by hard-edged domain boundaries and minor within-domain variation is neglected. The terracing operator was applied separately to aeromagnetic and gravity data from a 136 km X 123 km area in eastern Kansas. Results provide a reasonably good physical representation of both the gravity and the aeromagnetic data. Superposition of the results from the two data sets shows many areas of agreement that can be referenced to geologic features within the buried Precambrian crystalline basement. The emerging picture of basement geology is much better resolved than that obtained either from the scanty available drill data or from interpretation of the geophysical data by inspection.

Geophysics↗
USGS2002

Time‐lapse inversion of crosswell radar data

The combination of differential radar tomography with conventional tracer and/or hydraulic tests facilitates high‐resolution characterization of subsurface heterogeneity and enables the identification of preferential flow paths. In dynamic imaging, each tomogram is typically inverted independently, under the assumption that data sets are collected quickly relative to changes in the imaged property (e.g., attenuation or velocity); however, such “snapshot” tomograms may contain large errors if the imaged property changes significantly during data collection. Acquisition of less data over a shorter time interval could ameliorate the problem, but the resulting decrease in ray density and angular coverage could degrade model resolution. To address these problems, we propose a new sequential approach for time‐lapse tomographic inversion. The method uses space‐time parameterization and regularization to combine data collected at multiple times and to account for temporal variation. The inverse algorithm minimizes the sum of weighted squared residuals and a measure of solution complexity based on an a priori space‐time covariance function and a spatiotemporally variable mean. We demonstrate our approach using a synthetic 2‐D time‐lapse ( x , z , t ) data set based loosely on a field experiment in which difference‐attenuation radar tomography was used to monitor the migration of a saline tracer in fractured rock. We quantitatively show the benefits of space‐time inversion by comparing results for snapshot and time‐lapse inversion schemes. Inversion over both space and time results in superior estimation error, model resolution, and data reproduction compared to conventional snapshot inversion. Finally, we suggest strategies to improve time‐lapse cross‐hole inversions using ray‐based inversion constraints and a modified survey design in which different sets of rays are collected in alternating time steps.

Geophysics↗
USGS2006

Seismoelectric numerical modeling on a grid

Our finite-difference algorithm provides a new method for simulating how seismic waves in arbitrarily heterogeneous porous media generate electric fields through an electrokinetic mechanism called seismoelectric coupling. As the first step in our simulations, we calculate relative pore-fluid/grain-matrix displacement by using existing poroelastic theory. We then calculate the electric current resulting from the grain/fluid displacement by using seismoelectric coupling theory. This electrofiltration current acts as a source term in Poisson's equation, which then allows us to calculate the electric potential distribution. We can safely neglect induction effects in our simulations because the model area is within the electrostatic near field for the depth of investigation (tens to hundreds of meters) and the frequency ranges (10 Hz to 1 kHz) of interest for shallow seismoelectric surveys.We can independently calculate the electric-potential distribution for each time step in the poroelastic simulation without loss of accuracy because electro-osmotic feedback (fluid flow that is perturbed by generated electric fields) is at least 105 times smaller than flow that is driven by fluid-pressure gradients and matrix acceleration, and is therefore negligible. Our simulations demonstrate that, distinct from seismic reflections, the seismoelectric interface response from a thin layer (at least as thin as one-twentieth of the seismic wavelength) is considerably stronger than the response from a single interface. We find that the interface response amplitude decreases as the lateral extent of a layer decreases below the width of the first Fresnel zone. We conclude, on the basis of our modeling results and of field results published elsewhere, that downhole and/or crosswell survey geometries and time-lapse applications are particularly well suited to the seismoelectric method. ?? 2006 Society of Exploration Geophysicists.

Geophysics↗
USGS2007

Applications of the JARS method to study levee sites in southern Texas and southern New Mexico

We apply the joint analysis of refractions with surface waves (JARS) method to several sites and compare its results to traditional refraction-tomography methods in efforts of finding a more realistic solution to the inverse refraction-traveltime problem. The JARS method uses a reference model, derived from surface-wave shear-wave velocity estimates, as a constraint. In all of the cases JARS estimates appear more realistic than those from the conventional refraction-tomography methods. As a result, we consider, the JARS algorithm as the preferred method for finding solutions to the inverse refraction-tomography problems. ?? 2007 Society of Exploration Geophysicists.

Conference Paper↗
USGS2008

Using HEM surveys to evaluate disposal of by-product water from CBNG development in the Powder River Basin, Wyoming

Production of methane from thick, extensive coal beds in the Powder River Basin of Wyoming has created water management issues. Since development began in 1997, more than 650 billion liters of water have been produced from approximately 22,000 wells. Infiltration impoundments are used widely to dispose of by-product water from coal bed natural gas (CBNG) production, but their hydrogeologic effects are poorly understood. Helicopter electromagnetic surveys (HEM) were completed in July 2003 and July 2004 to characterize the hydrogeology of an alluvial aquifer along the Powder River. The aquifer is receiving CBNG produced water discharge from infiltration impoundments. HEM data were subjected to Occam's inversion algorithms to determine the aquifer bulk conductivity, which was then correlated towater salinity using site-specific sampling results. The HEM data provided high-resolution images of salinity levels in the aquifer, a result not attainable using traditional sampling methods. Interpretation of these images reveals clearly the produced water influence on aquifer water quality. Potential shortfalls to this method occur where there is no significant contrast in aquifer salinity and infiltrating produced water salinity and where there might be significant changes in aquifer lithology. Despite these limitations, airborne geophysical methods can provide a broadscale (watershed-scale) tool to evaluate CBNG water disposal, especially in areas where field-based investigations are logistically prohibitive. This research has implications for design and location strategies of future CBNG water surface disposal facilities within the Powder River Basin. ?? 2008 2008 Society of ExplorationGeophysicists. All rights reserved.

Geophysics↗
USGS2008

ALLTEM UXO detection and discrimination

ALLTEM is a multi-axis electromagnetic induction system designed for unexploded ordnance (UXO) applications. It uses a continuous triangle-wave excitation and provides good late-time signal-to-noise ratio (SNR) especially for ferrous targets. Multi-axis transmitter (Tx) and receiver (Rx) systems such as ALLTEM provide a richer data set from which to invert for the target parameters required to distinguish between clutter and UXO. Inversions of field data over the Army's UXO Calibration Grid and Blind Test Grid at the Yuma Proving Ground (YPG), Arizona in 2006 produced polarizability moment values for many buried UXO items that were reasonable and generally repeatable for targets of the same type buried at different orientations and depths. In 2007 a test stand was constructed that allows for collection of data with varying spatial data density and accurate automated position control. The behavior of inverted ALLTEM test stand data as a function of spatial data density, sensor SNR, and position error has been investigated. The results indicate that the ALLTEM inversion algorithm is more tolerant of sensor noise and position error than has been reported for single-axis systems. A high confidence level in inversion-derived target parameters is required when a target is declared to be harmless scrap metal that may safely be left in the ground. Unless high confidence can be demonstrated, state regulators will likely require that targets be dug regardless of any "no-dig" classifications produced from inversions, in which case remediation costs would not be decreased.

Conference Paper↗
USGS2010

Improved hydrogeophysical characterization and monitoring through parallel modeling and inversion of time-domain resistivity andinduced-polarization data

Electrical geophysical methods have found wide use in the growing discipline of hydrogeophysics for characterizing the electrical properties of the subsurface and for monitoring subsurface processes in terms of the spatiotemporal changes in subsurface conductivity, chargeability, and source currents they govern. Presently, multichannel and multielectrode data collections systems can collect large data sets in relatively short periods of time. Practitioners, however, often are unable to fully utilize these large data sets and the information they contain because of standard desktop-computer processing limitations. These limitations can be addressed by utilizing the storage and processing capabilities of parallel computing environments. We have developed a parallel distributed-memory forward and inverse modeling algorithm for analyzing resistivity and time-domain induced polar-ization (IP) data. The primary components of the parallel computations include distributed computation of the pole solutions in forward mode, distributed storage and computation of the Jacobian matrix in inverse mode, and parallel execution of the inverse equation solver. We have tested the corresponding parallel code in three efforts: (1) resistivity characterization of the Hanford 300 Area Integrated Field Research Challenge site in Hanford, Washington, U.S.A., (2) resistivity characterization of a volcanic island in the southern Tyrrhenian Sea in Italy, and (3) resistivity and IP monitoring of biostimulation at a Superfund site in Brandywine, Maryland, U.S.A. Inverse analysis of each of these data sets would be limited or impossible in a standard serial computing environment, which underscores the need for parallel high-performance computing to fully utilize the potential of electrical geophysical methods in hydrogeophysical applications.

Geophysics↗
USGS2010

Estimating structural dip from gravity and magnetic profile data

The anomalies of ideal sources, dipping magnetic contacts and dipping density sheets, provide the theoretical basis for estimating structural dip from gravity and magnetic profile data. The dip is always related to the local phase angle of a complex analytic signal evaluated directly over the source. For magnetic sheets, the complex analytic signal is constructed from the anomaly and its Hilbert transform. For magnetic contacts or density sheets, it is constructed from the first‐order horizontal and vertical derivatives of the anomaly. Dip estimation is implemented as an integral part of source‐location algorithms. The analytic signal and local wavenumber methods require interpolation of the local phase from sampled anomaly or derivative values. In the case of extended Euler deconvolution and multiple‐source Werner deconvolution the dip is derived as part of the least squares solution.

Conference Paper↗
USGS2011

Classification of MEC with the ALLTEM at Camp Stanley, Texas

The ALLTEM is a multi-axis electromagnetic induction system designed for unexploded ordnance UXO applications. It uses a continuous triangle-wave excitation and provides good late-time signal-to-noise ratio SNR especially for ferrous targets. Multi-axis transmitter Tx and receiver Rx systems such as ALLTEM provide a richer data set from which to invert for the target parameters required to distinguish between clutter and UXO. Inversions of field data acquired between 2006 and 2010 over the Army's UXO Standardized Test sites at the Yuma Proving Ground YPG in Arizona and at the Aberdeen Proving Ground APG in Maryland have produced reasonable and generally repeatable results for many UXO items buried at different orientations and depths. In February-March 2011 ALLTEM data was acquired at two locations on the Camp Stanley Storage Activity CSSA just north of San Antonio, Texas. Camp Stanley is used to store munitions as well as test, fire, and overhaul munitions components. Site B-20 is an open burn/open detonation OBOD area and Site B-27 consists of narrow trenches blasted into limestone containing buried range and munitions debris and possibly MEC. The processing, analysis, and classification techniques developed at the controlled environments of YPG and APG have been applied to these two "live" sites at Camp Stanley. ALLTEM data analysis includes both classical numerical inversion of data from each anomaly and clustering of the raw data by means of a self-organizing map SOM via generalized neural network algorithms. Final classification consists of an integration of both the numerical and SOM results. ?? 2011 Society of Exploration Geophysicists.

SEG Technical Program Expanded Abstracts↗
USGS2011

Imaging using cross-hole seismoelectric tomography

We propose a new cross-hole imaging approach based on seismoelectric conversions associated with the transmission of seismic waves from seismic sources located in a borehole to receivers electrodes located in a second borehole. The seismoelectric seismic-to-electric problem is solved using Biot theory coupled with a generalized Ohm's law with an electrokinetic coupling term. The components of the displacement of the solid phase, the fluid pressure, and the electrical potential are solved using a finite element approach with PML boundary conditions for the seismic waves and boundary conditions mimicking an infinite material for the electrostatic problem. We have developed an inversion algorithm using the electrical disturbances recorded in the second borehole to localize the position of the heterogeneities responsible for the seismoelectric conversions. Because of the ill-posed nature of the inverse problem, regularization is used to constrain the solution at each time in the seismoelectric time window comprised between the time of the seismic shot and the time of the first arrival of the seismic waves in the second borehole. All the inverted volumetric current source densities are stacked to produce an image of the position of the heterogeneities between the two boreholes. Two simple synthetic case studies are presented to test this concept. The first case study corresponds to a vertical discontinuity between two homogeneous sub-domains. The second case study corresponds to a poroelastic inclusion embedded into an homogenous poroelastic formation. In both cases, the position of the heterogeneity is fairly well-recovered using only the electrical disturbances associated with the seismoelectric conversions. ?? 2011 Society of Exploration Geophysicists.

SEG Technical Program Expanded Abstracts↗
USGS2014

Multielevation calibration of frequency-domain electromagnetic data

Systematic calibration errors must be taken into account because they can substantially impact the accuracy of inverted subsurface resistivity models derived from frequency-domain electromagnetic data, resulting in potentially misleading interpretations. We have developed an approach that uses data acquired at multiple elevations over the same location to assess calibration errors. A significant advantage is that this method does not require prior knowledge of subsurface properties from borehole or ground geophysical data (though these can be readily incorporated if available), and is, therefore, well suited to remote areas. The multielevation data were used to solve for calibration parameters and a single subsurface resistivity model that are self consistent over all elevations. The deterministic and Bayesian formulations of the multielevation approach illustrate parameter sensitivity and uncertainty using synthetic- and field-data examples. Multiplicative calibration errors (gain and phase) were found to be better resolved at high frequencies and when data were acquired over a relatively conductive area, whereas additive errors (bias) were reasonably resolved over conductive and resistive areas at all frequencies. The Bayesian approach outperformed the deterministic approach when estimating calibration parameters using multielevation data at a single location; however, joint analysis of multielevation data at multiple locations using the deterministic algorithm yielded the most accurate estimates of calibration parameters. Inversion results using calibration-corrected data revealed marked improvement in misfit, lending added confidence to the interpretation of these models.

Geophysics↗
USGS2022

Velocity modeling of supercritical pore fluids through porous media under reservoir conditions with applications for petroleum secondary migration and carbon sequestration plumes

Computational methods to characterize secondary migration in porous media traditionally rely on fluid transport equations with assumptions of time invariance, such as flowpath modeling of buoyancy vectors, statistical percolation algorithms, capillary pressure curves, or a form of Darcy’s Law which presumes instantaneous fluid transport. However, in petroleum systems modeling, the timeframe of secondary migration from source to reservoir is important to quantify in relation to other geologic factors such as timing of petroleum generation, fault movement, and seal formation. Additionally, quantifying migration velocities enables an estimation of the distance a plume of geologically sequestered carbon dioxide travels over time, as well as the identification of low-permeability strata appropriate for long-term containment. This study introduces a method to quantify transport velocities of supercritical fluids in low-permeability lithologies for a broad range of rock and fluid properties likely encountered in the sedimentary sequence. A time-dependent form of Darcy’s Law for pressure-driven viscous flow through homogeneous isotropic porous media was used to model flow velocities within a carrier bed. Thermodynamic equations of state were used to determine thermophysical properties of supercritical pore fluids under reservoir pressures ranging from 0–200 MPa (0–29,000 psi) to constrain the momentum equations. Three case studies were examined that (1) estimated fluid flow velocities of methane within the low-permeability Upper Jurassic Haynesville Formation, (2) defined permeability-based flow units to evaluate saline formations for long-term geologic carbon sequestration, and (3) calculated the migration distance of carbon dioxide plumes at the Decatur, Illinois injection and sequestration project.

SEG-AAPG Interpretation↗
USGS2022

Neural net detection of seismic features related to gas hydrates and free gas accumulations on the northern U.S. Atlantic margin

Bottom-simulating reflections (BSRs) that sometimes mark the base of the gas hydrate stability zone in marine sediments are often identified based on the reverse polarity reflections that cut across stratigraphic layering in seismic amplitude data. On the northern U.S. Atlantic margin (USAM) between Cape Hatteras and Hudson Canyon, legacy seismic data have revealed pronounced BSRs south of the deepwater extension of Hudson Canyon and more subtle ones from offshore Delaware south to Cape Hatteras, where the reflections sometimes follow stratigraphic layering. Using high-resolution seismic data acquired during the 2018 Mid-Atlantic Resource Imaging Experiment and a supervised neural net, we identify seismic features associated with gas hydrates and/or the top of gas between Hudson Canyon and Cape Hatteras. Using seismic attributes especially sensitive to the presence of gas, we train a neural network algorithm on seismic data from an area with strong BSRs and then apply the model to the rest of the data set. The results indicate that gas hydrate and/or shallow free gas are significantly more widespread on the northern part of the USAM than previously known. Seismic indicators of gas extend landward from the 2000 m isobath to the upper continental slope in sectors with (offshore Virginia) and, to a lesser extent, without (offshore New Jersey) pervasive upper slope methane seeps. Higher sand content and intermediate sediment thickness, factors related to the container size and gas charge in a petroleum systems framework, are associated with more robust gas indicators.

U.S. Atlantic margin↗
USGS2018

A flatfile of ground motion intensity measurements from induced earthquakes in Oklahoma and Kansas

We have produced a uniformly processed database of orientation-independent (RotD50, RotD100) ground motion intensity measurements containing peak horizontal ground motions (accelerations and velocities) and 5-percent-damped pseudospectral accelerations (0.1–10 s) from more than 3,800 M ≥ 3 earthquakes in Oklahoma and Kansas that occurred between January 2009 and December 2016. Ground motion time series were collected from regional, national, and temporary seismic arrays out to 500 km. We relocated the majority of the earthquake hypocenters using a multiple-event relocation algorithm to produce a set of near-uniformly processed hypocentral locations. Ground motion processing followed standard methods, with the primary objective of reducing the effects of noise on the measurements. Regional wave-propagation features and the high seismicity rate required careful selection of signal windows to ensure that we captured the entire ground motion record and that contaminating signals from extraneous earthquakes did not contribute to the database. Processing was carried out with an automated scheme and resulted in a database comprising more than 174,000 records ( https://dx.doi.org/10.5066/F73B5X8N ). We anticipate that these results will be useful for improved understanding of earthquake ground motions and for seismic hazard applications.

Kansas, Oklahoma↗
USGS2018

Correcting spacecraft jitter in HiRISE images

Mechanical oscillations or vibrations on spacecraft, also called pointing jitter, cause geometric distortions and/or smear in high-resolution digital images acquired from orbit. Geometric distortion is especially a problem with pushbroom sensors, such as the High Resolution Imaging Science Experiment (HiRISE) instrument on-board the Mars Reconnaissance Orbiter (MRO). Geometric distortions occur at a range of frequencies that may not be obvious in the image products, but can cause problems with stereo image correlation in the production of digital elevation models, and in measuring surface changes in time series with orthorectified images. The HiRISE focal plane comprises a staggered array of fourteen charge-coupled devices (CCDs) with pixel instantaneous field of view (IFOV) of 1 microradian. The high spatial resolution of HiRISE makes it both sensitive to, and an excellent recorder of jitter. We present an algorithm using Fourier analysis to resolve the jitter function for a HiRISE image that is then used to update instrument pointing information to remove geometric distortions from the image. Implementation of the jitter analysis and image correction is performed on selected HiRISE images made available to the public. Results show marked reduction of geometric distortions. This work has applications to similar cameras operating now (such as the Lunar Reconnaissance Orbiter Camera Narrow Angle Camera (LROC NAC) on-board the Lunar Reconnaissance Orbiter) and to the design of future instruments (such as the Europa Imaging System, planned for the Europa Clipper mission).

Book chapter↗
USGS2023

Prioritization of species status assessments for decision support

Species status assessments are used to inform U.S. Fish and Wildlife Service (USFWS) decision making for Endangered Species Act (ESA) classification decisions, recovery planning, and more. The large number of species that require assessment and uncertainty in the data available impede the process of assigning and completing the assessments, which makes creating a multiyear work plan extremely difficult. An optimized triaging system that maximizes the use of the best available information while managing the complex ESA workload and meeting deadlines is necessary. We used a structured decision-making framework to approach the problem with the goal of creating a prioritization tool that would be effective at scheduling assessments, given the best information available and priorities of the USFWS. We collected data on the species awaiting assessment and developed a value function that incorporates existing deadlines, taxonomic uncertainty, controversy of the species, and population and habitat data availability and quality. We used a constrained linear optimization algorithm to maximize the value function and ensure that workload capacity was not exceeded. A comparison of model scenarios indicates that imposed deadlines impact the model more than capacity constraints. Additionally, differential weighting of the metrics significantly affected the outcome of the model. In the future, elicitation of metric weights should be done routinely before the model is run for use in official planning to ensure alignment with current USFWS priorities. Output from this optimization can be used to inform a five-year work plan, allocate resources, and discuss workforce decisions.

Decision Analysis↗
USGS2010

Correction for the 17O interference in δ(13C) measurements when analyzing CO2 with stable isotope mass spectrometry

Measurements of δ( 13 C) determined on CO 2 with an isotope-ratio mass spectrometer (IRMS) must be corrected for the amount of 17 O in the CO 2 . For data consistency, this must be done using identical methods by different laboratories. This report aims at unifying data treatment for CO 2 IRMS by proposing (i) a unified set of numerical values, and (ii) a unified correction algorithm, based on a simple, linear approximation formula. Because the oxygen of natural CO 2 is derived mostly from the global water pool, it is recommended that a value of 0.528 be employed for the factor λ, which relates differences in 17 O and 18 O abundances. With the currently accepted N( 13 C)/N( 12 C) of 0.011 180(28) in VPDB (Vienna Peedee belemnite) reevaluation of data yields a value of 0.000 393(1) for the oxygen isotope ratio N( 17 O)/N( 16 O) of the evolved CO 2 . The ratio of these quantities, a ratio of isotope ratios, is essential for the 17 O abundance correction: [N( 17 O)/N( 16 O)]/[N( 13 C)/N( 12 C)] = 0.035 16(8). The equation [δ( 13 C) ≈ 45 δ VPDB-CO2 + 2 17 R/ 13 R ( 45 δ VPDB-CO2 – λ 46 δ VPDB-CO2 )] closely approximates δ( 13 C) values with less than 0.010 ‰ deviation for normal oxygen-bearing materials and no more than 0.026 ‰ in extreme cases. Other materials containing oxygen of non-mass-dependent isotope composition require a more specific data treatment. A similar linear approximation is also suggested for δ( 18 O). The linear approximations are easy to implement in a data spreadsheet, and also help in generating a simplified uncertainty budget.

Pure and Applied Chemistry↗
USGS2012

An accessible method for implementing hierarchical models with spatio-temporal abundance data

A common goal in ecology and wildlife management is to determine the causes of variation in population dynamics over long periods of time and across large spatial scales. Many assumptions must nevertheless be overcome to make appropriate inference about spatio-temporal variation in population dynamics, such as autocorrelation among data points, excess zeros, and observation error in count data. To address these issues, many scientists and statisticians have recommended the use of Bayesian hierarchical models. Unfortunately, hierarchical statistical models remain somewhat difficult to use because of the necessary quantitative background needed to implement them, or because of the computational demands of using Markov Chain Monte Carlo algorithms to estimate parameters. Fortunately, new tools have recently been developed that make it more feasible for wildlife biologists to fit sophisticated hierarchical Bayesian models (i.e., Integrated Nested Laplace Approximation, ‘INLA’). We present a case study using two important game species in North America, the lesser and greater scaup, to demonstrate how INLA can be used to estimate the parameters in a hierarchical model that decouples observation error from process variation, and accounts for unknown sources of excess zeros as well as spatial and temporal dependence in the data. Ultimately, our goal was to make unbiased inference about spatial variation in population trends over time.

PLoS ONE↗