Master of Science / Department of Geology / Matthew Totten / The Petersilie oil field in Ness County, Kansas produces out of the Mississippian System, a reservoir composed mainly of shallow water carbonates, at depths of around 4375 ft (1334 m). The lithology of the field ranges from limestone to dolomite, to interlaminated limestone-dolomite beds. Chert is commonly found throughout. Petersilie field lies to the west of the Central Kansas Uplift, and to the east of the Hugoton Embayment. The field saw much drilling activity in the 1960’s, when it reached a production peak of nearly 378,000 barrels of oil per year. Production declined swiftly after that until the late 1990’s, when waterflooding was successfully employed. In this study, a reservoir model was produced for the Mississippian as it occurs in Petersilie field using the Department of Energy’s EdBOAST reservoir modeling software, with the intent of providing a reference for future drilling activity in the Mississippian and determining reservoir characteristics that may have contributed to the effectiveness of waterflooding in this area. The reservoir model was checked by simulation with a companion reservoir simulator program, BOAST 98. Subsequent comparison of simulated and actual oil production curves demonstrates the reliability of well log and drill stem test data for the field and proves the reservoir model to be a good fit for the Mississippian in Petersilie. Production curve analysis of Petersilie indicates the field was an ideal candidate for waterflooding because it has a solution-gas drive mechanism. As the field approached depletion from primary recovery, oil saturations remained high. Petersilie also exhibits high porosity and good permeability. The BOAST software was found to be an effective and inexpensive means for understanding the Mississippian reservoir in central to south-central Kansas. It was determined that BOAST has potential for practical use by smaller independent oil companies targeting the Mississippian in Kansas.

Master of Science / Department of Geology / Abdelmoneam Raef / Amid increasing interest in geological sequestration of carbon dioxide (CO2), detailed rock formation characterization has emerged as priority to ensure successful sequestration. Utilizing recent advances in the field of 3D seismic attributes analysis, offers improved opportunities to provide more details when characterizing reservoir formations. In this study, several post-stack seismic attributes integrated with seismic modeling for highlighting critical structural elements and petrophysical facies variation of rock formations at Wellington and Anson-Bates fields, Sumner County, Kansas. A newly acquired 3D Seismic data set and several geophysical well logs are also used to achieve the objectives of this study. Results sought in this study are potentially important for understanding pathways for CO2 to migrate along. Seismic amplitude, coherency, and most negative curvature attributes were used to characterize the subsurface for structural effects on the rock formations of interest. These attributes detect multiple anomaly features that can be interpreted as small throw faults. However, in this study, there is a larger anomalous feature associated with the Mississippian formation that can be interpreted as a small throw fault or incised channel sand. Determining which of the two is very important for flow simulation models to be more exact. Modeling of the seismic was undertaken to help in the interpretation of the Mississippian amplitude anomaly. An artificial neural network, based on well log porosity cross-plots and three seismic attributes, was trained and implemented to yield a seismic petrophysical facies map. The neural network was trained using three volume seismic waveform attributes along with three wells with difference in well log porosity. A reworked lithofacies along small throw faults has been revealed based on comparing the seismic structural attributes and the seismic petrophysical facies. Arbuckle formation characterization was successful to a certain degree. Structural attributes showed multiple faults in the northern half of the survey. These faults are in agreement with known structure in the area associated with the Nemaha uplift. Further characterization of the Arbuckle was hindered by the lack of well data. This study emphasizes the need for greater attention to small-scale features when embarking upon characterization of a reservoir for CO2 based geosequestration.

Hydrocarbon deposits in the form of petroleum, natural gas, and natural gas hydrates occur offshore worldwide. Electromagnetic methods that measure the electrical resistivity of sediments can be used to map, assess, and monitor these resistive targets. In particular, quantitative assessment of hydrate content in marine deposits, which form within the upper few hundred meters of seafloor, is greatly facilitated by complementing conventional seismic methods with EM data. The North-East Pacific Time-series Undersea Network Experiment (NEPTUNE) is an underwater marine observatory that provides power and network connection to a host of instruments installed on the seafloor on the Cascadia Margin offshore Vancouver Island. The observatory’s aim is to provide a platform for very long-term studies in which access to data is available on a continuous basis. For this thesis project, a transient dipole-dipole time-domain EM system was constructed and deployed on the NEPTUNE network with the goal of long-term monitoring of a well-studied hydrate deposit in the area. The instrument includes a source transmitter of electrical current and individual receivers to measure small electric field variations. The instrument is powered by the NEPTUNE observatory and data can be collected remotely by connecting to the instrument through the web. Data collected so far from the instrument are consistent with a resistive structure. The best fitting model from 1D inversion is a 36 ± 3 m thick layer of 5.3 ± 0.3 Ωm resistivity, overlaying a less resistive 0.7 ± 0.1 Ωm halfspace. Average hydrate concentration, deduced with the aid of ODP-889 well-log derived Archie’s parameters, is approximately 72% of pore space in the resistive layer, consistent with the very high concentration of gas hydrates (~80%) recovered from seafloor cores. The weekly collection of data from the instrument shows that the resistive structure has changed little since monitoring began in October of 2010.

Hydrocarbon deposits in the form of petroleum, natural gas, and natural gas hydrates occur offshore worldwide. Electromagnetic methods that measure the electrical resistivity of sediments can be used to map, assess, and monitor these resistive targets. In particular, quantitative assessment of hydrate content in marine deposits, which form within the upper few hundred meters of seafloor, is greatly facilitated by complementing conventional seismic methods with EM data. The North-East Pacific Time-series Undersea Network Experiment (NEPTUNE) is an underwater marine observatory that provides power and network connection to a host of instruments installed on the seafloor on the Cascadia Margin offshore Vancouver Island. The observatory’s aim is to provide a platform for very long-term studies in which access to data is available on a continuous basis. For this thesis project, a transient dipole-dipole time-domain EM system was constructed and deployed on the NEPTUNE network with the goal of long-term monitoring of a well-studied hydrate deposit in the area. The instrument includes a source transmitter of electrical current and individual receivers to measure small electric field variations. The instrument is powered by the NEPTUNE observatory and data can be collected remotely by connecting to the instrument through the web. Data collected so far from the instrument are consistent with a resistive structure. The best fitting model from 1D inversion is a 36 ± 3 m thick layer of 5.3 ± 0.3 Ωm resistivity, overlaying a less resistive 0.7 ± 0.1 Ωm halfspace. Average hydrate concentration, deduced with the aid of ODP-889 well-log derived Archie’s parameters, is approximately 72% of pore space in the resistive layer, consistent with the very high concentration of gas hydrates (~80%) recovered from seafloor cores. The weekly collection of data from the instrument shows that the resistive structure has changed little since monitoring began in October of 2010.

Plane-layer geometry convection models remain useful for modelling planetary mantle dynamics however they yield significantly warmer mean temperatures than spherical shell models. For example, in a uniform property spherical shell with the same radius ratio, f, as the Earth's mantle; a bottom heating Rayleigh number, Ra, of 10^7 and a nondimensional internal heating rate, H, of 23 (arguably Earth-like values) are insufficient to heat the mean temperature, θ, above the mean of the non-dimensional boundary value temperatures (0.5), the temperature in a plane-layer model with no internal heating. This study investigates the impact of this geometrical effect in convection models featuring uniform and stratified viscosity. To address the effect of geometry, heat sinks are implemented to lower the mean temperature in 3D plane-layer isoviscous convection models. Over 100 models are analyzed, and their mean temperatures are used to derive a single equation for predicting θ, as a function of Ra, H and f in spherical and plane-layer systems featuring free-slip surfaces. The inclusion of first-order terrestrial characteristics is introduced to quantitatively assess the influence of system geometry on planetary scale simulations. Again, over 100 models are analyzed featuring a uniform upper mantle viscosity and a lower mantle viscosity that increases by a factor of 30 or 100. An effective Rayleigh number, Raη, is defined based on the average viscosity of the mantle. Equations for the relationship between θ, Raη, and H are derived for convection in a spherical shell with f = 0.547 and plane-layer geometries. These equations can be used to determine the appropriate heating rate for a plane-layer convection model to emulate spherical shell convection mean temperatures for effective Rayleigh numbers comparable to the Earth’s value and greater. Comparing cases with the same H and Raη, the increased lower mantle viscosity amplifies the mismatch in mean temperatures between spherical shell and plane-layer models. These findings emphasize the importance of adjusting heating rates in plane-layer geometry models and have important implications for studying convection with temperature-dependent parameters in plane-layer systems. The findings are particularly relevant to the study of convection in super-Earths where full spherical shell calculations remain intractable.

Flow instabilities occur at all scales in planetary systems. In this thesis we examine three cases of such instabilities, on three very different length scales. In the first part, we test the idea that Archean granite-greenstone belts (GGBs) form by crustal diapirism, or Rayleigh-Taylor instabilities. GGBs are characterized by large granitic domes (50-100 km in diameter) embedded in narrow keel-shaped greenstones. They are ubiquitous in Archean (> 2.5 Ga) terrains, but rare thereafter. We performed finite element calculations for a visco-elastic, temperature-dependent, non-Newtonian crust under conditions appropriate for the Archean, which show that dense low-viscosity volcanics overlying a felsic basement will overturn diapirically in as little as 10 Ma, displacing as much as 60 % of the volcanics to the lower crust. This surprisingly fast overturn rate suggests that diapiric overturn dominated crustal tectonics in the hot conditions of the Early Earth, becoming less important as the Earth cooled. Moreover, the deposition of large volumes of wet basaltic volcanics to the lower crust may provide the source for the formation of the distinctly Archean granitic rocks which dominate Earth's oldest continents. The second part examines the origin of Venusian coronae, circular volcanic features unique to Venus. Coronae are thought to result from small instabilities (diapirs) from the core-mantle boundary, which are typical of stagnant-lid convection. However, most young coronae are located in a region surrounded by long-lived hotspots, typical of a more active style of mantle convection. Using analogue experiments in corn syrup heated from below, we show that the co-existence of diapirs and long-lived mantle plumes are a direct consequence of the catastrophic overturn of the cold Venusian lithosphere thought to have occurred ~ 700 Ma ago. In the last part we analyze the thermal effect of fluid flow through a full-scale experiment testing clay and concrete tunnel seals in a Deep Geological Repository for nuclear was finite element software, we were able to show that the formation of fissures in the heated chamber between the two seals effectively limited heat flow, and could explain the discrepancy between the predicted and measured temperatures.

The primary objective of this study is to improve understanding of the mechanics of hydraulic fracturing in naturally-fractured reservoirs. The study focuses on enhancing the interpretation of hydraulic fracture-induced microseismic data using an S-wave Gaussian-beam method and numerical modeling techniques for interpretation. The S-wave Gaussian-beam method was comprehensively calibrated by synthetic and real data sets with different recording networks, and this showed the potential to retrieve additional microseismic data from hydraulic fracturing with linear receiver arrays. This approach could enhance current practice because a large number of induced events in these environments have very strong S-waves with P-wave amplitudes similar, or less than, background noise levels. The numerical study using the distinct element methods PFC2D and PFC3D was used to validate the understanding of the hydraulic fracturing mechanisms induced in laboratory and field fluid treatments in naturally-fractured reservoirs. This was achieved through direct comparison with the results of the geometry of hydraulic fractures and seismic source information (locations, magnitudes, and mechanisms) from both laboratory experiments and field observations. A suite of numerical models with fully-dynamic and hydro-mechanical coupling has been used to examine in detail the interaction between natural and induced fractures with the variations of the differential stresses and the orientations of the pre-fractures, and the relationship between the fluid front, the fracture tip, and the induced seismicity. The numerical results qualitatively agreed with the laboratory and field observations of the geometry of hydraulic fractures, confirmed the possible mechanics of new fracture development and their interactions with natural fractures, and illustrated the possible relationship between the fluid front and the fracture tip. The validated model could therefore help track the potential extent of induced fracturing in naturally-fractured reservoirs and the extent to which it can be detected by a microseismic monitoring array in order to assess the effectiveness of a hydraulic fracturing project.

Plane-layer geometry convection models remain useful for modelling planetary mantle dynamics however they yield significantly warmer mean temperatures than spherical shell models. For example, in a uniform property spherical shell with the same radius ratio, f, as the Earth's mantle; a bottom heating Rayleigh number, Ra, of 10^7 and a nondimensional internal heating rate, H, of 23 (arguably Earth-like values) are insufficient to heat the mean temperature, θ, above the mean of the non-dimensional boundary value temperatures (0.5), the temperature in a plane-layer model with no internal heating. This study investigates the impact of this geometrical effect in convection models featuring uniform and stratified viscosity. To address the effect of geometry, heat sinks are implemented to lower the mean temperature in 3D plane-layer isoviscous convection models. Over 100 models are analyzed, and their mean temperatures are used to derive a single equation for predicting θ, as a function of Ra, H and f in spherical and plane-layer systems featuring free-slip surfaces. The inclusion of first-order terrestrial characteristics is introduced to quantitatively assess the influence of system geometry on planetary scale simulations. Again, over 100 models are analyzed featuring a uniform upper mantle viscosity and a lower mantle viscosity that increases by a factor of 30 or 100. An effective Rayleigh number, Raη, is defined based on the average viscosity of the mantle. Equations for the relationship between θ, Raη, and H are derived for convection in a spherical shell with f = 0.547 and plane-layer geometries. These equations can be used to determine the appropriate heating rate for a plane-layer convection model to emulate spherical shell convection mean temperatures for effective Rayleigh numbers comparable to the Earth’s value and greater. Comparing cases with the same H and Raη, the increased lower mantle viscosity amplifies the mismatch in mean temperatures between spherical shell and plane-layer models. These findings emphasize the importance of adjusting heating rates in plane-layer geometry models and have important implications for studying convection with temperature-dependent parameters in plane-layer systems. The findings are particularly relevant to the study of convection in super-Earths where full spherical shell calculations remain intractable.

The primary objective of this study is to improve understanding of the mechanics of hydraulic fracturing in naturally-fractured reservoirs. The study focuses on enhancing the interpretation of hydraulic fracture-induced microseismic data using an S-wave Gaussian-beam method and numerical modeling techniques for interpretation. The S-wave Gaussian-beam method was comprehensively calibrated by synthetic and real data sets with different recording networks, and this showed the potential to retrieve additional microseismic data from hydraulic fracturing with linear receiver arrays. This approach could enhance current practice because a large number of induced events in these environments have very strong S-waves with P-wave amplitudes similar, or less than, background noise levels. The numerical study using the distinct element methods PFC2D and PFC3D was used to validate the understanding of the hydraulic fracturing mechanisms induced in laboratory and field fluid treatments in naturally-fractured reservoirs. This was achieved through direct comparison with the results of the geometry of hydraulic fractures and seismic source information (locations, magnitudes, and mechanisms) from both laboratory experiments and field observations. A suite of numerical models with fully-dynamic and hydro-mechanical coupling has been used to examine in detail the interaction between natural and induced fractures with the variations of the differential stresses and the orientations of the pre-fractures, and the relationship between the fluid front, the fracture tip, and the induced seismicity. The numerical results qualitatively agreed with the laboratory and field observations of the geometry of hydraulic fractures, confirmed the possible mechanics of new fracture development and their interactions with natural fractures, and illustrated the possible relationship between the fluid front and the fracture tip. The validated model could therefore help track the potential extent of induced fracturing in naturally-fractured reservoirs and the extent to which it can be detected by a microseismic monitoring array in order to assess the effectiveness of a hydraulic fracturing project.

Flow instabilities occur at all scales in planetary systems. In this thesis we examine three cases of such instabilities, on three very different length scales. In the first part, we test the idea that Archean granite-greenstone belts (GGBs) form by crustal diapirism, or Rayleigh-Taylor instabilities. GGBs are characterized by large granitic domes (50-100 km in diameter) embedded in narrow keel-shaped greenstones. They are ubiquitous in Archean (> 2.5 Ga) terrains, but rare thereafter. We performed finite element calculations for a visco-elastic, temperature-dependent, non-Newtonian crust under conditions appropriate for the Archean, which show that dense low-viscosity volcanics overlying a felsic basement will overturn diapirically in as little as 10 Ma, displacing as much as 60 % of the volcanics to the lower crust. This surprisingly fast overturn rate suggests that diapiric overturn dominated crustal tectonics in the hot conditions of the Early Earth, becoming less important as the Earth cooled. Moreover, the deposition of large volumes of wet basaltic volcanics to the lower crust may provide the source for the formation of the distinctly Archean granitic rocks which dominate Earth's oldest continents. The second part examines the origin of Venusian coronae, circular volcanic features unique to Venus. Coronae are thought to result from small instabilities (diapirs) from the core-mantle boundary, which are typical of stagnant-lid convection. However, most young coronae are located in a region surrounded by long-lived hotspots, typical of a more active style of mantle convection. Using analogue experiments in corn syrup heated from below, we show that the co-existence of diapirs and long-lived mantle plumes are a direct consequence of the catastrophic overturn of the cold Venusian lithosphere thought to have occurred ~ 700 Ma ago. In the last part we analyze the thermal effect of fluid flow through a full-scale experiment testing clay and concrete tunnel seals in a Deep Geological Repository for nuclear was finite element software, we were able to show that the formation of fissures in the heated chamber between the two seals effectively limited heat flow, and could explain the discrepancy between the predicted and measured temperatures.