AMPLITUDE-VERSUS-ANGLE ANALYSIS AND WIDE-ANGLE-INVERSION OF CROSSWELL SEISMIC DATA IN A CARBONATE RESERVOIR

Crosswell data set contains a range of angles limited only by the geometry of the source and receiver configuration, the separation of the boreholes and the depth to the target. However, the wide angles reflections present in crosswell imaging result in amplitude-versus-angle (AVA) features not usually observed in surface data. These features include reflections from angles that are near critical and beyond critical for many of the interfaces; some of these reflections are visible only for a small range of angles, presumably near their critical angle. High-resolution crosswell seismic surveys were conducted over a Silurian (Niagaran) reef at two fields in northern Michigan, Springdale and Coldspring. The Springdale wells extended to much greater depths than the reef, and imaging was conducted from above and from beneath the reef. Combining the results from images obtained from above with those from beneath provides additional information, by exhibiting ranges of angles that are different for the two images, especially for reflectors at shallow depths, and second, by providing additional constraints on the solutions for Zoeppritz equations. Inversion of seismic data for impedance has become a standard part of the workflow for quantitative reservoir characterization. Inversion of crosswell data using either deterministic or geostatistical methods can lead to poor results with phase change beyond the critical angle, however, the simultaneous pre-stack inversion of partial angle stacks may be best conducted with restrictions to angles less than critical. Deterministic inversion is designed to yield only a single model of elastic properties (best-fit), while the geostatistical inversion produces multiple models (realizations) of elastic properties, lithology and reservoir properties. Geostatistical inversion produces results with far more detail than deterministic inversion. The magnitude of difference in details between both types of inversion becomes increasingly pronounced for thinner reservoirs, particularly those beyond the vertical resolution of the seismic. For any interface imaged from above and from beneath, the results AVA characters must result from identical contrasts in elastic properties in the two sets of images, albeit in reverse order. An inversion approach to handle both datasets simultaneously, at pre-critical angles, is demonstrated in this work. The main exploration problem for carbonate reefs is determining the porosity distribution. Images of elastic properties, obtained from deterministic and geostatistical simultaneous inversion of a high-resolution crosswell seismic survey were used to obtain the internal structure and reservoir properties (porosity) of Niagaran Michigan reef. The images obtained are the best of any Niagaran pinnacle reef to date.

USE OF AN ELECTRICAL IMPEDANCE TOMOGRAPHY METHOD TO DETECT AND TRACK FRACTURES IN A GELATIN MEDIUM

Electrical impedance tomography is applied to the problem of detecting, locating, and tracking fractures in ballistics gelatin. The hardware developed is intended to be physically robust and based on off-the-shelf hardware. Fractures were created in two separate ways: by shooting a .22 caliber bullet into the gelatin and by injecting saline solution into the gelatin. The .22 caliber bullet created an air gap, which was seen as an increase in resistivity. The saline solution created a fluid filled gap, which was seen as a decrease in resistivity. A double linear array was used to take data for each of the fracture mechanisms and a two dimensional cross section was inverted from the data. The results were validated by visually inspecting the samples during the fracture event. It was found that although there were reconstruction errors present, it was possible to reconstruct a representation of the resistive cross section. Simulations were performed to better understand the reconstructed cross-sections and to demonstrate the ability of a ring array, which was not experimentally tested.

Exploration of Ground Penetrating Radar and Time Domain Reflectometry Methods for the Determination of Pavement Dielectric Constant

Traditionally, densities of newly built roadways are checked by direct sampling (cores) or by nuclear density gauge measurements. For roadway engineers, density of asphalt pavement surfaces is essential to determine pavement quality. Unfortunately, field measurements of density by direct sampling or by nuclear measurement are slow processes. Therefore, I have explored the use of rapidly-deployed ground penetrating radar (GPR) as an alternative means of determining pavement quality. The dielectric constant of pavement surface may be a substructure parameter that correlates with pavement density, and can be used as a proxy when density of asphalt is not known from nuclear or destructive methods. The dielectric constant of the asphalt can be determined using ground penetrating radar (GPR). In order to use GPR for evaluation of road surface quality, the relationship between dielectric constants of asphalt and their densities must be established. Field measurements of GPR were taken at four highway sites in Houghton and Keweenaw Counties, Michigan, where density values were also obtained using nuclear methods in the field. Laboratory studies involved asphalt samples taken from the field sites and samples created in the laboratory. These were tested in various ways, including, density, thickness, and time domain reflectometry (TDR). In the field, GPR data was acquired using a 1000 MHz air-launched unit and a ground-coupled unit at 200 and 500 MHz. The equipment used was owned and operated by the Michigan Department of Transportation (MDOT) and available for this study for a total of four days during summer 2005 and spring 2006. The analysis of the reflected waveforms included “routine” processing for velocity using commercial software and direct evaluation of reflection coefficients to determine a dielectric constant. The dielectric constants computed from velocities do not agree well with those obtained from reflection coefficients. Perhaps due to the limited range of asphalt types studied, no correlation between density and dielectric constant was evident. Laboratory measurements were taken with samples removed from the field and samples created for this study. Samples from the field were studied using TDR, in order to obtain dielectric constant directly, and these correlated well with the estimates made from reflection coefficients. Samples created in the laboratory were measured using 1000 MHz air-launched GPR, and 400 MHz ground-coupled GPR, each under both wet and dry conditions. On the basis of these observations, I conclude that dielectric constant of asphalt can be reliably measured from waveform amplitude analysis of GJPR data, based on the consistent agreement with that obtained in the laboratory using TDR. Because of the uniformity of asphalts studied here, any correlation between dielectric constant and density is not yet apparent.

DOING MORE WITH SHORT PERIOD DATA: DETERMINING MAGNITUDES FROM CLIPPED AND OVER-RUN SEISMIC DATA AT MOUNT ST. HELENS

How can we calculate earthquake magnitudes when the signal is clipped and over-run? When a volcano is very active, the seismic record may saturate (i.e., the full amplitude of the signal is not recorded) or be over-run (i.e., the end of one event is covered by the start of a new event). The duration, and sometimes the amplitude, of an earthquake signal are necessary for determining event magnitudes; thus, it may be impossible to calculate earthquake magnitudes when a volcano is very active. This problem is most likely to occur at volcanoes with limited networks of short period seismometers. This study outlines two methods for calculating earthquake magnitudes when events are clipped and over-run. The first method entails modeling the shape of earthquake codas as a power law function and extrapolating duration from the decay of the function. The second method draws relations between clipped duration (i.e., the length of time a signal is clipped) and the full duration. These methods allow for magnitudes to be determined within 0.2 to 0.4 units of magnitude. This error is within the range of analyst hand-picks and is within the acceptable limits of uncertainty when quickly quantifying volcanic energy release during volcanic crises. Most importantly, these estimates can be made when data are clipped or over-run. These methods were developed with data from the initial stages of the 2004-2008 eruption at Mount St. Helens. Mount St. Helens is a well-studied volcano with many instruments placed at varying distances from the vent. This fact makes the 2004-2008 eruption a good place to calibrate and refine methodologies that can be applied to volcanoes with limited networks.

Integrated Geophysical Investigation of the St. James Fault Complex: A Case Study

We noninvasively detected the characteristics and location of a regional fault in an area of poor bedrock exposure complicated by karst weathering features in the subsurface. Because this regional fault is associated with sinkhole formation, its location is important for hazard avoidance. The bedrock lithologies on either side of the fault trace are similar; hence, we chose an approach that capitalized on the complementary strengths of very low frequency (VLF) electromagnetic, resistivity, and gravity methods. VLF proved most useful as a first-order reconnaissance tool, allowing us to define a narrow target area for further geophysical exploration. Fault-related epikarst was delineated using resistivity. Ultimately, a high-resolution gravity survey and subsequent inverse modeling using the results of the resistivity survey helped to further constrain the location and approximate orientation of the fault. The combined results indicated that the location of the fault trace needed to be adjusted 53 m south of the current published location and was consistent with a north-dipping thrust fault. Additionally, a gravity low south of the fault trace agreed with the location of conductive material from the resistivity and VLF surveys. We interpreted these anomalies to represent enhanced epikarst in the fault footwall. We clearly found that a staged approach involving a progression of methods beginning with a reconnaissance VLF survey, followed by high-resolution gravity and electrical resistivity surveys, can be used to characterize a fault and fault-related karst in an area of poor bedrock surface exposure.

Triggering of volcanic activity by large earthquakes

Statistical analyses of temporal relationships between large earthquakes and volcanic eruptions suggest seismic waves may trigger eruptions even over great (>1000 km) distances, although the causative mechanism is not well constrained. In this study the relationship between large earthquakes and subtle changes in volcanic activity was investigated in order to gain greater insight into the relationship between dynamic stresses propagated by surface waves and volcanic response. Daily measurements from the Ozone Monitoring Instrument (OMI), onboard the Aura satellite, provide constraints on volcanic sulfur-dioxide (SO2) emission rates as a measure of subtle changes in activity. Time series of SO2 emission rates were produced from OMI data for thirteen persistently active volcanoes from 1 October 2004 to 30 September 2010. In order to quantify the affect of earthquakes at teleseismic distances, we modeled surface-wave amplitudes from the source mechanisms of moment magnitude (Mw) ≥7 earthquakes, and calculated the Peak Dynamic Stress (PDS). We assessed the influence of earthquakes on volcanic activity in two ways: 1) by identifying increases in the SO2 time series data and looking for causative earthquakes and 2) by examining the average emission rate before and after each earthquake. In the first, the SO2 time series for each volcano was used to calculate a baseline threshold for comparison with post-earthquake emission. Next, we generated a catalog of responses based on sustained SO2 emission increases above this baseline. Delay times between each SO2 response and each prior earthquake were analyzed using both the actual earthquake catalog, and a randomly generated catalog of earthquakes. This process was repeated for each volcano. Despite varying multiple parameters, this analysis did not demonstrate a clear relationship between earthquake-generated PDS and SO2 emission. However, the second analysis, which was based on the occurrence of large earthquakes indicated a response at most volcanoes. Using the PDS calculations as a filtering criterion for the earthquake catalog, the SO2 mass for each volcano was analyzed in 28-day windows centered on the earthquake origin time. If the average SO2 mass after the earthquake was greater than an arbitrary percentage of pre-earthquake mass, we identified the volcano as having a response to the event. This window analysis provided insight on what type of volcanic activity is more susceptible to triggering by dynamic stress. The volcanoes with very open systems included in this study, Ambrym, Gaua, Villarrica, Erta Ale and, Turrialba, showed a clear response to dynamic stress while the volcanoes with more closed systems, Merapi, Semeru, Fuego, Pacaya, and Bagana, showed no response.

Paleomagnetic Study of the Portage Lake Volcanics exposed in the Quincy Mine

A detailed paleomagnetic and rock-magnetic investigation was conducted on thirty six basaltic flows of the ~1095 Ma Portage Lake Volcanics. The flows were sampled along the East Adit of the Quincy Mine (Hancock, MI). Thirty two flows yielded well-defined primary magnetization directions carried by magnetite. A secondary magnetization component carried by hematite was also found in twenty nine flows. After correction for serial correlation between the flows, nineteen independent mean directions were calculated. The corresponding paleomagnetic pole is located at 25.5 °N, 182.1 °W (A95 = 3.5°). The new pole overlaps with the pole from the ~1087 Ma Lake Shore Traps suggesting a standstill of the North American plate during that time period. The low angular dispersion of virtual geomagnetic poles (S = 7.9°) suggests that the flows were erupted within a short time period, or that the strength of geomagnetic secular variation was lower than that of the recent field.

Seismic observations of regional gas blowdown at Teal South

A re-examination of seismic time-lapse data from the Teal South field provides support for a previously proposed model of regional pressure decline and the associated liberation of gas from nearby reservoirs due to the production from the only reservoir among them that is under production. The use of a specific attribute, instantaneous amplitude, and a series of time slices, however, provides increased detail in understanding fluid migration into or out of the reservoirs, and the path taken by pressure changes across faults. The regional decrease of pressure due to production in one reservoir has dramatic effects in nearby untapped reservoirs, one of which appears to exhibit evidence for the escape, and possible re-trapping nearby, of hydrocarbons from a spill point. The influx of water into the producing reservoir is also evidenced by a decrease in amplitude at one end of the oil-water contact.

Preliminary Construction of a Photo-Electric Cell Concentrator

The separation of the valuable portion from the waste portion of an ore is an individual problem for every ore. However, the various methods for accomplishing this end, more or less classify themselves by the physical properties of the constituents of the ore. Most of the properties of minerals have been utilized in some way or other to affect the separation of the valuable from the invaluable parts. Practically nothing has been done so far with color and luster to attain this purpose. It is believed that the photo—electric cell could also be used in concentrating a certain class of ores which are not well suited to other methods.

The Aplitic Body Bordering Roosevelt Drive South of Butte, Montana

The area studied consists of a large body of aplite situated 12 miles south of Butte. It extends eastward almost to the Butte-Pipestone road, southward to Lime Kiln Hill and westward a distance of 3000 feet. The petrographic work was supplemented by a map of the area. Special attention was given to those frequent "Limonite Specks" found in the aplite.

The Electrolytic Deposition of Chromium Upon Aluminum

The purpose or this investigation is primarily to determine the best conditions for plating chromium on aluminum. The work was carried out with the hope of obtaining coherent deposits, and of determin­ing the conditions under which such deposits may be duplicated.

The Effect of Hydrogen Overvoltage on the Electrolytic Deposition of Zinc

In the deposition of metallic zinc by electro­lysis from neutral or acid solution, little difficulty is experienced provided certain impurities are absent from the electrolyte. The use of the process has long been considered as a potential source, patents on the process having been issued as early as 1880. However, the early experimenters failed to realize the importance of impurities in the electrolyte, and for this reason, the process suffered several severe setbacks when commercial plants were built.

Joint Planes in the Rhyolites of the Butte District, Montana

The problem herein discussed deals with the pointing planes found in the area of rhyolite located in the northwestern portion of the Butte District. The question to be determined was whether or not the pointing planes in the rhyolites could be class­ified with the Butte Ore fissure systems.

A Bentonite Deposit near Warm Springs, Montana

This report deals with a bentonite deposit recently developed, approximately seven miles northeast of Warm Springs, Montana. A group of claims have been staked on the deposit and are owned by the Lincoln Mining Company of Anaconda, Montana. The company also has several claims prospected for silver one mile from its present site of operations, but the silver prospects have failed to produce. The bentonite deposit was discovered incidentally during the course of other development work, and at present two adits have been driven into the side of a mountain, each crosscutting a vein-like mass of bentonite varying from two to three feet in width.

Some Notes on the Composition of "Copper Pitch"

A supply of so-called "copper pitch" ore was received by the Montana Bureau of Mines and Geology in response to a request by them from a resident of Kalispell, who had previously sent a specimen to the Bureau for a mineralogical analysis. Since this material was little known and had apparently received but little study under a reflec­ting microscope, it was thought that such a study might throw some light on the mineralogical and chemical composition of the material.