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 73 records · Page 4Linked to original sources

Mathematical model of gamma-ray spectrometry borehole logging for quantitative analysis

A technique for analyzing gamma-ray spectral-logging data has been developed, in which a digital computer is used to calculate the effects of gamma-ray attentuation in a borehole environment. The computer model allows for the calculation of the effects of lithology, porosity, density, and the thickness of a horizontal layer of uniformly distributed radioactive material surrounding a centralized probe in a cylindrical borehole. The computer program also contains parameters for the calculation of the effects of well casing, drilling fluid, probe housing, and losses through the sodium-iodide crystal. Errors associated with the commonly used mathematical assumption of a point detector are eliminated in this model. (USGS)

Open-File Report

The physical hydrogeology of ore deposits

Hydrothermal ore deposits represent a convergence of fluid flow, thermal energy, and solute flux that is hydrogeologically unusual. From the hydrogeologic perspective, hydrothermal ore deposition represents a complex coupled-flow problem—sufficiently complex that physically rigorous description of the coupled thermal (T), hydraulic (H), mechanical (M), and chemical (C) processes (THMC modeling) continues to challenge our computational ability. Though research into these coupled behaviors has found only a limited subset to be quantitatively tractable, it has yielded valuable insights into the workings of hydrothermal systems in a wide range of geologic environments including sedimentary, metamorphic, and magmatic. Examples of these insights include the quantification of likely driving mechanisms, rates and paths of fluid flow, ore-mineral precipitation mechanisms, longevity of hydrothermal systems, mechanisms by which hydrothermal fluids acquire their temperature and composition, and the controlling influence of permeability and other rock properties on hydrothermal fluid behavior. In this communication we review some of the fundamental theory needed to characterize the physical hydrogeology of hydrothermal systems and discuss how this theory has been applied in studies of Mississippi Valley-type, tabular uranium, porphyry, epithermal, and mid-ocean ridge ore-forming systems. A key limitation in the computational state-of-the-art is the inability to describe fluid flow and transport fully in the many ore systems that show evidence of repeated shear or tensional failure with associated dynamic variations in permeability. However, we discuss global-scale compilations that suggest some numerical constraints on both mean and dynamically enhanced crustal permeability. Principles of physical hydrogeology can be powerful tools for investigating hydrothermal ore formation and are becoming increasingly accessible with ongoing advances in modeling software.

Economic Geology

Geodetic imaging: Reservoir monitoring using satellite interferometry

Fluid fluxes within subsurface reservoirs give rise to surface displacements, particularly over periods of a year or more. Observations of such deformation provide a powerful tool for mapping fluid migration within the Earth, providing new insights into reservoir dynamics. In this paper we use Interferometric Synthetic Aperture Radar (InSAR) range changes to infer subsurface fluid volume strain at the Coso geothermal field. Furthermore, we conduct a complete model assessment, using an iterative approach to compute model parameter resolution and covariance matrices. The method is a generalization of a Lanczos-based technique which allows us to include fairly general regularization, such as roughness penalties. We find that we can resolve quite detailed lateral variations in volume strain both within the reservoir depth range (0.4-2.5 km) and below the geothermal production zone (2.5-5.0 km). The fractional volume change in all three layers of the model exceeds the estimated model parameter uncertainly by a factor of two or more. In the reservoir depth interval (0.4-2.5 km), the predominant volume change is associated with northerly and westerly oriented faults and their intersections. However, below the geothermal production zone proper [the depth range 2.5-5.0 km], there is the suggestion that both north- and northeast-trending faults may act as conduits for fluid flow.

Geophysical Journal International

Low-frequency radiation from point sources in a fluid-filled borehole

Far-field displacement fields have been derived for an impulsive point force acting on a fluid-filled borehole wall under the assumption that the borehole diameter is small compared to the wavelength involved. The displacements due to an arbitrary source can be computed easily by combining the solutions for the impulsive sources. In general, the borehole source generates not only longitudinal and vertically polarized shear waves, but also horizontally polarized shear waves. This study also indicates that only the axisymmetric motion around the borehole due to normal stress is affected by the fluid in the borehole. In the long-wavelength limit, the presence of the fluid does not affect the radiation from tangential sources into the surrounding medium.

Geophysics

Semi-implicit finite difference methods for three-dimensional shallow water flow

A semi-implicit finite difference method for the numerical solution of three-dimensional shallow water flows is presented and discussed. The governing equations are the primitive three-dimensional turbulent mean flow equations where the pressure distribution in the vertical has been assumed to be hydrostatic. In the method of solution a minimal degree of implicitness has been adopted in such a fashion that the resulting algorithm is stable and gives a maximal computational efficiency at a minimal computational cost. At each time step the numerical method requires the solution of one large linear system which can be formally decomposed into a set of small three-diagonal systems coupled with one five-diagonal system. All these linear systems are symmetric and positive definite. Thus the existence and uniquencess of the numerical solution are assured. When only one vertical layer is specified, this method reduces as a special case to a semi-implicit scheme for solving the corresponding two-dimensional shallow water equations. The resulting two- and three-dimensional algorithm has been shown to be fast, accurate and mass-conservative and can also be applied to simulate flooding and drying of tidal mud-flats in conjunction with three-dimensional flows. Furthermore, the resulting algorithm is fully vectorizable for an efficient implementation on modern vector computers.

International Journal for Numerical Methods in Flu

Computer analysis of digital well logs

A comprehensive system of computer programs has been developed by the U.S. Geological Survey for analyzing digital well logs. The programs are operational on a minicomputer in a research well-logging truck, making it possible to analyze and replot the logs while at the field site. The minicomputer also serves as a controller of digitizers, counters, and recorders during acquisition of well logs. The analytical programs are coordinated with the data acquisition programs in a flexible system that allows the operator to make changes quickly and easily in program variables such as calibration coefficients, measurement units, and plotting scales. The programs are designed to analyze the following well-logging measurements: natural gamma-ray, neutron-neutron, dual-detector density with caliper, magnetic susceptibility, single-point resistance, self potential, resistivity (normal and Wenner configurations), induced polarization, temperature, sonic delta-t, and sonic amplitude. The computer programs are designed to make basic corrections for depth displacements, tool response characteristics, hole diameter, and borehole fluid effects (when applicable). Corrected well-log measurements are output to magnetic tape or plotter with measurement units transformed to petrophysical and chemical units of interest, such as grade of uranium mineralization in percent eU3O8, neutron porosity index in percent, and sonic velocity in kilometers per second.

Circular

Performance of ice meter and weight assemblies

The performance of three ice meters and weight assemblies used by the U.S. Geological Survey and Environment Canada were compared in a towing tank. Each meter was rated individually on a rod suspension and then rerated on a cable suspension, with the appropriate weight assembly. Vertical and veer cable angles were measured along with meter yaw angle. The effect of the weight assembly on the rod-suspension rating for each meter was illustrated by computing a correction coefficient which ranged between 0.88 and 1.10 depending on the weight system used and the fluid velocity. A sluch-n-all type weight assembly least affected the meter rating and was the most stable in all flow conditions. (USGS)

Water-Resources Investigations Report

Overpressure and hydrocarbon accumulations in Tertiary strata, Gulf Coast of Louisiana

Many oil and gas reservoirs in Tertiary strata of southern Louisiana are located close to the interface between a sand-rich, normally pressured sequence and an underlying sand-poor, overpressured sequence. This association, recognized for many years by Gulf Coast explorationists, is revisited here because of its relevance to an assessment of undiscovered oil and gas potential in the Gulf Coast of Louisiana. The transition from normally pressured to highly overpressured sediments is documented by converting mud weights to pressure, plotting all pressure data from an individual field as a function of depth, and selecting a top and base of the pressure transition zone. Vertical extents of pressure transition zones in 34 fields across southern onshore Louisiana range from 300 to 9000 ft and are greatest in younger strata and in the larger fields. Display of pressure transition zones on geologic cross sections illustrates the relative independence of the depth of the pressure transition zone and geologic age. Comparison of the depth distribution of pressure transition zones with production intervals confirms previous findings that production intervals generally overlap the pressure transition zone in depth and that the median production depth lies above the base of the pressure transition zone in most fields. However, in 11 of 55 fields with deep drilling, substantial amounts of oil and gas have been produced from depths deeper than 2000 ft below the base of the pressure transition zone. Mud-weight data in 7 fields show that "local" pressure gradients range from 0.91 to 1.26 psi/ft below the base of the pressure transition zone. Pressure gradients are higher and computed effective stress gradients are negative in younger strata in coastal areas, indicating that a greater potential for fluid and sediment movement exists there than in older Tertiary strata.

Louisiana

Computer simulation of reservoir depletion and oil flow from the Macondo well following the Deepwater Horizon blowout

This report describes the application of a computer model to simulate reservoir depletion and oil flow from the Macondo well following the Deepwater Horizon blowout. Reservoir and fluid data used for model development are based on (1) information released in BP's investigation report of the incident, (2) information provided by BP personnel during meetings in Houston, Texas, and (3) calibration by history matching to shut-in pressures measured in the capping stack during the Well Integrity Test. The model is able to closely match the measured shut-in pressures. In the simulation of the 86-day period from the blowout to shut in, the simulated reservoir pressure at the well face declines from the initial reservoir pressure of 11,850 pounds per square inch (psi) to 9,400 psi. After shut in, the simulated reservoir pressure recovers to a final value of 10,300 psi. The pressure does not recover back to the initial pressure owing to reservoir depletion caused by 86 days of oil discharge. The simulated oil flow rate declines from 63,600 stock tank barrels per day just after the Deepwater Horizon blowout to 52,600 stock tank barrels per day just prior to shut in. The simulated total volume of oil discharged is 4.92 million stock tank barrels. The overall uncertainty in the simulated flow rates and total volume of oil discharged is estimated to be + or - 10 percent.

Open-File Report

A three-dimensional method-of-characteristics solute-transport model (MOC3D)

This report presents a model, MOC3D, that simulates three-dimensional solute transport in flowing ground water. The model computes changes in concentration of a single dissolved chemical constituent over time that are caused by advective transport, hydrodynamic dispersion (including both mechanical dispersion and diffusion), mixing (or dilution) from fluid sources, and mathematically simple chemical reactions (including linear sorption, which is represented by a retardation factor, and decay). The transport model is integrated with MODFLOW, a three-dimensional ground-water flow model that uses implicit finite-difference methods to solve the transient flow equation. MOC3D uses the method of characteristics to solve the transport equation on the basis of the hydraulic gradients computed with MODFLOW for a given time step. This implementation of the method of characteristics uses particle tracking to represent advective transport and explicit finite-difference methods to calculate the effects of other processes. However, the explicit procedure has several stability criteria that may limit the size of time increments for solving the transport equation; these are automatically determined by the program. For improved efficiency, the user can apply MOC3D to a subgrid of the primary MODFLOW grid that is used to solve the flow equation. However, the transport subgrid must have uniform grid spacing along rows and columns. The report includes a description of the theoretical basis of the model, a detailed description of input requirements and output options, and the results of model testing and evaluation. The model was evaluated for several problems for which exact analytical solutions are available and by benchmarking against other numerical codes for selected complex problems for which no exact solutions are available. These test results indicate that the model is very accurate for a wide range of conditions and yields minimal numerical dispersion for advection-dominated problems. Mass-balance errors are generally less than 10 percent, and tend to decrease and stabilize with time.

Water-Resources Investigations Report

Probabilistic 3-D time-lapse inversion of magnetotelluric data: Application to an enhanced geothermal system

Surface-based monitoring of mass transfer caused by injections and extractions in deep boreholes is crucial to maximize oil, gas and geothermal production. Inductive electromagnetic methods, such as magnetotellurics, are appealing for these applications due to their large penetration depths and sensitivity to changes in fluid conductivity and fracture connectivity. In this work, we propose a 3-D Markov chain Monte Carlo inversion of time-lapse magnetotelluric data to image mass transfer following a saline fluid injection. The inversion estimates the posterior probability density function of the resulting plume, and thereby quantifies model uncertainty. To decrease computation times, we base the parametrization on a reduced Legendre moment decomposition of the plume. A synthetic test shows that our methodology is effective when the electrical resistivity structure prior to the injection is well known. The centre of mass and spread of the plume are well retrieved.We then apply our inversion strategy to an injection experiment in an enhanced geothermal system at Paralana, South Australia, and compare it to a 3-D deterministic time-lapse inversion. The latter retrieves resistivity changes that are more shallow than the actual injection interval, whereas the probabilistic inversion retrieves plumes that are located at the correct depths and oriented in a preferential north-south direction. To explain the time-lapse data, the inversion requires unrealistically large resistivity changes with respect to the base model. We suggest that this is partly explained by unaccounted subsurface heterogeneities in the base model from which time-lapse changes are inferred.

Geophysical Journal International

Rates of salt solution in the Permian basin

For safe, long-term storage of radioactive materials in salt beds, rates of solution of salt which might imperil such storage must be known. For solution to continue, fluid must move through the system. Major fluid discharge from the system is by surface streams. Using U.S. Geological Survey records of streamflow and chemical quality of water, computations of sodium chloride discharges of numerous subbasins have been made. A map has been constructed showing tons of sodium chloride discharged per day and tons of sodium chloride discharged per square mile per year. Also shown are locations of major saline springs and their rate of discharge of sodium chloride in tons per day.

Colorado, Kansas, New Mexico, Oklahoma, Texas

Using diurnal temperature signals to infer vertical groundwater-surface water exchange

Heat is a powerful tracer to quantify fluid exchange between surface water and groundwater. Temperature time series can be used to estimate pore water fluid flux, and techniques can be employed to extend these estimates to produce detailed plan-view flux maps. Key advantages of heat tracing include cost-effective sensors and ease of data collection and interpretation, without the need for expensive and time-consuming laboratory analyses or induced tracers. While the collection of temperature data in saturated sediments is relatively straightforward, several factors influence the reliability of flux estimates that are based on time series analysis (diurnal signals) of recorded temperatures. Sensor resolution and deployment are particularly important in obtaining robust flux estimates in upwelling conditions. Also, processing temperature time series data involves a sequence of complex steps, including filtering temperature signals, selection of appropriate thermal parameters, and selection of the optimal analytical solution for modeling. This review provides a synthesis of heat tracing using diurnal temperature oscillations, including details on optimal sensor selection and deployment, data processing, model parameterization, and an overview of computing tools available. Recent advances in diurnal temperature methods also provide the opportunity to determine local saturated thermal diffusivity, which can improve the accuracy of fluid flux modeling and sensor spacing, which is related to streambed scour and deposition. These parameters can also be used to determine the reliability of flux estimates from the use of heat as a tracer.

Groundwater

Groundwater

Introduction Groundwater represents the terrestrial subsurface component of the hydrologic cycle. As such, groundwater is generally in motion, moving from elevated areas of recharge to lower areas of discharge. Groundwater usually moves in accordance with Darcy’s law (Dalmont, Paris: Les Fontaines Publiques de la Ville de Dijon, 1856). Groundwater residence times can be under a day in small upland catchments to over a million years in subcontinental-sized desert basins. The broadest definition of groundwater includes water in the unsaturated zone, considered briefly here. Water chemically bound to minerals, as in gypsum (CaSO 4 • 2H 2 O) or hydrated clays, cannot flow in response to gradients in total hydraulic head (pressure head plus elevation head); such water is thus usually excluded from consideration as groundwater. In 1940, M. King Hubbert showed Darcy’s law to be a special case of thermodynamically based potential field equations governing fluid motion, thereby establishing groundwater hydraulics as a rigorous engineering science ( Journal of Geology 48, pp. 785–944). The development of computer-enabled numerical methods for solving the field equations with real-world approximating geometries and boundary conditions in the mid-1960s ushered in the era of digital groundwater modeling. An estimated 30 percent of global fresh water is groundwater, compared to 0.3 percent that is surface water, 0.04 percent atmospheric water, and 70 percent that exists as ice, including permafrost ( Shiklomanov and Rodda 2004 , cited under Groundwater Occurrence ). Groundwater thus constitutes the vast majority—over 98 percent—of the unfrozen fresh-water resources of the planet, excluding surface-water reservoirs. Environmental dimensions of groundwater are equally large, receiving attention on multiple disciplinary fronts. Riparian, streambed, and spring-pool habitats can be sensitively dependent on the amount and quality of groundwater inputs that modulate temperature and solutes, including nutrients and dissolved oxygen. Groundwater withdrawals can negatively impact riparian habitats by depriving ecosystems of adequate fresh water and fragmenting communities when streams go dry. Biochemical reactions in shallow groundwater can remove anthropogenically elevated nitrogen compounds and reduce—but only to a point—the greening of waterways and shorelines with periphyton and harmful algal blooms. Groundwater extraction for beneficial use is increasingly limited by water-quality constraints imposed by naturally occurring and introduced substances. Overdrafting can cause land-surface subsidence, damaging buildings and roads and disrupting canals, sewers, and other gravity-flow conveyances. Increases in groundwater levels can cause soil salinization in dry regions and erosive sapping and flooding in wet regions. Coastal saltwater intrusion, groundwater flooding, salinization associated with groundwater-irrigated agriculture, induced seismicity from injected wastes, and the detrimental impacts of groundwater depletion are among the major environmental challenges of our time.

Book chapter

Guide to the Revised Ground-Water Flow and Heat Transport Simulator: HYDROTHERM - Version 3

The HYDROTHERM computer program simulates multi-phase ground-water flow and associated thermal energy transport in three dimensions. It can handle high fluid pressures, up to 1 ? 109 pascals (104 atmospheres), and high temperatures, up to 1,200 degrees Celsius. This report documents the release of Version 3, which includes various additions, modifications, and corrections that have been made to the original simulator. Primary changes to the simulator include: (1) the ability to simulate unconfined ground-water flow, (2) a precipitation-recharge boundary condition, (3) a seepage-surface boundary condition at the land surface, (4) the removal of the limitation that a specified-pressure boundary also have a specified temperature, (5) a new iterative solver for the linear equations based on a generalized minimum-residual method, (6) the ability to use time- or depth-dependent functions for permeability, (7) the conversion of the program code to Fortran 90 to employ dynamic allocation of arrays, and (8) the incorporation of a graphical user interface (GUI) for input and output. The graphical user interface has been developed for defining a simulation, running the HYDROTHERM simulator interactively, and displaying the results. The combination of the graphical user interface and the HYDROTHERM simulator forms the HYDROTHERM INTERACTIVE (HTI) program. HTI can be used for two-dimensional simulations only. New features in Version 3 of the HYDROTHERM simulator have been verified using four test problems. Three problems come from the published literature and one problem was simulated by another partially saturated flow and thermal transport simulator. The test problems include: transient partially saturated vertical infiltration, transient one-dimensional horizontal infiltration, two-dimensional steady-state drainage with a seepage surface, and two-dimensional drainage with coupled heat transport. An example application to a hypothetical stratovolcano system with unconfined ground-water flow is presented in detail. It illustrates the use of HTI with the combination precipitation-recharge and seepage-surface boundary condition, and functions as a tutorial example problem for the new user.

Techniques and Methods

Assessment of the geoavailability of trace elements from minerals in mine wastes: analytical techniques and assessment of selected copper minerals

In this study, four randomly selected copper-bearing minerals were examined—azurite, malachite, bornite, and chalcopyrite. The objectives were to examine and enumerate the crystalline and chemical properties of each of the minerals, to determine which, if any, of the Cu-bearing minerals might adversely affect systems biota, and to provide a multi-procedure reference. Laboratory work included use of computational software for quantifying crystalline and amorphous material and optical and electron imaging instruments to model and project crystalline structures. Chemical weathering, human fluid, and enzyme simulation studies were also conducted. The analyses were conducted systematically: X-ray diffraction and microanalytical studies followed by a series of chemical, bio-leaching, and toxicity experiments.

Scientific Investigations Report

A compilation of rate parameters of water-mineral interaction kinetics for application to geochemical modeling

Geochemical reaction path modeling is useful for rapidly assessing the extent of water-aqueous-gas interactions both in natural systems and in industrial processes. Modeling of some systems, such as those at low temperature with relatively high hydrologic flow rates, or those perturbed by the subsurface injection of industrial waste such as CO2 or H2S, must account for the relatively slow kinetics of mineral-gas-water interactions. We have therefore compiled parameters conforming to a general Arrhenius-type rate equation, for over 70 minerals, including phases from all the major classes of silicates, most carbonates, and many other non-silicates. The compiled dissolution rate constants range from -0.21 log moles m-2 s-1 for halite, to -17.44 log moles m-2 s-1 for kyanite, for conditions far from equilibrium, at 25 ?C, and pH near neutral. These data have been added to a computer code that simulates an infinitely well-stirred batch reactor, allowing computation of mass transfer as a function of time. Actual equilibration rates are expected to be much slower than those predicted by the selected computer code, primarily because actual geochemical processes commonly involve flow through porous or fractured media, wherein the development of concentration gradients in the aqueous phase near mineral surfaces, which results in decreased absolute chemical affinity and slower reaction rates. Further differences between observed and computed reaction rates may occur because of variables beyond the scope of most geochemical simulators, such as variation in grain size, aquifer heterogeneity, preferred fluid flow paths, primary and secondary mineral coatings, and secondary minerals that may lead to decreased porosity and clogged pore throats.

Open-File Report

Modeling subsurface performance of a geothermal reservoir using machine learning

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells—increasing or decreasing the fluid flow rates across the wells—and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. In this study, we describe a new approach combining reservoir modeling and machine learning to produce models that enable such a strategy. Our computational approach allows us, first, to translate sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy, and second, to find optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an “open-source” reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 h, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 s. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs.

Energies