Search USGS⌕ Search

SEARCH · Search USGS

Results for “Mathematical Geology”

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 37 records · Page 2Linked to original sources
USGS1988

Experiments with microcomputer-based artificial intelligence environments

The U.S. Geological Survey (USGS) has been experimenting with the use of relatively inexpensive microcomputers as artificial intelligence (AI) development environments. Several AI languages are available that perform fairly well on desk-top personal computers, as are low-to-medium cost expert system packages. Although performance of these systems is respectable, their speed and capacity limitations are questionable for serious earth science applications foreseen by the USGS. The most capable artificial intelligence applications currently are concentrated on what is known as the "artificial intelligence computer," and include Xerox D-series, Tektronix 4400 series, Symbolics 3600, VAX, LMI, and Texas Instruments Explorer. The artificial intelligence computer runs expert system shells and Lisp, Prolog, and Smalltalk programming languages. However, these AI environments are expensive. Recently, inexpensive 32-bit hardware has become available for the IBM/AT microcomputer. USGS has acquired and recently completed Beta-testing of the Gold Hill Systems 80386 Hummingboard, which runs Common Lisp on an IBM/AT microcomputer. Hummingboard appears to have the potential to overcome many of the speed/capacity limitations observed with AI-applications on standard personal computers. USGS is a Beta-test site for the Gold Hill Systems GoldWorks expert system. GoldWorks combines some high-end expert system shell capabilities in a medium-cost package. This shell is developed in Common Lisp, runs on the 80386 Hummingboard, and provides some expert system features formerly available only on AI-computers including frame and rule-based reasoning, on-line tutorial, multiple inheritance, and object-programming. ?? 1988 International Association for Mathematical Geology.

Mathematical Geology↗
USGS1995

Optimization techniques for integrating spatial data

Two optimization techniques ta predict a spatial variable from any number of related spatial variables are presented. The applicability of the two different methods for petroleum-resource assessment is tested in a mature oil province of the Midcontinent (USA). The information on petroleum productivity, usually not directly accessible, is related indirectly to geological, geophysical, petrographical, and other observable data. This paper presents two approaches based on construction of a multivariate spatial model from the available data to determine a relationship for prediction. In the first approach, the variables are combined into a spatial model by an algebraic map-comparison/integration technique. Optimal weights for the map comparison function are determined by the Nelder-Mead downhill simplex algorithm in multidimensions. Geologic knowledge is necessary to provide a first guess of weights to start the automatization, because the solution is not unique. In the second approach, active set optimization for linear prediction of the target under positivity constraints is applied. Here, the procedure seems to select one variable from each data type (structure, isopachous, and petrophysical) eliminating data redundancy. Automating the determination of optimum combinations of different variables by applying optimization techniques is a valuable extension of the algebraic map-comparison/integration approach to analyzing spatial data. Because of the capability of handling multivariate data sets and partial retention of geographical information, the approaches can be useful in mineral-resource exploration. ?? 1995 International Association for Mathematical Geology.

Mathematical Geology↗
USGS1993

Boundary assessment under uncertainty: A case study

Estimating certain attributes within a geological body whose exact boundary is not known presents problems because of the lack of information. Estimation may result in values that are inadmissible from a geological point of view, especially with attributes which necessarily must be zero outside the boundary, such as the thickness of the oil column outside a reservoir. A simple but effective way to define the boundary is to use indicator kriging in two steps, the first for the purpose of extrapolating control points outside the body, the second to obtain a weighting function which expresses the uncertainty attached to estimations obtained in the boundary region. ?? 1993 International Association for Mathematical Geology.

Mathematical Geology↗
USGS2002

Comparison of two probability distributions used to model sizes of undiscovered oil and gas accumulations: Does the tail wag the assessment?

Undiscovered oil and gas assessments are commonly reported as aggregate estimates of hydrocarbon volumes. Potential commercial value and discovery costs are, however, determined by accumulation size, so engineers, economists, decision makers, and sometimes policy analysts are most interested in projected discovery sizes. The lognormal and Pareto distributions have been used to model exploration target sizes. This note contrasts the outcomes of applying these alternative distributions to the play level assessments of the U.S. Geological Survey's 1995 National Oil and Gas Assessment. Using the same numbers of undiscovered accumulations and the same minimum, medium, and maximum size estimates, substitution of the shifted truncated lognormal distribution for the shifted truncated Pareto distribution reduced assessed undiscovered oil by 16% and gas by 15%. Nearly all of the volume differences resulted because the lognormal had fewer larger fields relative to the Pareto. The lognormal also resulted in a smaller number of small fields relative to the Pareto. For the Permian Basin case study presented here, reserve addition costs were 20% higher with the lognormal size assumption. ?? 2002 International Association for Mathematical Geology.

Mathematical Geology↗
USGS1993

Relationship of geological and geothermal field properties: Midcontinent area, USA, an example

Quantitative approaches to data analysis in the last decade have become important in basin modeling and mineral-resource estimation. The interrelation of geological, geophysical, geochemical, and geohydrological variables is important in adjusting a model to a real-world situation. Revealing the interdependences of variables can contribute in understanding the processes interacting in sedimentary basins. It is reasonably simple to compare spatial data of the same type but more difficult if different properties are involved. Statistical techniques, such as cluster analysis or principal components analysis, or some algebraic approaches can be used to ascertain the relations of standardized spatial data. In this example, structural configuration on five different stratigraphic horizons, one total sediment thickness map, and four maps of geothermal data were copared. As expected, the structural maps are highly related because all had undergone about the same deformation with differing degrees of intensity. The temperature gradients derived (1) from shallow borehole logging measurements under equilibrium conditions with the surrounding rock, and (2) from non-equilibrium bottom-hole temperatures (BHT) from deeper depths are mainly independent of each other. This was expected and confirmed also for the two temperature maps at 1000 ft which were constructed using both types of gradient values. Thus, it is evident that the use of a 2-point (BHT and surface temperature) straightline calculation of a mean temperature gradient gives different information about the geothermal regime than using gradients from temperatures logged under equilibrium conditions. Nevertheless, it is useful to determine to what a degree the larger dataset of nonequilibrium temperatures could reflect quantitative relationships to geologic conditions. Comparing all maps of geothermal information vs. the structural and the sediment thickness maps, it was determined that all correlations are moderately negative or slightly positive. These results are clearly shown by the cluster analysis and the principal components. Considering a close relationship between temperature and thermal conductivity of the sediments as observed for most of the Midcontinent area and relatively homogeneous heat-flow density conditions for the study area these results support the following assumptions: (1) undifferentiated geothermal gradients, computed from temperatures of different depth intervals and differing sediment properties, cannot contribute to an improved understanding of the temperature structure and its controls within the sedimentary cover, and (2) the quantitative approach of revealing such relations needs refined datasets of temperature information valid for the different depth levels or stratigraphic units. ?? 1993 International Association for Mathematical Geology.

Mathematical Geology↗
USGS1988

A statistical method for estimating rates of soil development and ages of geologic deposits: A design for soil-chronosequence studies

A statistical method for estimating rates of soil development in a given region based on calibration from a series of dated soils is used to estimate ages of soils in the same region that are not dated directly. The method is designed specifically to account for sampling procedures and uncertainties that are inherent in soil studies. Soil variation and measurement error, uncertainties in calibration dates and their relation to the age of the soil, and the limited number of dated soils are all considered. Maximum likelihood (ML) is employed to estimate a parametric linear calibration curve, relating soil development to time or age on suitably transformed scales. Soil variation on a geomorphic surface of a certain age is characterized by replicate sampling of soils on each surface; such variation is assumed to have a Gaussian distribution. The age of a geomorphic surface is described by older and younger bounds. This technique allows age uncertainty to be characterized by either a Gaussian distribution or by a triangular distribution using minimum, best-estimate, and maximum ages. The calibration curve is taken to be linear after suitable (in certain cases logarithmic) transformations, if required, of the soil parameter and age variables. Soil variability, measurement error, and departures from linearity are described in a combined fashion using Gaussian distributions with variances particular to each sampled geomorphic surface and the number of sample replicates. Uncertainty in age of a geomorphic surface used for calibration is described using three parameters by one of two methods. In the first method, upper and lower ages are specified together with a coverage probability; this specification is converted to a Gaussian distribution with the appropriate mean and variance. In the second method, "absolute" older and younger ages are specified together with a most probable age; this specification is converted to an asymmetric triangular distribution with mode at the most probable age. The statistical variability of the ML-estimated calibration curve is assessed by a Monte Carlo method in which simulated data sets repeatedly are drawn from the distributional specification; calibration parameters are reestimated for each such simulation in order to assess their statistical variability. Several examples are used for illustration. The age of undated soils in a related setting may be estimated from the soil data using the fitted calibration curve. A second simulation to assess age estimate variability is described and applied to the examples. ?? 1988 International Association for Mathematical Geology.

Mathematical Geology↗
USGS2005

United States geological survey's reserve-growth models and their implementation

The USGS has developed several mathematical models to forecast reserve growth of fields both in the United States (U.S.) and the world. The models are based on historical reserve growth patterns of fields in the U.S. The patterns of past reserve growth are extrapolated to forecast future reserve growth. Changes of individual field sizes through time are extremely variable, therefore, the reserve growth models take on a statistical approach whereby volumetric changes for populations of fields are used in the models. Field age serves as a measure of the field-development effort that is applied to promote reserve growth. At the time of the USGS World Petroleum Assessment 2000, a reserve growth model for discovered fields of the world was not available. Reserve growth forecasts, therefore, were made based on a model of historical reserve growth of fields of the U.S. To test the feasibility of such an application, reserve growth forecasts were made of 186 giant oil fields of the world (excluding the U.S. and Canada). In addition, forecasts were made for these giant oil fields subdivided into those located in and outside of Organization of Petroleum Exporting Countries (OPEC). The model provided a reserve-growth forecast that closely matched the actual reserve growth that occurred from 1981 through 1996 for the 186 fields as a whole, as well as for both OPEC and non-OPEC subdivisions, despite the differences in reserves definition among the fields of the U.S. and the rest of the world. ?? 2005 International Association for Mathematical Geology.

Conference Paper↗
USGS2005

Guiding principles of USGS methodology for assessment of undiscovered conventional oil and gas resources

During the last 30 years, the methodology for assessment of undiscovered conventional oil and gas resources used by the Geological Survey has undergone considerable change. This evolution has been based on five major principles. First, the U.S. Geological Survey has responsibility for a wide range of U.S. and world assessments and requires a robust methodology suitable for immaturely explored as well as maturely explored areas. Second, the assessments should be based on as comprehensive a set of geological and exploration history data as possible. Third, the perils of methods that solely use statistical methods without geological analysis are recognized. Fourth, the methodology and course of the assessment should be documented as transparently as possible, within the limits imposed by the inevitable use of subjective judgement. Fifth, the multiple uses of the assessments require a continuing effort to provide the documentation in such ways as to increase utility to the many types of users. Undiscovered conventional oil and gas resources are those recoverable volumes in undiscovered, discrete, conventional structural or stratigraphic traps. The USGS 2000 methodology for these resources is based on a framework of assessing numbers and sizes of undiscovered oil and gas accumulations and the associated risks. The input is standardized on a form termed the Seventh Approximation Data Form for Conventional Assessment Units. Volumes of resource are then calculated using a Monte Carlo program named Emc2, but an alternative analytic (non-Monte Carlo) program named ASSESS also can be used. The resource assessment methodology continues to change. Accumulation-size distributions are being examined to determine how sensitive the results are to size-distribution assumptions. The resource assessment output is changing to provide better applicability for economic analysis. The separate methodology for assessing continuous (unconventional) resources also has been evolving. Further studies of the relationship between geologic models of conventional and continuous resources will likely impact the respective resource assessment methodologies. ?? 2005 International Association for Mathematical Geology.

Conference Paper↗
USGS2007

Porphyry copper deposit tract definition - A global analysis comparing geologic map scales

Geologic maps are a fundamental data source used to define mineral-resource potential tracts for the first step of a mineral resource assessment. Further, it is generally believed that the scale of the geologic map is a critical consideration. Previously published research has demonstrated that the U.S. Geological Survey porphyry tracts identified for the United States, which are based on 1:500,000-scale geology and larger scale data and published at 1:1,000,000 scale, can be approximated using a more generalized 1:2,500,000-scale geologic map. Comparison of the USGS porphyry tracts for the United States with weights-of-evidence models made using a 1:10,000,000-scale geologic map, which was made for petroleum applications, and a 1:35,000,000-scale geologic map, which was created as context for the distribution of porphyry deposits, demonstrates that, again, the USGS US porphyry tracts identified are similar to tracts defined on features from these small scale maps. In fact, the results using the 1:35,000,000-scale map show a slightly higher correlation with the USGS US tract definition, probably because the conceptual context for this small-scale map is more appropriate for porphyry tract definition than either of the other maps. This finding demonstrates that geologic maps are conceptual maps. The map information shown in each map is selected and generalized for the map to display the concepts deemed important for the map maker's purpose. Some geologic maps of small scale prove to be useful for regional mineral-resource tract definition, despite the decrease in spatial accuracy with decreasing scale. The utility of a particular geologic map for a particular application is critically dependent on the alignment of the intention of the map maker with the application. ?? International Association for Mathematical Geology 2007.

Natural Resources Research↗
USGS2008

Are fractal dimensions of the spatial distribution of mineral deposits meaningful?

It has been proposed that the spatial distribution of mineral deposits is bifractal. An implication of this property is that the number of deposits in a permissive area is a function of the shape of the area. This is because the fractal density functions of deposits are dependent on the distance from known deposits. A long thin permissive area with most of the deposits in one end, such as the Alaskan porphyry permissive area, has a major portion of the area far from known deposits and consequently a low density of deposits associated with most of the permissive area. On the other hand, a more equi-dimensioned permissive area, such as the Arizona porphyry permissive area, has a more uniform density of deposits. Another implication of the fractal distribution is that the Poisson assumption typically used for estimating deposit numbers is invalid. Based on datasets of mineral deposits classified by type as inputs, the distributions of many different deposit types are found to have characteristically two fractal dimensions over separate non-overlapping spatial scales in the range of 5-1000 km. In particular, one typically observes a local dimension at spatial scales less than 30-60 km, and a regional dimension at larger spatial scales. The deposit type, geologic setting, and sample size influence the fractal dimensions. The consequence of the geologic setting can be diminished by using deposits classified by type. The crossover point between the two fractal domains is proportional to the median size of the deposit type. A plot of the crossover points for porphyry copper deposits from different geologic domains against median deposit sizes defines linear relationships and identifies regions that are significantly underexplored. Plots of the fractal dimension can also be used to define density functions from which the number of undiscovered deposits can be estimated. This density function is only dependent on the distribution of deposits and is independent of the definition of the permissive area. Density functions for porphyry copper deposits appear to be significantly different for regions in the Andes, Mexico, United States, and western Canada. Consequently, depending on which regional density function is used, quite different estimates of numbers of undiscovered deposits can be obtained. These fractal properties suggest that geologic studies based on mapping at scales of 1:24,000 to 1:100,000 may not recognize processes that are important in the formation of mineral deposits at scales larger than the crossover points at 30-60 km. ?? 2008 International Association for Mathematical Geology.

Natural Resources Research↗
USGS2009

Evaluation of procedures for prediction of unconventional gas in the presence of geologic trends

This study extends the application of local spatial nonparametric prediction models to the estimation of recoverable gas volumes in continuous-type gas plays to regimes where there is a single geologic trend. A transformation is presented, originally proposed by Tomczak, that offsets the distortions caused by the trend. This article reports on numerical experiments that compare predictive and classification performance of the local nonparametric prediction models based on the transformation with models based on Euclidean distance. The transformation offers improvement in average root mean square error when the trend is not severely misspecified. Because of the local nature of the models, even those based on Euclidean distance in the presence of trends are reasonably robust. The tests based on other model performance metrics such as prediction error associated with the high-grade tracts and the ability of the models to identify sites with the largest gas volumes also demonstrate the robustness of both local modeling approaches. ?? International Association for Mathematical Geology 2009.

Natural Resources Research↗
USGS2008

Map scale effects on estimating the number of undiscovered mineral deposits

Estimates of numbers of undiscovered mineral deposits, fundamental to assessing mineral resources, are affected by map scale. Where consistently defined deposits of a particular type are estimated, spatial and frequency distributions of deposits are linked in that some frequency distributions can be generated by processes randomly in space whereas others are generated by processes suggesting clustering in space. Possible spatial distributions of mineral deposits and their related frequency distributions are affected by map scale and associated inclusions of non-permissive or covered geological settings. More generalized map scales are more likely to cause inclusion of geologic settings that are not really permissive for the deposit type, or that include unreported cover over permissive areas, resulting in the appearance of deposit clustering. Thus, overly generalized map scales can cause deposits to appear clustered. We propose a model that captures the effects of map scale and the related inclusion of non-permissive geologic settings on numbers of deposits estimates, the zero-inflated Poisson distribution. Effects of map scale as represented by the zero-inflated Poisson distribution suggest that the appearance of deposit clustering should diminish as mapping becomes more detailed because the number of inflated zeros would decrease with more detailed maps. Based on observed worldwide relationships between map scale and areas permissive for deposit types, mapping at a scale with twice the detail should cut permissive area size of a porphyry copper tract to 29% and a volcanic-hosted massive sulfide tract to 50% of their original sizes. Thus some direct benefits of mapping an area at a more detailed scale are indicated by significant reductions in areas permissive for deposit types, increased deposit density and, as a consequence, reduced uncertainty in the estimate of number of undiscovered deposits. Exploration enterprises benefit from reduced areas requiring detailed and expensive exploration, and land-use planners benefit from reduced areas of concern. ?? 2008 International Association for Mathematical Geology.

Natural Resources Research↗
USGS2010

Quality of the log-geometric distribution extrapolation for smaller undiscovered oil and gas pool size

The U.S. Geological Survey procedure for the estimation of the general form of the parent distribution requires that the parameters of the log-geometric distribution be calculated and analyzed for the sensitivity of these parameters to different conditions. In this study, we derive the shape factor of a log-geometric distribution from the ratio of frequencies between adjacent bins. The shape factor has a log straight-line relationship with the ratio of frequencies. Additionally, the calculation equations of a ratio of the mean size to the lower size-class boundary are deduced. For a specific log-geometric distribution, we find that the ratio of the mean size to the lower size-class boundary is the same. We apply our analysis to simulations based on oil and gas pool distributions from four petroleum systems of Alberta, Canada and four generated distributions. Each petroleum system in Alberta has a different shape factor. Generally, the shape factors in the four petroleum systems stabilize with the increase of discovered pool numbers. For a log-geometric distribution, the shape factor becomes stable when discovered pool numbers exceed 50 and the shape factor is influenced by the exploration efficiency when the exploration efficiency is less than 1. The simulation results show that calculated shape factors increase with those of the parent distributions, and undiscovered oil and gas resources estimated through the log-geometric distribution extrapolation are smaller than the actual values. ?? 2010 International Association for Mathematical Geology.

Natural Resources Research↗
USGS2021

Realizations

In statistics, a realization is an observed value of a random variable (Gubner 2006 ). In mathematical geology, the most important realizations are those in the form of maps of spatially correlated regionalized variables. Spatial description of random variables within complex domains and making certain decisions about those require complete knowledge of the attribute of interest at each point in space. However, it is virtually impossible to sample from every location within the domain to gain a complete spatial understanding of the random variables with certainty at different scales. Therefore, limited sampling leaves us with incomplete information, which is the source of uncertainty. Understanding the uncertainty and quantifying it are essential to minimize the risks of decision making. Geostatistical simulation techniques aim to quantify spatial uncertainty of random variables by numerically reproducing the reality, which we have limited knowledge of, in a discretized...

Book chapter↗
USGS2008

Volatility of bitumen prices and implications for the industry

Sustained crude oil price increases have led to increased investment in and production of Canadian bitumen to supplement North American oil supplies. For new projects, the evaluation of profitability is based on a prediction of the future price path of bitumen and ultimately light/medium crude oil. This article examines the relationship between the bitumen and light crude oil prices in the context of a simple error-correction economic-adjustment model. The analysis shows bitumen prices to be significantly more volatile than light crude prices. Also, the dominant effect of an oil price shock on bitumen prices is immediate and is amplified, both in absolute terms and percentage price changes. It is argued that the bitumen industry response to such market risks will likely be a realignment toward vertical integration via new downstream construction, mergers, or on a de facto basis by the establishment of alliances. ?? 2008 International Association for Mathematical Geology.

Natural Resources Research↗
USGS2009

A test and re-estimation of Taylor's empirical capacity-reserve relationship

In 1977, Taylor proposed a constant elasticity model relating capacity choice in mines to reserves. A test of this model using a very large (n = 1,195) dataset confirms its validity but obtains significantly different estimated values for the model coefficients. Capacity is somewhat inelastic with respect to reserves, with an elasticity of 0.65 estimated for open-pit plus block-cave underground mines and 0.56 for all other underground mines. These new estimates should be useful for capacity determinations as scoping studies and as a starting point for feasibility studies. The results are robust over a wide range of deposit types, deposit sizes, and time, consistent with physical constraints on mine capacity that are largely independent of technology. ?? 2009 International Association for Mathematical Geology.

Natural Resources Research↗
USGS2000

Analyzing slug tests in wells screened across the watertable: A field assessment

The slug test is the most widely used technique for the in situ estimation of hydraulic conductivity in confined and unconfined formations. Currently, there are no generally accepted methods in the groundwater literature for the analysis of response data from slug tests performed in wells screened across the watertable. A field study was undertaken in an attempt to develop a set of practical guidelines for tests conducted in such wells. Three wells, screened within unconsolidated material exhibiting a range of hydraulic conductivities (.05-30.0 m/day), were installed to depths of up to 9 m (30 ft) in Kansas River alluvium that ranges in thickness from 15 m to 21 m (50 ft to 70 ft) near Lawrence, Kansas. Intensive well-development efforts removed any drilling debris that could interfere with well-formation hydraulics. Once the wells were developed properly, a series of slug tests was performed at each well. The tests were designed to assess the role of the unsaturated zone and the appropriateness of assuming a fixed hydraulic head upper boundary. The results of this investigation can be summarized as follows: (1) the sufficiency of well development should be based on repeat slug tests and not the clarity of pumped water; (2) the effective screen radius for best model analysis should be based on a mass balance and not nominal screen dimensions; (3) the watertable can be represented as a constant head boundary and flow in the unsaturated zone can be ignored in most situations; (4) conventional techniques for the analysis of slug-test data seem to be reasonable for slug tests conducted in wells screened across the watertable, when used with the appropriate effective screen radius and normalized head range; and (5) fluctuations in the watertable elevation through time can be exploited to obtain some insight into the nature of vertical variation in hydraulic conductivity at a well. The results of this investigation indicate that multiple slug tests should be performed at wells screened across the watertable in order to reliably assess the sufficiency of well development and the appropriateness of conventional theory. ?? 2000 International Association for Mathematical Geology.

Natural Resources Research↗
USGS1996

Application of a feedforward neural network in the search for kuroko deposits in the Hokuroku District, Japan

A feedforward neural network with one hidden layer and five neurons was trained to recognize the distance to kuroko mineral deposits. Average amounts per hole of pyrite, sericite, and gypsum plus anhydrite as measured by X-rays in 69 drillholes were used to train the net. Drillholes near and between the Fukazawa, Furutobe, and Shakanai mines were used. The training data were selected carefully to represent well-explored areas where some confidence of the distance to ore was assured. A logarithmic transform was applied to remove the skewness of distance and each variable was scaled and centered by subtracting the median and dividing by the interquartile range. The learning algorithm of annealing plus conjugate gradients was used to minimize the mean squared error of the scaled distance to ore. The trained network then was applied to all of the 152 drillholes that had measured gypsum, sericite, and pyrite. A contour plot of the neural net predicted distance to ore shows fairly wide areas of 1 km or less to ore; each of the known deposit groups is within the 1 km contour. The high and low distances on the margins of the contoured distance plot are in part the result of boundary effects of the contouring algorithm. For example, the short distances to ore predicted west of the Shakanai (Hanaoka) deposits are in basement. However, the short distances to ore predicted northeast of Furotobe, just off the figure, coincide with the location of the Nurukawa kuroko deposit and the Omaki deposit, south of the Shakanai-Hanaoka deposits, seems to be on an extension of short distance to ore contour, but is beyond the 3 km limit from drillholes. Also of interest are some areas only a few kilometers from the Fukazawa and Shakanai groups of deposits that are estimated to be many kilometers from ore, apparently reflecting the network's recognition of the extreme local variability of the geology near some deposits.

Mathematical Geology↗