Optimization of grid-drilling using computer simulation

A computer simulation method has been developed to find efficient drilling grids for mineral deposits. A well-known ore deposit is used as a model to develop an efficient pattern for undiscovered ore bodies in the same area or in other prospects where similar geometry is suspected. The model for this study is the Austinville, Virginia deposit, a Mississippi Valley-type deposit composed of 17 ore bodies totaling 34 million short tons (30 million metric tons). The method employs a computer program that simulates drilling the model deposit with different patterns, including various levels of follow-up drilling. Follow-up holes are drilled in fences at one half the original spacing around holes in the grid that show ore-grade mineralization. Each pattern is drilled 100 times from random starting locations to provide a range of outcomes of drilling, including the best, worst, and most likely. For this study, patterns of 100 drill holes were composed of 10 fences spaced 1000–5000 feet (305–1524 m) apart, each with 10 holes spaced 200–1000 feet (61–305 m) apart. In all, 25 grids were used with zero to three levels of follow-up drilling. The 600/2000 grid, with drill holes spaced 600 feet (183 m) apart in fences spaced 2000 feet (610 m) apart, was compared with the 200/5000 grid because they represented contrasting outcomes. The 600/2000 grid penetrated many ore bodies consistently but with few multiple hits to individual ore bodies; whereas the 200/5000 grid inconsistently penetrated few ore bodies with many multiple hits. The 600/2000 grid was more efficient than the 200/5000 grid at hitting large ore bodies of 1,000,000 short tons or greater (900,000 metric tons or greater) and was made more effective by adding one cycle of follow-up drilling. The 600/2000 grid had a 97% chance of hitting one or more large ore bodies with at least one drill hole per ore body, and the 200/5000 grid had a 64% chance. Once hit, there was an 82% chance that the largest ore body would be penetrated by three or more holes when using the 600/2000 grid and an 88% chance using the 200/5000 grid.     

Planetary cratering 1. The question of multiple impactor populations: Lunar evidence

Abstract— This paper addresses several current issues related to use of craters in interpreting planetary surface histories. The primary goal is to test the widely adopted hypothesis of multiple populations of impactors at different times or places in the Solar System. New data presented here revise a “lunar highland” crater diameter distribution that has been widely used as evidence of an early distinct population of impactors. This curve, which has a depression of the size distribution at mid-sizes, does not, in fact, represent the lunar highlands generally. I show that it is associated with regions of intercrater plains. The more extensive the obliteration by intercrater plains, the deeper the depression. Modeling indicates that the depression of the curve is caused by the obliteration process itself. The oldest, most cratered regions of lunar highlands do not show the depression. These findings call into question earlier interpretations of multiple populations of impactors in the Solar System and of a distinctive primordial population. The present work is consistent, instead, with (1) a relatively uniform size distribution of interplanetary impactors, of mixed origins, back to 4 Ga ago and throughout the sampled Solar System; (2) fragmentation as the process that produced that size distribution; (3) saturation equilibrium on the most heavily cratered surfaces; and (4) differences in structure in the size distribution caused not by distinct impactor populations but by episodes of endogenic obliteration. If accepted, these results would modify some studies of solar system evolution, including assertions of two to five distinct populations of impactors, assumptions of lack of saturation equilibrium, and identifications of specific heliocentric or planetocentric sources for impactors within outer planet satellite systems.     

Structure and dynamics of the Bok globule B335

CO maps of the Bok globule B335 are presented and used to derive its density profile, mass distribution, and rotational velocity structure. It is found that the cloud is in nearly hydrostatic equilibrium with a density profile that varies roughly as r^?1 in the core and r^?3 in the envelope. The observed rotation is unimportant in the force balance at the present stage of evolution.     

On the practice of estimating fractal dimension

Coastlines epitomize deterministic fractals and fractal (Hausdorff-Besicovitch) dimensions; a divider [compass] method can be used to calculate fractal dimensions for these features. Noise models are used to develop another notion of fractals, a stochastic one. Spectral and variogram methods are used to estimate fractal dimensions for stochastic fractals. When estimating fractal dimension, the objective of the analysis must be consistent with the method chosen for fractal dimension calculation. Spectal and variogram methods yield fractal dimensions which indicate the similarity of the feature under study to noise (e.g., Brownian noise). A divider measurement method yields a fractal dimension which is a measure of complexity of shape.     

SHRIMP zircon dating and Sm/Nd isotopic investigations of Neoproterozoic granitoids, Eastern Desert, Egypt

There is an increasing evidence for the involvement of pre-Neoproterozoic zircons in the Arabian–Nubian Shield, a Neoproterozoic crustal tract that is generally regarded to be juvenile. The source and significance of these xenocrystic zircons are not clear. In an effort to better understand this problem, older and younger granitoids from the Egyptian basement complex were analyzed for chemical composition, SHRIMP U–Pb zircon ages, and Sm–Nd isotopic compositions. Geochemically, the older granitoids are metaluminous and exhibit characteristics of I-type granites and most likely formed in a convergent margin (arc) tectonic environment. On the other hand, the younger granites are peraluminous and exhibit the characteristics of A-type granites; these are post-collisional granites. The U–Pb SHRIMP dating of zircons revealed the ages of magmatic crystallization as well as the presence of slightly older, presumably inherited zircon grains. The age determined for the older granodiorite is 652.5 ± 2.6 Ma, whereas the younger granitoids are 595–605 Ma. Xenocrystic zircons are found in most of the younger granitoid samples; the xenocrystic grains are all Neoproterozoic, but fall into three age ranges that correspond to the ages of other Eastern Desert igneous rocks, viz. 710–690, 675–650 and 635–610 Ma. The analyzed granitoids have (+3.8 to +6.5) and crystallization ages, which confirm previous indications that the Arabian–Nubian Shield is juvenile Neoproterozoic crust. These results nevertheless indicate that older Neoproterozoic crust contributed to the formation of especially the younger granite magmas.     

Quick-look assessments to identify optimal CO2 EOR storage sites

A newly developed, multistage quick-look methodology allows for the efficient screening of an unmanageably large number of reservoirs to generate a workable set of sites that closely match the requirements for optimal CO₂ enhanced oil recovery (EOR) storage. The objective of the study is to quickly identify miscible CO₂ EOR candidates in areas that contain thousands of reservoirs and to estimate additional oil recovery and sequestration capacities of selected top options through dimensionless modeling and reservoir characterization. Quick-look assessments indicate that the CO₂ EOR resource potential along the US Gulf Coast is 4.7 billion barrels, and CO₂ sequestration capacity is 2.6 billion metric tons. In the first stage, oil reservoirs are screened and ranked in terms of technical and practical feasibility for miscible CO₂ EOR. The second stage provides quick estimates of CO₂ EOR potential and sequestration capacities. In the third stage, a dimensionless group model is applied to a selected set of sites to improve the estimates of oil recovery and storage potential using appropriate inputs for rock and fluid properties, disregarding reservoir architecture and sweep design. The fourth stage validates and refines the results by simulating flow in a model that describes the internal architecture and fluid distribution in the reservoir. The stated approach both saves time and allows more resources to be applied to the best candidate sites.     

Li-Sr-Nd isotope signatures of the plume and cratonic lithospheric mantle beneath the margin of the rifted Tanzanian craton (Labait)

Lithium concentrations and isotopic compositions of olivine and ⁸⁷Sr/⁸⁶Sr and ¹⁴³Nd/¹⁴⁴Nd of coexisting clinopyroxene from peridotite xenoliths from the Quaternary Labait volcano, Tanzania, document the influence of rift-related metasomatism on the ancient cratonic mantle. Olivines show negative correlations between Fo content and both δ⁷Li and Li concentrations. Olivines in iron-rich peridotites (Fo_85–87) have high Li concentrations (3.2–4.8 ppm) and heavy δ⁷Li (+5.2 to +6.6). In contrast, olivines in ancient, refractory peridotites have lower Li concentrations (∼2 ppm) and relatively light δ⁷Li (+2.6 to +3.5). This reflects mixing between ancient, refractory cratonic lithosphere and asthenosphere-derived rift magmas. A uniquely fertile, deformed, high-temperature garnet lherzolite, interpreted to be from the base of the lithosphere, has a ⁸⁷Sr/⁸⁶Sr of 0.7029 and ¹⁴³Nd/¹⁴⁴Nd of 0.51286, similar to HIMU oceanic basalts. It provides the best estimate of the Sr–Nd isotope composition of the upwelling mantle (i.e., plume, sensu lato) underlying this portion of the East African Rift, and is slightly less radiogenic compared to previous estimates of the plume that were based on rift basalts. Although elevated δ⁷Li are not exclusive to HIMU source regions, the data collectively indicate that the plume beneath Labait has HIMU characteristics in Sr, Nd and Li isotope composition. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Exurbia from the bottom-up: Confronting empirical challenges to characterizing a complex system

We describe empirical results from a multi-disciplinary project that support modeling complex processes of land-use and land-cover change in exurban parts of Southeastern Michigan. Based on two different conceptual models, one describing the evolution of urban form as a consequence of residential preferences and the other describing land-cover changes in an exurban township as a consequence of residential preferences, local policies, and a diversity of development types, we describe a variety of empirical data collected to support the mechanisms that we encoded in computational agent-based models. We used multiple methods, including social surveys, remote sensing, and statistical analysis of spatial data, to collect data that could be used to validate the structure of our models, calibrate their specific parameters, and evaluate their output. The data were used to investigate this system in the context of several themes from complexity science, including have (a) macro-level patterns; (b) autonomous decision making entities (i.e., agents); (c) heterogeneity among those entities; (d) social and spatial interactions that operate across multiple scales and (e) nonlinear feedback mechanisms. The results point to the importance of collecting data on agents and their interactions when producing agent-based models, the general validity of our conceptual models, and some changes that we needed to make to these models following data analysis. The calibrated models have been and are being used to evaluate landscape dynamics and the effects of various policy interventions on urban land-cover patterns.