Search USGSSearch

SEARCH · Search USGS

Results for “Computers and Fluids”

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 55 records · Page 3Linked to original sources

Euler-Lagrangian computation for estuarine hydrodynamics

The transport of conservative and suspended matter in fluid flows is a phenomenon of Lagrangian nature because the process is usually convection dominant. Nearly all numerical investigations of such problems use an Eulerian formulation for the convenience that the computational grids are fixed in space and because the vast majority of field data are collected in an Eulerian reference frame. Several examples are given in this paper to illustrate a modeling approach which combines the advantages of both the Eulerian and Lagrangian computational techniques.

Conference Paper

A finite-element simulation model for saturated-unsaturated, fluid-density-dependent ground-water flow with energy transport or chemically- reactive single-species solute transport

SUTRA (Saturated-Unsaturated Transport) is a computer program which can be used to simulate the movement of fluid and the transport of either energy or dissolved substances in a subsurface environment. The model employs a two-dimensional hybrid finite-element and integrated-finite-difference method to approximate the governing equations that describe the two interdependent processes that are simulated by SUTRA: (1) fluid density-dependent saturated or unsaturated groundwater flow, and either (2a) transport of a solute in the groundwater, in which the solute may be subject to: equilibrium adsorption on the porous matrix, and both first-order and zero-order production or decay, or, (2b) transport of thermal energy in the groundwater and solid matrix of the aquifer. SUTRA provides, as the primary calculated results, fluid pressures and either solute concentrations or temperatures, as they vary with time, everywhere in the simulated subsurface system. SUTRA may also be used to simulate simpler subsets of the above process. SUTRA may be employed for areal and cross-sectional models of saturated groundwater flow systems, and for cross-sectional models of unsaturated zone flow. Solute transport simulation using SUTRA may be used to simulate natural or man-induced chemical transport, solute sorption, production and decay. SUTRA may be used for simulation of variable density leachate movement, and for cross-sectional simulation of salt-water intrusion in aquifers at near-well or regional scales, with either dispersed or relatively sharp transition zones between fresh water and salt water. SUTRA energy transport simulation may be employed to model thermal regimes in aquifers, subsurface heat conduction, aquifer thermal energy storage systems, geothermal reservoirs, thermal pollution of aquifers, and natural hydrogeologic convection systems. (USGS)

Water-Resources Investigations Report

Hawaiian fissure fountains: Quantifying vent and shallow conduit geometry, episode 1 of the 1969-1974 Mauna Ulu eruption

Geometries of shallow magmatic pathways feeding volcanic eruptions are poorly constrained, yet many key interpretations about eruption dynamics depend on knowledge of these geometries. Direct quantification is difficult because vents typically become blocked with lava at the end of eruptions. Indirect geophysical techniques have shed light on some volcanic conduit geometries, but the scales are too coarse to resolve narrow fissures (widths typically 1 m). Kīlauea's Mauna Ulu eruption, which started with <50 m high Hawaiian fountains along a 4.5 km fissure on 24 May 1969, provides a unique opportunity to measure the detailed geometry of a shallow magmatic pathway, as the western vents remain unobstructed to depths >30 m. Direct measurements at the ground surface were augmented by tripod-mounted lidar measurements to quantify the shallow conduit geometry for three vents at a resolution <4 cm. We define the form of the fissure in terms of aspect ratio, flaring ratio, irregularity, sinuosity, and segmentation and discuss the factors influencing these parameters. In the past, simplified first-order fissure geometries have been used in computational modeling. Our data can provide more accurate conduit shapes for better understanding of shallow fissure fluid dynamics and how it controls eruptive behavior, especially if incorporated into computer models.

Hawaii

Comment on “The reduction of friction in long-runout landslides as an emergent phenomenon” by Brandon C. Johnson et al.

Results from a highly idealized, 2-D computational model indicate that dynamic normal-stress rarefactions might cause friction reduction in long-runout landslides, but the physical relevance of the idealized dynamics has not been confirmed by experimental tests. More importantly, the model results provide no evidence that refutes alternative hypotheses about friction reduction mechanisms. One alternative hypothesis, which is strongly supported by field evidence, experimental data, and the predictions of a well-constrained computational model, involves development of high pore fluid pressures in deforming landslide material or overridden bed material. However, no scientific basis exists for concluding that a universal mechanism is responsible for friction reduction in all long-runout landslides.

Journal of Geophysical Research F: Earth Surface

Hydrologic and geochemical monitoring in Long Valley Caldera, Mono County, California, 1985

Hydrologic and geochemical monitoring, to detect changes caused by magmatic and tectonic processes in the Long Valley caldera has continued through 1985. The monitoring included the collection of the following types of data: chemical and isotopic composition of water and gases from springs, wells, and steam vents; temperatures in wells, springs, and steam vents; flow rates of springs and streams; water levels in wells; and barometric pressure and precipitation at several sites. In addition, reservoir temperatures for the geothermal system were estimated from computations based on chemical geothermometers applied to fluid samples from wells and springs. Estimates of thermal water discharged from springs were made on the basis of boron and chloride fluxes in surface waters for selected sites in the Casa Diablo area and along the Mammoth-Hot Creek drainage. These data are presented in tables and graphs. The Long Valley area was relatively quiescent throughout 1985 in terms of geodetic changes and seismic activity. As a consequence , the hydrologic system varied mainly in response to seasonal influences of temperature, atmospheric pressure, and precipitation. However, spring flows near Casa Diablo were influenced by pumping at the geothermal production well field nearby. (Author 's abstract)

Water-Resources Investigations Report

Modeling transport in transient ground-water flow: An unacknowledged approximation

During unsteady or transient ground-water flow, the fluid mass per unit volume of aquifer changes as the potentiometric head changes, and solute transport is affected by this change in fluid storage. Three widely applied numerical models of two-dimensional transport partially account for the effects of transient flow by removing terms corresponding to the fluid continuity equation from the transport equation, resulting in a simpler governing equation. However, fluid-storage terms remaining in the transport equation that change during transient flow are, in certain cases, held constant in time in these models. For the case of increasing heads, this approximation, which is unacknowledged in these models' documentation, leads to transport velocities that are too high, and increased concentration at fluid and solute sources. If heads are dropping in time, computed transport velocities are too low. Using parameters that somewhat exaggerate the effects of this approximation, an example numerical simulation indicates solute travel time error of about 14 percent but only minor errors due to incorrect dilution volume. For horizontal flow and transport models that assume fluid density is constant, the product of porosity and aquifer thickness changes in time: initial porosity times initial thickness plus the change in head times the storage coefficient. This formula reduces to the saturated thickness in unconfined aquifers if porosity is assumed to be constant and equal to specific yield. The computational cost of this more accurate representation is insignificant and is easily incorporated in numerical models of solute transport.

Ground Water

2022 Crustal Deformation Modeling Workshop Report

The 2022 Crustal Deformation Modeling Workshop was held June 20–24 at the Colorado School of Mines in Golden, Colorado. The workshop included two days of tutorials on the use of the open-source software PyLith for crustal deformation modeling followed by three days of science talks and discussions. The workshop focused on three primary themes: (1) Earthquake cycle modeling; (2) Inversions for fault slip; and (3) Faulting, fluids, and surface loading. The talks highlighted how computational models provide insight into intriguing observations of Earth and planetary behavior. These included (1) earthquake synchronization of rupture patches due to their close proximity to each other, (2) the influence of fault geometry and damage zones on hypocenter depth and rupture propagation, (3) a lack of steady-state faulting behavior due to long time scales for grain size evolution in the mid-crust, (4) crustal deformation due to tidal, hydrological, and atmospheric loads, and (5) plumes of gas and icy particles due to tidal driven faulting on Enceladus (one of Saturn’s moons). The talks also described new computational modeling capabilities for incorporating complex geologic structure into Bayesian inversions for fault slip and efficient implementation of earthquake cycle models using a symmetric interior discontinuous Galerkin method. The complete agenda is available on the Computational Infrastructure for Geodynamic (CIG) website.

Report

Mudflow rheology in a vertically rotating flume

Joint research by the U.S. Geological Survey and the University of Missouri-Rolla currently (1990) is being conducted on a 3.05 meters in diameter vertically rotating flume used to simulate mudflows under steady-state conditions. Observed mudflow simulations indicate flow patterns in the flume are similar to those occurring in natural mudflows. Variables such as mean and surface velocity, depth, and average boundary shear stress can be measured in this flume more easily than in the field or in a traditional tilting flume. Sensitive variables such as sediment concentration, grain-size distribution, and Atterberg limits also can be precisely and easily controlled. A known Newtonian fluid, SAE 30 motor oil, was tested in the flume and the computed value for viscosity was within 12.5 percent of the stated viscosity. This provided support that the data from the flume can be used to determine the rheological properties of fluids such as mud. Measurements on mud slurries indicate that flows with sediment concentrations ranging from 81 to 87 percent sediment by weight can be approximated as Bingham plastic for strain rates greater than 1 per second. In this approximation, the yield stress and Bingham viscosity were extremely sensitive to sediment concentration. Generally, the magnitude of the yield stress was large relative to the change in shear stress with increasing mudflow velocity.

Conference Paper

Environmental contaminants and biomarker responses in fish from the Rio Grande and its U.S. tributaries: Spatial and temporal trends

We collected, examined, and analyzed 368 fish of seven species from 10 sites on rivers of the Rio Grande Basin (RGB) during late 1997 and early 1998 to document temporal and geographic trends in the concentrations of accumulative contaminants and to assess contaminant effects on the fish. Sites were located on the mainstem of the Rio Grande and on the Arroyo Colorado and Pecos River in Texas (TX), New Mexico (NM), and Colorado. Common carp (Cyprinus carpio) and largemouth bass (Micropterus salmoides) were the targeted species. Fish were examined in the field for internal and external visible gross lesions, selected organs were weighed to compute ponderal and organosomatic indices, and samples of tissues and fluids were obtained and preserved for analysis of fish health and reproductive biomarkers. Whole fish from each station were composited by species and gender and analyzed for organochlorine chemical residues and elemental contaminants using instrumental methods, and for 2,3,7,8-tetrachloro dibenzo-p-dioxin-like activity (TCDD-EQ) using the H4IIE rat hepatoma cell bioassay. Overall, fish from lower RGB stations contained greater concentrations of organochlorine pesticide residues and appeared to be less healthy than those from sites in the central and upper parts of the basin, as indicated by a general gradient of residue concentrations and biomarker responses. A minimal number of altered biomarkers and few or no elevated contaminant concentrations were noted in fish from the upper RGB. The exception was elevated concentrations [up to 0.46 ??g/g wet-weight (ww)] of total mercury (Hg) in predatory species from the Rio Grande at Elephant Butte Reservoir, NM, a condition documented in previous studies. Arsenic (As) and selenium (Se) concentrations were greatest in fish from sites in the central RGB; Se concentrations in fish from the Pecos River at Red Bluff Lake, TX and from the Rio Grande at Langtry, TX and Amistad International Reservoir, TX exceeded published fish and wildlife toxicity thresholds. In the lower RGB, residues of p,p???-DDT metabolites (???1.69 ??g/g ww), chlordane-related compounds (???0.21 ??g/g ww), dieldrin (???0.0.05 ??g/g ww), and toxaphene (???2.4 ??g/g ww) were detected in fish from most sites; maximum concentrations were in channel catfish (Ictalurus punctatus) from the Arroyo Colorado at Harlingen, TX. Concentrations of one or more residues exceeded toxicity thresholds for fish and wildlife in fish from this site and from the Rio Grande at Mission, TX and Brownsville, TX; however, concentrations were lower than those reported by previous studies. In addition, the proportional concentrations of p,p???-DDT at all sites were low, indicating weathered DDT rather than the influx of new material. Concentrations of total PCBs (<0.05 ??g/g ww) and TCDD-EQ (???6 pg/g ww) were comparatively low in all samples. Hepatic ethoxyresorufin O-deethylase (EROD) activity in some fish was elevated relative to reference rates at most sites, but was generally lower than previously reported activity in fish from heavily contaminated locations. The comparatively low PCB and TCDD-EQ concentrations together with elevated EROD activity may reflect exposure to polycyclic aromatic hydrocarbons. Reproductive biomarkers were consistent with chronic contaminant exposure at lower RGB sites; comparatively large percentages of intersex male largemouth bass, relatively low gonadosomatic indices, and elevated plasma vitellogenin concentrations in male fish were noted at three of the four stations. Large percentages of atretic eggs were also observed in the ovaries of female common carp from the Rio Grande at Brownsville, TX. Although many of the conditions noted may have other causes in addition to contaminant exposure, the biomarker results for the lower RGB sites are consistent with subtle responses of fish to contaminants, an interpretation supported by the chemical data of this and other investigations.

Science of the Total Environment

Poroelastic properties of the Arbuckle Group in Oklahoma derived from well fluid level response to the 3 September 2016 M w 5.8 Pawnee and 7 November 2016 M w 5.0 Cushing earthquakes

The Arbuckle Group (Arbuckle) is a basal sedimentary unit that is the primary target for saltwater disposal in Oklahoma. Thus, the reservoir characteristics of the Arbuckle, including how the poroelastic properties change laterally and over time are of significant interest. We report observations of fluid level changes in two monitoring wells in response to the 3 September 2016 M w 5.8 Pawnee and the 7 November 2016 M w 5.0 Cushing earthquakes. We investigate the relationship between static strain resulting from these events and the fluid level changes observed in the wells. We model the fluid level response by estimating static strains from a set of earthquake source parameters and spatiotemporal poroelastic properties of the Arbuckle in the neighborhood of the monitoring wells. Results suggest that both the direction of the observed fluid level step and the amplitude can be predicted from the computed volumetric strain change and a reasonable set of poroelastic parameters. Modeling results indicate that poroelastic parameters differ at the time of the Pawnee and Cushing earthquakes, with a moderately higher Skempton’s coefficient required to fit the response to the Cushing earthquake. This may indicate that dynamic shaking resulted in physical alteration of the Arbuckle at distances up to ∼50 km from the Pawnee earthquake.

Oklahoma

Biomonitoring of Environmental Status and Trends (BEST) Program: Environmental contaminants and their effects on fish in the Rio Grande Basin

We collected, examined, and analyzed 368 fish of seven species from 10 sites in the Rio Grande Basin (RGB) during late 1997 and early 1998. Four sites were National Contaminant Biomonitoring Program (NCBP) stations where organochlorine ad elemental contaminants in fish had been monitored from 1969 and through 1986. The other six were USGS-National Stream Quantity Accounting Network (NASQAN) stations where water quality is monitored. The objectives were to document temporal and geographic trends in the concentrations of accumulative organic and inorganic contaminants in RGB fish and the effects of contaminants on the fish; to continue testing the feasibility of incorporating biomarkers (that is, biochemical, histopathological, and other biological indicators of contaminant exposure or effects) into a monitoring program for large U.S. rivers; and to evaluate the compatibility of monitoring methods based on the analysis of fish with those used to monitor water by NASQAN. Common carp ( Cyprinus carpio ; carp) and black basses ( Micropterus sp.; bass) were the targeted species; together, they represented 77% of the fish collected. Each fish was examined in the field for externally and internally visible gross lesions, selected organs were weighed to compute various ponderal and organosomatic indices, and samples of tissues and fluids were obtained and preserved for analysis of fish health and reproductive biomarkers. Composite samples of whole fish from each station were grouped by species and gender and analyzed by instrumental methods for persistent organic and inorganic contaminants and for dioxin-like activity (TCDD-EQ) using the H4IIE rat hepatoma cell bioassay.

Rio Grande basin

Borehole geophysics applied to ground-water investigations

The purpose of this manual is to provide hydrologists, geologists, and others who have the necessary training with the basic information needed to apply the most useful borehole-geophysical-logging techniques to the solution of problems in ground-water hydrology. Geophysical logs can provide information on the construction of wells and on the character of the rocks and fluids penetrated by those wells, in addition to changes in the character of these factors with time. The response of well logs is caused by: petrophysical factors; the quality; temperature, and pressure of interstitial fluids; and ground-water flow. Qualitative and quantitative analysis of the analog records and computer analysis of digitized logs are used to derive geohydrologic information. This information can then be extrapolated vertically within a well and laterally to other wells using logs. The physical principles by which the mechanical and electronic components of a logging system measure properties of rocks, fluids and wells, and the principles of measurement need to be understood to correctly interpret geophysical logs. Planning the logging operation involves selecting the equipment and the logs most likely to provide the needed information. Information on well construction and geohydrology are needed to guide this selection. Quality control of logs is an important responsibility of both the equipment operator and log analyst and requires both calibration and well-site standardization of equipment. Logging techniques that are widely used in ground-water hydrology or that have significant potential for application to this field include: spontaneous potential, resistance, resistivity, gamma, gamma spectrometry, gamma-gamma, neutron, acoustic velocity, acoustic televiewer, caliper, and fluid temperature, conductivity, and flow. The following topics are discussed for each of these techniques: principles and instrumentation, calibration and standardization, volume of investigation, extraneous effects, and interpretation and applications.

Open-File Report

Borehole geophysics applied to ground-water investigations

The purpose of this manual is to provide hydrologists, geologists, and others who have the necessary background in hydrogeology with the basic information needed to apply the most useful borehole-geophysical-logging techniques to the solution of problems in ground-water hydrology. Geophysical logs can provide information on the construction of wells and on the character of the rocks and fluids penetrated by those wells, as well as on changes in the character of these factors over time. The response of well logs is caused by petrophysical factors, by the quality, temperature, and pressure of interstitial fluids, and by ground-water flow. Qualitative and quantitative analysis of analog records and computer analysis of digitized logs are used to derive geohydrologic information. This information can then be extrapolated vertically within a well and laterally to other wells using logs. The physical principles by which the mechanical and electronic components of a logging system measure properties of rocks, fluids, and wells, as well as the principles of measurement, must be understood if geophysical logs are to be interpreted correctly. Plating a logging operation involves selecting the equipment and the logs most likely to provide the needed information. Information on well construction and geohydrology is needed to guide this selection. Quality control of logs is an important responsibility of both the equipment operator and the log analyst and requires both calibration and well-site standardization of equipment. Logging techniques that are widely used in ground-water hydrology or that have significant potential for application to this field include spontaneous potential, resistance, resistivity, gamma, gamma spectrometry, gamma-gamma, neutron, acoustic velocity, acoustic televiewer, caliper, and fluid temperature, conductivity, and flow. The following topics are discussed for each of these techniques: principles and instrumentation, calibration and standardization, volume of investigation, extraneous effects, and interpretation and applications.

Techniques of Water-Resources Investigations

Hydrogeophysical methods for analyzing aquifer storage and recovery systems

Hydrogeophysical methods are presented that support the siting and monitoring of aquifer storage and recovery (ASR) systems. These methods are presented as numerical simulations in the context of a proposed ASR experiment in Kuwait, although the techniques are applicable to numerous ASR projects. Bulk geophysical properties are calculated directly from ASR flow and solute transport simulations using standard petrophysical relationships and are used to simulate the dynamic geophysical response to ASR. This strategy provides a quantitative framework for determining site‐specific geophysical methods and data acquisition geometries that can provide the most useful information about the ASR implementation. An axisymmetric, coupled fluid flow and solute transport model simulates injection, storage, and withdrawal of fresh water (salinity ∼500 ppm) into the Dammam aquifer, a tertiary carbonate formation with native salinity approximately 6000 ppm. Sensitivity of the flow simulations to the correlation length of aquifer heterogeneity, aquifer dispersivity, and hydraulic permeability of the confining layer are investigated. The geophysical response using electrical resistivity, time‐domain electromagnetic (TEM), and seismic methods is computed at regular intervals during the ASR simulation to investigate the sensitivity of these different techniques to changes in subsurface properties. For the electrical and electromagnetic methods, fluid electric conductivity is derived from the modeled salinity and is combined with an assumed porosity model to compute a bulk electrical resistivity structure. The seismic response is computed from the porosity model and changes in effective stress due to fluid pressure variations during injection/recovery, while changes in fluid properties are introduced through Gassmann fluid substitution.

Groundwater

Biomonitoring of Environmental Status and Trends (BEST) Program: Environmental contaminants and their effects on fish in the Mississippi River Basin

We collected, examined, and analyzed 1378 fish of 22 species from 47 sites in the Mississippi River basin (MRB) during 1995 and from a reference site in 1996. The sampling sites in the MRB represented National Contaminant Biomonitoring Program (NCBP) stations situated at key points on major rivers and National Water- Quality Assessment Program (NAWQA) stations located on lower-order rivers and streams in the Eastern Iowa Basins (EIB) and Mississippi Embayment (MSE) Study Units. The reference site was the water supply system of the USGS-Leetown Science Center in rural Jefferson County, WV. Common carp (Cyprinus carpio; carp) and black basses (Micropterus spp.; bass), the targeted species, together represented 82% of the fish collected. Each fish was examined in the field for externally and internally visible gross lesions, selected organs were weighed to compute various ponderal and organo-somatic indices, and selected tissues and fluids were obtained and preserved for analysis of biomarkers. Fish health indicators included splenic macrophage aggregates, lysozyme activity, and hispathological analysis of liver, kidney, and other tissues. Reproductive biomarkers included analysis of plasma concentrations of vitellogenin (vtg) and the sex steroid hormones 17-estradiol (E2) and 11-ketotestosterone (11- kt); and the histological determination of percent oocyte atresia (in female fish) and gonadal stage. Hepatic ethoxyresorufin O-deethylase (EROD) activity was also measured. Composite samples of whole fish from each station were grouped by species and gender and analyzed for persistent organochlorine and elemental contaminants and for dioxin-like activity (TCDD-EQ) using the H4IIE rat hepatoma cell bioassay. Organochlorine and inorganic contaminant concentrations in fish were generally low relative to historical levels at most sites, but remained present at concentrations representing threats to piscivorous wildlife in some locations. Toxaphene and DDT (mostly as p,p?-DDE) concentrations remained elevated in fish from the cottongrowing regions of the lower Mississippi valley, and were generally greater in the smaller streams draining agricultural areas (that is, in the MSE Study Unit) than at large river sites. Cyclodiene pesticide concentrations were also greatest in the EIB Study Unit and elsewhere in the corn-growing regions of the mid-MRB. Former point-sources of organochlorine pesticides also remained evident, especially in the Mississippi River near Memphis, TN. Consistent with previous findings, total PCB concentrations tended to be greatest (1-3 g/g) in the industrialized and urbanized Ohio River and Upper Mississippi sub-basins and at Memphis, TN, and were generally correlated with TCDD-EQ and EROD activity. Conversely, PCB concentrations were low (<0.1 g/g) in the more agricultural parts of the MRB. Concentrations of inorganic contaminants were also relatively low and stable or declining relative to past levels at most sites. Exceptions were Hg and Se; Hg concentrations were slightly elevated (>0.3 g/g) in bass from the Mississippi River at Memphis and several other sites and in carp from one MSE site. Concentrations of Se were also great enough to constitute a hazard to piscivorous wildlife (>0.6 g/g) at several MRB sites in the western parts of the MRB and were especially high (4-5 g/g) in fish from John Martin Reservoir, CO, where elevated concentrations were reported previously. Biomarker results indicated that fish from many stations had been exposed to contaminants, but at no sites did findings indicate exposure to high concentrations of toxic chemicals. Noteworthy among biomarker findings was that 73% of the male smallmouth bass (Micropterus dolomieui) from the Mississippi River at Lake City, MN (Lake Pepin) were intersex as indicated by the histological detection of ovotestes; and the combined EROD and H4IIE results indicated that fish from several rural sites in the

Mississippi River basin

Machine learning provides reconnaissance-type estimates of carbon dioxide storage resources in oil and gas reservoirs

Oil and gas reservoirs represent suitable containers to sequester carbon dioxide (CO 2 ) in a supercritical state because they are accessible, reservoir properties are known, and they previously contained stored buoyant fluids. However, planners must quantify the relative magnitude of the CO 2 storage resource in these reservoirs to formulate a comprehensive strategy for CO 2 mitigation. Even reconnaissance-type estimates of CO 2 storage resources of known oil and gas reservoirs may require complicated calculations involving 1) estimates of recoverable oil and gas, 2) reservoir properties (depth, temperature, pressure, etc.), and 3) the physical qualities of the retained fluids. We demonstrate the application of machine learning (ML) algorithms to bypass these computations to yield more rapid estimates of CO 2 storage resources in reservoirs capable of hosting CO 2 in a supercritical state. ML algorithms are computationally efficient because they do not impose the strong assumptions on the data-generating process that standard statistical or engineering procedures require. Further, ML algorithms can capture highly complex, particularly nonlinear, relationships among predictor variables. We demonstrate the application of four different ML algorithms using data from onshore and offshore oil and gas reservoirs in Europe, and show they perform well when predictions are compared to engineering estimates. The proposed methods and models provide an effective and novel way to more rapidly and directly determine the subsurface CO 2 storage capacity of oil and gas reservoirs around the world, information that operators, researchers, and policymakers alike require to meet energy transition and decarbonization goals.

Frontiers in Enviornmental Science

Nonlinear inversion of borehole-radar tomography data to reconstruct velocity and attenuation distribution in earth materials

A nonlinear tomographic inversion method that uses first-arrival travel-time and amplitude-spectra information from cross-hole radar measurements was developed to simultaneously reconstruct electromagnetic velocity and attenuation distribution in earth materials. Inversion methods were developed to analyze single cross-hole tomography surveys and differential tomography surveys. Assuming the earth behaves as a linear system, the inversion methods do not require estimation of source radiation pattern, receiver coupling, or geometrical spreading. The data analysis and tomographic inversion algorithm were applied to synthetic test data and to cross-hole radar field data provided by the US Geological Survey (USGS). The cross-hole radar field data were acquired at the USGS fractured-rock field research site at Mirror Lake near Thornton, New Hampshire, before and after injection of a saline tracer, to monitor the transport of electrically conductive fluids in the image plane. Results from the synthetic data test demonstrate the algorithm computational efficiency and indicate that the method robustly can reconstruct electromagnetic (EM) wave velocity and attenuation distribution in earth materials. The field test results outline zones of velocity and attenuation anomalies consistent with the finding of previous investigators; however, the tomograms appear to be quite smooth. Further work is needed to effectively find the optimal smoothness criterion in applying the Tikhonov regularization in the nonlinear inversion algorithms for cross-hole radar tomography.

New Hampshire

Profiles of reservoir properties of oil-bearing plays for selected petroleum provinces in the United States

Profiles of reservoir properties of oil-bearing plays for selected petroleum provinces in the United States were developed to characterize the database to be used for a potential assessment by the U.S. Geological Survey (USGS) of oil that would be technically recoverable by the application of enhanced oil recovery methods using injection of carbon dioxide (CO 2 -EOR). The USGS assessment methodology may require reservoir-level data for the purposes of screening conventional oil reservoirs and projecting CO 2 -EOR performance in terms of the incremental recoverable oil. The information used in this report is based on reservoir properties from the &ldquo;Significant Oil and Gas Fields of the United States Database&rdquo; prepared by Nehring Associates, Inc. (2012). As described by Nehring Associates, Inc., the database &ldquo;covers all producing provinces (basins) in the United States except the Appalachian Basin and the Cincinnati Arch.&rdquo; Under contract to the USGS, INTEK, Inc., developed and applied algorithms to estimate variables useful in projecting EOR performance at the reservoir level and to complete some partial reservoir records of the &ldquo;Significant Oil and Gas Fields of the United States Database&rdquo; (Nehring Associates, Inc., 2012). The augmented database is referred to here as the &ldquo;Comprehensive Resource Database&rdquo; (CRD). The CRD play and province classification scheme corresponds to the definitions used in the 1995 USGS National Oil and Gas Assessment (NOGA). The profiles in this report consist of a resource table and a six-part figure showing the variation of reservoir parameters selected because of their importance in the choice of a miscible or immiscible method for CO 2 -EOR and in the assessment of potential oil recovery using the EOR processes. A subset of these reservoirs may be available for either miscible- or immiscible-type flooding for CO 2 -EOR. Plays with fewer than 10 oil reservoirs were not graphed and were omitted from the province profiles. For this report and for the purposes of screening reservoirs as candidates for the application of CO 2 -EOR methods, oil reservoirs must have no more than 10,000 standard cubic feet of natural gas per barrel of oil at surface conditions. Oil-bearing plays presented in this report must contain at least one oil reservoir so defined. The profile plots allow geologists to evaluate the range of empirical and default values of the oil reservoir characteristics within a play and across plays that belong to the same province in the CRD. For most plays, the default estimates can be identified by the stacking of points at a single value on strip charts in the profiles. Reasonable default values should be within the range of the reservoir parameter values assigned by Nehring Associates, Inc. (2012), to reservoirs of that particular play. Each province profile figure consists of five strip charts and a boxplot. The five strip charts display for individual plays the following reservoir-fluid and reservoir properties: A , oil density (American Petroleum Institute [API] gravity in degrees); B , computed pseudo-Dykstra-Parsons coefficient; C , reservoir porosity (in percent); D , reservoir permeability (in millidarcies); and E , estimates of the original oil in place (OOIP) per unit volume of reservoir rock (in barrels per acre-foot). The OOIP per unit volume of reservoir rock is an indicator of the relative richness of the oil reservoir and is derived from estimates in the CRD of OOIP, reservoir acreage, and net pay. The net pay is the interval of productive reservoir rock. The same data for OOIP per unit volume are graphed as a strip chart ( E ) and a boxplot ( F ).

Open-File Report