Search USGS⌕ Search

SEARCH · Search USGS

Results for “Computers & Geosciences”

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
USGS1993

Volumetric calculations in an oil field: The basis method

The basis method for estimating oil reserves in place is compared to a traditional procedure that uses ordinary kriging. In the basis method, auxiliary variables that sum to the net thickness of pay are estimated by cokriging. In theory, the procedure should be more powerful because it makes full use of the cross-correlation between variables and forces the original variables to honor interval constraints. However, at least in our case study, the practical advantages of cokriging for estimating oil in place are marginal. ?? 1993.

Computers & Geosciences↗
USGS1993

A fortran program for Monte Carlo simulation of oil-field discovery sequences

We have developed a program for performing Monte Carlo simulation of oil-field discovery histories. A synthetic parent population of fields is generated as a finite sample from a distribution of specified form. The discovery sequence then is simulated by sampling without replacement from this parent population in accordance with a probabilistic discovery process model. The program computes a chi-squared deviation between synthetic and actual discovery sequences as a function of the parameters of the discovery process model, the number of fields in the parent population, and the distributional parameters of the parent population. The program employs the three-parameter log gamma model for the distribution of field sizes and employs a two-parameter discovery process model, allowing the simulation of a wide range of scenarios. ?? 1993.

Computers & Geosciences↗
USGS1994

PETRO.CALC.PLOT, Microsoft Excel macros to aid petrologic interpretation

PETRO.CALC.PLOT is a package of macros which normalizes whole-rock oxide data to 100%, calculates the cation percentages and molecular proportions used for normative mineral calculations, computes the apices for ternary diagrams, determines sums and ratios of specific elements of petrologic interest, and plots 33 X-Y graphs and five ternary diagrams. PETRO.CALC.PLOT also may be used to create other diagrams as desired by the user. The macros run in Microsoft Excel 3.0 and 4.0 for Macintosh computers and in Microsoft Excel 3.0 and 4.0 for Windows. Macros provided in PETRO.CALC.PLOT minimize repetition and time required to recalculate and plot whole-rock oxide data for petrologic analysis. ?? 1994.

Computers & Geosciences↗
USGS1994

JKPLOT version 2.00: A device-independent plotting system written in QuickBasic for an IBM PC

Version 2.0 of JKPLOT is a simple, device-independent plotting system that can provide a base for building plotting facilities tailored to the needs of geologists. The programs were written in Microsoft QuickBasic for use on equipment compatible with an IBM personal computer. The system provides subroutines and functions for plotting in both 2 and 3 dimensions with the capability for making scaled plots with annotated axes. Plot output can be directed to the PC screen, an Epson printer, any HPGL pen plotter, or an HP laser printer. The user's documentation and program diskettes of the JKPLOT system are available from U.S. Geological Survey Open File Services. ?? 1994.

Computers & Geosciences↗
USGS1994

Digital hum filtering

Data may be overprinted by a steady-state cyclical noise (hum). Steady-state indicates that the noise is invariant with time; its attributes, frequency, amplitude, and phase, do not change with time. Hum recorded on seismic data usually is powerline noise and associated higher harmonics; leakage from full-waveform rectified cathodic protection devices that contain the odd higher harmonics of powerline frequencies; or vibrational noise from mechanical devices. The fundamental frequency of powerline hum may be removed during data acquisition with the use of notch filters. Unfortunately, notch filters do not discriminate signal and noise, attenuating both. They also distort adjacent frequencies by phase shifting. Finally, they attenuate only the fundamental mode of the powerline noise; higher harmonics and frequencies other than that of powerlines are not removed. Digital notch filters, applied during processing, have many of the same problems as analog filters applied in the field. The method described here removes hum of a particular frequency. Hum attributes are measured by discrete Fourier analysis, and the hum is canceled from the data by subtraction. Errors are slight and the result of the presence of (random) noise in the window or asynchrony of the hum and data sampling. Error is minimized by increasing window size or by resampling to a finer interval. Errors affect the degree of hum attenuation, not the signal. The residual is steady-state hum of the same frequency. ?? 1994.

Computers & Geosciences↗
USGS1994

Dissolution of bedded rock salt: A seismic profile across the active eastern margin of the Hutchinson Salt Member, central Kansas

Since late Tertiary, bedded rock salt of the Permian Hutchinson Salt Member has been dissolved more-or-less continuously along its active eastern margin in central Kansas as a result of sustained contact with unconfined, undersaturated groundwater. The associated westward migration of the eastern margin has resulted in surface subsidence and the contemporaneous sedimentation of predominantly valley-filling Quarternary alluvium. In places, these alluvium deposits extend more than 25 km to the east of the present-day edge of the main body of contiguous rock salt. The margin could have receded this distance during the past several million years. From an environmental perspective, the continued leaching of the Hutchinson Salt is a major concern. This predominantly natural dissolution occurs in a broad zone across the central part of the State and adversely affects groundwater and surface-water quality as nonpoint source pollution. Significant surface subsidence occurs as well. Most of these subsidence features have formed gradually; others developed in a more catastrophic manner. The latter in particular pose real threats to roadways, railways, and buried oil and gas pipelines. In an effort to further clarify the process of natural salt dissolution in central Kansas and with the long-term goal of mitigating the adverse environmental affects of such leaching, the Kansas Geological Survey acquired a 4-km seismic profile across the eastern margin of the Hutchinson Salt in the Punkin Center area of central Kansas. The interpretation of these seismic data (and supporting surficial and borehole geologic control) is consistent with several hypotheses regarding the process and mechanisms of dissolution. More specifically these data support the theses that: 1. (1) Dissolution along the active eastern margin of the Hutchinson Salt Member was initiated during late Tertiary. Leaching has resulted in the steady westward migration of the eastern margin, surface subsidence, and the contemporaneous deposition of predominantly valley-filling Quarternary alluvium. 2. (2) Along the active eastern margin, the rock salt has been leached vertically from the top down, and horizontally along the uppermost remnant bedded soluble layer(s). As a result, the eastern margin thickens gradually (up to 90 m) and in a stepwise manner from east to west for distances on the order 5-15 km. 3. (3) In places, the Hutchinson Salt Member has been leached locally along NNE-trending paleoshear zones situated to the west of the present-day edge of the main body of contiguous rock salt. Leaching at these sites initiated when the main dissolution front impinged upon preexisting shear zones. ?? 1994.

Computers & Geosciences↗
USGS1995

Environmental probabilistic quantitative assessment methodologies

In this paper, four petroleum resource assessment methodologies are presented as possible pollution assessment methodologies, even though petroleum as a resource is desirable, whereas pollution is undesirable. A methodology is defined in this paper to consist of a probability model and a probabilistic method, where the method is used to solve the model. The following four basic types of probability models are considered: 1) direct assessment, 2) accumulation size, 3) volumetric yield, and 4) reservoir engineering. Three of the four petroleum resource assessment methodologies were written as microcomputer systems, viz. TRIAGG for direct assessment, APRAS for accumulation size, and FASPU for reservoir engineering. A fourth microcomputer system termed PROBDIST supports the three assessment systems. The three assessment systems have different probability models but the same type of probabilistic method. The type of advantages of the analytic method are in computational speed and flexibility, making it ideal for a microcomputer. -from Author

Computers & Geosciences↗
USGS1996

Semivariogram modeling by weighted least squares

Permissible semivariogram models are fundamental for geostatistical estimation and simulation of attributes having a continuous spatiotemporal variation. The usual practice is to fit those models manually to experimental semivariograms. Fitting by weighted least squares produces comparable results to fitting manually in less time, systematically, and provides an Akaike information criterion for the proper comparison of alternative models. We illustrate the application of a computer program with examples showing the fitting of simple and nested models. Copyright ?? 1996 Elsevier Science Ltd.

Computers & Geosciences↗
USGS1996

A digital photogrammetric method for measuring horizontal surficial movements on the slumgullion earthflow, Hinsdale county, Colorado

The traditional approach to making aerial photographic measurements uses analog or analytic photogrammetric equipment. We have developed a digital method for making measurements from aerial photographs which uses geographic information system (GIS) software, and primarily DOS-based personal computers. This method, which is based on the concept that a direct visual comparison can be made between images derived from two sets of aerial photographs taken at different times, was applied to the surface of the active portion of the Slumgullion earthflow in Colorado to determine horizontal displacement vectors from the movements of visually identifiable objects, such as trees and large rocks. Using this method, more of the slide surface can be mapped in a shorter period of time than using the standard photogrammetric approach. More than 800 horizontal displacement vectors were determined on the active earthflow surface using images produced by our digital photogrammetric technique and 1985 (1:12,000-scale) and 1990 (1:6,000-scale) aerial photographs. The resulting displacement field shows, with a 2-m measurement error (??? 10%), that the fastest moving portion of the landslide underwent 15-29 m of horizontal displacement between 1985 and 1990. Copyright ?? 1996 Elsevier Science Ltd.

Computers & Geosciences↗
USGS2004

CORRELATOR 5.2 - A program for interactive lithostratigraphic correlation of wireline logs

The limited radius of investigation of petrophysical measurements made in boreholes and the relatively large distances between wells result in an incomplete sensing of the subsurface through well logging. CORRELATOR is a program for estimating geological properties between logged boreholes. An initial and fundamental step is the lithostratigraphic correlation of logs in different wells. The method employed by the program closely emulates the process of visual inspection used by experienced subsurface geologists in manual correlation. Mathematically, the determination of lithostratigraphical equivalence is based on the simultaneous assessment of similarity in shale content, similarity in the patterns of vertical variation in a petrophysical property that is measured with high vertical resolution, and spatial consistency of stratigraphic relationships as determined by an expert system. Multiple additional options for processing log readings allow maximization in the extraction of information from pairs of logs per well and great flexibility in the final display of results in the form of cross sections and dip diagrams. ?? 2004 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗
USGS2004

A visual basic program to generate sediment grain-size statistics and to extrapolate particle distributions

Measures that describe and summarize sediment grain-size distributions are important to geologists because of the large amount of information contained in textural data sets. Statistical methods are usually employed to simplify the necessary comparisons among samples and quantify the observed differences. The two statistical methods most commonly used by sedimentologists to describe particle distributions are mathematical moments (Krumbein and Pettijohn, 1938) and inclusive graphics (Folk, 1974). The choice of which of these statistical measures to use is typically governed by the amount of data available (Royse, 1970). If the entire distribution is known, the method of moments may be used; if the next to last accumulated percent is greater than 95, inclusive graphics statistics can be generated. Unfortunately, earlier programs designed to describe sediment grain-size distributions statistically do not run in a Windows environment, do not allow extrapolation of the distribution's tails, or do not generate both moment and graphic statistics (Kane and Hubert, 1963; Collias et al., 1963; Schlee and Webster, 1967; Poppe et al., 2000) 1 . Owing to analytical limitations, electro-resistance multichannel particle-size analyzers, such as Coulter Counters, commonly truncate the tails of the fine-fraction part of grain-size distributions. These devices do not detect fine clay in the 0.6–0.1 μm range (part of the 11-phi and all of the 12-phi and 13-phi fractions). Although size analyses performed down to 0.6 μm microns are adequate for most freshwater and near shore marine sediments, samples from many deeper water marine environments (e.g. rise and abyssal plain) may contain significant material in the fine clay fraction, and these analyses benefit from extrapolation. The program (GSSTAT) described herein generates statistics to characterize sediment grain-size distributions and can extrapolate the fine-grained end of the particle distribution. It is written in Microsoft Visual Basic 6.0 and provides a window to facilitate program execution. The input for the sediment fractions is weight percentages in whole-phi notation (Krumbein, 1934; Inman, 1952), and the program permits the user to select output in either method of moments or inclusive graphics statistics (Fig. 1). Users select options primarily with mouse-click events, or through interactive dialogue boxes.

Computers & Geosciences↗
USGS2005

Use of XML and Java for collaborative petroleum reservoir modeling on the Internet

The GEMINI (Geo-Engineering Modeling through INternet Informatics) is a public-domain, web-based freeware that is made up of an integrated suite of 14 Java-based software tools to accomplish on-line, real-time geologic and engineering reservoir modeling. GEMINI facilitates distant collaborations for small company and academic clients, negotiating analyses of both single and multiple wells. The system operates on a single server and an enterprise database. External data sets must be uploaded into this database. Feedback from GEMINI users provided the impetus to develop Stand Alone Web Start Applications of GEMINI modules that reside in and operate from the user's PC. In this version, the GEMINI modules run as applets, which may reside in local user PCs, on the server, or Java Web Start. In this enhanced version, XML-based data handling procedures are used to access data from remote and local databases and save results for later access and analyses. The XML data handling process also integrates different stand-alone GEMINI modules enabling the user(s) to access multiple databases. It provides flexibility to the user to customize analytical approach, database location, and level of collaboration. An example integrated field-study using GEMINI modules and Stand Alone Web Start Applications is provided to demonstrate the versatile applicability of this freeware for cost-effective reservoir modeling. ?? 2005 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗
USGS2005

Dynamic modeling of Tampa Bay urban development using parallel computing

Urban land use and land cover has changed significantly in the environs of Tampa Bay, Florida, over the past 50 years. Extensive urbanization has created substantial change to the region's landscape and ecosystems. This paper uses a dynamic urban-growth model, SLEUTH, which applies six geospatial data themes (slope, land use, exclusion, urban extent, transportation, hillside), to study the process of urbanization and associated land use and land cover change in the Tampa Bay area. To reduce processing time and complete the modeling process within an acceptable period, the model is recoded and ported to a Beowulf cluster. The parallel-processing computer system accomplishes the massive amount of computation the modeling simulation requires. SLEUTH calibration process for the Tampa Bay urban growth simulation spends only 10 h CPU time. The model predicts future land use/cover change trends for Tampa Bay from 1992 to 2025. Urban extent is predicted to double in the Tampa Bay watershed between 1992 and 2025. Results show an upward trend of urbanization at the expense of a decline of 58% and 80% in agriculture and forested lands, respectively.

Computers & Geosciences↗
USGS2005

Statistical analysis of water-quality data containing multiple detection limits: S-language software for regression on order statistics

Trace contaminants in water, including metals and organics, often are measured at sufficiently low concentrations to be reported only as values below the instrument detection limit. Interpretation of these "less thans" is complicated when multiple detection limits occur. Statistical methods for multiply censored, or multiple-detection limit, datasets have been developed for medical and industrial statistics, and can be employed to estimate summary statistics or model the distributions of trace-level environmental data. We describe S-language-based software tools that perform robust linear regression on order statistics (ROS). The ROS method has been evaluated as one of the most reliable procedures for developing summary statistics of multiply censored data. It is applicable to any dataset that has 0 to 80% of its values censored. These tools are a part of a software library, or add-on package, for the R environment for statistical computing. This library can be used to generate ROS models and associated summary statistics, plot modeled distributions, and predict exceedance probabilities of water-quality standards. ?? 2005 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗
USGS2006

Nonlinear inversion of potential-field data using a hybrid-encoding genetic algorithm

Using a genetic algorithm to solve an inverse problem of complex nonlinear geophysical equations is advantageous because it does not require computer gradients of models or "good" initial models. The multi-point search of a genetic algorithm makes it easier to find the globally optimal solution while avoiding falling into a local extremum. As is the case in other optimization approaches, the search efficiency for a genetic algorithm is vital in finding desired solutions successfully in a multi-dimensional model space. A binary-encoding genetic algorithm is hardly ever used to resolve an optimization problem such as a simple geophysical inversion with only three unknowns. The encoding mechanism, genetic operators, and population size of the genetic algorithm greatly affect search processes in the evolution. It is clear that improved operators and proper population size promote the convergence. Nevertheless, not all genetic operations perform perfectly while searching under either a uniform binary or a decimal encoding system. With the binary encoding mechanism, the crossover scheme may produce more new individuals than with the decimal encoding. On the other hand, the mutation scheme in a decimal encoding system will create new genes larger in scope than those in the binary encoding. This paper discusses approaches of exploiting the search potential of genetic operations in the two encoding systems and presents an approach with a hybrid-encoding mechanism, multi-point crossover, and dynamic population size for geophysical inversion. We present a method that is based on the routine in which the mutation operation is conducted in the decimal code and multi-point crossover operation in the binary code. The mix-encoding algorithm is called the hybrid-encoding genetic algorithm (HEGA). HEGA provides better genes with a higher probability by a mutation operator and improves genetic algorithms in resolving complicated geophysical inverse problems. Another significant result is that final solution is determined by the average model derived from multiple trials instead of one computation due to the randomness in a genetic algorithm procedure. These advantages were demonstrated by synthetic and real-world examples of inversion of potential-field data. ?? 2005 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗
USGS2006

An efficient Matlab script to calculate heterogeneous anisotropically elastic wave propagation in three dimensions

We have created a second-order finite-difference solution to the anisotropic elastic wave equation in three dimensions and implemented the solution as an efficient Matlab script. This program allows the user to generate synthetic seismograms for three-dimensional anisotropic earth structure. The code was written for teleseismic wave propagation in the 1-0.1 Hz frequency range but is of general utility and can be used at all scales of space and time. This program was created to help distinguish among various types of lithospheric structure given the uneven distribution of sources and receivers commonly utilized in passive source seismology. Several successful implementations have resulted in a better appreciation for subduction zone structure, the fate of a transform fault with depth, lithospheric delamination, and the effects of wavefield focusing and defocusing on attenuation. Companion scripts are provided which help the user prepare input to the finite-difference solution. Boundary conditions including specification of the initial wavefield, absorption and two types of reflection are available. ?? 2005 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗
USGS2006

KGS-HighK: A Fortran 90 program for simulation of hydraulic tests in highly permeable aquifers

Slug and pumping tests (hydraulic tests) are frequently used by hydrogeologists to obtain in-situ estimates of the transmissive and storage properties of a formation (Streltsova, 1988; Kruseman and de Ridder, 1990; Butler, 1998). In aquifers of high hydraulic conductivity, hydraulic tests are affected by mechanisms that are not considered in the analysis of tests in less permeable media (Bredehoeft et al., 1966). Inertia-induced oscillations in hydraulic head are the most common manifestation of such mechanisms. Over the last three decades, a number of analytical solutions that incorporate these mechanisms have been developed for the analysis of hydraulic tests in highly permeable aquifers (see Butler and Zhan (2004) for a review of this previous work). These solutions, however, are restricted to a subset of the conditions commonly encountered in the field. Recently, a more general solution has been developed that builds on this previous work to remove many of the limitations imposed by these earlier approaches (Butler and Zhan, 2004). The purpose of this note is to present a Fortran 90 program, KGS-HighK, for the evaluation of this new solution. This note begins with a brief overview of the conceptual model that motivated the development of the solution of Butler and Zhan (2004) for pumping- and slug-induced flow to/from a central well. The major steps in the derivation of that solution are described, but no details are given. Instead, a Mathematica notebook is provided for those interested in the derivation details. The key algorithms used in KGS-HighK are then described and the program structure is briefly outlined. A field example is provided to demonstrate program performance. The note concludes with a short summary section. ?? 2005 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗