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

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↗

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↗

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↗

Discovery sequence and the nature of low permeability gas accumulations

There is an ongoing discussion regarding the geologic nature of accumulations that host gas in low-permeability sandstone environments. This note examines the discovery sequence of the accumulations in low permeability sandstone plays that were classified as continuous-type by the U.S. Geological Survey for the 1995 National Oil and Gas Assessment. It compares the statistical character of historical discovery sequences of accumulations associated with continuous-type sandstone gas plays to those of conventional plays. The seven sandstone plays with sufficient data exhibit declining size with sequence order, on average, and in three of the seven the trend is statistically significant. Simulation experiments show that both a skewed endowment size distribution and a discovery process that mimics sampling proportional to size are necessary to generate a discovery sequence that consistently produces a statistically significant negative size order relationship. The empirical findings suggest that discovery sequence could be used to constrain assessed gas in untested areas. The plays examined represent 134 of the 265 trillion cubic feet of recoverable gas assessed in undeveloped areas of continuous-type gas plays in low permeability sandstone environments reported in the 1995 National Assessment. ?? 2005 International Association for Mathematical Geology.

Natural Resources Research↗

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↗

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↗

Probabilistic prediction models for aggregate quarry siting

Weights-of-evidence (WofE) and logistic regression techniques were used in a GIS framework to predict the spatial likelihood (prospectivity) of crushed-stone aggregate quarry development. The joint conditional probability models, based on geology, transportation network, and population density variables, were defined using quarry location and time of development data for the New England States, North Carolina, and South Carolina, USA. The Quarry Operation models describe the distribution of active aggregate quarries, independent of the date of opening. The New Quarry models describe the distribution of aggregate quarries when they open. Because of the small number of new quarries developed in the study areas during the last decade, independent New Quarry models have low parameter estimate reliability. The performance of parameter estimates derived for Quarry Operation models, defined by a larger number of active quarries in the study areas, were tested and evaluated to predict the spatial likelihood of new quarry development. Population density conditions at the time of new quarry development were used to modify the population density variable in the Quarry Operation models to apply to new quarry development sites. The Quarry Operation parameters derived for the New England study area, Carolina study area, and the combined New England and Carolina study areas were all similar in magnitude and relative strength. The Quarry Operation model parameters, using the modified population density variables, were found to be a good predictor of new quarry locations. Both the aggregate industry and the land management community can use the model approach to target areas for more detailed site evaluation for quarry location. The models can be revised easily to reflect actual or anticipated changes in transportation and population features. ?? International Association for Mathematical Geology 2007.

Natural Resources Research↗

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↗

A new reserve growth model for United States oil and gas fields

Reserve (or field) growth, which is an appreciation of total ultimate reserves through time, is a well-recognized phenomenon, particularly in mature petroleum provinces. The importance of forecasting reserve growth accurately in a mature petroleum province made it necessary to develop improved growth functions, and a critical review of the original Arrington method was undertaken. During a five-year (1992-1996), the original Arrington method gave 1.03% higher than the actual oil reserve growth, whereas the proposed modified method gave a value within 0.3% of the actual growth, and therefore it was accepted for the development for reserve growth models. During a five-year (1992-1996), the USGS 1995 National Assessment gave 39.3% higher oil and 33.6% lower gas than the actual growths, whereas the new model based on Modified Arrington method gave 11.9% higher oil and 29.8% lower gas than the actual growths. The new models forecast predict reserve growths of 4.2 billion barrels of oil (2.7%) and 30.2 trillion cubic feet of gas (5.4%) for the conterminous U.S. for the next five years (1997-2001). ?? 2005 International Association for Mathematical Geology.

Natural Resources Research↗

Wellsite, laboratory, and mathematical techniques for determining sorbed gas content of coals and gas shales utilizing well cuttings

Drill cuttings can be used for desorption analyses but with more uncertainty than desorption analyses done with cores. Drill cuttings are not recommended to take the place of core, but in some circumstances, desorption work with cuttings can provide a timely and economic supplement to that of cores. The mixed lithologic nature of drill cuttings is primarily the source of uncertainty in their analysis for gas content, for it is unclear how to apportion the gas generated from both the coal and the dark-colored shale that is mixed in usually with the coal. In the Western Interior Basin Coal Basin in eastern Kansas (Pennsylvanian-age coals), dark-colored shales with normal (??? 100 API units) gamma-ray levels seem to give off minimal amounts of gas on the order of less than five standard cubic feet per ton (scf/ton). In some cuttings analyses this rule of thumb for gas content of the shale is adequate for inferring the gas content of coals, but shales with high-gamma-ray values (>150 API units) may yield several times this amount of gas. The uncertainty in desorption analysis of drill cuttings can be depicted graphically on a diagram identified as a "lithologic component sensitivity analysis diagram." Comparison of cuttings desorption results from nearby wells on this diagram, can sometimes yield an unique solution for the gas content of both a dark shale and coal mixed in a cuttings sample. A mathematical solution, based on equating the dry, ash-free gas-contents of the admixed coal and dark-colored shale, also yields results that are correlative to data from nearby cores. ?? 2007 International Association for Mathematical Geology.

Natural Resources Research↗

Assessment of exploration bias in data-driven predictive models and the estimation of undiscovered resources

The spatial distribution of discovered resources may not fully mimic the distribution of all such resources, discovered and undiscovered, because the process of discovery is biased by accessibility factors (e.g., outcrops, roads, and lakes) and by exploration criteria. In data-driven predictive models, the use of training sites (resource occurrences) biased by exploration criteria and accessibility does not necessarily translate to a biased predictive map. However, problems occur when evidence layers correlate with these same exploration factors. These biases then can produce a data-driven model that predicts known occurrences well, but poorly predicts undiscovered resources. Statistical assessment of correlation between evidence layers and map-based exploration factors is difficult because it is difficult to quantify the "degree of exploration." However, if such a degree-of-exploration map can be produced, the benefits can be enormous. Not only does it become possible to assess this correlation, but it becomes possible to predict undiscovered, instead of discovered, resources. Using geothermal systems in Nevada, USA, as an example, a degree-of-exploration model is created, which then is resolved into purely explored and unexplored equivalents, each occurring within coextensive study areas. A weights-of-evidence (WofE) model is built first without regard to the degree of exploration, and then a revised WofE model is calculated for the "explored fraction" only. Differences in the weights between the two models provide a correlation measure between the evidence and the degree of exploration. The data used to build the geothermal evidence layers are perceived to be independent of degree of exploration. Nevertheless, the evidence layers correlate with exploration because exploration has preferred the same favorable areas identified by the evidence patterns. In this circumstance, however, the weights for the "explored" WofE model minimize this bias. Using these revised weights, posterior probability is extrapolated into unexplored areas to estimate undiscovered deposits. ?? International Association for Mathematical Geology 2007.

Natural Resources Research↗

Development and implementation of a Bayesian-based aquifer vulnerability assessment in Florida

The Florida Aquifer Vulnerability Assessment (FAVA) was designed to provide a tool for environmental, regulatory, resource management, and planning professionals to facilitate protection of groundwater resources from surface sources of contamination. The FAVA project implements weights-of-evidence (WofE), a data-driven, Bayesian-probabilistic model to generate a series of maps reflecting relative aquifer vulnerability of Florida's principal aquifer systems. The vulnerability assessment process, from project design to map implementation is described herein in reference to the Floridan aquifer system (FAS). The WofE model calculates weighted relationships between hydrogeologic data layers that influence aquifer vulnerability and ambient groundwater parameters in wells that reflect relative degrees of vulnerability. Statewide model input data layers (evidential themes) include soil hydraulic conductivity, density of karst features, thickness of aquifer confinement, and hydraulic head difference between the FAS and the watertable. Wells with median dissolved nitrogen concentrations exceeding statistically established thresholds serve as training points in the WofE model. The resulting vulnerability map (response theme) reflects classified posterior probabilities based on spatial relationships between the evidential themes and training points. The response theme is subjected to extensive sensitivity and validation testing. Among the model validation techniques is calculation of a response theme based on a different water-quality indicator of relative recharge or vulnerability: dissolved oxygen. Successful implementation of the FAVA maps was facilitated by the overall project design, which included a needs assessment and iterative technical advisory committee input and review. Ongoing programs to protect Florida's springsheds have led to development of larger-scale WofE-based vulnerability assessments. Additional applications of the maps include land-use planning amendments and prioritization of land purchases to protect groundwater resources. ?? International Association for Mathematical Geology 2007.

Natural Resources Research↗

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↗

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↗

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↗

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↗

Geological applications of multi-response permutation procedures

The multi-purpose permutation procedures (MRPP) test statistic is designed to analyze multivariate data at the ordinal or higher levels. It is based on the weighted averages of symmetric distance functions over all paired objects within a priori disjoint groups of objects from a finite population of objects where each object's response is a point in an r-dimensional space. Thus, the arguments of each r-dimensional point correspond to the r measured responses of each object in the finite population of objects. The null hypothesis underlying MRPP is that the observed sample of objects within groups of a specified size structure is randomly obtained from the pooled collection of objects comprising the finite population. The procedure is used to test the presumed geomorphic differences among three reaches of the Lower Mississippi River. A combination of the proposed reaches is favored and several variable-based reach configurations are proposed. ?? 1993 International Association for Mathernatical Geology.

Mathematical Geology↗

U.S. Geological Survey probabilistic methodology for oil and gas resource appraisal of the United States

Probabilistic methodology used by the U.S. Geological Survey is described for estimating the quantity of undiscovered recoverable conventional resources of oil and gas in the United States. A judgmental probability distribution of the "quantity of resource" and its properties is determined for a geologic province or basin. From this distribution, point and interval estimates of the quantity of undiscovered resource are obtained. Distributions and their properties are established for each of the following resources: (1) oil and nonassociated gas from estimates of the probability of the resource being present and the conditional probability distribution of the quantity of resource given that the resource is present, (2) associated-dissolved gas from its corresponding oil distribution, (3) total gas, (4) oil and total gas in two or more provinces. Computer graphics routines are illustrated with examples from the U.S. Geological Survey Circular 860. ?? 1984 Plenum Publishing Corporation.

Journal of the International Association for Mathe↗