In this research project the author introduces the use of geometrical modelling techniques alongside geostatistical methods to model a stratigraphic orebody and to present a graphics display system developed as a fIrst step towards a general integrated system for computer aided design and planning in mining. Geometrical modelling techniques and geostatistical methods are combined to carry out the process of modelling a stratigraphic orebody. From a mining point of view, there are two main features of interest in a stratigraphic ore-body: a-The modelling of the geometry of the orebody. b- The modelling (estimation) of the physical properties (grades, etc ... ) of the orebody. The first feature is the subject of this research project. Modelling methods and techniques developed elsewhere and for different applications, such as Computer Aided Design, have been applied successfully to model the geometry of stratigraphic orebodies. The modelling process consists of the applications of surface modelling techniques to represent the hangingwall and the footwall of the stratigraphic orebody and thereby to produce the space where the physical properties are geostatistically to be estimated. The graphics display system is presented to highlight the use of computer graphics techniques to communicate graphically all sorts of information concerning the modelling of stratigraphic orebodies and also to display the end product of the modelling process, such as cross-sections, plane-sections, wireframe and solid models of the orebody. The graphics system itself is part of a computer based system for mine design and planning similar to computer aided design systems used mainly in the manufacturing industry. The presentation of the research project in this thesis started by the review of the literature of some existing ore reserves estimation methods in the mineral industry, particularly geostatistical methods. Then an overview and the scope of this research project have been given. The second chapter describes the type of data which could be encountered while building a geometrical model of a stratigraphic orebody and a description of data from a nickel vein deposit used as a case study for this research project. The accumulations have been estimated geostatistically subject to geometrical control. The geometrical control concept and surface modelling techniques are presented in chapter three together with the numerical application of modelling a nickel vein deposit using two different surface modelling techniques. Chapter four describes the graphics display system developed to display several geometric features of stratigraphic orebodies in two and three dimensions. The summary of this research project and some concluding remarks are given in chapter five.

Thesis (Ph. D.)--University of Washington, 2001. / Vita. Includes bibliographical references (leaves 110-132).

The modern plate tectonic paradigm provides a predictive model to understand what mediates dynamic processes at both plate margins and intraplate settings. At some locations on the Earth, the geological record provides evidence of apparent violations of this theoretical framework. In this dissertation, I examine a region on the rifted continental margin of Eastern North America, where at least four distinct episodes of magmatism occurred (in the Late Triassic, Early Jurassic, Early Cretaceous) since the onset of rifting and ultimate breakup of the supercontinent Pangea. It also coincides with a present-day low seismic velocity anomaly in the upper mantle. No other region on the Eastern North American Margin has a record of such anomalous dynamic processes occurring and persisting for more than 200 Myr. In this dissertation, I primarily use zircon U-Pb geochronology to establish the basic chronological framework in which magmatic and magmatic-hydrothermal systems in this region existed and persisted, establishing the temporal parameter space in which it will be possible to test geodynamic mechanisms for their formation. In Chapter 2, I use ultra-high precision zircon U-Pb geochronology via Chemical Abrasion-Isotope Dilution-Thermal Ionization Mass Spectrometry (CA-ID-TIMS) to test whether the onset of magmatism in the largest igneous body in this region (the White Mountain Batholith) is linked to the eruption of the vast flood basalts within the Central Atlantic Magmatic Province (CAMP) and whether its apparent duration is indeed > 50 Myr, as previous workers using whole rock K-Ar and mineral Ar-Ar geochronologic methods suggest. My work uncovered a previously unknown episode of rift-related magmatism in the region that precedes the both the CAMP and the emplacement of the White Mountain Batholith by 3 – 5 Myr. In Chapter 3, I use a combination of high-precision zircon U-Pb geochronology and absolute plate motion models to test whether the Cretaceous igneous province in this region resulted from hot spot magmatism as North America moved over the purported Great Meteor Hotspot. These results cannot falsify the hotspot hypothesis and the new zircon U-Pb ages therefore provide the best available chronological constraints for one of the longest-lived hot spot tracks on the Earth. In Chapter 4, motivated by the confirmation of age discrepancies between low- and intermediate-temperature chronometers and the zircon U-Pb ages presented in Chapter 2, I use a combination of both CA-ID-TIMS and LA-ICP-MS zircon U-Pb geochronologic techniques to place constraints on the timing and duration of magmatism for the Late Triassic-Early Jurassic province. The results of this chapter demonstrate that the magmatism of the White Mountain Magma Series occurs in discrete pulses through much of the Jurassic. Together with zircon Hf isotopic analyses from select samples, I synthesize these age results and construct a hypothesis testing framework in which it will be possible for future investigators to unravel the geodynamic complexities in this region. I provide recommendations for future work and emphasize the need for unified approaches coupling geochronology, geochemistry, and geophysics, to test the range of possible mechanisms responsible for these episodes of anomalous tectono-magmatic activity.

Leakage of stored bulk phase CO₂ is one risk for sequestration in deep saline aquifers. As the less dense CO₂ migrates upward within a storage formation or in layers above the formation, the security of its storage depends upon the trapping mechanisms that counteract the migration. The trapping mechanism motivating this research is local capillary trapping (LCT), which occurs during buoyancy-driven migration of bulk phase CO₂ within a saline aquifer with spatially heterogeneous petrophysical properties. When a CO₂ plume rising by buoyancy encounters a region where capillary entry pressure is locally larger than average, CO₂ accumulates beneath the region. One benefit of LCT, applied specifically to CO₂ sequestration and storage, is that saturation of stored CO₂ phase is larger than the saturation for other permanent trapping mechanisms. Another potential benefit is security: CO₂ that occupies local capillary traps remains there, even if the overlying formation that provides primary containment were to be compromised and allow leakage. Most work on LCT has involved numerical simulation (Saadatpoor 2010, Ganesh 2012); the research work presented here is a step toward understanding local capillary trapping at the bench scale. An apparatus and set of fluids are described which allow examining the extent of local capillary trapping, i.e. buoyant nonwetting phase immobilization beneath small-scale capillary barriers, which can be expected in typical heterogeneous storage formation. The bench scale environment analogous to CO₂ and brine in a saline aquifer is created in a quasi-two dimensional experimental apparatus with dimension of 63 cm by 63 cm by 5 cm, which allows for observation of plume migration with physically representative properties but at experimentally convenient ambient conditions. A surrogate fluid pair is developed to mimic the density, viscosity and interfacial tension relationship found at pressure and temperature typical of storage aquifers. Porous media heterogeneity, pressure boundary conditions, migration modes of uprising nonwetting phase, and presence of fracture/breach in the capillary barrier are studied in series of experiments for their influences on LCT. A variety of heterogeneous porous media made of a range of sizes of loosely packed silica beads are used to validate and test the persistence of local capillary trapping mechanism. By adjusting the boundary conditions (fluid levels in reservoirs attached to top and to bottom ports of the apparatus), the capillary pressure gradient across the domain was manipulated. Experiments were conducted with and without the presence of fracture/potential leakage pathway in the capillary seal. The trapped buoyant phase remained secure beneath the local capillary barriers, as long as the effective capillary pressure exerted by the trapped phase (proportional to column height of the phase) is smaller than the capillary entry pressure of the barrier. The local capillary trapping mechanism remained persistent even under forced imbibition, in which a significantly higher hydraulic potential gradient, and therefore a larger gradient in capillary pressure, was applied to the system. The column height of buoyant fluid that remained beneath the local capillary barrier was smaller by a factor corresponding to the increase in capillary pressure gradient. Mimicking a breach of the caprock by opening valves at the top of the apparatus allowed buoyant mobile phase held beneath the valves to escape, but buoyant phase held in local traps at saturations above residual, and therefore potentially mobile, was undisturbed. This work provides systematic validation of a novel concept, namely the long-term security of CO₂ that fills local (small-scale) capillary traps in heterogeneous storage formations. Results from this work reveal the first ever unequivocal experimental evidence on persistence of local capillary trapping mechanism. Attempts to quantify the nonwetting phase saturation and extent of LCT persistence serve as the initial steps to potentially reduce the risks associated with long-term storage security. / text

Master of Science / Department of Geology / Matthew W. Totten / The area of study for this project is the Leach Field, which is located in Jackson County, Kansas. Production in the Leach Field has historically been disappointing, with 388,787 barrels of oil being produced since the field’s discovery in 1963 (KGS, 2015). Production of the field has been highly variable, with only 20,568 barrels of oil being produced in the last 20 years. Economic and other concerns that have impacted production and production rates of the field include: low oil prices soon after its discovery, numerous changes of ownership, and lack of significant production infrastructure in the area. Stroke of Luck Energy & Exploration, LLC. has recently purchased the majority of the leases and wells in the Leach Field, and is reestablishing the field as a productive oil field. Plans include: washing down several plugged and abandoned wells, and drill new wells to increase production in the field. The goal of this study was to determine the major geologic factors controlling reservoir quality in the Hunton and Viola Limestone Formations in the Leach Field, so that a future exploration model can be developed to help increase and stabilize the field's overall production. This model was created by applying several testing methods including: well logging analysis, microscope analysis, and subsurface mapping. Based on these results it was determined that the quality of the reservoir rocks is controlled by the degree of dolomitizaiton in both formations. Reservoir quality is as important as structure in determining well productivity in the Leach Field.

The injection operator for a carbon dioxide sequestration project must control the reservoir and associated pore space within the project boundaries to allow for orderly development of the storage facility. A large number of interest owners within a project area is likely to make reaching unanimous agreement among all owners of pore space unlikely, and thus control of the reservoir difficult. In order to facilitate geologic sequestration of carbon dioxide on privately owned land in the United States, or on land for which the minerals or pore space are privately owned, a scheme for aggregating the ownership of pore space is needed. To allow geologic sequestration projects to move forward with less than unanimous consent of interest owners, states can employ various methods of aggregating pore space ownership. This paper examines oil and gas unitization statues and statutes creating groundwater districts to find legislative regimes useful for achieving pore space ownership aggregation. Among the approaches discussed, aggregation of pore space ownership through a unitization model is the most likely choice. Taking that one step further and setting up new unit operating agreements for enhanced oil recovery to serve as a repository for incremental geologic sequestration, and eventual full sequestration activities, provides a firm path toward reducing carbon dioxide emissions while respecting property rights. This paper also compares the few existing pore space aggregation statutes in the United States, which achieve aggregation of pore space ownership through either unitization or eminent domain. The state that appears to be the best equipped to deal with aggregation of pore space ownership is Wyoming. Wyoming has been a leader in developing legislation to deal with pore space ownership before other states. North Dakota and Utah are also very well situated to move forward with carbon sequestration activities. / text

No description available.

<P> This study investigates ways to improve the accuracy of 3D geologic models by assessing the impact of data quality, grid complexity, data quantity and distribution, interpolation algorithm and program selection on model accuracy. The first component of this research examines the impact of variable quality data on 3D model outputs and presents a new methodology to optimize the impact of high quality data, while minimizing the impact of low quality data on the model results. This 'Quality Weighted' modelling approach greatly improves model accuracy when compared with un-weighted models. </p> <p> The second component of the research assesses the variability and influence of data quantity, data distribution, algorithm selection, and program selection on the accuracy of 3D geologic models. A series of synthetic grids representing environments of varying complexity were created from which data subsets were extracted using specially developed MA TLAB scripts. The modelled data were compared back to the actual synthetic values and statistical tests were conducted to quantify the impact of each variable on the accuracy of the model predictions. The results indicate that grid complexity is the predominant control on model accuracy, more data do not necessarily produce more accurate models, and data distribution is particularly important when relatively simple environments are modelled. A major finding of this study is that in some situations, the software program selected for modelling can have a greater influence on model accuracy than the algorithm used for interpolation. When modelling spatial data there is always a high level of uncertainty, especially in subsurface environments where the unit(s) of interest are defined by data only available in select locations. The research presented in this thesis can be used to guide the selection of modelling parameters used in 3D subsurface investigations and will allow the more effective and efficient creation of accurate 3D models. </p> / Thesis / Doctor of Philosophy (PhD)

Previamente à cada estimativa ou simulação geoestatística os domínios geológicos do depósito devem ser modelados, o que tradicionalmente é feito de forma manual por um geomodelador, em um processo laborioso, demorado e subjetivo. Por essa razão novas técnicas conhecidas como métodos implícitos veem surgindo. Essas técnicas fornecem algoritmos que substituem o processo de digitalização manual dos métodos explícitos por alguma forma de procedimento automático. Essa dissertação visita alguns métodos implícitos bem estabelecidos com atenção especial à modelagem geológica implícita com funções distância assinalada. Um estudo de caso em um banco de dados real é apresentado e a aplicabilidade do método discutida. Embora não substitua por completo um geomodelador experiente, o método provou ser capaz de gerar modelos geológicos semi-automáticos realistas a partir dos dados amostrais, e se mostra útil principalmente nas fases iniciais da pesquisa mineral. / Prior to every geostatistical estimation or simulation study there is a need for delimiting the geologic domains of the deposit, which is traditionally done manually by a geomodeler in a laborious, time consuming and subjective process. For this reason, novel techniques referred to as implicit modelling have appeared. These techniques provide algorithms that replace the manual digitization process of the traditional methods by some form of automatic procedure. This dissertation covers a few well established implicit methods currently available with special attention to the signed distance function methodology. A case study based on a real dataset was performed and its applicability discussed. Although it did not replace an experienced geomodeler, the method proved to be capable in creating semi-automatic geological models from the sampling data, especially in the early stages of exploration.

Previamente à cada estimativa ou simulação geoestatística os domínios geológicos do depósito devem ser modelados, o que tradicionalmente é feito de forma manual por um geomodelador, em um processo laborioso, demorado e subjetivo. Por essa razão novas técnicas conhecidas como métodos implícitos veem surgindo. Essas técnicas fornecem algoritmos que substituem o processo de digitalização manual dos métodos explícitos por alguma forma de procedimento automático. Essa dissertação visita alguns métodos implícitos bem estabelecidos com atenção especial à modelagem geológica implícita com funções distância assinalada. Um estudo de caso em um banco de dados real é apresentado e a aplicabilidade do método discutida. Embora não substitua por completo um geomodelador experiente, o método provou ser capaz de gerar modelos geológicos semi-automáticos realistas a partir dos dados amostrais, e se mostra útil principalmente nas fases iniciais da pesquisa mineral. / Prior to every geostatistical estimation or simulation study there is a need for delimiting the geologic domains of the deposit, which is traditionally done manually by a geomodeler in a laborious, time consuming and subjective process. For this reason, novel techniques referred to as implicit modelling have appeared. These techniques provide algorithms that replace the manual digitization process of the traditional methods by some form of automatic procedure. This dissertation covers a few well established implicit methods currently available with special attention to the signed distance function methodology. A case study based on a real dataset was performed and its applicability discussed. Although it did not replace an experienced geomodeler, the method proved to be capable in creating semi-automatic geological models from the sampling data, especially in the early stages of exploration.